Dear Forum, I have a question about interaction estimate in the Cox model: why the hazard ratio of the interaction is not produced in the summary of the model? (Instead, the estimate of the coefficient is given in the print of the model.)
# Example: modINT <-cph( Surv(T_BASE, T_FIN,STATUS)~ NYHA + ASINI + RFP + FE_REC + XX_PR*XX_DISF) print(modINT) coef se(coef) z p NYHA=2 1.254 0.584 2.15 0.031767 ASINI 0.665 0.409 1.62 0.104247 RFP=2 0.725 0.704 1.03 0.302578 FE_REC=2 -1.637 0.810 -2.02 0.043331 XX_PR 2.189 0.649 3.37 0.000748 XX_DISF 3.233 1.000 3.23 0.001222 XX_PR * XX_DISF -2.852 1.280 -2.23 0.025852 summary(modINT) Effects Response : Surv(T_BASE, T_FIN, STATUS) Factor Low High Diff. Effect S.E. Lower 0.95 Upper 0.95 ASINI 2.0725 2.85 0.7775 0.52 0.32 -0.11 1.14 Hazard Ratio 2.0725 2.85 0.7775 1.68 NA 0.90 3.13 XX_PR 0.0000 1.00 1.0000 2.19 0.65 0.92 3.46 Hazard Ratio 0.0000 1.00 1.0000 8.92 NA 2.50 31.86 XX_DISF 0.0000 1.00 1.0000 3.23 1.00 1.27 5.19 Hazard Ratio 0.0000 1.00 1.0000 25.35 NA 3.57 179.88 NYHA - 2:1 1.0000 2.00 NA 1.25 0.58 0.11 2.40 Hazard Ratio 1.0000 2.00 NA 3.50 NA 1.12 11.00 RFP - 2:1 1.0000 2.00 NA 0.73 0.70 -0.65 2.10 Hazard Ratio 1.0000 2.00 NA 2.07 NA 0.52 8.20 FE_REC - 2:1 1.0000 2.00 NA -1.64 0.81 -3.23 -0.05 Hazard Ratio 1.0000 2.00 NA 0.19 NA 0.04 0.95 Adjusted to: XX_PR=0 XX_DISF=0 ____________________________________________________________________________________ Be a better friend, newshound, and [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.