powered by
Simulate data as in Kong et al. (2004).
sim.kong(gamma, beta, design = "FULL", rho = 0.9, seed = 1, impute = FALSE, ppi)