## Aufgabe 1
daten <- read.table("C:\\Dokumente und Einstellungen\\Administrator\\Desktop\\Wolfgang\\Öko\\Blatt 5\\swisslabor.dat", header=TRUE)
attach(daten)

## (a)
glm.labor1 <- glm(Teilnahme ~ 1+Alter+Bildung+JKinder+AKinder+Herkunft, family=binomial("logit"))
summary(glm.labor1)

## (b)
glm.labor2 <- glm(Teilnahme ~ 1+Alter+Bildung+JKinder+AKinder+Herkunft+I(Alter^2), family=binomial("logit"))
summary(glm.labor2)

## (e)
glm.swiss1 <- glm(Teilnahme ~ 1+Alter, family=binomial("logit"))
glm.swiss2 <- glm(Teilnahme ~ 1+Alter+Bildung+JKinder+AKinder+Herkunft, family=binomial("logit"))
anova(glm.swiss1, glm.swiss2, test="Chisq")

## (f)
glm.swiss1 <- glm(Teilnahme ~ 1+Alter+Bildung+JKinder+AKinder, family=binomial("logit"))
glm.swiss2 <- glm(Teilnahme ~ 1+Alter+Bildung+JKinder+AKinder+Herkunft, family=binomial("logit"))
anova(glm.swiss1, glm.swiss2, test="Chisq")


## Aufgabe 2
detach(daten)
daten <- read.table("C:\\Dokumente und Einstellungen\\Administrator\\Desktop\\Wolfgang\\Öko\\Blatt 5\\motorins.dat", header=TRUE)
attach(daten)

## (a)
glm.ins <- glm(Schaeden ~ 1+Kilometer+Zone+Bonus+Marke+Versicherte, family=poisson(link="log"))
summary(glm.ins)

## (b)
glm.ins1 <- glm(Schaeden ~ 1, family=poisson(link="log"))
glm.ins2 <- glm(Schaeden ~ 1+Kilometer+Zone+Bonus+Marke+Versicherte, family=poisson(link="log"))
anova(glm.ins1, glm.ins2, test="Chisq")

## (c)
glm.ins1 <- glm(Schaeden ~ 1+Kilometer+Bonus+Marke+Versicherte, family=poisson(link="log"))
glm.ins2 <- glm(Schaeden ~ 1+Kilometer+Zone+Bonus+Marke+Versicherte, family=poisson(link="log"))
anova(glm.ins1, glm.ins2, test="Chisq")