楼主: baoaibaobao
6554 5

[问答] 求助,关于glmnet提示错误! [推广有奖]

  • 0关注
  • 25粉丝

副教授

38%

还不是VIP/贵宾

-

威望
0
论坛币
9669 个
通用积分
13.1440
学术水平
80 点
热心指数
128 点
信用等级
38 点
经验
12060 点
帖子
548
精华
0
在线时间
1008 小时
注册时间
2009-5-21
最后登录
2024-4-30

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

感谢您参与论坛问题回答

经管之家送您两个论坛币!

+2 论坛币
最近在学习用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),哪位大侠能帮帮忙!
二维码

扫码加我 拉你入群

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

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

关键词:glmnet NET MNE response object 求助 提示 glmnet

沙发
qoiqpwqr 发表于 2011-5-30 07:31:59 |只看作者 |坛友微信交流群
没用过,帮你顶一下。

使用道具

藤椅
ltx5151 发表于 2011-5-30 08:59:42 |只看作者 |坛友微信交流群
问题太抽象了,没法帮忙,lz至少也要把code写出来呀。

使用道具

板凳
baoaibaobao 发表于 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)

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

使用道具

报纸
wahahagu 发表于 2015-4-13 20:00:17 |只看作者 |坛友微信交流群
ss要是0-1变量哦~

使用道具

地板
wahahagu 发表于 2015-4-13 20:01:29 |只看作者 |坛友微信交流群
ss要是0-1变量哦~

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-5-14 03:40