This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem imputing (large) file into GenABEL

I have a very large imputation file (363GB in total) so I need to convert to a filevector format with GenABEL. My imputation file is in IMPUTE2 format and looks like so standard .imputed file:

   21 rs142159069    45000079   T T 1 0 0  etc
   etc

In other words, each row is a SNP with columns: chr, rsid, position, alleleA, alleleB, and then probabilities of being AA,AB or BB for each sample. There are 80 million SNPs in this file for 766 samples.

I used the following command in GenABEL to run the command:

data<-impute2databel(gen,pheno,outfile,makeprob = TRUE, old = FALSE, dataOutType = "FLOAT")

I get the following errors (and also if I remove the makeprob = TRUE, old = FALSE, dataOutType = "FLOAT" options)

Options in effect: 
     --infile    = c(21, 21, 21, 21, 21, 21, 21, 21, 21, 21)
     --outfile   = c(7, 5, 10, 6, 9, 1, 2, 8, 3, 4)
     --skiprows  = OFF
     --skipcols  = 5
     --cnrow     = ON, using data from file 'c(5,'
     --rncol     = ON, using data from file 'c(45000079,'
     --transpose = ON
     --Rmatrix   = OFF
     --nanString = c(5,
Can not open file 'c(21,' for reading
ERROR in Rstuff:failed in text2fvf_R
Index file not exists: ./tmp528216.fvi
Error in (function (cl, name, valueClass)  :

Clearly something is wrong with my impute gen file. However, I keep checking it looks normal. Has anyone had this issue before? I tried this on a much smaller file of just 10 SNPs and I also get the same error. Many thanks.

genabel impute2 impute2databel

0 answers

No answers yet.

Log in to answer this question.