1
我正在学习如何使用和解释R中包nnet
的multinom()
的输出。我用于练习的数据集包含一些缺少值的变量。我相信在使用这个函数时,包含缺失值的观察值将被排除在多项式回归之外。我想知道如何获得模型中包含的观测值的数量。如何获取使用R中的函数multinom创建的模型中包含的观察值的数量?
一年前有人问过类似的问题(How to change the output after running multinom() in R),但问题的这个特定方面没有回答。我希望这次有人会有答案。
我的数据集包含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
由此看来,什么是我需要写分析得到的观测数的命令?有什么建议么?
谢谢
有时它也是'长度()',因为结果是一个向量。这是例如用于'logitmfx()',大概是从'mfx'包中进行的其他回归。 – MERose 2015-11-26 15:24:36
@MERose你可以使用'length()'来获得'mfx'包中'logitmfx()'或其他模型的观察次数吗?如果我拟合一个称为模型的模型:'model < - logitmfx(dv〜iv1 + iv2,data = df,robust = FALSE,clustervar1 =“iv3”,clustervar2 =“iv4”)',则输出不会提供N和“长度(模型)”不提供正确的N。 – coip 2017-07-06 17:52:59