Dear all, I want to determine the standard error or the mean squared error for the parameter estimate for beta and eta base on the real data. Any help on how to obtain these estimated errors.
library(survival) d <- data.frame(ob=c(149971, 70808, 133518, 145658, 175701, 50960, 126606, 82329), state=1) s <- Surv(d$ob,d$state) sr <- survreg(s~1,dist="weibull") beta<-1/sr$scale p1=(beta) p1 eta<-exp(sr$coefficients[1]) b=(eta) b Thank you Chris Guure Researcher Institute for Mathematical Research UPM ______________________________________________ R-help@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.