Normalized values in oligo::crlmm and oligo::snprma don't match?
Entering edit mode
Denise Mak ▴ 10
Last seen 7.1 years ago
Dear List, I have two questions regarding oligo::crlmm for normalizing/genotyping Affy SNP 6.0 arrays. 1) I want to extract normalized values and I'm using the readSummaries() function after calling crlmm(). I found that the results are quite different from oligo::snprma. I've pasted the first bit of results from my calls to crlmm and snprma below (for extracting normalized A allele info). Maybe I'm using the wrong function from crlmm? But I also can't tell if the normalization function inside crlmm uses the normalizeToHapmap = TRUE setting (I also tested the false setting as well). 2) Ideally, what I would like is to extract normalized SNP and CN probe info. As you can see below, I'm only extracting the 906,600 SNP probes. Is there a way to use oligo functions to retrieve CN probe normalized info? Unfortunately, I cannot upgrade my R version or oligo packages because we're in the final stages of our analysis (which uses genotype info from oligo::crlmm). Thank you kindly, Denise I have checked to make sure that the sample order is the same. > crlmm(files, outputDir, verbose = TRUE); > alleleA <- readSummaries("alleleA", outputDir); > alleleA <- alleleA[order(rownames(alleleA)),]; > str(alleleA) num [1:906600, 1:76] 11.19 10.03 9.99 11.57 10.06 ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:906600] "SNP_A-1780270" "SNP_A-1780271" "SNP_A-1780272" "SNP_A-1780274" ... > snprma.results <- snprma(files); > str(A(snprma.results)) num [1:906600, 1:76] 6.29 5.98 6.62 6.19 6.13 ... - attr(*, "dimnames")=List of 2 ..$ : chr [1:906600] "SNP_A-1780270" "SNP_A-1780271" "SNP_A-1780272" "SNP_A-1780274" ... sessionInfo() R version 2.12.1 (2010-12-16) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] pd.genomewidesnp.6_1.1.0 RSQLite_0.9-4 [3] DBI_0.2-5 genomewidesnp6Crlmm_1.0.2 [5] oligo_1.14.0 oligoClasses_1.12.1 [7] Biobase_2.12.1 getopt_1.17 loaded via a namespace (and not attached): [1] AnnotationDbi_1.12.0 Biostrings_2.20.1 IRanges_1.10.4 [4] affxparser_1.24.0 affyio_1.20.0 annotate_1.28.0 [7] bit_1.1-6 crlmm_1.8.7 ellipse_0.3-5 [10] ff_2.2-1 genefilter_1.30.0 mvtnorm_0.9-96 [13] preprocessCore_1.12.0 splines_2.12.1 survival_2.36-2 [16] tools_2.12.1 xtable_1.5-6 [[alternative HTML version deleted]]
SNP Normalization oligo crlmm SNP Normalization oligo crlmm • 636 views

Login before adding your answer.

Traffic: 251 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6