Plot survreg and intcox
On Aug 26, 2014, at 4:17 AM, Silong Liao wrote:
Dear R users, I'm trying to plot survival probability against time(in years) using "survreg" and "intcox". Please can you help me with this problem?
That line at the bottom of the message delivered to the list subscribers saying that the HTL version was deleted explains why this is such a mess. Please read the Posting Guide and post in plain text:
David. > My code shows below:>mod.reg1=survreg(s_new~type+sex+eye+preopiop+preopva,dist="weibull")>summary(mod.reg1) > Call:survreg(formula = s_new ~ type + sex + eye + preopiop + preopva, dist = "weibull") Value Std. Error z p(Intercept) 40.539 20.582 1.970 4.89e-02typeTrab -6.606 4.279 -1.544 1.23e-01sexM -1.055 3.765 -0.280 7.79e-01eyeR -2.112 3.587 -0.589 5.56e-01preopiop -0.308 0.269 -1.147 2.52e-01preopva -0.461 1.771 -0.260 7.95e-01Log(scale) 2.058 0.285 7.222 5.12e-13 > Scale= 7.83 > Weibull distributionLoglik(model)= -78.7 Loglik(intercept only)= -81.4 Chisq= 5.37 on 5 degrees of freedom, p= 0.37 Number of Newton-Raphson Iterations: 10 n= 339 > ------------------------------------------------------------------------------------------------------------->cox.fit=intcox(s_new~type+eye+sex+age+preopiop+preopva,data=glaucoma_new)>summary(cox.fit)Call:intcox(formula = s_new ~ type + eye + sex + age + preopiop + preopva, data = glaucoma_new) > n= 339 > coef exp(coef) se(coef) z Pr(>|z|)typeTrab 0.59391 1.81106 NA NA NAeyeR 0.28419 1.32868 NA NA NAsexM -0.11597 0.89050 NA NA NAage -0.06556 0.93655 NA NA NApreopiop 0.03903 1.03980 NA NA NApreopva -0.05517 0.94632 NA NA NA > exp(coef) exp(-coef) lower .95 upper .95typeTrab 1.8111 0.5522 NA NAeyeR 1.3287 0.7526 NA NAsexM 0.8905 1.1230 NA NAage 0.9365 1.0678 NA NApreopiop 1.0398 0.9617 NA NApreopva 0.9463 1.0567 NA NA > Rsquare= NA (max possible= 0.327 )Likelihood ratio test= NA on 6 df, p=NAWald test = NA on 6 df, p=NAScore (logrank) test = NA on 6 df, p=NA > > > Kind regards,Sid > > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. David Winsemius Alameda, CA, USA