On May 13, 2011, at 10:49 AM, viostorm wrote:
I would like to get the probabilities
.... of what?
from cph into a vector so I can put
them programatically on a graph.
Does anyone know how to get a list of the probabilities?
require(rms)
?rms
Perhaps:
?Predict
Or
?survest
(You do need to read more about the rms/Hmisc system. It is an
interdependent set of functions. For instance, you probably also need
to set up a datadist and assign it to the datadist option. There are
both plot and bplot methds for Predict objects. All this _is_ in the
help pages.)
All I can seem to
do is $coef which only contains coefficients.
----------------------------
cph(formula = srv.comb ~ collaterals + age)
Model Tests Discrimination
Indexes
Obs 555 LR chi2 7.46 R2 0.014
Events 136 d.f. 2 g 0.263
Center -0.2202 Pr(> chi2) 0.0240 gr 1.301
Score chi2 7.56
Pr(> chi2) 0.0229
Coef S.E. Wald Z Pr(>|Z|)
collaterals 0.4531 0.1721 2.63 0.0085
age -0.0067 0.0082 -0.82 0.4123
------------------------------
David Winsemius, MD
West Hartford, CT
______________________________________________
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.
