Description: Simulate joint distribution given nonparametric marginals and their covariance structure characterized by a correlation matrix. The simulator engages the problem from a purely computational perspective. It assumes no statistical models such as copulas and parametric distributions, and can approximate the target correlations regardless of theoretical feasibility. The algorithm integrates and advances the Iman-Conover (1982) approach (<doi:10.1080/03610918208812265>) and the Ruscio-Kaczetow (2008) iteration (<doi:10.1080/00273170802285693>). The package is implemented in C++ and carefully designed toward computing speed, suitable for large input in a manycore environment. Precision of the approximation and computing speed both outperform various CRAN packages to date by a substantial margin. Benchmarks are detailed in function examples.

Description: Simulate multivariate correlated data given nonparametric marginal distributions and their covariance structure characterized by a correlation matrix. The simulator engages the problem from a purely computational perspective. It assumes no statistical models such as copulas or parametric distributions, and can approximate the target correlations regardless of theoretical feasibility. The algorithm integrates and advances the Iman-Conover approach (1982) <doi:10.1080/03610918208812265> and the Ruscio-Kaczetow iteration (2008) <doi:10.1080/00273170802285693>. Package functions are carefully implemented in C++ for pushing computing speed, suitable for large input in a manycore environment. Precision of the approximation and computing speed both outperform various CRAN packages to date by a substantial margin. Benchmarks are detailed in function examples. Additionally, a simple heuristic algorithm initially designed to optimize the simulated joint distribution demonstrated not only strong error reduction capability, but also the potential of achieving the same level of precision of approximation without the enhanced Iman-Conover-Ruscio-Kaczetow as a primer, especially when the functional relationships between marginal distributions are highly nonlinear.