[R] How to extract se(coef) from cph?

Frank Harrell f.harrell at Vanderbilt.Edu
Fri Aug 6 08:20:18 CEST 2010


In an upcoming release of the rms package, all fit objects can be 
printed using LaTeX if putting LaTeX code directly to the console 
(this is optimized for Sweave).  You will be able to say print(fit, 
latex=TRUE).

Frank E Harrell Jr   Professor and Chairman        School of Medicine
                      Department of Biostatistics   Vanderbilt University

On Thu, 5 Aug 2010, David Winsemius wrote:

>
> On Aug 5, 2010, at 4:03 PM, Biau David wrote:
>
>> Hello,
>>
>> I am modeling some survival data wih cph (Design). I have modeled a
>> predictor
>> which showed non linear effect with restricted cubic splines. I
>> would like to
>> retrieve the se(coef) for other, linear, predictors.
>
> The cph object has a "var". The vcov function is an extractor
> function. You would probably be using something like:
>
> diag(vcov(fit))^(1/2)
>
>> This is just to make nice
>> LateX tables automatically.
>
> Are you sure Frank has not already programed that for you somewhere?
> Perhaps latex.cph?
>
>> I have the coefficients with coef().
>>
>> How do I do that?
>>
>> Thanks,
>>
>> David Biau.
>>
>
> -- 
>
> David Winsemius, MD
> West Hartford, CT
>
> ______________________________________________
> R-help at 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.
>



More information about the R-help mailing list