Skip to content

values from a linear model

5 messages · ONKELINX, Thierry, Dimitris Rizopoulos, Vladimir Eremeev +1 more

#
Dear R users,
how can I extrapolate values listed in the summary of an lm model but not 
directly available between object values such as the the standard errors of 
the calculated parameters?

for example I got a model:

mod <- lm(Crd ~ 1 + Week, data=data)

and its summary:
Call:
lm(formula = Crd ~ 1 + Week, data = data, model = TRUE, y = TRUE)

Residuals:
       Min         1Q     Median         3Q        Max 
-4.299e-03 -1.653e-03  2.628e-05  1.291e-03  5.130e-03 

Coefficients:
             Estimate Std. Error  t value Pr(>|t|)    
(Intercept) 1.000e+01  3.962e-04 25238.73   <2e-16 ***
Week        5.038e-04  6.812e-06    73.96   <2e-16 ***
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 

Residual standard error: 0.001966 on 98 degrees of freedom
Multiple R-Squared: 0.9824,     Adjusted R-squared: 0.9822 
F-statistic:  5469 on 1 and 98 DF,  p-value: < 2.2e-16

I'm interested in values of Std. Error of coefficients...

thank you very much

Best regards
	Manuele
#
It's not very clear to me but I think summary(mod)$coef[, "Std. Error"]
is wat you need?

Cheers,

Thierry

------------------------------------------------------------------------
----
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics,
methodology and quality assurance
Gaverstraat 4
9500 Geraardsbergen
Belgium
tel. + 32 54/436 185
Thierry.Onkelinx op inbo.be
www.inbo.be 

Do not put your faith in what statistics say until you have carefully
considered what they do not say.  ~William W. Watt
A statistical analysis, properly conducted, is a delicate dissection of
uncertainties, a surgery of suppositions. ~M.J.Moroney
#
try this:

coef(summary(mod))[, "Std. Error"]


I hope it helps.

Best,
Dimitris

----
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven

Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
     http://www.student.kuleuven.be/~m0390867/dimitris.htm


----- Original Message ----- 
From: "Manuele Pesenti" <amicogodzilla at bruttocarattere.org>
To: <r-help at stat.math.ethz.ch>
Sent: Tuesday, July 24, 2007 12:02 PM
Subject: [R] values from a linear model
Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
#
Manuele Pesenti wrote:
If you want to assign these values to some other variables, try assigning
the result of the summary() to a variable and working with its components
(the result is a list, use $ or [[]] to get its members)

mod.sum<-summary(mod)
then
coef(mod.sum)[,2]
or
mod.sum$coefficients[,2]

will give you those Std. Errors
#
On Tuesday 24 July 2007 12:02:58 Manuele Pesenti wrote:
thank you very much for all interesting answer
:)

	Manuele