[R] Getting Pr from Summary(lm)

[Peter Dalgaard]

"David J. Netherway" <david.netherway at adelaide.edu.au> writes:

> Hello,
> 
> I am trying to get the P values from the output of a summary for lm.
> 
> lm <- lm(y ~ age + sex)
> s <- summary(lm)
> 
> I thought that I might be able to get them using a combination of
> scan, grep and sub.
> But I got stuck on the first step - being able to process "s" as a
> text string.
> I could perhaps write it to file than scan it back but there is
> probably an easier
> way to do the whole thing.
> 
> Help would be welcome, David

Summaries are (usually) R objects which can be manipulated
programmatically. They should be confused with what their print
methods print. Consider

example(summary.lm)
coef(sld90)[,4,drop=F]


-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45) 35327907