An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20100107/4acb4256/attachment.pl>
Extract p-value from linear model
3 messages · Trafim Vanishek, Peter Ehlers, (Ted Harding)
coef(summary(lm(..stuff..)))[,4] or coef(summary(lm(..stuff..)))[,4,drop=FALSE] -Peter Ehlers
Trafim Vanishek wrote:
Dear all, I have the following question. Is it possible in R to call for the p-value directly in a model like this a <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) b <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) summary(lm(a~(b+0))) or even in this a <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) b <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) summary(lm(a~b)) to call the p-values directly for every estimated parameter. Thanks a lot for the help! [[alternative HTML version deleted]]
______________________________________________ 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.
Peter Ehlers University of Calgary 403.202.3921
On 07-Jan-10 14:15:07, Trafim Vanishek wrote:
Dear all, I have the following question. Is it possible in R to call for the p-value directly in a model like this a <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) b <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) summary(lm(a~(b+0)))
summary(lm(a~(b+0))) # [...] # Coefficients: # Estimate Std. Error t value Pr(>|t|) # b 0.90777 0.07134 12.72 4.67e-07 *** # --- summary(lm(a~(b+0)))$coef[,4] # [1] 4.665795e-07
or even in this a <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) b <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) summary(lm(a~b)) to call the p-values directly for every estimated parameter. Thanks a lot for the help!
summary(lm(a~b)) # [...] # Coefficients: # Estimate Std. Error t value Pr(>|t|) # (Intercept) 7.7957 2.1661 3.599 0.00699 ** # b -0.6230 0.4279 -1.456 0.18351 # --- summary(lm(a~b))$coef[,4] # (Intercept) b # 0.006992306 0.183511560 Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <Ted.Harding at manchester.ac.uk> Fax-to-email: +44 (0)870 094 0861 Date: 07-Jan-10 Time: 14:39:46 ------------------------------ XFMail ------------------------------