Nothing
#######################################################################################
###############
############### AntMAN Package : Tests and Examples
###############
###############
#######################################################################################
##############################################
### Load the AntMan package
##############################################
library("AntMAN")
set.seed(123)
##############################################
### Load the IRIS package
##############################################
data("iris")
mat <- matrix(0, 150, 4)
mat[,1] <- iris$Sepal.Length
mat[,2] <- iris$Sepal.Width
mat[,3] <- iris$Petal.Length
mat[,4] <- iris$Petal.Width
mixture_mvn_params <- AM_mix_hyperparams_multinorm(mu0 = rep(0, 4),
ka0 = 0.4,
nu0 = 40,
Lam0=1*diag(4) )
mcmc_params<- AM_mcmc_parameters(niter=10000,
burnin=2000,
thin=1,
verbose=1)
components_prior <- AM_mix_components_prior_pois (a=10,b=2,init = 1)
weights_prior <- AM_mix_weights_prior_gamma(a=1, b=1, init= 1)
fit <- AM_mcmc_fit(
y = mat,
mix_kernel_hyperparams = mixture_mvn_params,
mix_components_prior =components_prior,
mix_weight_prior = weights_prior,
mcmc_parameters = mcmc_params)
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.