全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
7156 5
2011-05-29
最近在学习用glmnet,在用glmnet实现cox惩罚回归时,对照例子改了代码,但是老提示错误,错误内容是:
提交这句后
> cv.fit <- cv.glmnet(x,y, family = "cox", maxit = 1000)
会提示错误:
错误: Cox model requires a matrix with columns 'time' (>0) and 'status' (binary) as a response; a 'Surv' object suffices
但是我明明已经定义了y为y=Surv(t,SS),哪位大侠能帮帮忙!
二维码

扫码加我 拉你入群

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

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

全部回复
2011-5-30 07:31:59
没用过,帮你顶一下。
二维码

扫码加我 拉你入群

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

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

2011-5-30 08:59:42
问题太抽象了,没法帮忙,lz至少也要把code写出来呀。
二维码

扫码加我 拉你入群

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

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

2011-5-30 10:13:56
3# ltx5151 set.seed(10101)
N=1000;p=30
nzc=p/3
x=matrix(rnorm(N*p),N,p)
beta=rnorm(nzc)
fx=x[,seq(nzc)]%*%beta/3
hx=exp(fx)
ty=rexp(N,hx)
tcens=rbinom(n=N,prob=.3,size=1)# censoring indicator
y=cbind(time=ty,status=1-tcens) # y=Surv(ty,1-tcens) with library(survival)
fit=glmnet(x,y,family="cox")
plot(fit)


这是自带的例子,可以运行,一切正常。
我套用自己的数据,是模拟数据
y=Surv(t,SS)
y
x<-matrix(X,100,1000)
x
cv.fit <- cv.glmnet(x,y, family = "cox", maxit = 1000)
fit<-glmnet(x,y,family = "cox", maxit = 1000)

plot(cv.fit)
plot(fit)

运行之后出现上述错误提示!!求解啊!!
二维码

扫码加我 拉你入群

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

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

2015-4-13 20:00:17
ss要是0-1变量哦~
二维码

扫码加我 拉你入群

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

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

2015-4-13 20:01:29
ss要是0-1变量哦~
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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