library(faraway)
data(chicago)
chicago
chicago[1:10,]
plot(chicago)
summary(ch)
summary(chicago)
ch <- data.frame(chicago[,c(1,2,3,4,6)],income=chicago[,7]/1000)
ch
boxplot(ch)
boxplot(ch$involunact)
boxplot(ch$involact)
summary(lm(involact ~ race))
summary(lm(involact ~ race,ch))
g <- lm(involact~race+fire+theft+age+income,ch)
summary(g)
summary(ch)
g <- lm(involact~race+fire+theft+age+log(income),ch)
summary(g)
plot(g$fit,g$res)
qqnorm(g$res)
plot(cooks.distance(g))
text(1:47,cooks.distance(g))
ch(c(6,24))
ch[c(6,24),]
summary(ch$fire)
ch$fire
sort(ch$fire)
ch[c(6,24),]
sort(ch$fire)
sort(ch$theft)
ch[6,]
x <- model.matrix(g)
lev <- hat(x)
plot(lev)
text(1:47,lev)
plot(g$res)
g <- lm(involact~rtace+fire+theft+age+log(income),data=ch,subset=(1:47)[-c(6,2
4)])
g <- lm(involact~race+fire+theft+age+log(income),data=ch,subset=(1:47)[-c(6,24
)])
summary(g)
g <- lm(involact~race+fire,data=ch,subset=(1:47)[-c(6,24)])
summary(g)
library(leaps)
y <- ch$involact
a <- leaps(x,y)
x
a <- leaps(x,y)
x
a <- leaps(x[,-1],y)
a
plot(a$size,a$Cp)
abline(0,1)
plot(a$size,a$Cp,ylim=c(0,20))
abline(0,1)
x <- x[-c(6,24),]
x <- model.matrix(g)
x
length(x[,1])
g <- lm(involact~race+fire+theft+age,data=ch,subset=(1:47)[-c(6,24)])
summary(g)
galt <- lm(involact~race+fire+log(income),data=ch,subset=(1:47)[-c(6,24)])
summary(galt)
galt <- lm(involact~race+fire+log(income),data=ch)
summary(galt)
galt <- lm(involact~race+fire+theft+age,data=ch)
summary(galt)
data(chiczip)
chiczip
galt <- lm(involact~race+fire+theft+age+chiczip,data=ch)
ch <- data.frame(ch,loca=chiczip)
galt <- lm(involact~race+fire+theft+age+chiczip,data=ch)
galt <- lm(involact~race+fire+theft+age+loca,data=ch)
summary(galt)
g <- lm(involact~race+fire+theft+age,subset=(chiczip=="s"),ch)
summary(g)
g <- lm(involact~race+fire+theft+age,subset=(chiczip=="n"),ch)
summary(g)