Using Thomas Lumley's solution with the lung dataset, I obtained:
curve(pweibull(x, scale=1/coef(lung.wbs), shape=1/lung.wbs$scale,lower.tail=FALSE),from=0, to=max(lung$time))
That's because you read the part of the message where I got the reparametrization wrong, rather than the part where I got right.
It should be
curve(pweibull(x, scale=exp(coef(lung.wbs)), shape=1/lung.wbs$scale,
lower.tail=FALSE),from=0, to=max(lung$time))
-thomas
______________________________________________ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
