R/oppg10.4.R

In SiljeHermansen/laerdegR: Stottepakke Til "Laer Deg R"

oppg10.4<-function(){
  data(kap10, envir = environment())
  kap10 <- kap10

  #Partivalg modellert som funksjon av inntekt.
  mod2<-glm(formula=FrP~Skepsis+Utjevn+Inntekt,
            data=kap10,
            family = binomial(link="logit"))

  #Definere hypotetisk datasett
  nyedata<-data.frame(Skepsis=mean(kap10$Skepsis, na.rm=TRUE),
                      Utjevn=mean(kap10$Utjevn, na.rm=TRUE),
                      Inntekt=c(9, 10)
  )

  #Predikere sannsynligheter
  preds<-predict(mod2, nyedata, type = "response")
  preds

  #Svar
  paste("Sannsynligheten for a stemme FrP framfor Hoyre synker med",
        round(preds[1]-preds[2],4)*100 ,
        "prosent (fra", round(preds[1],4)*100, "til", round(preds[2],4)*100,
        ") nar inntekt oker fra 9 til 10 mens andre variabler holdes konstant.")

}