[R] coxph | anova results
Daniel Malter
daniel at umd.edu
Thu Oct 18 08:41:37 CEST 2007
Hi,
I ran a coxph model and in the summary the independent variables are
reported highly significant. However, if I run anova() on that model, some
of the variables are reported to explain basically no deviance and other
variables which are reported insignificant in the model summary are reported
to explain some deviance here. Can somebody give advice how to make sense of
this output and interpret the results if the two outputs are somewhat
contradictory? Is it due to anova() including the variables sequentially? Do
I rely on the model or on the anova for interpretation?
Thanks much,
Daniel
Call:
coxph(formula = Surv(DurBeg, Closed) ~ SumInvested + numberofgainstocksheld
+
numberoflossstocksheld + TotalInvested + WorseThanMarket +
Loss + posilow + posihigh + negalow + negahigh + Loss:posilow +
Loss:posihigh + Loss:negalow + Loss:negahigh)
n=11232 (28 observations deleted due to missingness)
coef exp(coef) se(coef) z p
SumInvested -2.24e-05 1.000 1.12e-05 -1.999 4.6e-02
numberofgainstocksheld -5.37e-02 0.948 1.20e-02 -4.493 7.0e-06
numberoflossstocksheld -1.72e-01 0.842 1.22e-02 -14.032 0.0e+00
TotalInvested 3.18e-05 1.000 1.65e-05 1.919 5.5e-02
WorseThanMarket1 -9.83e-02 0.906 4.60e-02 -2.139 3.2e-02
Loss1 3.45e-01 1.412 9.67e-02 3.570 3.6e-04
posilow -6.72e-02 0.935 3.73e-02 -1.803 7.1e-02
posihigh 6.93e-02 1.072 3.62e-02 1.914 5.6e-02
negalow 1.87e-01 1.206 4.28e-02 4.367 1.3e-05
negahigh -2.42e-01 0.785 4.78e-02 -5.060 4.2e-07
Loss1:posilow 9.06e-03 1.009 5.25e-02 0.172 8.6e-01
Loss1:posihigh -1.04e-01 0.901 5.07e-02 -2.054 4.0e-02
Loss1:negalow -2.98e-01 0.742 6.10e-02 -4.891 1.0e-06
Loss1:negahigh 2.66e-01 1.305 6.70e-02 3.971 7.2e-05
exp(coef) exp(-coef) lower .95 upper .95
SumInvested 1.000 1.000 1.000 1.000
numberofgainstocksheld 0.948 1.055 0.926 0.970
numberoflossstocksheld 0.842 1.187 0.822 0.863
TotalInvested 1.000 1.000 1.000 1.000
WorseThanMarket1 0.906 1.103 0.828 0.992
Loss1 1.412 0.708 1.168 1.707
posilow 0.935 1.069 0.869 1.006
posihigh 1.072 0.933 0.998 1.151
negalow 1.206 0.829 1.109 1.311
negahigh 0.785 1.274 0.715 0.862
Loss1:posilow 1.009 0.991 0.910 1.118
Loss1:posihigh 0.901 1.110 0.816 0.995
Loss1:negalow 0.742 1.348 0.658 0.836
Loss1:negahigh 1.305 0.766 1.144 1.488
Rsquare= 0.035 (max possible= 0.967 )
Likelihood ratio test= 399 on 14 df, p=0
Wald test = 400 on 14 df, p=0
Score (logrank) test = 412 on 14 df, p=0
> anova(cluster4,test="Chisq")
Analysis of Deviance Table
Cox model: response is Surv(DurBeg, Closed)
Terms added sequentially (first to last)
Df Deviance Resid. Df Resid. Dev P(>|Chi|)
NULL 11232 38469
SumInvested 1 0 11231 38526 1
numberofgainstocksheld 1 55 11230 38471 1.032e-13
numberoflossstocksheld 1 271 11229 38200 7.454e-61
TotalInvested 1 4 11228 38196 3.964e-02
WorseThanMarket 1 7 11227 38189 8.439e-03
Loss 1 13 11226 38176 4.021e-04
posilow 1 66 11225 38110 4.306e-16
posihigh 1 1.502e-01 11224 38110 1
negalow 1 3.036e-01 11223 38110 1
negahigh 1 10 11222 38100 1.321e-03
Loss:posilow 1 1 11221 38099 3.385e-01
Loss:posihigh 1 2 11220 38096 1.208e-01
Loss:negalow 1 11 11219 38085 9.707e-04
Loss:negahigh 1 16 11218 38070 7.283e-05
More information about the R-help
mailing list