q()
library(faraway)
data()
data(strongx)
ls()
strongx
plot(crossx,energy)
plot(crossx,energy,data=strongx)
plot(strongx$crossx,strongx$energy)
g <- lm(crossx~energy,data=strongx)
plot(strongx$energy,strongx$crossx)
abline(g)
plot(strongx$crossx,g$res)
abline(h=0)
plot(g$fit,g$res)
abline(h=0)
plot(strongx$crossx,g$res)
data(savings)
savings
g <- lm(sr~pop15+pop75+dpi+ddpi,savings)
plot(g$fit,g$res)
plot(savings$sr,g$res)
plot(savings$pop15,g$res)
plot(savings$pop75,g$res)
plot(savings$pop15,g$res)
par(mfrow=c(3,3))
for(i in 1:9) plot(1:50,rnorm(50))
for(i in 1:9) plot(1:50,(1:50)*rnorm(50))
for(i in 1:9) plot(1:50,sqrt(1:50)*rnorm(50))
for(i in 1:9) plot(1:50,cos((1:50)*pi/25)*rnorm(50))
plot(savings$pop15,g$res)
par(mfrow=c(1,1))
plot(savings$pop15,g$res)
var(g$res[savings$pop15 < 35])
var(g$res[savings$pop15 > 35])
ls()
savings[savings$pop15 > 35,])
savings[(savings$pop15 > 35),])
savings[(savings$pop15 > 35),]
savings[(savings$pop15 < 35),]
savings[(savings$pop15 > 35),]
plot(g$fit,abs(g$res))
data(gala)
gala
g <- lm(Species ~ Area + Elevation + Scruz + Nearest + Adjacent, gala)
plot(g$fit,g$res)
gs <- lm(sqrt(Species ~ Area + Elevation + Scruz + Nearest + Adjacent, gala)
)
gs <- lm(sqrt(Species) ~ Area + Elevation + Scruz + Nearest + Adjacent, gala)
plot(gs$fit,gs$res)
qqnorm(gs$res)
par(mfrow=c(3,3))
for(i in 1:9) qqnorm(rnorm(50))
for(i in 1:9) qqnorm(exp(rnorm(50)))
for(i in 1:9) qqnorm(runif(50)))
for(i in 1:9) qqnorm(runif(50))
for(i in 1:9) qqnorm(rnorm(50))
for(i in 1:9) qqnorm(rt(50,5))
for(i in 1:9) qqnorm(rt(50,20))
for(i in 1:9) qqnorm(rt(50,2))
for(i in 1:9) qqnorm(rcauchy(50))
halfnorm(cooks.distance(g))
par(mfrow=c(1,1))
g <- lm(savings~pop15+pop75+dpi+ddpi,data=savings)
g <- lm(sr~pop15+pop75+dpi+ddpi,data=savings)
cooks.distance(g)
plot(cooks.distance(g))
halfnorm(cooks.distance)
halfnorm(cooks.distance(g))
halfnorm(sqrt(cooks.distance(g)))
q()