Ударная гамма-регрессия в R

Я знаю пакеты http://cran.r-project.org/web/packages/glmnet/index.html и http://cran.r-project.org/web/packages/penalized/index.html, но ни один из них, похоже, не поддерживает Gamma GLM.Ударная гамма-регрессия в R

Я бы хотел использовать эластичную сетку для гамма-GLM в R, что это самый простой способ сделать это?

(мета: также обсуждает, следует ли это пойти на крест ПРОВЕРЯЕМОГО лучшие ответов?)

источник

2015-01-07 kevinykuo

Что касается мета: если это конкретно «как его закодировать в R?», Тогда SO отлично (но IMO это не будет offtopic на stats.SE либо); если это «какую статистическую процедуру я должен использовать?», то это прочная статистика. – NPE

Я думаю, что это для SO - ОП указал, что они хотят использовать эластичную сеть (они не спрашивают, должны ли они или нет). Несколько мыслей (не ответ на самом деле) (1) вы могли бы установить штраф за эластичную сетку (или, по крайней мере, я могу видеть, как вы создадите хребет или лассо - я не совсем уверен, как реализовать эластично- net) с помощью грубой силы в 'optim' или' bbmle :: mle2', но, вероятно, это будет довольно неэффективно; (2) у вас есть веская причина предпочитать Gamma над журнальным нормальным анализом (т. Е. Ответом на преобразование логарифма и штрафованием наименьших квадратов)? Latter был бы намного проще и удобнее. –

также, вроде как «ищет пакет» ... –

Вы можете использовать пакет HDtweedie. Гамма является частным случаем распределения Tweedie с p = 2. Это относительно новый пакет, поэтому ожидайте некоторые проблемы с прорезями.

источник

2015-01-20 16:50:32

Ударная гамма-регрессия в R

ответ

Смежные вопросы