[R] help with probabilites in cph

[David Winsemius]

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