exp <- TA_cell_Young@assays[["RNA"]]@counts
fdata <- data.frame(gene_short_name = row.names(TA_cell_Young), row.names = row.names(TA_cell_Young))
pdata <- TA_cell_Young@meta.data
fd <- new("AnnotatedDataFrame", data = fdata)
pd <- new("AnnotatedDataFrame", data = pdata)
CDS <- newCellDataSet(cellData = exp,phenoData = pd,featureData = fd)
报错:
> CDS <- newCellDataSet(cellData = exp,phenoData = pd,featureData = fd)
Error in validObject(.Object) :
invalid class “CellDataSet” object: 1: feature numbers differ between assayData and featureData
invalid class “CellDataSet” object: 2: featureNames differ between assayData and featureData
报错解决:
identical(rownames(fd),rownames(exp))
exp<- exp[rownames(fd), ]