Recall genotypes from Illumina idat files using crlmm package
0
1
Entering edit mode
nadne ▴ 10
@nadne-9806
Last seen 7.8 years ago

Hello,

I want to get the genotype recall from idat files. I have few thousands samples, and I'm not sure how to do that efficiently. I used the next set of commands:

library(crlmm)
library(ff)
library(snow)
library(doSNOW)
cl <- makeCluster(90, type="SOCK")
registerDoSNOW(cl)
options("krlmm.cores" = 90)
ocProbesets(150e3)
ocSamples(5000)
dir_name = ‘set1’ # Directory contains 500 samples
g = genotype.Illumina(cdfName='humanomni25quadv1b', path=dir_name, call.method='krlmm', verbose=TRUE)

Can someone please help me in fine tuning the different parameters: 

ocProbesets, ocSamples, krlmm.cores, snow, ff

Are all needed or are they masking each-other?

Can someone give me rough estimation of runtime? Should I expect hours, days or weeks?

I have access to different machine with up to 100 cpus and up to 1T RAM, how can I take advantage of such?

crlmm genotype • 1.4k views
ADD COMMENT

Login before adding your answer.

Traffic: 1111 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6