R kode for ekstraoppgave 7

 

 

# Punkt a)

# Vi leser inn dataene tilpasser en line?r regresjonsmodell med alle forklaringsvariablene:

leiligheter=read.table("http://www.uio.no/studier/emner/matnat/math/STK1110/h10/leiligheter.txt",header=T)

attach(leiligheter)

fit.a=lm(pris~areal+rom+leie+balkong+garasje+x+y)

summary(fit.a)

 

 

# Punkt b)

# Vi fjerner etter tur de forklaringsvariblene som har h?yest P-verdi (?n om gangen, og s? lenge P-verdien er over 5%):

 

#Fjerner rom

fit.b1=lm(pris~areal+leie+balkong+garasje+x+y)?

summary(fit.b1)

 

#Fjerner y

fit.b2=lm(pris~areal+leie+balkong+garasje+x)?

summary(fit.b2)

 

#Fjerner balkong

fit.b3=lm(pris~areal+leie+garasje+x)?

summary(fit.b3)

 

 

# Punkt c)

# Vi lager ulike plot av residualene

 

# Normalfordelingplott av residualene (for sjekk av normalitet):

qqnorm(fit.b3$res)

 

# Plott av residualene mot de tilpassede verdiene (for sjekk av lik varians)

plot(fit.b3$fit,fit.b3$res)

 

# Plott av residualene mot de numeriske forklaringsvariablene (for sjekk av linearitet)

plot(leiligheter$areal,fit.b3$res)

plot(leiligheter$leie,fit.b3$res)

plot(leiligheter$x,fit.b3$res)