用R语言进行生存分析，并假设数据符合指数分布，调用的是内部的数据包，程序如下> library(survival)
> data(ovarian)
> s1<-Surv(ovarian$futime,ovarian$fustat)    #判断完全数据和截尾数据
> s1
[1]   59   115   156   421+  431   448+  464   475   477+  563   638   744+
[13]  769+  770+  803+  855+ 1040+ 1106+ 1129+ 1206+ 1227+  268   329   353
[25]  365   377+
> fit<-survreg(s1~1,dist='exponential')  #按照指数分布拟合函数，
到这里就已经拟合完了，下面是结果，可以看到参数是7.169
> fit

      Call:
      survreg(formula = s1 ~ 1, dist = "exponential")

      Coefficients:
     (Intercept)
               7.16935

       Scale fixed at 1

       Loglik(model)= -98   Loglik(intercept only)= -98
       n= 26

但是最后的拟合图不会画了，想知道怎么才能画出如下所示的这样的图（这是我截取的例子里面的图）