Search
Question: Error from scater package
0
gravatar for huiyunwu
11 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 11 months ago by Aaron Lun21k • written 11 months ago by huiyunwu0

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

ADD REPLYlink modified 11 months ago • written 11 months ago by Peter Hickey380
0
gravatar for Aaron Lun
11 months ago by
Aaron Lun21k
Cambridge, United Kingdom
Aaron Lun21k 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 11 months ago by Aaron Lun21k
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: 154 users visited in the last hour