library("ez") ezANOVA(data=Schule, dv=pangst, between=geschlecht, wid=vpnr, white.adjust=TRUE) library("car") library("lmtest") LinMod <- lm (pangst ~ geschlecht, data=Schule) coeftest (LinMod, vcov = hccm)