rho_M1M2: Approximate Correlation between Two Continuous Mixture...

In AFialkowski/SimCorrMix: Simulation of Correlated Data with Multiple Variable Types Including Continuous and Count Mixture Distributions

Description

This function approximates the expected correlation between two continuous mixture variables M1 and M2 based on their mixing proportions, component means, component standard deviations, and correlations between components across variables. The equations can be found in the Expected Cumulants and Correlations for Continuous Mixture Variables vignette. This function can be used to see what combination of component correlations gives a desired correlation between M1 and M2.

Usage

rho_M1M2(mix_pis = list(), mix_mus = list(), mix_sigmas = list(),
  p_M1M2 = NULL)

Arguments

mix_pis	a list of length 2 with 1st component a vector of mixing probabilities that sum to 1 for component distributions of M1 and likewise for 2nd component and M2
mix_mus	a list of length 2 with 1st component a vector of means for component distributions of M1 and likewise for 2nd component and M2
mix_sigmas	a list of length 2 with 1st component a vector of standard deviations for component distributions of M1 and likewise for 2nd component and M2
p_M1M2	a matrix of correlations with rows corresponding to M1 and columns corresponding to M2; i.e., p_M1M2[1, 2] is the correlation between the 1st component of M1 and the 2nd component of M2

Value

the expected correlation between M1 and M2

References

Davenport JW, Bezder JC, & Hathaway RJ (1988). Parameter Estimation for Finite Mixture Distributions. Computers & Mathematics with Applications, 15(10):819-28.

Pearson RK (2011). Exploring Data in Engineering, the Sciences, and Medicine. In. New York: Oxford University Press.

See Also

rho_M1Y

Examples

# M1 is mixture of N(-2, 1) and N(2, 1);
# M2 is mixture of Logistic(0, 1), Chisq(4), and Beta(4, 1.5)
# pairwise correlation between components across M1 and M2 set to 0.35
L <- calc_theory("Logistic", c(0, 1))
C <- calc_theory("Chisq", 4)
B <- calc_theory("Beta", c(4, 1.5))
rho_M1M2(mix_pis = list(c(0.4, 0.6), c(0.3, 0.2, 0.5)),
  mix_mus = list(c(-2, 2), c(L[1], C[1], B[1])),
  mix_sigmas = list(c(1, 1), c(L[2], C[2], B[2])),
  p_M1M2 = matrix(0.35, 2, 3))

AFialkowski/SimCorrMix documentation built on May 30, 2019, 3:47 p.m.