R/validate.target.cormat.R

In OrdNor: Concurrent Generation of Ordinal and Normal Data with Given Correlation Matrix and Marginal Distributions

validate.target.cormat <-
function(plist, CorrMat, no.ord, no.norm){

if (is.positive.definite(CorrMat) == FALSE) {
stop("Specified correlation matrix is not positive definite! \n")
}
if (isSymmetric(CorrMat) == FALSE) {
stop("Specified correlation matrix is not symmetric! \n")
}

if ( ncol(CorrMat) != (no.ord + no.norm) ){
stop("Dimension of correlation matrix does not match the number of variables!\n")
} 

Limits= valid.limits(plist, no.ord, no.norm) ## plist is validated here.

minmat = Limits$lower
maxmat = Limits$upper
rangemat = minmat<=CorrMat & CorrMat<=maxmat

if (sum(!rangemat)>0) 
{

d<-ncol(rangemat)
cat("Target matrix is not valid. The following values are invalid.\n")
for (i in 1:d){
for (j in i:d){ 
if(rangemat[i,j]==FALSE){ cat("Corr[",i,",",j,"] must be between", round(minmat[i,j],3), "and",round(maxmat[i,j],3),"\n")}
}
}
stop ("Range violation occurred in the target correlation matrix.\n")
}
return(TRUE)

}