Skip to content
Prev 32940 / 398503 Next

Regression slopes

Martin Biuw wrote:
There might be an easier way, but the brute force method would be:

R> set.seed(1)
R> x <- data.frame(x = 1:10, y = 1:10 + rnorm(10))
R> lm.x <- lm(y ~ x, data = x)
R> coef.x <- summary(lm.x)$coef
R> # t-statistic comparing slope to 1
R> t.x <- (coef.x["x","Estimate"]-1)/coef.x["x","Std. Error"]
R> # p-value
R> 2 * pt(abs(t.x), lm.x$df, lower = FALSE)
[1] 0.5559868

Hope this helps,

Sundar