multinom：如何获得观察数量

Question

我正在使用multinom中的nnet从大量临床数据库中创建逻辑模型。我使用的语法是：

library(nnet)
fit=multinom(group ~ sex + age + var3 + var4,
            data=d, na.action = na.omit)

现在，每一行（患者）都有不同数量的NA，因为并非所有患者都记录了所有临床数据，而且我不清楚该模型是否仅使用所有变量所做的那些行不包含任何NA。 更一般地说，获得模型所基于的观察值的Ns是有用的，我怀疑这个数据帧的行数（N个患者）小。 我到处都看，但我似乎无法找到如何做到这一点。

Answer 1

我相信length(residuals(fit))应该有用。

如果要浏览用于拟合模型的数据集，可以执行以下操作：

mf <- model.frame(group ~ sex + age + var3 + var4, data=d,
                  na.action=na.omit)

（这是multinom处理您的数据的方式），然后计算行数，将不同类别的观察数量列表为不同的变量（例如table(mf$var3)）等。{{1应按类别列出每个变量的观察数量。

您可能会发现lapply(mf,table)包中的describe函数有用：

Hmisc