Entering edit mode
Hi, Can anyone help me resolve the following, thanks. I'm following The DMRcate Package User's Guide and get the following error -
library(DMRcate)
> data(dmrcatedata)
> myMs<-logit2(myBetas)
> nrow(illuminaSNPs)
[1] 153113
> nrow(myMs)
[1] 10042
> myMs.noSNPs<-rmSNPandCH(myMs,dist=2,mafcut=0.05)
> nrow(myMs.noSNPs)
[1] 9403
> patient<-factor(sub("-.*","",colnames(myMs)))
> type<-factor(sub(".*-","",colnames(myMs)))
> design<-model.matrix(~patient+type)
> myannotation<-cpg.annotate(myMs.noSNPs,analysis.type="differential",design=design,coeff=39)
Error: datatype must be one of 'array' or 'sequencing'
Error in cpg.annotate(myMs.noSNPs, analysis.type = "differential", design = design, :
object 'annotated' not found
In addition: Warning messages:
1: In if (datatype == "array") { :
the condition has length > 1 and only the first element will be used
2: In if (datatype == "sequencing") { :
the condition has length > 1 and only the first element will be used
sessionInfo()
R version 3.2.3 (2015-12-10)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.11.2 (El Capitan)
locale:
[1] en_AU.UTF-8/en_AU.UTF-8/en_AU.UTF-8/C/en_AU.UTF-8/en_AU.UTF-8
attached base packages:
[1] splines stats4 parallel stats graphics grDevices utils datasets methods base
other attached packages:
[1] DMRcate_1.6.1 DMRcatedata_1.6.1 DSS_2.10.0 bsseq_1.6.0 minfi_1.16.0 bumphunter_1.10.0 locfit_1.5-9.1
[8] iterators_1.0.8 foreach_1.4.3 Biostrings_2.38.2 XVector_0.10.0 SummarizedExperiment_1.0.1 GenomicRanges_1.22.2 GenomeInfoDb_1.6.1
[15] IRanges_2.4.6 S4Vectors_0.8.5 lattice_0.20-33 Biobase_2.30.0 BiocGenerics_0.16.1 BiocInstaller_1.20.1
loaded via a namespace (and not attached): [1] nlme_3.1-122 bitops_1.0-6 matrixStats_0.50.1 RColorBrewer_1.1-2 tools_3.2.3 doRNG_1.6 nor1mix_1.2-1 rpart_4.1-10 [9] Hmisc_3.17-1 DBI_0.3.1 Gviz_1.14.0 colorspace_1.2-6 nnet_7.3-11 gridExtra_2.0.0 base64_1.1 preprocessCore_1.32.0 [17] chron_2.3-47 pkgmaker_0.22 rtracklayer_1.30.1 scales_0.3.0 genefilter_1.52.0 quadprog_1.5-5 stringr_1.0.0 digest_0.6.8 [25] Rsamtools_1.22.0 foreign_0.8-66 illuminaio_0.12.0 siggenes_1.44.0 R.utils_2.2.0 GEOquery_2.36.0 dichromat_2.0-0 BSgenome_1.38.0 [33] limma_3.26.3 RSQLite_1.0.0 mclust_5.1 BiocParallel_1.4.3 gtools_3.5.0 acepack_1.3-3.3 R.oo_1.19.0 VariantAnnotation_1.16.4 [41] RCurl_1.95-4.7 magrittr_1.5 Formula_1.2-1 futile.logger_1.4.1 Rcpp_0.12.2 munsell_0.4.2 proto_0.3-10 R.methodsS3_1.7.0 [49] stringi_1.0-1 MASS_7.3-45 zlibbioc_1.16.0 plyr_1.8.3 grid_3.2.3 multtest_2.26.0 GenomicFeatures_1.22.7 annotate_1.48.0 [57] beanplot_1.2 igraph_1.0.1 rngtools_1.2.4 corpcor_1.6.8 codetools_0.2-14 biomaRt_2.26.1 mixOmics_5.2.0 futile.options_1.0.0 [65] XML_3.98-1.3 biovizBase_1.18.0 latticeExtra_0.6-26 lambda.r_1.1.7 data.table_1.9.6 gtable_0.1.2 reshape_0.8.5 ggplot2_2.0.0 [73] xtable_1.8-0 survival_2.38-3 GenomicAlignments_1.6.1 AnnotationDbi_1.32.2 registry_0.3 ellipse_0.3-8 cluster_2.0.3 rgl_0.95.1429

Excellent, thanks Tim.