Hello!
I am trying to fit confidence intervals for my server object as well, but I
can only get an interval estimation for the intercept coefficient and am
getting NAs for my scale parameter. Did you also have this issue?
Below is a rough outline of my code
>newsurvobj <-Surv(studytimenew, censors)
> newserv <- survfit(newsurvobj~1, type="kaplan-meier")
> survreg(formula=newsurvobj~1, dist="weibull")
Call:
survreg(formula = newsurvobj ~ 1, dist = "weibull")
Coefficients:
(Intercept)
2.726151
Scale= 0.4861742
Loglik(model)= -160.8 Loglik(intercept only)= -160.8
n= 100
> ??confint
> survival <- survreg(formula=newsurvobj~1, dist="weibull")
> confint(survival, level=.95)
2.5 % 97.5 %
(Intercept) 2.540459 2.911843
confint(survival, "coefficients$scale",level=.95)
2.5 % 97.5 %
coefficients$scale NA NA
Thank you so much for any help you can offer!
-Ann
--
View this message in context:
http://r.789695.n4.nabble.com/Confidence-Intervals-for-survreg-survival-tp4674137p4704283.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
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.
