Hoe doe je op de juiste manier planned contrasts in R?
Ik heb het bestand van Field (puppies_contrast.sav) in R geladen. Daarna kan ik volgens deze stappen inderdaad resultaten krijgen, zie hieronder.
Ik krijg per contrast de between groups variantie + F waarde, dat is mooi. Maar ik krijg niet de t-waarde zoals je die in SPSS krijgt. En de output van R bevat niet de output die je bij SPSS ziet onder "does not assume equal variances".
Mijn vraag: kan je in R ook de "Contrast Tests"-output krijgen zoals in SPSS?
-----
data <- rosetta::getData("puppies_contrast.sav")
levels(data$Dose)
# [1] "Control" "15 mins" "30 mins"
# vertellen welke groepen vergeleken moeten worden
c1 <- c(-2, 1, 1)
c2 <- c(0, -1, 1)
# matrix van maken
matrix = cbind(c1, c2)
# contrasten instellen voor Dose
contrasts(data$Dose) <- matrix
# en analyse doen
model_met_c <- aov(Happiness ~ Dose, data=data)
# omnibusresultaten
summary(model_met_c)
# Df Sum Sq Mean Sq F value Pr(>F)
# Dose 2 20.13 10.067 5.119 0.0247 *
# Residuals 12 23.60 1.967
# contrasten laten zien
summary.aov(model_met_c, split=list(Dose=list("Controle vs Experimenteel"=1,
"Kort vs Lang"=2)))
# Df Sum Sq Mean Sq F value Pr(>F)
# Dose 2 20.13 10.067 5.119 0.0247 *
# Dose: Controle vs Experimenteel 1 12.03 12.033 6.119 0.0293 *
# Dose: Kort vs Lang 1 8.10 8.100 4.119 0.0652 .
# Residuals 12 23.60 1.967
# ---
# Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1