Search
Question: Error from scater package
0
gravatar for huiyunwu
8 months ago by
huiyunwu0
huiyunwu0 wrote:

Hi there,

While playing with the example provided online, I got a question and need your help. Below is the workflow:

Source:

https://www.bioconductor.org/help/workflows/simpleSingleCell/

Code:

all.counts <- as.data.frame(read_excel("D:\\HDexamples\\GSE61533_HTSEQ_count_results.xls", sheet=1))
rownames(all.counts) <- all.counts$ID
all.counts <- as.matrix(all.counts[,-1])

library(SingleCellExperiment) # installed from local file download
sce <- SingleCellExperiment(assays=list(counts=all.counts))
is.spike <- grepl("^ERCC", rownames(sce))
isSpike(sce, "ERCC") <- is.spike
is.mito <- grepl("^mt-", rownames(sce))

library(scater)
sce <- calculateQCMetrics(sce, feature_controls=list(ERCC=is.spike, Mt=is.mito))

Error in calculateQCMetrics(sce, feature_controls = list(ERCC = is.spike,  :

  object must be an SCESet object.

 

Package "scater" was installed a few days ago.

 

Thanks,

William Wu

ADD COMMENTlink modified 8 months ago by Aaron Lun19k • written 8 months ago by huiyunwu0

Please include the output of sessionInfo() in your question (http://www.bioconductor.org/help/support/posting-guide/)

ADD REPLYlink modified 8 months ago • written 8 months ago by Peter Hickey370
0
gravatar for Aaron Lun
8 months ago by
Aaron Lun19k
Cambridge, United Kingdom
Aaron Lun19k wrote:

The workflows have been updated to Bioconductor 3.6 in preparation for the upcoming release. If you want to use the code on the workflow page, you need to upgrade to BioC-devel for the time being.

library(BiocInstaller)
useDevel()
biocLite()
ADD COMMENTlink written 8 months ago by Aaron Lun19k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 130 users visited in the last hour