# GMM model # // library(mvtnorm) ) n1 = n2 = mu1 = c(,) mu2 = c(-,-) sigma1 = matrix(c(,.,.,),nrow=) sigma2 = matrix(c(,.,.,),nrow=) y1 = rep(,n1) y2 = rep(,n2) x1 = rmvnorm(n1, mean=mu1, sigma=sigma1) x2 = rmvnorm(n2, mean=mu2, sigma=sigma2) x = rb…