Question: Importing huge data into GenABEL
gravatar for Paula Sanchez
5.3 years ago by
Paula Sanchez0 wrote:

I am trying to import a huge data file into GenABEL, but I am having trouble. I have 15,000 bulls (30 Chromosomes) genotyped with 700,000 SNPs.

My genotype files is coded as:


 Name Chr   Pos 238570 238632 238757 238591

BT_17218   1 17218     AB     AB     AB     AB

BT_18426   1 18426     AA     --     AA     AB

BT_20658   1 20658     AB     AA     AB     BB

BT_28296   1 28296     BB     AB     BB     --

BT_31152   1 31152     AB     AB     BB     AB

When I try to convert into GenABEL object I get the following error:

data <- = "pheno.txt", gen = "geno.raw", force=T)
ids loaded...
marker names loaded...
chromosome data loaded...
map data loaded...
allele coding data loaded...
strand data loaded...
Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings,  :
  too many items

What would be the best option for me to convert this data? Thanks in advance.

gwas snp R genome • 2.3k views
ADD COMMENTlink modified 4.9 years ago by Biostar ♦♦ 20 • written 5.3 years ago by Paula Sanchez0

Can you try splitting the data by chromosome, read it in, and merge the resulting gwaa-data objects?

ADD REPLYlink written 5.3 years ago by komal.rathi3.7k

This may be linked to hitting R's maximum number of rows. See here and here.

ADD REPLYlink modified 11 months ago by _r_am31k • written 4.9 years ago by Jean-Karim Heriche23k

What version of R? What version of GenABEL?

ADD REPLYlink modified 11 months ago by _r_am31k • written 4.9 years ago by Sean Davis26k

Version: R/R-3.2.2 GenABEL: 1.8.0 Did someone find an answer for this? Total data size is 4000 samples and 6 Millions variants !!

ADD REPLYlink written 3.7 years ago by always_learning1.1k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 974 users visited in the last hour