Mam dane, które dopasowałem za pomocą modelu LOESS w R, co daje mi to:
Dane mają jeden predyktor i jedną odpowiedź i są heteroscedastyczne.
Dodałem również przedziały ufności. Problem polega na tym, że przedziały są przedziałami ufności dla linii, podczas gdy mnie interesują przedziały prognozowania. Na przykład dolny panel jest bardziej zmienny niż górny, ale nie jest to rejestrowane w odstępach.
To pytanie jest nieco powiązane:
Zrozumienie przedziału ufności z regresji wielomianowej , szczególnie odpowiedź @AndyW, jednak w swoim przykładzie używa stosunkowo prostego interval="predict"
argumentu, który istnieje predict.lm
, ale go nie ma predict.loess
.
Mam więc dwa bardzo powiązane pytania:
- Jak uzyskać przedziały prognozowania punktowego dla LOESS?
- Jak mogę przewidzieć wartości, które będą przechwytywać ten przedział, tj. Wygenerować wiązkę liczb losowych, które ostatecznie będą wyglądać trochę jak oryginalne dane?
Możliwe, że nie potrzebuję MIŁOŚCI i powinienem użyć czegoś innego, ale nie znam moich opcji. Zasadniczo powinien pasować do linii za pomocą regresji lokalnej lub wielokrotnej regresji liniowej, dając mi oszacowania błędów dla linii, a ponadto różne wariancje dla różnych zmiennych objaśniających, dzięki czemu mogę przewidzieć rozkład zmiennej odpowiedzi (y) przy pewnych wartościach x .
źródło
Odpowiedzi:
Nie wiem, jak zrobić pasma predykcyjne z oryginalną
loess
funkcją, aleloess.sd
wmsir
pakiecie jest taka funkcja ! Prawie dosłownie zmsir
dokumentacji:Twoje drugie pytanie jest nieco trudniejsze, ponieważ
loess.sd
nie ma funkcji przewidywania, ale możesz zhakować je razem, interpolując liniowo przewidywane środki i SD, z których się wydostajeszloess.sd
(używaszapprox
). Można je z kolei wykorzystać do symulacji danych przy użyciu rozkładu normalnego z przewidywanymi średnimi i SD:źródło
loess.sd
, nie różni się zbytnio od tego, co @rnso zasugerował w komentarzu do innego mojego pytania . Dzięki!