[BioC] crlmm problem

Aaerp al3n70rn at gmail.com
Sat Feb 1 11:27:23 CET 2014


Hi all,

I am trying to analyze some Illumina beadarray SNP data 
(humanomniexpress12v1b) using crlmm. Do you have any thoughts on what 
might be happening here?

Kind regards,
Agusti

library(crlmm)
library(ff)

ocProbesets(150e3)

datadir <- "/Volumes/My Passport/SNP_LGG/6182351100"
samplesheet = read.csv(file.path(datadir, 
"SampleSheet_PJ1111109_11S46v2.csv"), header=TRUE, as.is=TRUE,sep="\t")

arrayNames <- file.path(datadir, unique(samplesheet[, "SentrixPosition_A"]))

all(file.exists(paste(arrayNames, "_Grn.idat", sep="")))
#[1] TRUE
all(file.exists(paste(arrayNames, "_Red.idat", sep="")))
#[1] TRUE
cdfName <- "humanomniexpress12v1b"
batch <- rep("1", nrow(samplesheet))
arrayInfo <- list(barcode=NULL, position="SentrixPosition_A")

cnSet <- genotype.Illumina(sampleSheet=samplesheet,
arrayNames=arrayNames,
arrayInfoColNames=arrayInfo,
cdfName="humanomniexpress12v1b",
batch=batch)

Instantiate CNSet container.
path arg not set.  Assuming files are in local directory, or that 
complete path is provided
Initializing container for genotyping and copy number estimation
Loading required package: humanomniexpress12v1bCrlmm
Welcome to humanomniexpress12v1bCrlmm version 1.0.1
reading /Volumes/My Passport/SNP_LGG/6182351100/R01C01_Grn.idat
Erreur dans G$RunInfo[1, 1] : indice hors limites

sessionInfo()
R version 3.0.2 (2013-09-25)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] fr_FR.UTF-8/fr_FR.UTF-8/fr_FR.UTF-8/C/fr_FR.UTF-8/fr_FR.UTF-8

attached base packages:
[1] tools     parallel  stats     graphics  grDevices utils datasets  
methods   base

other attached packages:
[1] ff_2.2-12 bit_1.1-10                       
humanomniexpress12v1bCrlmm_1.0.1
[4] crlmm_1.20.1 preprocessCore_1.22.0            oligoClasses_1.22.0
[7] BiocGenerics_0.6.0

loaded via a namespace (and not attached):
  [1] affyio_1.28.0        Biobase_2.20.1       BiocInstaller_1.10.4 
Biostrings_2.28.0    codetools_0.2-8
  [6] ellipse_0.3-8        foreach_1.4.1        GenomicRanges_1.12.5 
grid_3.0.2           illuminaio_0.2.0
[11] IRanges_1.18.4       iterators_1.0.6      lattice_0.20-23 
Matrix_1.1-1.1       matrixStats_0.8.12
[16] mvtnorm_0.9-9996     R.methodsS3_1.5.2    RcppEigen_0.3.2.0 
stats4_3.0.2         VGAM_0.9-3
[21] zlibbioc_1.6.0



More information about the Bioconductor mailing list