Crlmm Package And Genotyping With Snp 5 Chips (Affy)
Hi
I'm trying to perform some analysis on SNP5 arrays (Affymetrix) and am running problems for performing genotyping commands.
My ultimate aim is to use the genotyping information to perform a subsequent crlmmCopyNumber step in my R code.
The analysis progresses smoothly until determining gender estimations and then it shows as error message.
Please advise!!
library(crlmm)
cdfName <- "genomewidesnp5"
celPath <- "input"
outDir <- "output"
Listfiles <- list.celfiles()
ocProbesets(50e3)
batch <- as.factor(rep("A",length(Listfiles)))
plate<-substr(basename(Listfiles),1,8)
library(ff)
cnSet <- genotype(filenames = Listfiles,
cdfName = "genomewidesnp5", batch = plate)
Loading required package:genomewidesnp5Crlmm
Welcome to genomewidesnp5Crlmmversion1.0.4
Loading annotations and mixtliure model parameters.
Initializing ff objects.
Instantiating CNSet container
Preprocessing 169 files.
Cloning A and B matrices to store
genotype calls and confidence scores.
Quantile normalizing nonpolymorphicmarkers
Processing nonpolymorphic probes for 169 files.
Calling 605453 SNPs for recalibration...
Loading annotations.
Determining gender.
Error in nrow(theCenters) : subscript out of bounds
sessionInfo()
R version 2.13.0 (2011-04-13)
Platform: x86_64-redhat-linux-gnu (64-bit)
locale:
[1] LC_CTYPE=en_US.utf8 LC_NUMERIC=C
[3] LC_TIME=en_US.utf8 LC_COLLATE=en_US.utf8
[5] LC_MONETARY=C LC_MESSAGES=en_US.utf8
[7] LC_PAPER=en_US.utf8 LC_NAME=C
[9] LC_ADDRESS=C LC_TELEPHONE=C
[11] LC_MEASUREMENT=en_US.utf8 LC_IDENTIFICATION=C
attached base packages:
[1] tools stats graphics grDevices utils datasets methods
[8] base
other attached packages:
[1] genomewidesnp5Crlmm_1.0.4 ff_2.2-2
[3] bit_1.1-7 crlmm_1.10.0
[5] oligoClasses_1.14.0 Biobase_2.12.2
loaded via a namespace (and not attached):
[1] affyio_1.20.0 annotate_1.30.0 AnnotationDbi_1.14.1
[4] Biostrings_2.20.1 DBI_0.2-5 ellipse_0.3-5
[7] genefilter_1.34.0 IRanges_1.10.4 mvtnorm_0.9-9991
[10] preprocessCore_1.14.0 RSQLite_0.9-4 splines_2.13.0
[13] survival_2.36-9 xtable_1.5-6
• 3,686 views
•
link
0 answers
No answers yet.
Log in to answer this question.
Edited for better formatting. Please indent code and output with 4 spaces.
might want to try here: http://www.bioconductor.org/help/mailing-list/#bioconductor and update to the latest R (2.14) since there would be a newer crlmm 1.12.1 that might have fixed this