0

伙计们。

我刚开始学习如何分析 RNAseq 数据。

在最开始的步骤中,我正在努力解决一个非常基本的问题。我试图在 R 上制作一个 Seurat 对象,但它一直失败。

HC2.data <- Read10X(data.dir = "/My/Path/scRNAseq/multi_HC/HC2/sample_feature_bc_matrix/")

10X 数据包含多个类型,并作为包含每种类型矩阵的列表返回。

HC2 <- CreateSeuratObject(counts = HC2.data, project = "HC2", min.cells = 2, min.features = 200)

CreateAssayObject(counts = counts, min.cells = min.cells, min.features = min.features) 中的错误:输入矩阵中不存在单元名称 (colnames) 名称

我认为我们必须在我的数据框中提供列名。我有如下两列。

名称(HC2.data)

[1]“基因表达”
[2]“多路捕获”

但是,没有如下列的名称。

列名(HC2.data)

无效的

接下来,我尝试根据网站(https://stat.ethz.ch/R-manual/R-devel/library/base/html/colnames.html)编写此命令

colnames(HC2.data,do.NULL = FALSE)

[1] “col1”

我也尝试使用这个命令,正如与我有同样问题的人之前推荐的( https://github.com/satijalab/seurat/issues/1650 )。

colnames(counts) <- paste0(colnames(counts), 1:ncol(counts))

is.data.frame(x) 中的错误:找不到对象“计数”

这是我所做的。我还不熟悉R。你们能帮我创建Seurat对象吗?

4

0 回答 0