#Innlesing av data
d = read.table("exmp12-36.txt",header=TRUE,sep=",")
#Liste opp navnene p? de ulike variable
names(d)
#Tar vekk leverage variable
d = d[,-4]
#Lage interaksjonsplott
pairs(d)
#Korrelasjoner mellom variable
cor(d)

#Bruk av lm kommando
fit = lm(wingspan~height+foot,data=d)
#Finne ut hvilken summary kommando som brukes
class(fit)
help(summary.lm)
#Summary, tar med korrelasjonsmatrise for beta
summary(fit,correlation=TRUE)
anova(fit)
plot(fit)

#Leverage
print(cbind(d,lm.influence(fit)$hat))
n = dim(d)[1]
tt = round(lm.influence(fit)$hat,2)
plot(d$foot,d$height);text(d$foot,d$height+0.3,tt)