#######################
library(MASS)
pfad = "/Volumes/Data_1/d/down"
queen = read.table(paste(pfad, "queen.txt", sep="/"))
attach(queen)
plot(f0, Alter)

(b)
r = lm(Alter ~ f0 + I(f0^2))
summary(r)
stepAIC(r)

(c)
shapiro.test(resid(r))
plot(resid(r))
acf(resid(r))

(d)
k = coef(r)
plot(f0, Alter)
curve(k[1] + k[2]*x + k[3]*x^2, 200, 300, add=T)

(e)
plot(Alter, f0)
r = lm(f0 ~ Alter + I(Alter^2))
k = coef(r)
curve(k[1] + k[2]*x + k[3]*x^2, add=T)
neuerwert = data.frame(Alter = 100)
ergebnis = predict(r, neuerwert, se.fit=T)
ergebnis$fit