r语言 - 为多项分布中的样本构建多项式 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")

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

Answer 1

如统计包的文档中所述，没有 glm 族 = 多项式。 GLM 不适合这种类型的模型。也许 glmnet 或 nnet 软件包将满足您的需求。

有几个关于这个主题的资源。例如，以下教程是使用以下条件进行 Google 搜索的第一个结果：r 多项式回归。