Как получить количество наблюдений, включенных в модель, созданную с помощью функции multinom в R?

Question

Я учусь использовать и интерпретировать вывод multinom() пакета nnetв R. Набор данных, который я использую для упражнения, содержит переменные с некоторыми пропущенными значениями. Я полагаю, что при использовании этой функции наблюдения, содержащие пропущенные значения, исключаются из полиномиальной регрессии. Я хотел бы знать, как я могу получить количество наблюдений, которые были включены в модель.

Подобный вопрос был задан около года назад ( Как изменить вывод после запуска multinom() в 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            

Исходя из этого, какую команду мне нужно написать, чтобы проанализировать количество наблюдений? Какие-либо предложения?

Спасибо

Answer 1

fitted() извлекает установленные значения модели:

nrow(fitted(model))