Обратный прогноз
X0 = invpred(X,Y,Y0)
[X0,DXLO,DXUP] = invpred(X,Y,Y0)
[X0,DXLO,DXUP] = invpred(X,Y,Y0,name1,val1,name2,val2,...)
X0 = invpred(X,Y,Y0) принимает векторы X и Y той же длины, соответствует простой регрессии и возвращает ориентировочное значение X0, для которого высота строки равна Y0. Вывод, X0, имеет тот же размер как Y0, и Y0 может быть массивом любого размера.
[X0,DXLO,DXUP] = invpred(X,Y,Y0) также вычисляет 95%-е обратные интервалы прогноза. DXLO и DXUP задают интервалы с нижней границей X0–DXLO и верхняя граница X0+DXUP. И DXLO и DXUP имеют тот же размер как Y0.
Интервалы не одновременны и не обязательно конечны. Некоторые интервалы могут расширить от конечного значения до -Inf или +Inf, и некоторые могут расширить по целой действительной строке.
[X0,DXLO,DXUP] = invpred(X,Y,Y0, задает пары имя/значение дополнительного аргумента, выбранные из следующего списка. Имена аргумента являются нечувствительными к регистру, и частичные соответствия позволены.name1,val1,name2,val2,...)
| Имя | Значение |
|---|---|
'alpha' | Значение между 0 и 1 определением доверительного уровня % |
'predopt' | Или |
Сгенерируйте выборочные данные.
x = 4*rand(25,1); y = 10 + 5*x + randn(size(x));
Сделайте scatterplot данных.
scatter(x,y)
Предскажите x значение для данного y значения 20.
x0 = invpred(x,y,20)
x0 = 1.9967