Mam wykres wartości resztkowych modelu liniowego w funkcji dopasowanych wartości, w których heteroscedastyczność jest bardzo wyraźna. Jednak nie jestem pewien, jak powinienem postępować teraz, ponieważ o ile rozumiem ta heteroscedastyczność powoduje, że mój model liniowy jest nieważny. (Czy to prawda?)
Użyj solidnego dopasowania liniowego za pomocą
rlm()
funkcjiMASS
pakietu, ponieważ jest najwyraźniej odporny na heteroscedastyczność.Ponieważ standardowe błędy moich współczynników są błędne z powodu heteroscedastyczności, mogę po prostu dostosować standardowe błędy, aby były odporne na heteroscedastyczność? Korzystając z metody opublikowanej w sekcji Przepełnienie stosu tutaj: regresja z heterometryczną korekcją błędów standardowych
Której metody najlepiej użyć do rozwiązania mojego problemu? Jeśli użyję rozwiązania 2, to czy moje możliwości przewidywania mojego modelu są całkowicie bezużyteczne?
Test Breuscha-Pagana potwierdził, że wariancja nie jest stała.
Moje resztki w funkcji dopasowanych wartości wyglądają następująco:
(większa wersja)
gls
i jednej ze struktur wariancji z pakietu nlme.