genes <- c('MT.CO1', 'MT.CO2', 'MT.CO3', 'MT.ND1', 'MT.ND2', 'MT.ND3', 'MT.ND4', 'MT.ND5', 'MT.ND6', 'MT.ND4L', 'MT.ATP6', 'MT.ATP8', 'MT.CYB')
genelevels <- c('MT.CO1_level', 'MT.CO2_level', 'MT.CO3_level', 'MT.ND1_level', 'MT.ND2_level', 'MT.ND3_level', 'MT.ND4_level', 'MT.ND5_level', 'MT.ND6_level', 'MT.ND4L_level', 'MT.ATP6_level', 'MT.ATP8_level', 'MT.CYB_level')
#画图
PAADsurv$SurvObj <- with(PAADsurv,Surv(time = time,event = status))
for (i in 1:length(genelevels)){
p <- survfit(Surv(time, status)~get(genelevels[[i]]), data = PAADsurv)
ggsurvplot(p,
palette = "lancet",
size = 1,
conf.int = TRUE
)
pdf(file=paste0("PAADmtsurvival.",genes[[i]],".pdf"),onefile = FALSE,
width = 5.5,
height =5)
print(p)
dev.off()
}
遇到了问题
Error in (function (..., row.names = NULL, check.rows = FALSE, check.names = TRUE, :
参数值意味着不同的行数: 139, 13, 1
求助怎么解决