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:> summary(mod)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 -- Manuele Pesenti manuele a inventati.org amicogodzilla a jabber.linux.it http://mpesenti.polito.it
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> -----Oorspronkelijk bericht----- > Van: r-help-bounces op stat.math.ethz.ch > [mailto:r-help-bounces op stat.math.ethz.ch] Namens Manuele Pesenti > Verzonden: dinsdag 24 juli 2007 12:03 > Aan: r-help op stat.math.ethz.ch > Onderwerp: [R] values from a linear model > > 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: > > > summary(mod) > > 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 > > -- > Manuele Pesenti > manuele op inventati.org > amicogodzilla op jabber.linux.it > http://mpesenti.polito.it > > ______________________________________________ > R-help op stat.math.ethz.ch 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. >
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
> 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:
>
>> summary(mod)
>
> 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
>
> --
> Manuele Pesenti
> manuele at inventati.org
> amicogodzilla at jabber.linux.it
> http://mpesenti.polito.it
>
> ______________________________________________
> R-help at stat.math.ethz.ch 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.
>
Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
Manuele Pesenti wrote:> > 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: > >> summary(mod) > > 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 >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 -- View this message in context: http://www.nabble.com/values-from-a-linear-model-tf4134911.html#a11760459 Sent from the R help mailing list archive at Nabble.com.
On Tuesday 24 July 2007 12:02:58 Manuele Pesenti wrote:> 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?thank you very much for all interesting answer :) Manuele -- Manuele Pesenti manuele a inventati.org amicogodzilla a jabber.linux.it http://mpesenti.polito.it