EmSkewMOD: Calculate modes

In EMMIXskew: The EM Algorithm and Skew Mixture Distribution

Description

Calculate the mode points for each component of skew mixture models.

Usage

EmSkewMOD(p,g,distr,mu,sigma,dof,delta,nrand=10000)

Arguments

p	The dimension of the data
g	The number of components to be fit
distr	A three letter string of distribution id
mu	A numeric matrix with each column corresponding to the mean
sigma	An array of dimension (p,p,g) with first two dimensions corresponding covariance matrix of each component
dof	A vector of degrees of freedom for each component
delta	A matrix with each column as skew parameter vector
nrand	The number of random numbers.

Value

A p by g matrix of mode points for each component of the skew mixture model.

Examples

p=2
g=3

sigma<-array(0,c(2,2,3))
for(h in 2:3) sigma[,,h]<-diag(2)
sigma[,,1]<-cbind( c(1,0),c(0,1))
mu  <- cbind(c(4,-4),c(3.5,4),c( 0, 0))

delta <- cbind(c(3,3),c(1,5),c(-3,1))
dof    <- c(3,5,5)

distr="mst"

EmSkewMOD(p,g,distr,mu,sigma,dof,delta,nrand=10000)

EMMIXskew documentation built on May 2, 2019, 11:07 a.m.