# NOT RUN {
if(require(mvtnorm)){
## Generate two random samples of size 50 from two multivariate normal distributions
# sample size
n<-50
# true parameters of class 1 and class 2
param.class1<-simulation$condition1
param.class2<-simulation$condition2$`5`$`2`
# simulated dataset
data.class1<-rmvnorm(n = n,mean =param.class1$mu ,sigma =param.class1$S)
data.class2<-rmvnorm(n = n,mean =param.class2$mu ,sigma=param.class2$S)
data<-rbind(data.class1,data.class2)
classes<-c(rep(1,nrow(data.class1)),rep(2,nrow(data.class2)))
s<-cov(data)
s1<-cov(data.class1)
s2<-cov(data.class2)
testMeanVariance(S = s,S1 =s1, S2 = s2, n1 = n, n2 = n)
## equivalently...
# estimated parameters: maximum likelihood estimate
est.param<-parameters(data = data,classes =classes ,shrink = FALSE)
testMeanVariance(est.param$S,est.param$S1,est.param$S2,est.param$n1,est.param$n2)
}
# }
Run the code above in your browser using DataLab