Weibull survival regression

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