[R] Obtaining p-values for coefficients from LRM function (package Design) - plaintext

Frank E Harrell Jr f.harrell at vanderbilt.edu
Sat Dec 13 20:08:57 CET 2008

joris meys wrote:
> Sent this mail in rich text format before. Excuse me for this.
> ------------------------
> Dear all,
> I'm using the lrm function from the package "Design", and I want to
> extract the p-values from the results of that function. Given an lrm
> object constructed as follows :
> fit <- lrm(Y~(X1+X2+X3+X4+X5+X6+X7)^2, data=dataset)
> I need the p-values for the coefficients printed by calling "fit".

The individual p-values are not very meaningful.  anova(fit) returns a 
matrix of meaningful tests with P-values.


> fit$coef (gives a list of only the coefficients)
> fit$pval, fit$p, fit$pvalue, fit$p.value,... : nothing works
> str(fit) : no hints there
> fit[1,4] : gives dimension errors
> help files don't seem to give me a function that extracts them. Yet,
> they are calculated and printed, based on the Wald statistics. So they
> must be reachable.
> Anybody knows how?
> Thank you in advance
> Kind regards
> Joris
> ______________________________________________
> 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.

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

More information about the R-help mailing list