Остатки и графики в упорядоченной полиномиальной регрессии

Question

Мне нужно построить график остаточного бинарного графика с подгоночными и остаточными значениями из упорядоченной мультиноминальной логит-регрессии.

Как извлечь остатки при использовании polr?Есть ли какая-либо другая функция, выполняющая полиноминальный логит, в котором можно извлечь остатки?

Это код, который я использовал

options(contrasts = c("contr.treatment", "contr.poly"))
mod1 <- polr(as.ordered(y) ~ x1 + x2 + x3, data, method='logistic')

fit <- mod1$fitted.values
res <- residuals(mod1)
binnedplot(fit, res)

Проблема в том, что объект 'res' равен 'null'.

Спасибо

Answer 1

В polr () нет функции, которая возвращает остаток.Вы должны вручную рассчитать его, используя его определение.

Answer 2

Для начала, можете ли вы сказать нам, как будут определены остатки в принципе для модели с категориальными ответами?fitted.values - это матрица вероятностей.Вы можете определить остатки в терминах правильного прогноза (определяя наиболее вероятный результат как прогноз, как в методе predict по умолчанию для polr объектов) - или вы можете вычислить таблицу истинных значений n-by-nпрогнозные значения.В качестве альтернативы вы можете уменьшить порядковые данные до целочисленной шкалы и вычислить среднее значение результата в качестве прогноза ... но я не вижу, чтобы вообще был какой-то уникальный способ определения остатков.