1

我正在使用multinomfromnnet从一个庞大的临床数据库创建一个逻辑模型。我使用的语法是：

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

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

4

1 回答 1

1

我相信length(residuals(fit))应该工作。

如果您想探索用于拟合模型的数据集，您可以执行以下操作：

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

（这是multinom处理您的数据的方式）然后计算行数，将不同变量（例如table(mf$var3)）的不同类别中的观察数量制成表格，等等。lapply(mf,table)应该按类别将每个变量的观察数量制成表格。

您可能会发现包中的describe函数Hmisc很有用：

Hmisc::describe(group ~ sex + age + var3 + var4, data=d)
于 2017-12-07T13:34:51.217 回答