From: Charles Annis, P.E. <Charles.Annis_at_statisticalengineering.com>

Date: Tue 03 May 2005 - 13:37:46 EST

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue May 03 13:43:03 2005

Date: Tue 03 May 2005 - 13:37:46 EST

Dear R-Helpers:

I have tried everything I can think of and hope not to appear too foolish when my error is pointed out to me.

sN.lm<- lm(log(Cycles) ~ log10(stress), data=survreg.df)

vvvvvvvvvvvgaussian.survreg<- survreg(formula=Surv(time=log(Cycles), event) ~ log10(stress), dist="gaussian", data=survreg.df)

produce identical parameter estimates and differ slightly in the residual standard error and scale, which is accounted for by scale being the MLE and thus biased. Correcting by sqrt(18/16) produces agreement. Using predict() for the lm, and predict.survreg() for the survreg model and correcting for the differences in stdev, produces identical plots of the fit and the upper and lower confidence intervals. All of this is as it should be.

And,

vvvvvvlognormal.survreg<- survreg(formula=Surv(time=(Cycles), event) ~ log10(stress), dist="lognormal", data=survreg.df)

produces summary() results that are identical to the earlier call to survreg(), except for the call, of course. The parameter estimates and SE are identical. Again this is as I would expect it.

But since the call uses Cycles, rather than log(Cycles) predict.survreg() returns $fit in Cycles units, rather than logs, and of course the fits are identical when plotted on a log-log grid and also agree with lm()

Here is the fly in the ointment: The upper and lower confidence intervals, based on the $se.fit for the dist="lognormal" are quite obviously different from the other two methods, and although I have tried everything I could imagine I cannot reconcile the differences.

So I have missed something. Would some kind soul please point out my error?

Thanks.

Charles Annis, P.E.

Charles.Annis@StatisticalEngineering.com
phone: 561-352-9699

eFax: 614-455-3265

http://www.StatisticalEngineering.com

R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Tue May 03 13:43:03 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:31:32 EST
*