我正在尝试在MATLAB中实现Polychoric和Polyserial相关Stata解决方案，但我无法完全匹配Stas Kolenikov的数字。 Stas Kolenikov为这些相关矩阵编写了Stata命令。
从我在Stata命令中找到的文档中可以看出，Stas正在使用两步估计，其中使用边际分布计算阈值，然后将其用作已知并且双变量正态分布的MLE最大化相关性。
我使用了Drasgow论文作为我的参考，它使用了其他技术中的两步估计。我很乐意在这里分享我的MATLAB代码。我附上了Drasgow论文供您参考。我的MATLAB代码能够匹配Drasgow数字。
对于Polychoric，我使用第71页的Drasgow中列出的dl / dsigma一阶条件，然后使用fzero MATLAB函数。
对于Polyserial，Drasgow在第69页的方程式1中产生Polyserial与乘积矩相关之间的关系。我使用该方程计算Polyserial相关性。
在这两种情况下，都使用边际分布来识别阈值。
我想知道这里是否有人或Stas本人可以告诉我Stata程序究竟是什么，所以我可以完全匹配数字。如果某些论文中列出了Stas的假设和第一顺序条件，那么如果有人可以链接到它，我真的很感激。我正在复制的Drasgow论文可以找到here.