Difference between revisions of "Least squares modelling (2013)"

From Statistics for Engineering
Jump to navigation Jump to search
m
Line 6: Line 6:
| vimeoID3 = 60592349
| vimeoID3 = 60592349
| vimeoID4 = 60688909
| vimeoID4 = 60688909
| vimeoID5 =
| vimeoID5 = 60862565
| vimeoID6 =
| vimeoID6 =
| vimeoID7 =
| vimeoID7 =
Line 27: Line 27:
| video_download_link4_MP4_size = 290 M
| video_download_link4_MP4_size = 290 M
| video_notes4 =
| video_notes4 =
| video_download_link5_MP4 =
| video_download_link5_MP4 =http://learnche.mcmaster.ca/media/4C3-2013-Class-07C.mp4
| video_download_link5_MP4_size =  
| video_download_link5_MP4_size = 375 M
| video_notes5 =
| video_notes5 =
| video_download_link6_MP4 =  
| video_download_link6_MP4 =  

Revision as of 20:32, 1 March 2013

Class date(s): 08 February 2013
Nuvola mimetypes pdf.png (PDF) Course slides
Download video: Link (plays in Google Chrome) [352 M]

Download video: Link (plays in Google Chrome) [347 M]

Download video: Link (plays in Google Chrome) [352 M]

Download video: Link (plays in Google Chrome) [290 M]

Download video: Link (plays in Google Chrome) [375 M]

Course notes and slides


Software source code

Take a look at the software tutorial.


Code used in class

Least squares demo

x <- c(10, 8, 13, 9, 11, 14, 6, 4, 12, 7, 5)
y <- c(8.04, 6.95, 7.58, 8.81, 8.33, 9.96, 7.24, 4.26, 10.84, 4.82, 5.68)
plot(x,y)
model.ls <- lm(y ~ x)
summary(model.ls)

coef(model.ls)
confinf(model.ls)

names(model.ls)
model$resisduals
resid(model.ls)

plot(x, y)
abline(model.ls)


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
library(car)
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")