pfad = "Verzeichnis wo Sie kj.txt und lneu.txt gespeichert haben"
dd = read.delim(paste(pfad, "kj.txt", sep="/"))
lneu = read.table(paste(pfad, "lneu.txt", sep="/"))

head(dd)
names(dd)
tab = with(dd, table(str, age))

barplot(tab, beside=T, legend=T)
g = glm(str ~ age, binomial, data=dd)
anova(g, test="Chisq")


head(lneu)
names(lneu)
tab = with(lneu, table(Vowel, Jahr))
p = tab[1,]/apply(tab, 1, sum)
zeit = as.numeric(colnames(tab))
plot(zeit, p, type="b", xlab="Jahr (Re. 1950)", ylab="Proportion HIGH")
g = glm(Vowel ~ Jahr, binomial, data=lneu)
anova(g, test="Chisq")

# zwei Faktoren
attach(lneu)
beide = paste(G, Jahr, sep=".")
tab = table(beide, Vowel)

temp = Vowel == "high"
interaction.plot(Jahr, G, temp)
mehrg = glm(Vowel ~ Jahr + G, binomial, data=lneu)
anova(mehrg, test="Chisq")
g = glm(Vowel ~ Jahr + G + Jahr:G, binomial, data=lneu)
g = glm(Vowel ~ Jahr * G, binomial, data=lneu)
anova(g, test="Chisq")
library(MASS)
stepAIC(g)