[R] Preferable contrasts?

Grathwohl,Dominik,LAUSANNE,NRC/NT dominik.grathwohl at rdls.nestle.com
Thu Nov 7 10:36:50 CET 2002


Dear all,

I'm working with Cox-regression, because data could be censored. 
But in this particular case not. 
Now I have a simple example: PRO and PRE are (0,1) coded.
The response is not normal distributed. 
We are interested in a model which could describe interaction.
But my results are depending strongly in the choose of the contrast option.
It is clear that there is some dependence in the contrasts, but in this
simple case 
I could get the vice versa effect.
My R output:

> options(contrasts = c(unordered = "contr.treatment", ordered =
"contr.poly"))
> summary(coxph(Surv(ILOG, alive) ~ factor(PRO)*factor(PRE)))
...
                             coef exp(coef) se(coef)      z     p
factor(PRO)1               0.6576     1.930    0.302  2.177 0.029
factor(PRE)1               0.0681     1.070    0.304  0.224 0.820
factor(PRO)1:factor(PRE)1 -0.7703     0.463    0.431 -1.789 0.074
...
> options(contrasts = c(unordered = "contr.SAS", ordered = "contr.poly"))
> summary(coxph(Surv(ILOG, alive) ~ factor(PRO)*factor(PRE)))
...
                            coef exp(coef) se(coef)      z     p
factor(PRO)0               0.113     1.119    0.304  0.370 0.710
factor(PRE)0               0.702     2.018    0.299  2.350 0.019
factor(PRO)0:factor(PRE)0 -0.770     0.463    0.431 -1.789 0.074
...

What would the experts recommend?

Kind regards,

Dominik

Dominik Grathwohl 
Biostatistician 
Nestlé Research Center 
PO Box 44, CH-1000 Lausanne 26 
Phone: + 41 21 785 8034 
Fax: + 41 21 785 8556 
e-mail: dominik.grathwohl at rdls.nestle.com 


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list