It computes the tetrachoric correlation matrix using the algorithm described in Emrich and Piedmonte (1991). If the resulting matrix is non-positive definite, a nearest positive definite matrix is returned and the warning message will be printed.
compute.sigma.star(prop.vec.bin, corr.mat)
Tetrachoric correlation matrix
Vector of marginal probabilities
Correlation matrix of the binary data
phi2tetra
and nearPD