Pode umraster Objeto (em R) tem camadas de modo diferente (tipo de dados)?

Em face disso parece que somos sempre forçados a um tipo:

library(raster)
## create a SpatialPixelsDataFrame with (trivially) two different "layer" types
d <- data.frame(expand.grid(x = 1:10, y = 2:11), z = 1:100, a = sample(letters, 100, replace = TRUE), stringsAsFactors = FALSE)
coordinates(d) <- 1:2
gridded(d) <- TRUE

## now coerce this to a raster brick or stack and our "a" is crushed to numeric NA
all(is.na(getValues(brick(d)[[2]])))
[1] TRUE

Existe algo como um rasterDataFrame?

Além disso, observe que, presumivelmente, não podemos usar os fatores de R, já que os dados raster @ são uma matriz ou, de outro modo, coagidos para numeric / integer. Estou esquecendo de algo?