full2meta: Convert between parameter formats

Description

These functions converts the parameters between the general Gaussian mixture (copula) model and the special GMCM. Most functions of the GMCM packages use the theta format described in rtheta.

Usage

full2meta(theta)
meta2full(par, d)

Arguments

theta

A list of parameters for the full model. Formatted as described in rtheta.

par

A vector of length 4 where par[1] is the probability of coming from the first component, par[2] is the mean value, par[3] is the standard deviation, and par[4] is the correlation of the reproducible comp

The dimension of the mixture distribution.

Value

full2meta returns a numeric vector of length 4 formatted as par.
meta2full returns a formatted list of parameters as described by rtheta.

Details

If a theta is supplied which is not on the form of Li et. al. (2011) the output is coerced by simply picking the first elements of the first mean vector and first covariance matrix as mean and standard deviation, respectively.

References

Li, Q., Brown, J. B. J. B., Huang, H., & Bickel, P. J. (2011). Measuring reproducibility of high-throughput experiments. The Annals of Applied Statistics, 5(3), 1752-1779. doi:10.1214/11-AOAS466

Tewari, A., Giering, M., & Raghunathan, A. (2011). Parametric Characterization of Multimodal Distributions with Non-gaussian Modes. IEEE 11th International Conference on Data Mining Workshops, 2011, 286-292. doi:10.1109/ICDMW.2011.135

Examples

Run this code

theta <- GMCM:::rtheta(m = 2, d = 2)
print(par <- full2meta(theta))
print(theta.special.case <- meta2full(par, d = 2))

Run the code above in your browser using DataLab