Problem with newSCESet(countData=all.counts) function
1
1
Entering edit mode
5.7 years ago
bk11 ★ 3.0k

Hi, I m running through this for single cell RNAseq data analysis. And I got this error after running a few codes. If some one could figure out that what might have I done wrong, I will appreciate. Thanks

library(R.utils)
gunzip("GSE61533_HTSEQ_count_results.xls.gz", remove=FALSE, overwrite=TRUE)
library(gdata)
all.counts <- read.xls(’GSE61533_HTSEQ_count_results.xls’, sheet=1, header=TRUE, row.names=1)
library(scater)

if (!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")
BiocManager::install("SingleCellExperiment", version = "3.8")

library(SingleCellExperiment)
sce <- newSCESet(countData=all.counts)

Error in newSCESet(countData = all.counts) : 
  could not find function "newSCESet"
single cell experiment single cell RNA seq • 3.9k views
ADD COMMENT
0
Entering edit mode

Very strange! Check again in a fresh R-session, especially double check that scater library is loaded

library(scater)

ADD REPLY
0
Entering edit mode

Thank you Santosh. library(scater) was loaded and still could not find function "newSECSet"

ADD REPLY
0
Entering edit mode
5.7 years ago

Scater has switched to SingleCellExperiment object and SCESet has been removed from new scater releases.

https://support.bioconductor.org/p/101956/

See the latest doc on how SingleCellExperiment works

https://bioconductor.org/packages/release/bioc/vignettes/scater/inst/doc/vignette-intro.html

See also the workflow for single cell analysis

http://master.bioconductor.org/packages/release/workflows/html/simpleSingleCell.html

ADD COMMENT

Login before adding your answer.

Traffic: 855 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6