为多项分布中的样本构建多项logit模型

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")

我会对拦截感兴趣，因为我的期望是他们将大部分时间花在例如第二季度。

Answer 1

正如stats包的文档中所述，没有glm family = multinomial。 glm不适合这种类型的模型。也许glmnet或nnet软件包可以满足您的需求。

此主题有几个资源。例如，以下tutorial是使用条件进行Google搜索的第一个结果：r multinomial regression。