CorrToolBox-package {CorrToolBox}R Documentation

Modeling Correlational Magnitude Transformations in Discretization Contexts


This package implements the computational algorithms for modeling the correlation transitions under specified distributional assumptions within the realm of discretization in the context of the latency and threshold concepts. Functions that compute the correlational magnitude changes in both directions (identification of the pre-discretization correlation value in order to attain a specified post-discretization magnitude, and the other way around) are provided.

This package consists of eight main functions. Computing the tetrachoric correlation from the phi coefficient and vice versa are done in phi2tet and tet2phi, respectively. Computing the polychoric correlation from the ordinal phi coefficient and vice versa are done in ophi2poly and poly2ophi, respectively. Computing the biserial correlation from the point-biserial correlation and vice versa are done in pbs2bs and bs2pbs, respectively. Computing the polyserial correlation from the point-polyserial correlation and vice versa are done in pps2ps and ps2pps, respectively.

Auxiliary functions are also provided. corrY2corrZ, corrZ2corrY, corrZ2ophi, corrZ2phi, and ophi2corrZ are intermediate functions utilized within the main functions but can be used as stand-alone functions. ordY discretizes a continuous variable, and mps2cps provides cumulative probabilities for each set of marginal probabilities in a list. Additional intermediate functions from imported packages include phi2tetra from the psych package, ordcont and contord from the GenOrd package, skewness and kurtosis from the moments package, validation.skewness.kurtosis from the BinNonNor package, and pmvnorm from the mvtnorm package.

Within each correlation transition function, the correlation boundaries for the given marginal distributions are compared to the specified input correlation to ensure there are no violations according to Demirtas and Hedeker (2011). The function valid.limits.BinOrdNN in the package BinOrdNonNor is utilized for this step. Additionally, Fleishman.coef.NN in the package BinOrdNonNor is used wherever Fleishman coefficients need to be calculated for a continuous variable.


Package: CorrToolBox
Type: Package
Version: 1.6.4
Date: 2022-02-21
License: GPL-2 | GPL-3


Rawan Allozi, Hakan Demirtas, Ran Gao

Maintainer: Ran Gao <>


Demirtas, H. (2016). A note on the relationship between the phi coefficient and the tetrachoric correlation under nonnormal underlying distributions. The American Statistician, 70(2), 143-148.

Demirtas, H., Ahmadian, R., Atis, S., Can, F.E., and Ercan, I. (2016). A nonnormal look at polychoric correlations: modeling the change in correlations before and after discretization. Computational Statistics, 31(4), 1385-1401.

Demirtas, H. and Hedeker, D. (2011). A practical way for computing approximate lower and upper correlation bounds. The American Statistician, 65(2), 104-109.

Demirtas, H. and Hedeker, D. (2016). Computing the point-biserial correlation under any underlying continuous distribution. Communications in Statistics-Simulation and Computation, 45(8), 2744-2751.

Demirtas, H., Hedeker, D., and Mermelstein, R. J. (2012). Simulation of massive public health data by power polynomials. Statistics in Medicine, 31(27), 3337-3346.

Demirtas, H. and Vardar-Acar, C. (2017). Anatomy of correlational magnitude transformations in latency and discretization contexts in Monte-Carlo studies. In ICSA Book Series in Statistics, John Dean Chen and Ding-Geng (Din) Chen (Eds): Monte-Carlo Simulation-Based Statistical Modeling. Singapore: Springer, 59-84.

Ferrari, P.A. and Barbiero, A. (2012). Simulating ordinal data. Multivariate Behavioral Research, 47(4), 566-589.

Fleishman A.I. (1978). A method for simulating non-normal distributions. Psychometrika, 43(4), 521-532.

Vale, C.D. and Maurelli, V.A. (1983). Simulating multivariate nonnormal distributions. Psychometrika, 48(3), 465-471.

[Package CorrToolBox version 1.6.4 Index]