[R] multi-dimension array of raw

[Gerald Jansen]

I would like to store and manipulate large sets of marker genotypes 
compactly using "raw" data arrays. This works fine for vectors or
matrices, but I run into the error shown in the example below as soon
as I try to use 3 dimensional arrays (eg. animal x marker x allele).

> a <- array(as.raw(1:6),c(2,3))
> a
     [,1] [,2] [,3]
[1,]   01   03   05
[2,]   02   04   06
> a[1,] <- raw(3)
> a
     [,1] [,2] [,3]
[1,]   00   00   00
[2,]   02   04   06
> b <- array(as.raw(1:6),c(1,2,3))
> b[1,,]
     [,1] [,2] [,3]
[1,]   01   03   05
[2,]   02   04   06
> b[1,1,] <- raw(3)
Error: incompatible types (from raw to raw) in array subset assignment

I can work around this with computed indices, but I wonder if this is
expected behaviour. 

Gerald Jansen