På dette punktet har du ikke beskrevet «heteroscedasticity» -strukturen i gruppen din i modellen din. Prøv weights=varPower() som vist i eksemplet i ?gls . Det blir kvitt heteroscedastisiteten i ditt tilfelle.
Sammenlign:
m1 <- gls(salary ~ age*sex) plot(m1)
m2 <- gls(salary ~ age*sex, weights=varPower()) plot(m2)
Også hvis du ser i kapittel 5.2.1 (side 208) i modeller for blandede effekter i S og S -Plus av Pinheiro og Bates 2000, det er mye informasjon om variasjonsfunksjonene i nlme. Dette svaret kan også være nyttig: Regresjonsmodellering med ulik varians .