Computes polychoric and polyserial correlations by quick "two-step" methods or
ML, optionally with standard errors; tetrachoric and biserial correlations are
special cases.

See also: https://polycor.r-forge.r-project.org/