我使用Stata vers。 14.2并且想要知道如何使用polychoric关联来测量在5点序数上测量的26个项目，并使用mi impute chained添加16个项目。 polychoric是一个用户编写的命令，可以运行findit polychoric
我曾尝试使用Steve Samuels的代码：Polychoric correlation (Stata) using multiple imputations and a complex sample design
但是，这种方法不起作用，因为Stata报告没有观察结果。
怎么了？
local pcvars IC08Q01 IC08Q02 IC08Q03 IC08Q04 IC08Q05 IC08Q06 IC08Q07 IC08Q08 IC08Q09 IC09Q01 IC09Q02 IC09Q03 IC09Q04 IC09Q05 IC09Q06 IC09Q07 IC10Q01 IC10Q02 IC10Q03 IC10Q04 IC10Q05 IC10Q06 IC10Q07 IC10Q08 IC10Q09
local pwtvar W_FSTUWT
mata:
stata("qui mi query")
M= st_numscalar("r(M)")
vnames = tokens(st_local("pcvars"))
nv = cols(vnames)
sum_r = J(nv,nv,0)
sum_atr = J(nv,nv,0)
for (j = 1; j<=M; j++) {
st_numscalar("j",j)
stata("qui polychoric `pcvars' [pw = `pwtvar'] if _mi_m==j")
r = st_matrix("r(R)")
sum_r = sum_r + r
sum_atr = sum_atr + atanh(r)
}
.no observations
.no observations
.no observations
.no observations
sum IC08Q01 IC08Q02 IC08Q03 IC08Q04 IC08Q05 IC08Q06 IC08Q07 IC08Q08 IC08Q09
IC09Q01 IC09Q02 IC09Q03 IC09Q04 IC09Q05 IC09Q06 IC09Q07 IC10Q01 IC10Q02
IC10Q03IC10Q04 IC10Q05 IC10Q06 IC10Q07 IC10Q08 IC10Q09
Variable | Obs Mean Std. Dev. Min Max
-------------+---------------------------------------------------------
IC08Q01 | 6,524 2.132894 1.37049 1 5
IC08Q02 | 6,538 2.558581 1.4809 1 5
IC08Q03 | 6,525 3.130268 1.256819 1 5
IC08Q04 | 6,515 3.726477 1.373942 1 5
IC08Q05 | 6,516 4.394721 1.098578 1 5
-------------+---------------------------------------------------------
IC08Q06 | 6,514 4.289377 .9817563 1 5
IC08Q07 | 6,506 3.418229 1.281393 1 5
IC08Q08 | 6,500 3.277077 1.206995 1 5
IC08Q09 | 6,499 3.255732 1.372518 1 5
IC09Q01 | 6,423 3.078468 1.048085 1 5
-------------+---------------------------------------------------------
IC09Q02 | 6,415 2.209509 1.178517 1 5
IC09Q03 | 6,379 2.002038 1.012696 1 5
IC09Q04 | 6,394 2.207226 1.22071 1 5
IC09Q05 | 6,408 2.87937 1.355711 1 5
IC09Q06 | 6,380 3.486364 1.075735 1 5
-------------+---------------------------------------------------------
IC09Q07 | 6,392 2.193367 1.177622 1 5
IC10Q01 | 6,375 2.795294 1.482978 1 5
IC10Q02 | 6,341 2.163381 1.26046 1 5
IC10Q03 | 6,343 3.271007 1.060689 1 5
IC10Q04 | 6,348 2.254411 1.275265 1 5
-------------+---------------------------------------------------------
IC10Q05 | 6,354 1.919893 1.185344 1 5
IC10Q06 | 6,310 1.605071 1.066042 1 5
IC10Q07 | 6,333 2.044687 1.168547 1 5
IC10Q08 | 6,305 2.843616 1.266879 1 5
IC10Q09 | 6,326 2.486405 1.222085 1 5
插补
sum _1_IC08Q01-_16_IC10Q09
Variable | Obs Mean Std. Dev. Min Max
-------------+---------------------------------------------------------
_1_IC08Q01 | 7,315 2.152153 1.381924 1 5
_2_IC08Q01 | 7,315 2.154067 1.38092 1 5
_3_IC08Q01 | 7,315 2.1635 1.387099 1 5
_4_IC08Q01 | 7,315 2.15188 1.379875 1 5
_5_IC08Q01 | 7,315 2.152563 1.383807 1 5
-------------+---------------------------------------------------------
_6_IC08Q01 | 7,315 2.163363 1.386474 1 5
_7_IC08Q01 | 7,315 2.155571 1.382484 1 5
_8_IC08Q01 | 7,315 2.155707 1.384 1 5
_9_IC08Q01 | 7,315 2.148462 1.378017 1 5
_10_IC08Q01 | 7,315 2.159672 1.386756 1 5
-------------+---------------------------------------------------------
_11_IC08Q01 | 7,315 2.150239 1.380847 1 5
_12_IC08Q01 | 7,315 2.152153 1.380836 1 5
_13_IC08Q01 | 7,315 2.144771 1.377367 1 5
_14_IC08Q01 | 7,315 2.158578 1.384514 1 5
_15_IC08Q01 | 7,315 2.148735 1.380663 1 5
-------------+---------------------------------------------------------
_16_IC08Q01 | 7,315 2.155981 1.380309 1 5
_1_IC08Q02 | 7,318 2.570101 1.484904 1 5
_2_IC08Q02 | 7,318 2.578163 1.491298 1 5
_3_IC08Q02 | 7,318 2.58199 1.485583 1 5
_4_IC08Q02 | 7,318 2.570101 1.48343 1 5
-------------+---------------------------------------------------------
_5_IC08Q02 | 7,318 2.579667 1.489109 1 5
_6_IC08Q02 | 7,318 2.580623 1.492908 1 5
_7_IC08Q02 | 7,318 2.5783 1.485138 1 5
_8_IC08Q02 | 7,318 2.578573 1.491185 1 5
_9_IC08Q02 | 7,318 2.570374 1.484246 1 5
-------------+---------------------------------------------------------
_10_IC08Q02 | 7,318 2.589779 1.485317 1 5
_11_IC08Q02 | 7,318 2.573791 1.490146 1 5
_12_IC08Q02 | 7,318 2.574201 1.485073 1 5
_13_IC08Q02 | 7,318 2.567095 1.481265 1 5
_14_IC08Q02 | 7,318 2.5783 1.487621 1 5
-------------+---------------------------------------------------------
_15_IC08Q02 | 7,318 2.564635 1.488093 1 5
_16_IC08Q02 | 7,318 2.581033 1.486004 1 5
_1_IC08Q03 | 7,319 3.134718 1.258433 1 5
_2_IC08Q03 | 7,319 3.125973 1.253574 1 5
_3_IC08Q03 | 7,319 3.135128 1.260504 1 5
-------------+---------------------------------------------------------
_4_IC08Q03 | 7,319 3.142369 1.261495 1 5
_5_IC08Q03 | 7,319 3.146878 1.260816 1 5
_6_IC08Q03 | 7,319 3.13909 1.260453 1 5
_7_IC08Q03 | 7,319 3.128979 1.258708 1 5
_8_IC08Q03 | 7,319 3.140456 1.2565 1 5
-------------+---------------------------------------------------------
_9_IC08Q03 | 7,319 3.123104 1.260218 1 5
_10_IC08Q03 | 7,319 3.146331 1.263153 1 5
_11_IC08Q03 | 7,319 3.14073 1.260054 1 5
_12_IC08Q03 | 7,319 3.127886 1.259362 1 5
_13_IC08Q03 | 7,319 3.135674 1.257515 1 5
-------------+---------------------------------------------------------
_14_IC08Q03 | 7,319 3.131439 1.262465 1 5
_15_IC08Q03 | 7,319 3.131712 1.256795 1 5
_16_IC08Q03 | 7,319 3.146058 1.259827 1 5
....
-------------+---------------------------------------------------------
_1_IC10Q09 | 7,306 2.506844 1.228242 1 5
_2_IC10Q09 | 7,306 2.496441 1.225131 1 5
_3_IC10Q09 | 7,306 2.504654 1.22736 1 5
_4_IC10Q09 | 7,306 2.507528 1.228683 1 5
_5_IC10Q09 | 7,306 2.497673 1.227701 1 5
-------------+---------------------------------------------------------
_6_IC10Q09 | 7,306 2.505886 1.226908 1 5
_7_IC10Q09 | 7,306 2.492883 1.231913 1 5
_8_IC10Q09 | 7,306 2.51095 1.233217 1 5
_9_IC10Q09 | 7,306 2.499316 1.231711 1 5
_10_IC10Q09 | 7,306 2.495209 1.225797 1 5
-------------+---------------------------------------------------------
_11_IC10Q09 | 7,306 2.514782 1.232844 1 5
_12_IC10Q09 | 7,306 2.495209 1.229365 1 5
_13_IC10Q09 | 7,306 2.506981 1.223663 1 5
_14_IC10Q09 | 7,306 2.518478 1.230571 1 5
_15_IC10Q09 | 7,306 2.498221 1.230375 1 5
-------------+---------------------------------------------------------
_16_IC10Q09 | 7,306 2.507802 1.23024 1 5