5 鲁棒性检验代码

您好我用的csv文件 文件格式如下 ,这个是包里自带的数据集我转成csv了

源数据是cel我该怎么办


> library(rbsurv)
>  setwd("C:/Users/Administrator/Desktop")
>  Mydata<-read.table(file = "C:\\Users\\Administrator\\Desktop\\gliomaSet.txt",header=TRUE)
> head(Mydata)[0:6,0:6]
           g1      g2      g3      g4      g5      g6
Chip1 2504.19 2229.21 3024.23 3509.83 1938.45 735.089
Chip2 3029.66 2878.98 2412.68 4291.43  902.45 771.696
Chip3 4280.25 3702.75 4053.38 5346.73 1251.64 714.212
Chip4 3763.27 2968.65 3326.00 3758.37 1630.62 792.536
Chip5 2433.61 2130.91 2466.12 3634.65 1444.51 794.392
Chip6 5229.69 2891.58 4725.19 5572.60 1926.45 562.855
> time <- gliomaSet$Time
> status <- gliomaSet$Status
> z <- cbind(gliomaSet$Age, gliomaSet$Gender)
> fit <- rbsurv(time=time, status=status, x=n, method="efron", max.n.genes=20, n.iter=10, n.seq=1)
Please wait...Error in model.frame.default(formula = Surv(time, status) ~ x[, i]) :
  变数的长度不一样('x[, i]')
> fit$model
错误: 找不到对象'fit'


请先 登录 后评论

1 个回答

microRNA

1. 你这代码中的变量n 是什么?

2. timestatus : 这些是啥?能取到值吗? 建议你先将课程中的代码跑通,看看每步的输入,输出的格式是什么样的。

请先 登录 后评论