W tym momencie nie opisałeś jeszcze struktury heteroskedastyczności wewnątrz grupy w swoim modelu. Spróbuj weights=varPower(), jak pokazano w przykładzie w ?gls . To eliminuje heteroskedastyczność w Twoim przypadku.
Porównaj:
m1 <- gls(salary ~ age*sex) plot(m1)
m2 <- gls(salary ~ age*sex, weights=varPower()) plot(m2)
Także jeśli zajrzysz do rozdziału 5.2.1 (strona 208) w modelach z efektami mieszanymi w S i S – Oprócz tego Pinheiro i Bates 2000, jest dużo informacji o funkcjach wariancji w nlme. Ta odpowiedź może być również pomocna: Modelowanie regresji z nierówną wariancją .