R语言新新新手，用指数分布进行生存分析拟合后，图形怎么不对呢

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

正在看一个例子，是用指数分布进行生存分析数据的拟合，前面都没问题
>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，当然不对了。

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

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：生存分析 指数分布 R语言 distribution Exponential survival library Error

求大神帮帮我这个给导师打工的苦逼学生吧