Can someone please suggest me a way to remove unwanted probes (XY probes, SNP associated probes, etc.) from my 450k dataset prior to normalization. I don't want them to screw up the downstream data analysis.
I tried removing these probes from rgSet object of minfi but it didn't help. Moreover, after this I was not able to convert it to grset object. Following is the error message:
> RGsetEx <- read.450k.exp(targets = targets, extended = TRUE)
> dim(RGsetEx)
Features Samples
622399 10
> detP <- detectionP(RGsetEx)
Loading required package: IlluminaHumanMethylation450kmanifest
> keep <- rowSums(detP < 0.01) == ncol(RGsetEx)
> RGsetEx <- RGsetEx[keep,]
> dim(RGsetEx)
Features Samples
619508 10
> grset <- preprocessFunnorm(RGsetEx, nPCs=8, sex = NULL, bgCorr = TRUE, dyeCorr = TRUE, verbose = TRUE)
[preprocessFunnorm] Background and dye bias correction with noob
Error in getGreen(object)[IRed$AddressA, ] : subscript out of bounds
> sessionInfo()
R version 3.2.3 (2015-12-10)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)
locale:
[1] LC_COLLATE=English_India.1252 LC_CTYPE=English_India.1252 LC_MONETARY=English_India.1252
[4] LC_NUMERIC=C LC_TIME=English_India.1252
attached base packages:
[1] stats4 parallel stats graphics grDevices utils datasets methods base
other attached packages:
[1] IlluminaHumanMethylation450kmanifest_0.4.0 missMethyl_1.4.0
[3] RSQLite_1.0.0 DBI_0.3.1
[5] ENmix_1.4.1 doParallel_1.0.10
[7] minfi_1.16.1 bumphunter_1.10.0
[9] locfit_1.5-9.1 iterators_1.0.8
[11] foreach_1.4.3 Biostrings_2.38.4
[13] XVector_0.10.0 SummarizedExperiment_1.0.2
[15] GenomicRanges_1.22.4 GenomeInfoDb_1.6.3
[17] IRanges_2.4.8 S4Vectors_0.8.11
[19] lattice_0.20-33 Biobase_2.30.0
[21] BiocGenerics_0.16.1
loaded via a namespace (and not attached):
[1] nor1mix_1.2-1
[2] splines_3.2.3
[3] ellipse_0.3-8
[4] statmod_1.4.24
[5] doRNG_1.6
[6] Rsamtools_1.22.0
[7] methylumi_2.16.0
[8] impute_1.44.0
[9] limma_3.26.8
[10] quadprog_1.5-5
[11] digest_0.6.9
[12] RColorBrewer_1.1-2
[13] colorspace_1.2-6
[14] preprocessCore_1.32.0
[15] Matrix_1.2-4
[16] plyr_1.8.3
[17] GEOquery_2.36.0
[18] siggenes_1.44.0
[19] XML_3.98-1.4
[20] mixOmics_5.2.0
[21] biomaRt_2.26.1
[22] genefilter_1.52.1
[23] zlibbioc_1.16.0
[24] xtable_1.8-2
[25] corpcor_1.6.8
[26] scales_0.4.0
[27] BiocParallel_1.4.3
[28] annotate_1.48.0
[29] beanplot_1.2
[30] pkgmaker_0.22
[31] mgcv_1.8-12
[32] ggplot2_2.1.0
[33] GenomicFeatures_1.22.13
[34] survival_2.38-3
[35] magrittr_1.5
[36] mclust_5.1
[37] nlme_3.1-125
[38] MASS_7.3-45
[39] tools_3.2.3
[40] registry_0.3
[41] org.Hs.eg.db_3.2.3
[42] matrixStats_0.50.1
[43] stringr_1.0.0
[44] munsell_0.4.3
[45] rngtools_1.2.4
[46] AnnotationDbi_1.32.3
[47] lambda.r_1.1.7
[48] base64_1.1
[49] futile.logger_1.4.1
[50] grid_3.2.3
[51] RCurl_1.95-4.8
[52] igraph_1.0.1
[53] bitops_1.0-6
[54] gtable_0.2.0
[55] codetools_0.2-14
[56] multtest_2.26.0
[57] reshape_0.8.5
[58] IlluminaHumanMethylation450kanno.ilmn12.hg19_0.2.1
[59] ruv_0.9.6
[60] illuminaio_0.12.0
[61] GenomicAlignments_1.6.3
[62] rtracklayer_1.30.2
[63] wateRmelon_1.10.0
[64] futile.options_1.0.0
[65] stringi_1.0-1
[66] sva_3.18.0
[67] Rcpp_0.12.3
[68] geneplotter_1.48.0
[69] rgl_0.95.1441
Can some one please suggest me how to remove unwanted probes before normalization.
If you want help, you need to be very explicit. Nobody but you knows what you mean by
I tried removing these probes from rgSet object of minfi but it didn't help. Moreover, after this I was not able to convert it to mset object.
Instead of saying what you did, it's better if you show a very limited amount of code that isn't doing what you expect. In addition you should indicate what type of object you are dealing with (using the
class
function) and also show what versions of R/BioC you are using (by showing the results of runningsessionInfo()
after you have run all your code).Hi James,
I have included the script and the error code here.