第一次发帖,请大大家帮我看看怎么解决这个问题,以下是我具体的代码和数据:
library(Seurat)
library(ggplot2)
library(psych)
data = read.csv('./各个表达矩阵/第1个维数是_59_的表达矩阵.csv',header = T)
rownames(data) = data[,1]
data = data[,-1]
colnames(data) = c(1:dim(data)[2])
cui <- CreateSeuratObject(counts = data)
table(
cui@meta.data$orig.ident)
cui <- subset(cui, subset = nFeature_RNA > 0)
dim(cui)
cui <- NormalizeData(cui, normalization.method = "LogNormalize", scale.factor = 10000)
cui <- FindVariableFeatures(cui, selection.method = "vst", nfeatures = 50)
all.genes <- rownames(cui)
cui <- ScaleData(cui, features = all.genes)
cui <- RunPCA(cui,npcs = 10,features = VariableFeatures(object = cui),verbose = FALSE )
出现的错误:Error in irlba(A = t(x = object), nv = npcs, ...) :
max(nu, nv) must be positive 应该怎么修改呢?