[R] extracting the standard error in lrm

[David Winsemius]

On Aug 11, 2010, at 6:56 AM, david dav wrote:

> Hi,
> I would like to extract the coefficients of a logistic regression
> (estimates and standard error as well) in lrm as in glm with
>
> summary(fit.glm)$coef
>

?coef

Try instead:

coef(fit.glm)

?vcov

The accessor function for the covariance matrix is vcov, so se's for  
Wald-type tests can be derived:

diag( vcov(fit.glm) )^0.5  # or
sqrt(diag( vcov(fit.glm))

--
David Winsemius, MD
West Hartford, CT