全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
2321 1
2016-12-12
正在看一个例子,是用指数分布进行生存分析数据的拟合,前面都没问题
>library(survival)
>data(ovarian)
> s2=survreg(Surv(futime, fustat)~1 , ovarian, dist=’exponential’)
> summary(s2)
            Call:
             survreg(formula = Surv(futime, fustat) ~ 1, data = ovarian, dist = "exponential")
             Value Std. Error z p
             (Intercept) 7.17 0.289 24.8 3.72e-136
             Scale fixed at 1
               Exponential distribution
             Loglik(model)= -98 Loglik(intercept only)= -98
              Number of Newton-Raphson Iterations: 4
               n= 26
>plot(T,1-pexp(T,exp(-7.169)),xlab="t",ylab=expression(hat(S)*"(t)"))

就这么几行程序,我完全按照这个输的,但是图拟合出来不对,左边是例子里附的图,右边是我做出来的,因为plot里用了“T”,代表Ture,所有x和y都是1,当然不对了。

求大神指导该怎么修改这个程序才能画出左边的图
QQ图片20161212214942.png


二维码

扫码加我 拉你入群

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

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

全部回复
2016-12-13 18:07:11
求大神帮帮我这个给导师打工的苦逼学生吧
二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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