Построить полиномиальную модель логита для выборки из полиномиального распределения

Question

У меня есть данные от нескольких субъектов о том, какую долю своего времени они проводили в каждой четверти комнаты квадратной формы. Вот пример данных:

a <- rmultinom(12, 100, c(0.1, 0.4, 0.3, 0.2))
a / apply(a,2,sum)

### not running, but I would like to have something like this
# glm(a ~ 1, family=multinomial("logit")

Я был бы заинтересован в перехвате, так как я ожидаю, что они проведут большую часть своего времени, например, во 2 квартале.

Answer 1

Как указано в документации к пакету статистики, glm family = multinomial отсутствует. GLM не подходит для этого типа модели. Возможно, пакеты glmnet или nnet подойдут вам.

Есть несколько ресурсов на эту тему. Например, следующий учебник является первым результатом поиска в Google с использованием критериев: r полиномиальная регрессия.