R/identifyTheEquivalent.ma.R

In MXM: Feature Selection (Including Multiple Solutions) and Bayesian Networks

identifyTheEquivalent.ma = function(queues, target, dataset, cvar, z, test, threshold, statistic, 
                                    univariateModels, pvalues, hash, stat_hash, pvalue_hash) {
  z = t(z);
  #case 3
  #if we have more than one equivalent vars in z , we select the first one
  #for loop if we dont use the below lapply function
  #equalsY = 0;
  for ( i in 1:ncol(z) )  {
    w = z[,i];
    w = t( t(w) );
    zPrime = c(setdiff(z , w) , cvar);
    cur_results = test(target = target, dataset = dataset, xIndex = w, csIndex = zPrime, statistic = statistic, univariateModels = univariateModels, hash = hash, stat_hash = stat_hash, pvalue_hash = pvalue_hash);
    
    if ( cur_results$pvalue > threshold ) {  
      queues[[w]] = as.matrix( c(queues[[w]] , queues[[cvar]]) );
      break;
      #equalsY = equalsY+1;
    }
  }
  #cat("\nEquals Ys in the current z subset = %d",equalsY);
  return(queues);
}