[R] weibull in R

[Ronaldo Reis Jr.]

Hi all
I try to make a weibull survival analysis on R.
I know make this on GLIM, and now I try to make the GLIM exercice GLEX8 on R 
to learning and compare the test.
The variables are:
time censor group bodymass

In GLIM I make:
$calc %s=1 $ to fit weibull rather than exponential
$input %pcl weibull $
$macro model group*bodymass $endmac$
$use weibull t w %s $

Then, GLIM estimate an alpha for the model group*bodymass, the value on GLIM 
is 0.736.

I try to make this on R.

> survreg(Surv(time,censor)~group*bodymass,dist="weibull",scale=0)

The result is an alpha (Scale on R) different.

Scale= 1.36 

Where is the error????

Exist any reference to learning survival in R?

Thanks
ROnaldo

 
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._