伙计们。
我刚开始学习如何分析 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对象吗?