[R] bootstrapped cox regression (rms package)

Eric Claus clausericd at gmail.com
Thu Nov 29 23:02:48 CET 2012


Hi Frank,
Below is the actual output from the anova(out) command.  I had copied in the
p-values and from the previous output from anova(out) and the confidence
intervals from print(quantile(out$boot.Coef[,i], c(.025, .975))) to
illustrate that the confidence intervals were similar to SPSS while the
p-values were not.

Actual output from anova.rms(out):

 Wald Statistics          Response: Surv(months, recidivate) 

 Factor   Chi-Square d.f. P     
fac1  0.27       1    0.6055
fac2  0.20       1    0.6514
fac3  0.01       1    0.9338
fac4  0.05       1    0.8311
fac5  1.06       1    0.3036
fac6  0.33       1    0.5647
fac7  0.81       1    0.3670
fac8  0.30       1    0.5832
 TOTAL   1.48       8    0.9930

Regarding your second question, it looks like SPSS is using the original
estimate of Cox beta coefficients in the test (i.e. a new point estimate is
not generated for the statistical test)

Thanks again,
Eric




--
View this message in context: http://r.789695.n4.nabble.com/bootstrapped-cox-regression-rms-package-tp4651306p4651363.html
Sent from the R help mailing list archive at Nabble.com.




More information about the R-help mailing list