我正在学习如何使用和解释Rmultinom()
中包的输出nnet
。我用于练习的数据集包含一些缺失值的变量。我相信在使用这个函数时,包含缺失值的观测值会被排除在多项式回归之外。我想知道如何获得模型中包含的观察数。
大约一年前提出了一个类似的问题(如何在 R 中运行 multinom() 后更改输出),但问题的这个特定方面没有得到回答。我希望这次有人会给出答案。
我的数据集包含 42 个变量和 593 个观察值。为了说明我的问题,我运行了一个模型,其中包含感兴趣的结果(4 个组)和 3 个协变量(2 个带有一些缺失值)。
输出是:
model=multinom(GROUP ~ female + BMI + born_Canada, data)
# weights: 20 (12 variable)
initial value 679.284237
iter 10 value 637.285180
final value 637.129747
converged
summary(model)
Call:
multinom(formula = GROUP ~ female + BMI + born_Canada, data = data)
Coefficients:
(Intercept) female BMI born_Canada
2 0.2827215 -0.1446951 0.007395264 0.1870738
3 0.5513568 -0.3101297 -0.007632127 0.1948097
4 -0.3698947 -0.5663595 -0.024331298 0.6115428
Std. Errors:
(Intercept) female BMI born_Canada
2 0.6753277 0.2694006 0.0197283 0.4484126
3 0.7161380 0.2806695 0.0215023 0.4632981
4 1.0415322 0.3829313 0.0312360 0.6920262
Residual Deviance: 1274.259
AIC: 1298.259
由此,我需要编写什么命令来获取分析的观察次数?有什么建议么?
谢谢