#Der Datensatz steht leider nicht zur VerfŸgung
sgs <- read.table("sgs.forMondrian2",header=T,sep="\t", quote="")
attach(sgs)
names(sgs)

plot(THQ,GEWICHT,pch=19)

iplot(THQ,GEWICHT)


M01<-glm(GEWICHT~THQ+FL+BIP)
summary(M01)
M01

#Korrelationen von ParameterschŠtzern
vcov(M01)
v<-vcov(M01)
vv<-1/sqrt(diag(v))
c1<-diag(vv)
corr<-c1%*%v%*%c1
corr

anova(M01)

#FŸr BoxCox braucht man Y>0
min(GEWICHT)
sg2<-subset(sgs,GEWICHT>0)
dim(sg2)
dim(sgs)
attach(sg2)

#Modelle
M01<-glm(GEWICHT~THQ+FL+BIP)
summary(M01)
boxcox(M01)

gw<-sqrt(GEWICHT)
M01s<-glm(gw~THQ+FL+BIP)
summary(M01s)
plot(M01s)