如何从GSE数据集创建Seurat对象?

mccptt67  于 2023-09-27  发布在  其他
关注(0)|答案(1)|浏览(104)

我知道,我们可以创建Seurat对象从细胞游侠输出文件(条形码,功能,矩阵)。但我有一个csv格式的GSE单细胞数据集。我该如何用它来创建修拉对象呢?
谢谢.
a = read.csv(<csv file >, header = TRUE)
b = t(a)
c = CreateSeuratObject(counts = b, project = "my_single_cell", min.cells = 3, min.features = 200)
Error in CreateAssayObject(counts = counts, min.cells = min.cells, min.features = min.features, :
输入矩阵中不存在单元格名称(列名)
dput(a[1:10, 1:10])

2.0293824672699, 0, 0, 1.67916893959045, 1.25670552253723, 0, 
0, 2.35019993782043, 0), Lypla1 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 
0), Gm37988 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0), Tcea1 = c(0, 0, 
0, 0, 0, 0, 0, 0, 0, 0), Atp6v1h = c(0, 0, 0, 1.69750428199768, 
0, 1.25670552253723, 0, 0, 0, 0), Rb1cc1 = c(0, 0, 0, 0, 0, 0, 
0, 0, 0, 0), X4732440D04Rik = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0), 
    Pcmtd1 = c(0, 0, 0, 0, 0, 1.25670552253723, 0, 2.10469627380371, 
    0, 0), Gm26901 = c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0)), row.names = c("AAACCCATCTGGAAGG-WT1", 
"AAACGAACAACTCGTA-WT1", "AAACGAACACAGACGA-WT1", "AAACGAATCACGAACT-WT1", 
"AAACGAATCAGGGTAG-WT1", "AAACGCTAGCTGACTT-WT1", "AAACGCTTCGACCTAA-WT1", 
"AAAGAACAGAATCCCT-WT1", "AAAGTCCAGAGTGACC-WT1", "AAAGTCCAGTAAACAC-WT1"
), class = "data.frame")```
>
idfiyjo8

idfiyjo81#

原始表格似乎行中有细胞,列中有基因,第一列中有细胞名称。要获得正确的行名称,请尝试以下操作:

a = read.csv(data.csv, row.names = 1)

b = t(a)

c = CreateSeuratObject(counts = b, project = "my_single_cell", min.cells = 3, min.features = 200)

编辑:该表似乎包含标准化的表达式值,而不是原始计数。修拉可能不是处理这种数据的合适工具。

相关问题