Обратное предсказание
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));
Сделайте разброс данных.
scatter(x,y)
Спрогнозировать значение x для заданного значения y, равного 20.
x0 = invpred(x,y,20)
x0 = 1.9967