savings
g1 <- lm(sr~pop75+dpi+ddpi,data=savings)
gx1 <- lm(pop15~pop75+dpi+ddpi,data=savings)
plot(gx1$res,g1$res)
plot(gx1$res,savings$pop15)
plot(savings$pop15,g1$res)
plot(gx1$res,savings$pop15)
plot(gx1$res,g1$res)
g <- lm(sr~.,data=savings)
summary(g)
prplot(g)
prplots(g)
cr.plots(g)
cr.plot(g)
plot(g)
plot(g)
help(prplot)
g
pr <- 
savings$sr-g$coef[1]-g$coef[3]*savings$pop75-g$coef[4]*savings$dpi-g$c
oef[5]*savings$ddpi
plot(savings$pop15,pr)
sb35 <- savings[savings$pop15<35,]
sb35
sa35 <- savings[savings$pop15>35,]
ga35 <- lm(savings~.,sa35)
ga35 <- lm(sr~.,sa35)
gb35 <- lm(sr~.,sb35)
summary(ga)
summary(ga35)
summary(gb35)
source("boxcox.R")
source("boxcox.S")
boxcox
x <- model.matrix(g)
x
y <- savings$sr
lambda <- seq=(0.5,1.5,by=0.1)
lambda <- seq(0.5,1.5,by=0.1)
boxcox(x,y,lambda)
bc <- boxcox(x,y,lambda)
warnings()
x <- x[,-1]
bc <- boxcox(x,y,lambda)
plot(bc$lambda,bc$loglik)
plot(bc$lambda,bc$loglik,type="l")