Assumpties toetsen------------------------------------------------------------------------------------------------------------
Variantie: leveneTest(reken$score,group= reken$group) #levenetest kan voor meerdere
aantallen
Normaliteit: shapiro.test(reken$score)
Extra’s -----------------------------------------------------------------------------------------------------------------------------
Van numeriek naar categorisch as.factor
reken$group= as.factor(reken$group) -> maak de numerieke variabele "group" categorisch
attach(meren)
detach(meren)
om 6 decimalen te laten zien zodat met afronden beter kan bij de summary van fit
print(summary(fit),digits=6)
MS groups -----------------------------------------------------------------------------------------------------------------------
MSgroups met de hand
groep1= reken$score[reken$group==1]
groep2= reken$score[reken$group==2]
groep3= reken$score[reken$group==3]
gem= mean(reken$score)
gem1= mean(groep1)
gem2= mean(groep2)
gem3= mean(groep3)
sum1= sum(length(groep1)*((gem1-gem)^2))
sum2= sum(length(groep2)*((gem2-gem)^2))
sum3= sum(length(groep3)*((gem3-gem)^2))
teller= sum1+sum2+sum3
df= 2
MSgroep= teller/df
MS error--------------------------------------------------------------------------------------------------------------------------
Met de hand
df= length(reken$score)-3
s1= var(groep1)
s2= var(groep2)
s3= var(groep3)
sum1= s1*(length(groep1)-1)
sum2= s2*(length(groep2)-1)
sum3= s3*(length(groep3)-1)
teller= sum1+sum2+sum3
MSerror= teller/df
ANOVA ---------------------------------------------------------------------------------------------------------