[GenABEL-dev] databel float type

[Frank, Alvaro Jesus]

Hi all,

what is wrong with this? R:
---------------------------------------------------------------------------------------------


library(DatABEL)

n = 100 # number of individuals
l = 2    # number of covariates
m = 10 # number of snps
t = 5  # number of traits

set.seed(1001)
runif(3)

XL <- matrix(rnorm((l+1)*n),ncol=(l+1)) # first column should be ones (intercept)
for(i in 1:n*l){ if(sample(1:10,1) > 5) XL[i]=0/0}

colnames(XL) <- c("intercept", paste("cov",1:l,sep=""))
rownames(XL) <- paste("ind",1:n,sep="")
XL_db <- matrix2databel(XL,filename="XL",type="FLOAT")



---------------------------------------------------------------------------------------------
The error that I get is this:

file contains data of unknown type 5
Error in !.Call("assignDoubleMatrix", x at data, newi, newj, as.double(value),  :
  invalid argument type

I cannot create any more databel flaot files....


A Frank
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/genabel-devel/attachments/20140728/c43c330d/attachment.html>