全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2173 0
2022-06-25
第一次发帖,请大大家帮我看看怎么解决这个问题,以下是我具体的代码和数据:
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 应该怎么修改呢?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群