Я пытаюсь запустить полиномиальную обратную регрессию, как описано здесь .Пакет R называется textir.Это работает просто отлично, но я не могу отобразить параметры в пакете R (называемые «nlambda» и «gamma») на (как я понимаю) гипер-параметры модели, как описано в статье (параметры формы и масштаба rи s гамма-гиперприора для каждой лямбды).
install.packages("textir")
library("textir")
data(we8there)
cl <- NULL
fits <- dmr(cl, we8thereRatings[,'Overall',drop=FALSE],
we8thereCounts, bins=5, gamma=1, nlambda=100)
Используя это в качестве отправной точки, я хотел бы воспроизвести Приложение 3 «A.3 Детали исследования прогноза вне выборки», третья панель («Рейтинг ресторана») на странице 31бумага, где авторы изменяют параметр формы "s".Я не вижу, как это реализовано в пакете R.
Интересно, что в другой статье Мэтта Тэдди он использует
fit <- mnlm(cl=NULL,counts=X, covars=V, bins=3,verb=TRUE, penalty=c(4,1))
и передает параметр penalty для предоставления двух чисел (mnlm() - это просто другой способдля вызова dmr(), согласно textir документации).Я не могу найти документацию по этому параметру ни в одном из связанных пакетов textir, dmr и gamlr.Кроме того, если вы отбросите их или измените их на абсурдные значения (даже на строки типа «bla»), вывод практически не изменится.