problem-10.9
problem-10.9
For illustration, we type the data into a text file with two columns,
the first being the elevation. Then we use read.table() to read it in.
> x = read.table(file="tmp.txt")
> names(x) = c("elev","Temp")
> res = lm(Temp ~ elev, x)
> summary(res)
...
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 59.290662 1.717329 34.52 3.9e-08 ***
elev -0.005115 0.000921 -5.55 0.0014 **
---
...
> T = (-0.005115 - (-.00534)) / .000921
> T
[1] 0.2443
> 2*pt(T,df=6,lower.tail=FALSE) # two-sided
[1] 0.8151
The p-value is not statistically significant.