2014-09-29 53 views
1

我正在学习如何使用和解释R中包nnetmultinom()的输出。我用于练习的数据集包含一些缺少值的变量。我相信在使用这个函数时,包含缺失值的观察值将被排除在多项式回归之外。我想知道如何获得模型中包含的观测值的数量。如何获取使用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    

由此看来,什么是我需要写分析得到的观测数的命令?有什么建议么?

谢谢

回答

2

fitted()提取模型的拟合值:

nrow(fitted(model))

+0

有时它也是'长度()',因为结果是一个向量。这是例如用于'logitmfx()',大概是从'mfx'包中进行的其他回归。 – MERose 2015-11-26 15:24:36

+0

@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

相关问题