EMMIXskew (version 1.0.3)

EmSkewMOD: Calculate modes

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

# NOT RUN {
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)

# }