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