double* discICAP(uint k, uint noOfSamples, uint noOfFeatures, double **featureMatrix, double *classColumn, double *outputFeatures, double *featureScores) {
uint *intFeatures = (uint *) checkedCalloc(noOfSamples*noOfFeatures,sizeof(uint));
uint *intClass = (uint *) checkedCalloc(noOfSamples,sizeof(uint));
uint *intOutputs = (uint *) checkedCalloc(k,sizeof(uint));
uint **intFeature2D = (uint**) checkedCalloc(noOfFeatures,sizeof(uint*));
int i;
for (i = 0; i < noOfFeatures; i++) {
intFeature2D[i] = intFeatures + i*noOfSamples;
normaliseArray(featureMatrix[i],intFeature2D[i],noOfSamples);
}
normaliseArray(classColumn,intClass,noOfSamples);
ICAP(k, noOfSamples, noOfFeatures, intFeature2D, intClass, intOutputs, featureScores);
for (i = 0; i < k; i++) {
outputFeatures[i] = intOutputs[i];
}
FREE_FUNC(intFeatures);
FREE_FUNC(intClass);
FREE_FUNC(intOutputs);
FREE_FUNC(intFeature2D);
intFeatures = NULL;
intClass = NULL;
intOutputs = NULL;
intFeature2D = NULL;
return outputFeatures;
}/*discICAP(int,int,int,double[][],double[],double[],double[])*/
/******************************************************************************
** entry point for the mex call
** nlhs - number of outputs
** plhs - pointer to array of outputs
** nrhs - number of inputs
** prhs - pointer to array of inputs
******************************************************************************/
void mexFunction(int nlhs, mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
/*************************************************************
** this function takes 4-6 arguments:
** flag = which algorithm to use,
** k = number of features to select,
** featureMatrix[][] = matrix of features,
** classColumn[] = targets,
** optionalParam1 = (path angle or beta value),
** optionalParam2 = (gamma value),
** the arguments should all be discrete integers.
** and has one output:
** selectedFeatures[] of size k
*************************************************************/
int flag, k;
double optionalParam1, optionalParam2;
int numberOfFeatures, numberOfSamples, numberOfTargets;
double *featureMatrix, *targets, *output, *outputFeatures;
double entropyTest;
int i,j;
/************************************************************
** number to function map
** 1 = MIFS
** 2 = mRMR
** 3 = CMIM
** 4 = JMI
** 5 = DISR
** 6 = CIFE
** 7 = ICAP
** 8 = CondRed
** 9 = BetaGamma
** 10 = CMI
*************************************************************/
if (nlhs > 1)
{
printf("Incorrect number of output arguments");
}/*if not 1 output*/
if ((nrhs < 4) || (nrhs > 6))
{
printf("Incorrect number of input arguments");
return;
}/*if not 4-6 inputs*/
/*get the flag for which algorithm*/
flag = (int) mxGetScalar(prhs[0]);
/*get the number of features to select, cast out as it is a double*/
k = (int) mxGetScalar(prhs[1]);
numberOfFeatures = mxGetN(prhs[2]);
numberOfSamples = mxGetM(prhs[2]);
numberOfTargets = mxGetM(prhs[3]);
if (nrhs == 6)
{
optionalParam1 = (double) mxGetScalar(prhs[4]);
optionalParam2 = (double) mxGetScalar(prhs[5]);
}
else if (nrhs == 5)
{
optionalParam1 = (double) mxGetScalar(prhs[4]);
optionalParam2 = 0.0;
}
if (numberOfTargets != numberOfSamples)
{
printf("Number of targets must match number of samples\n");
printf("Number of targets = %d, Number of Samples = %d, Number of Features = %d\n",numberOfTargets,numberOfSamples,numberOfFeatures);
plhs[0] = mxCreateDoubleMatrix(0,0,mxREAL);
return;
}/*if size mismatch*/
else if ((k < 1) || (k > numberOfFeatures))
{
printf("You have requested k = %d features, which is not possible\n",k);
plhs[0] = mxCreateDoubleMatrix(0,0,mxREAL);
return;
}
else
{
featureMatrix = mxGetPr(prhs[2]);
targets = mxGetPr(prhs[3]);
/*double calculateEntropy(double *dataVector, int vectorLength)*/
entropyTest = calculateEntropy(targets,numberOfSamples);
if (entropyTest < 0.0000001)
{
printf("The class label Y has entropy of 0, therefore all mutual informations containing Y will be 0. No feature selection is performed\n");
plhs[0] = mxCreateDoubleMatrix(0,0,mxREAL);
return;
}
else
{
/*printf("Flag = %d, k = %d, numFeatures = %d, numSamples = %d\n",flag,k,numberOfFeatures,numberOfSamples);*/
switch (flag)
{
case 1: /* MIFS */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
if (nrhs == 4)
{
/* MIFS is Beta = 1, Gamma = 0 */
optionalParam1 = 1.0;
optionalParam2 = 0.0;
}
/*void BetaGamma(int k, long noOfSamples, long noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures, double beta, double gamma)*/
BetaGamma(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output,optionalParam1,optionalParam2);
incrementVector(output,k);
break;
}
case 2: /* mRMR */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/*void mRMR_D(int k, int noOfSamples, int noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures)*/
mRMR_D(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output);
incrementVector(output,k);
break;
}
case 3: /* CMIM */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/*void CMIM(int k, int noOfSamples, int noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures)*/
CMIM(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output);
incrementVector(output,k);
break;
}
case 4: /* JMI */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/*void JMI(int k, int noOfSamples, int noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures)*/
JMI(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output);
incrementVector(output,k);
break;
}
case 5: /* DISR */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/*void DISR(int k, int noOfSamples, int noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures)*/
DISR(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output);
incrementVector(output,k);
break;
}
case 6: /* CIFE */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/* CIFE is Beta = 1, Gamma = 1 */
optionalParam1 = 1.0;
optionalParam2 = 1.0;
/*void BetaGamma(int k, long noOfSamples, long noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures, double beta, double gamma)*/
BetaGamma(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output,optionalParam1,optionalParam2);
incrementVector(output,k);
break;
}
case 7: /* ICAP */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/*void ICAP(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output);*/
ICAP(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output);
incrementVector(output,k);
break;
}
case 8: /* CondRed */
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/* CondRed is Beta = 0, Gamma = 1 */
optionalParam1 = 0.0;
optionalParam2 = 1.0;
/*void BetaGamma(int k, long noOfSamples, long noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures, double beta, double gamma)*/
BetaGamma(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output,optionalParam1,optionalParam2);
incrementVector(output,k);
break;
}
case 9: /* BetaGamma */
{
if (nrhs != 6)
{
printf("Insufficient arguments specified for Beta Gamma FS\n");
plhs[0] = mxCreateDoubleMatrix(0,0,mxREAL);
return;
}
else
{
plhs[0] = mxCreateDoubleMatrix(k,1,mxREAL);
output = (double *)mxGetPr(plhs[0]);
/*void BetaGamma(int k, long noOfSamples, long noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures, double beta, double gamma)*/
BetaGamma(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output,optionalParam1,optionalParam2);
incrementVector(output,k);
}
break;
}
case 10: /* CMI */
{
output = (double *)mxCalloc(k,sizeof(double));
/*void CondMI(int k, int noOfSamples, int noOfFeatures,double *featureMatrix, double *classColumn, double *outputFeatures)*/
CondMI(k,numberOfSamples,numberOfFeatures,featureMatrix,targets,output);
i = 0;
while((output[i] != -1) && (i < k))
{
i++;
}
plhs[0] = mxCreateDoubleMatrix(i,1,mxREAL);
outputFeatures = (double *)mxGetPr(plhs[0]);
for (j = 0; j < i; j++)
{
outputFeatures[j] = output[j] + 1; /*C indexes from 0 not 1*/
}/*for number of selected features*/
mxFree(output);
output = NULL;
break;
}
}/*switch on flag*/
return;
}
}
}/*mex function entry*/