tests/testthat/test-all.R

In MVNBayesian: Bayesian Analysis Framework for MVN (Mixture) Distribution

ipak <- function(pkg) {
  new.pkg <- pkg[!(pkg %in% installed.packages()[, "Package"])]
  if (length(new.pkg))
    install.packages(new.pkg, dependencies = TRUE)
  sapply(pkg, require, character.only = TRUE)
}

#libraries
packages <- c("mvtnorm","stats","plyr","Rfast","rgl")
ipak(packages)

# Test for MND framework:
# step.1
head(dataset1)
# step.2
BP <- MND_BayesianPosteriori(dataset1)
# step.3
BP_Gibbs <- MND_GibbsSampler(5000, BP)
# step.4
result_MND <- MND_MCMC(BP, 5000, c(1), c(76.53))

# Test for MND_BayesianIteratior:
result_Iter <- MND_BayesianIterator(dataset1,c(80,16,3))

# Test for MixMND framework:
# step.1
head(dataset2)
dataset2_test <- dataset2[,1:4]
# step.2
MixBP <- MixMND_BayesianPosteriori(dataset2_test,3)
# step.3
MixBP_Gibbs <- MixMND_GibbsSampler(8000, MixBP, random_method = "Fast")
# step.4
result_MixMND <- MixMND_MCMC(MixBP, 8000, c(1), c(1))