为来自多项分布的样本构建多项 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

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

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