R/rMSCEC.R

In skewMLRM: Estimation for Scale-Shape Mixtures of Skew-Normal Distributions

rMSCEC <-
function(n, mu, Sigma, lambda, nu=0.5, gamma=0.5){
if(is.null(n))
    stop("sample size must be specified")
  if(is.null(mu))
    stop("mu must be specified")
  if(is.null(Sigma))
    stop("Sigma must be specified")
    if(is.vector(Sigma)) Sigma=matrix(Sigma)
  if(is.null(lambda))
    stop("lambda must be specified")
  if(round(n)!=n | n<=0)
    stop("sample size must be a positive integer")      
  if(!is.matrix(mu)) mu<-matrix(mu, ncol=length(mu),nrow=1)
  if (!isSymmetric(Sigma))
        stop("Sigma is not symmetric")
  if (min(eigen(Sigma)$values)<=0)
        stop("Sigma is not a positive-definite matrix")
  if (ncol(Sigma) != ncol(mu))
        stop("The dimension of mu does not agree with the dimension of Sigma")
  if (ncol(Sigma) != length(lambda))
        stop("The dimension of lambda does not agree with the dimension of Sigma")
  if (ncol(mu) != length(lambda))
        stop("The dimension of mu does not agree with the dimension of lambda")
  if (nu<0 | nu>1)
        stop("nu must be between 0 and 1")
  if (gamma<0 | gamma>1)
        stop("gamma must be between 0 and 1")
  p<-ncol(Sigma)
    lambda=as.vector(lambda)
    delta <- lambda/c(sqrt(1 + t(lambda)%*%lambda))
  y <- matrix(0,n,p)
    for (i in 1:n)
   {
        uu<-rbinom(1,1,nu)
        u<-gamma*uu+1-uu
        Z<-rMSNC(1, mu=rep(0,p), Sigma=diag(p), lambda=lambda)
        y[i,]=c(mu)+u^(-1/2)*matrix.sqrt(Sigma)%*%t(Z) 
  }
    return(y)
  }