Interwał przewidywania regresji liniowej

24

Jeśli najlepszym przybliżeniem liniowym (przy użyciu najmniejszych kwadratów) moich punktów danych jest linia y=mx+b , jak mogę obliczyć błąd przybliżenia? Jeśli obliczę odchylenie standardowe różnic między obserwacjami i prognozami ei=real(xi)(mxi+b) , czy mogę później powiedzieć, że rzeczywista (ale nie zaobserwowana) wartość yr=real(x0) należy do przedziału ( y p = m x 0 + b ) z prawdopodobieństwem ~ 68%, przy założeniu rozkładu normalnego?[ypσ,yp+σ]yp=mx0+b

Wyjaśnić:

Dokonałem obserwacji dotyczących funkcji , oceniając ją w kilku punktach x i . Dopasowuję te obserwacje do linii l ( x ) = m x + b . Dla x 0 , którego nie obserwowałem, chciałbym wiedzieć, jak duże może być f ( x 0 ) - l ( x 0 ) . Stosując powyższą metodę, czy poprawne jest stwierdzenie, że f ( x 0 ) [ l ( x 0f(x)xil(x)=mx+bx0f(x0)l(x0)f(x0)[l(x0)σ,l(x0)+σ] z prob. ~ 68%?

BMX
źródło
1
Myślę, że pytasz o przedziały prognozowania. Pamiętaj jednak, że używasz „ ” zamiast „ y i ”. Czy to literówka? My nie przewidują x s. xiyix
gung - Przywróć Monikę
@gung: Używam do oznaczenia na przykład czasu, ay wartości jakiejś zmiennej w tym czasie, więc y = f ( x ) oznacza, że ​​dokonałem obserwacji y w czasie x . Chcę wiedzieć, jak daleko mogą być przewidywania funkcji dopasowania od rzeczywistych wartości y. Czy to ma sens? Funkcja r e a l ( x i ) zwraca „poprawną” wartość y dla x i , a moje punkty danych składają się z ( x i , r e a (xyy=f(x)yxreal(xi)yxi . (xi,real(xi))
bmx
1
To wydaje się całkowicie rozsądne. Części, na których się koncentruję, to np. „ ”, zwykle myślimy o błędach / resztkach w modelu reg jako „ e i = y i - ( m x i + b ) ”. SD z reszt ma odgrywać rolę w obliczaniu przedziały predykcji. Chodzi o to, że „ x iei=real(xi)(mxi+b)ei=yi(mxi+b)xi„to dla mnie dziwne; zastanawiam się, czy to literówka, czy pytasz o coś, czego nie rozpoznaję.
Gung - Przywróć Monikę
Myślę, że widzę; Brakowało mi twojej edycji. Sugeruje to, że system jest w pełni deterministyczny i jeśli miał dostęp do rzeczywistej funkcji podstawowej, zawsze można przewidzieć, doskonale w / o błędzie. Nie tak zazwyczaj myślimy o modelach reg. yi
gung - Przywróć Monikę
4
bmx, Wydaje mi się, że masz jasne pojęcie o swoim pytaniu i dobrą świadomość niektórych problemów. Być może zainteresuje Cię przegląd trzech ściśle powiązanych wątków. stats.stackexchange.com/questions/17773 opisuje przedziały prognozowania w terminach nietechnicznych; stats.stackexchange.com/questions/26702 daje bardziej matematyczny opis; a na stronie stats.stackexchange.com/questions/9131 Rob Hyndman zapewnia formułę, której szukasz. Jeśli nie odpowiedzą w pełni na twoje pytanie, przynajmniej mogą dać ci standardową notację i słownictwo, aby je wyjaśnić.
whuber

Odpowiedzi:

30

@whuber wskazał ci trzy dobre odpowiedzi, ale być może nadal mogę napisać coś wartościowego. Wasze wyraźne pytanie, jak rozumiem, brzmi:

Biorąc pod uwagę, my wyposażona y^i=m^xi+b^ (zauważ, że dodana 'kapeluszy') , i przy założeniu, że reszty moi zazwyczaj rozproszone, , można przewidzieć, że dotychczas zauważony reakcji, r n e w , o znanej wartości czynnikiem, x n e w , będzie mieścić się w przedziale ( Y - σ e , y + σN(0,σ^e2)ynewxnew , z prawdopodobieństwem 68%?(y^σe,y^+σe)

Intuitively, the answer seems like it should be 'yes', but the true answer is maybe. This will be the case when the parameters (i.e., m,b, & σ) are known and without error. Since you estimated these parameters, we need to take their uncertainty into account.

tdf errort

y^new±t(1α/2, df error)s, instead of y^new±z(1α/2)s, and go about our merry way? Unfortunately, no. The bigger issue is that there is uncertainty about your estimate of the conditional mean of the response at that location due to the uncertainty in your estimates m^ & b^. Thus, the standard deviation of your predictions needs to incorporate more than just serror. Because variances add, the estimated variance of the predictions will be:

spredictions(new)2=serror2+Var(m^xnew+b^)
Notice that the "x" is subscripted to represent the specific value for the new observation, and that the "s2" is correspondingly subscripted. That is, your prediction interval is contingent on the location of the new observation along the x axis. The standard deviation of your predictions can be more conveniently estimated with the following formula:
spredictions(new)=serror2(1+1N+(xnewx¯)2(xix¯)2)
As an interesting side note, we can infer a few facts about prediction intervals from this equation. First, prediction intervals will be narrower the more data we had when we built the prediction model (this is because there's less uncertainty in m^ & b^). Second, predictions will be most precise if they are made at the mean of the x values you used to develop your model, as the numerator for the third term will be 0. The reason is that under normal circumstances, there is no uncertainty about the estimated slope at the mean of x, only some uncertainty about the true vertical position of the regression line. Thus, some lessons to be learned for building prediction models are: that more data is helpful, not with finding 'significance', but with improving the precision of future predictions; and that you should center your data collection efforts on the interval where you will need to be making predictions in the future (to minimize that numerator), but spread the observations as widely from that center as you can (to maximize that denominator).

Having calculated the correct value in this manner, we can then use it with the appropriate t distribution as noted above.

gung - Reinstate Monica
źródło