Hi all,
This is a new error I haven't experienced in eight months working with minfi, until last week. When attempting to load a GenomicMethylSet, either as part of an environment or as an individual .rda file, I receive the following error message:
Error in nrow(x@elementMetadata) : error in evaluating the argument 'x' in selecting a method for function 'nrow': Error: no slot of name "elementMetadata" for this object of class "GenomicMethylSet"
In the following instance, I attempted to load an .Rdata file containing a few additional data frames and a numeric matrix. I am able to load these, but then I receive the error message and the gset does not appear in my working environment. I realize saved .rda and .Rdata files can become corrupted, but it just seems odd to have this error specifically pertaining to the GenomicMethylSet. Haven't seen this issue posted elsewhere, so hoping it is just an end user error :)
Here is my session code as an example ("env_cluster_110515.RData" contains the GenomicMethylSet):
> load("~/dnaarray/env_cluster_110515.RData")
Loading required package: minfi
Loading required package: BiocGenerics
Loading required package: parallel
Attaching package: ‘BiocGenerics’
The following objects are masked from ‘package:parallel’:
clusterApply, clusterApplyLB, clusterCall, clusterEvalQ, clusterExport,
clusterMap, parApply, parCapply, parLapply, parLapplyLB, parRapply, parSapply,
parSapplyLB
The following objects are masked from ‘package:stats’:
IQR, mad, xtabs
The following objects are masked from ‘package:base’:
anyDuplicated, append, as.data.frame, as.vector, cbind, colnames, do.call,
duplicated, eval, evalq, Filter, Find, get, grep, grepl, intersect,
is.unsorted, lapply, lengths, Map, mapply, match, mget, order, paste, pmax,
pmax.int, pmin, pmin.int, Position, rank, rbind, Reduce, rownames, sapply,
setdiff, sort, table, tapply, union, unique, unlist, unsplit
Loading required package: Biobase
Welcome to Bioconductor
Vignettes contain introductory material; view with 'browseVignettes()'. To cite
Bioconductor, see 'citation("Biobase")', and for packages
'citation("pkgname")'.
Loading required package: lattice
Loading required package: GenomicRanges
Loading required package: S4Vectors
Loading required package: stats4
Loading required package: IRanges
Loading required package: GenomeInfoDb
Loading required package: SummarizedExperiment
Loading required package: Biostrings
Loading required package: XVector
Loading required package: bumphunter
Loading required package: foreach
foreach: simple, scalable parallel programming from Revolution Analytics
Use Revolution R for scalability, fault tolerance and more.
http://www.revolutionanalytics.com
Loading required package: iterators
Loading required package: locfit
locfit 1.5-9.1 2013-03-22
Setting options('download.file.method.GEOquery'='auto')
Setting options('GEOquery.inmemory.gpl'=FALSE)
Error in nrow(x@elementMetadata) :
error in evaluating the argument 'x' in selecting a method for function 'nrow': Error: no slot of name "elementMetadata" for this object of class "GenomicMethylSet"
> sessionInfo()
R version 3.2.2 (2015-08-14)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 14.04.2 LTS
locale:
[1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C LC_TIME=en_US.UTF-8
[4] LC_COLLATE=en_US.UTF-8 LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
[7] LC_PAPER=en_US.UTF-8 LC_NAME=C LC_ADDRESS=C
[10] LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
attached base packages:
[1] stats4 parallel stats graphics grDevices utils datasets methods base
other attached packages:
[1] minfi_1.16.0 bumphunter_1.10.0 locfit_1.5-9.1
[4] iterators_1.0.8 foreach_1.4.3 Biostrings_2.38.0
[7] XVector_0.10.0 SummarizedExperiment_1.0.1 GenomicRanges_1.22.1
[10] GenomeInfoDb_1.6.1 IRanges_2.4.1 S4Vectors_0.8.1
[13] lattice_0.20-33 Biobase_2.30.0 BiocGenerics_0.16.1
loaded via a namespace (and not attached):
[1] mclust_5.1 rgl_0.95.1367 base64_1.1
[4] Rcpp_0.12.1 Rsamtools_1.22.0 digest_0.6.8
[7] plyr_1.8.3 futile.options_1.0.0 ellipse_0.3-8
[10] RSQLite_1.0.0 ggplot2_1.0.1 zlibbioc_1.16.0
[13] GenomicFeatures_1.22.4 annotate_1.48.0 preprocessCore_1.32.0
[16] proto_0.3-10 splines_3.2.2 BiocParallel_1.4.0
[19] stringr_1.0.0 munsell_0.4.2 igraph_1.0.1
[22] pheatmap_1.0.7 RCurl_1.95-4.7 biomaRt_2.26.0
[25] rtracklayer_1.30.1 multtest_2.26.0 pkgmaker_0.22
[28] GEOquery_2.36.0 quadprog_1.5-5 codetools_0.2-14
[31] matrixStats_0.15.0 XML_3.98-1.3 reshape_0.8.5
[34] GenomicAlignments_1.6.1 MASS_7.3-44 bitops_1.0-6
[37] grid_3.2.2 gtable_0.1.2 nlme_3.1-122
[40] xtable_1.8-0 registry_0.3 DBI_0.3.1
[43] magrittr_1.5 scales_0.3.0 stringi_1.0-1
[46] reshape2_1.4.1 genefilter_1.52.0 doRNG_1.6
[49] RPMM_1.20 limma_3.26.1 futile.logger_1.4.1
[52] nor1mix_1.2-1 lambda.r_1.1.7 RColorBrewer_1.1-2
[55] mixOmics_5.1.2 siggenes_1.44.0 tools_3.2.2
[58] illuminaio_0.12.0 rngtools_1.2.4 survival_2.38-3
[61] colorspace_1.2-6 AnnotationDbi_1.32.0 cluster_2.0.3
[64] beanplot_1.2

Thanks for the reply, Kasper. I will add this step to my workflow.
cheers,
Sean
I have the same issue.
I created the object on the server and am now loading it in RStudio.
`load(paste0(QC_loc,"/aems450k.RGChannelSetQC.Rdata"))`
But apparently your trick won't work, as it won't be loaded at all... Any alternative, other then making the object from scratch...?
Please include the result of running
sessionInfo()to help diagnose the problem. Does it work if you load the object from command line R rather than from Rstudio?