R: Удаление переменной из регрессии из одной скрытой переменной вводит значения NA в другие скрытые переменные - PullRequest
Новая
       11

R: Удаление переменной из регрессии из одной скрытой переменной вводит значения NA в другие скрытые переменные

0 голосов
/ 19 сентября 2019

Для анализа набора данных мне передали код, написанный кем-то еще в R, который использует моделирование структурных уравнений для анализа скрытых переменных набора данных.

При передаче инструкций по привязке переменных к скрытым переменным для регрессии (к пакету Lavaan в R) мы сначала инструктируем Lavaan регрессировать двоичную переменную «пол» (мужской или женский) с определенной скрытой переменной.При повторном запуске программы без передачи «sex» мы вводим новые значения NA для стандартных ошибок в других скрытых переменных.

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

Я проверил документацию по Lavaan, но, возможно, я упустил любые решения моей проблемы.

Вот строка, в которой мы передаем отношения регрессии 

'''r
# regressions  
ab_com ~ word_id_time + word_id_acc + naming_acc + naming_speed + phon 
+age.tidy+sex  
word_id_time ~ naming_speed  
word_id_acc ~ naming_acc + phon  
'''

Мое изменение просто удаляет "секс"

Вот результаты, когда "секс" включен: 

               Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
  ab_com ~                                                              
    word_id_time      0.015    0.048    0.320    0.749    0.017    0.018
    word_id_acc      -0.496    8.080   -0.061    0.951   -6.295   -6.547
    naming_acc        1.775   58.464    0.030    0.976    1.775    1.846
    naming_speed     -0.064    0.071   -0.898    0.369   -0.064   -0.066
    phon              5.671  161.071    0.035    0.972    5.671    5.899
    age.tidy         -0.100    0.051   -1.953    0.051   -0.100   -0.092
    sex               0.101    0.089    1.132    0.258    0.101    0.052
  word_id_time ~                                                        
    naming_speed      0.473    0.065    7.292    0.000    0.428    0.428
  word_id_acc ~                                                         
    naming_acc        3.832   56.202    0.068    0.946    0.302    0.302
    phon             10.451  155.749    0.067    0.947    0.824    0.824

Вот результаты, когда "секс" удален: 

               Estimate  Std.Err  z-value  P(>|z|)   Std.lv  Std.all
  ab_com ~                                                              
    word_id_time      0.015    0.048    0.301    0.764    0.016    0.017
    word_id_acc      -0.683       NA                    -10.601  -11.010
    naming_acc        3.078       NA                      3.078    3.197
    naming_speed     -0.060    0.071   -0.853    0.394   -0.060   -0.063
    phon              9.234       NA                      9.234    9.591
    age.tidy         -0.098    0.051   -1.923    0.054   -0.098   -0.090
  word_id_time ~                                                        
    naming_speed      0.473    0.065    7.293    0.000    0.428    0.428
  word_id_acc ~                                                         
    naming_acc        4.703       NA                      0.303    0.303
    phon             12.793       NA                      0.824    0.824
...