1

我正在学习如何使用和解释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            

由此,我需要编写什么命令来获取分析的观察次数?有什么建议么?

谢谢

4

1 回答 1

2

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

nrow(fitted(model))

于 2014-09-30T02:00:58.717 回答