Re: [R] Scale parameter in Weibull distribution

2012-03-07 Thread David Winsemius
On Mar 6, 2012, at 6:19 PM, FU-WEN LIANG wrote: Thanks for your advise, David. I did read the help for survreg and using the followings to calculate. survreg's scale =1/(rweibull shape) survreg's intercept = log(rweibull scale) However, the scale in rweibull has been transformed by exp(

Re: [R] Scale parameter in Weibull distribution

2012-03-06 Thread FU-WEN LIANG
Thanks for your advise, David. I did read the help for survreg and using the followings to calculate. survreg's scale =1/(rweibull shape) survreg's intercept = log(rweibull scale) However, the scale in rweibull has been transformed by exp(betaX's). In my case, the baseline hazard for T is 0

Re: [R] Scale parameter in Weibull distribution

2012-03-06 Thread David Winsemius
On Mar 6, 2012, at 5:53 PM, FU-WEN LIANG wrote: Hi all, I'm trying to generate a Weibull distribution including four covariates in the model. Here is the code I used: T = rweibull(200, shape=1.3, scale=0.004*exp(-(-2.5*b1+2.5*b2+0.9*x1-1.3*x2)/1.3)) C = rweibull(n, shape=1.5, scale=0.008)

[R] Scale parameter in Weibull distribution

2012-03-06 Thread FU-WEN LIANG
Hi all, I'm trying to generate a Weibull distribution including four covariates in the model. Here is the code I used: T = rweibull(200, shape=1.3, scale=0.004*exp(-(-2.5*b1+2.5*b2+0.9*x1-1.3*x2)/1.3)) C = rweibull(n, shape=1.5, scale=0.008) #censoring time time = pmin(T,C) #observed time is m