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

______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html

Reply via email to