生物导体单细胞RNA seq错误，assayData()函数不起作用，为什么？

Question

我遵循了“单细胞RNA测序的生物导体工作流程:规范化、降维、聚类和谱系推断”的代码。但是，在预处理步骤中，函数assayData()出现了错误。以下是论文的链接：https://www.bioconductor.org/help/course-materials/2017/BioC2017/Day2/Workshops/singleCell/doc/workshop.html#introduction

我按照以下方式对代码进行了研究：

# Bioconductor
library(BiocParallel)
library(clusterExperiment)
library(scone)
library(zinbwave)

# GitHub
library(slingshot)

# CRAN
library(doParallel)
library(gam)
library(RColorBrewer)

set.seed(20)

##Parallel computing
register(SerialParam())
NCORES <- 2
mysystem = Sys.info()[["sysname"]]
if (mysystem == "Darwin"){
  registerDoParallel(NCORES)
  register(DoparParam())
}else if (mysystem == "Linux"){
  register(bpstart(MulticoreParam(workers=NCORES)))
}else{
  print("Please change this to allow parallel computing on your computer.")
  register(SerialParam())
}

## Pre-processing
data_dir <- "C:/Users/kuosh/Documents/"
if (!dir.exists(data_dir)) system(sprintf('mkdir %s', data_dir))

urls = c("https://www.ncbi.nlm.nih.gov/geo/download/?acc=GSE95601&format=file&file=GSE95601%5FoeHBCdiff%5FCufflinks%5FeSet%2ERda%2Egz",
         "https://raw.githubusercontent.com/rufletch/p63-HBC-diff/master/ref/oeHBCdiff_clusterLabels.txt")

if(!file.exists(paste0(data_dir, "GSE95601_oeHBCdiff_Cufflinks_eSet.Rda"))) {
  if (!dir.exists(data_dir)) system(sprintf('mkdir %s', data_dir))
  download.file(urls[1], paste0(data_dir, "GSE95601_oeHBCdiff_Cufflinks_eSet.Rda.gz"))
  R.utils::gunzip(paste0(data_dir, "GSE95601_oeHBCdiff_Cufflinks_eSet.Rda.gz"))
  assayData(Cufflinks_eSet)$exprs = NULL
  assayData(Cufflinks_eSet)$fpkm_table = NULL
  assayData(Cufflinks_eSet)$tpm_table = NULL
  save(Cufflinks_eSet, file='data/GSE95601_oeHBCdiff_Cufflinks_eSet_reduced.Rda')
}

该文件已成功下载和解压缩，但我遇到了一个错误，如：Error in assayData(CufflinkseSet)$exprs = NULL : object 'CufflinkseSet‘not。有人能帮我吗？非常感谢。

Answer 1

如果我理解得很好，您已经将数据存储在对象"Cufflinks_eSet“中。但是您正在从"CufflinkseSet“中获取数据。所以这应该是漏掉下划线"_“的错误。请尝试：

assayData(Cufflinks_eSet)$exprs

而不是。