Обратное предсказание
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