Difference between revisions of "Least squares modelling (2013)"
Jump to navigation
Jump to search
Kevin Dunn (talk | contribs) m |
Kevin Dunn (talk | contribs) |
||
| Line 56: | Line 56: | ||
Thermocouple data | Thermocouple data | ||
<syntaxhighlight lang=" | <syntaxhighlight lang="sas"> | ||
V <- c(0.01, 0.12, 0.24, 0.38, 0.51, 0.67, 0.84, 1.01, 1.15, 1.31) | V <- c(0.01, 0.12, 0.24, 0.38, 0.51, 0.67, 0.84, 1.01, 1.15, 1.31) | ||
T <- c(273, 293, 313, 333, 353, 373, 393, 413, 433, 453) | T <- c(273, 293, 313, 333, 353, 373, 393, 413, 433, 453) | ||
Revision as of 16:00, 27 February 2013
| Class date(s): | 08 February 2013 | ||||
 (PDF)
|
Course slides | ||||
| |||||
| |||||
| |||||
Course notes and slides
- Course textbook (print chapter 4)
- Slides for class
Software source code
Take a look at the software tutorial.
Code used in class
Thermocouple data
V <- c(0.01, 0.12, 0.24, 0.38, 0.51, 0.67, 0.84, 1.01, 1.15, 1.31)
T <- c(273, 293, 313, 333, 353, 373, 393, 413, 433, 453)
plot(V, T)
model <- lm(T ~ V)
summary(model)
coef(model)
confint(model) # get the coefficient confidence intervals
resid(model) # model residuals
qqPlot(resid(model)) # q-q plot of the residuals to check normality
plot(V, T)
v.new <- seq(0, 1.5, 0.1)
t.pred <- coef(model)[1] + coef(model)[2]*v.new
lines(v.new, t.pred, type="l", col="blue")
# Plot x against the residuals to check for non-linearity
plot(V, resid(model))
abline(h=0)
# Plot the raw data and the regression line in red
plot(V, T)
abline(model, col="red")