Vraag 9 Test: Voorwaarden

Het is vaak niet eenvoudig om te beoordelen of de varianties sterk van mekaar verschillen of niet. Om een beter idee te krijgen, kunnen we eens een aantal boxplots simuleren met dezelfde steekproefgrootte en in de veronderstelling dat de varianties gelijk zijn. Daarnaast bestaat er ook een test in R die test of de varianties al dan niet gelijk zijn. De nulhypothese van deze test is dat de varianties gelijk zijn, de alternatieve hypothese is dat ze verschillend zijn.

set.seed(52)
par(mfrow=c(3,3), mar=c(3,2,1,1))
sd1<-sd(stokstaartjes_wide$verschil) 
mean1<-mean(stokstaartjes_wide %>% filter(verbannen == 1) %>% pull(verschil))
mean2<-mean(stokstaartjes_wide %>% filter(verbannen == 0) %>% pull(verschil))
for(i in 1:9){
  x<-rnorm(15, mean=mean1, sd=sd1)
  y<-rnorm(15, mean=mean2, sd=sd1)
  boxplot(x,y)
}

var.test(stokstaartjes_wide$verschil~stokstaartjes_wide$verbannen)

#>   F test to compare two variances
#> data:  stokstaartjes_wide$verschil by stokstaartjes_wide$verbannen
#> F = 0.55608, num df = 14, denom df = 14, p-value = 0.2842
#> alternative hypothesis: true ratio of variances is not equal to 1
#> 95 percent confidence interval:
#>   0.1866918 1.6563266
#> sample estimates:
#> ratio of variances 
#>          0.5560779