################################################################################
#
# Beispiel: Wirksamkeit einer Blutdrucktherapie
#
# Vorlesung Lineare Modelle (Helmut Küchenhoff) 28.05.2014
# Thema: Stetige Einflussgrößen 
#
# data: Daten von 44 Patienten, denen randomisiert eine Blutdrucktherapie 
#       zugeordnet wurde
#       file: blutdruck.sav (http://www.statistik.lmu.de/~abender/daten/Blutdruck.sav) 
#
# contents: (1) Einlesen der Daten
#		        (2) Modellbildung
#
# last.modified: AB 30.05.2014
#
################################################################################
setwd("C:/Users/Küchenhoff/Desktop/RLimo")
#### (1) Einlesen der Daten ####################################################
library(foreign)
blutdruck <- read.spss('Blutdruck.sav', 
  use.value.labels=TRUE, max.value.labels=Inf, to.data.frame=TRUE)

show(blutdruck)
#### (2) Modellbildung #########################################################
library(car)
library(effects)
library(lattice)
library(nlme)
## umbenennung der Levels
levels(blutdruck$Gruppe)<-c("K1","K2","E2S","E3S")
### (a) einfaches Modell mit Haupteffekten
lm1 <- lm(Ende ~ Gruppe + Anfang , data=blutdruck)
summary(lm1)
plot(lm1)
Anova(lm1)
plot(residuals(lm1)~blutdruck$Gruppe)
gls1<-gls(Ende ~ Gruppe + Anfang ,
    weights=varIdent(form=~1|Gruppe),method="REML", data=blutdruck)
summary(gls1)

# Darstellung der Effekt
plot(allEffects(gls1), ask=FALSE)
dev.off()