이 시점에서 아직 모델에서 그룹 내이 분산 구조를 설명하지 않았습니다. ?gls
의 예에 표시된대로 weights=varPower()
를 사용해보세요. . 그러면 귀하의 경우 이분산성이 제거됩니다.
비교 :
m1 <- gls(salary ~ age*sex) plot(m1)
m2 <- gls(salary ~ age*sex, weights=varPower()) plot(m2)
S와 S의 혼합 효과 모델에서 5.2.1 장 (208 페이지)을 보면 -Pinheiro 및 Bates 2000에 의해 nlme
의 분산 함수에 대한 많은 정보가 있습니다.이 답변도 도움이 될 수 있습니다. 불균등 분산을 사용한 회귀 모델링 .