Bonferroni correction

//*********************************************************************
// statsBC
//
// Shawn Driscoll
// Salk Institute for Biological Studies
//
// Implements Bonferroni correction. With a Family-wise error rate of
// 0.05, the expected number of false positives will be 0.05.
//
// Parameters:
// w - input wave of p-values
//
// Optional Parameters:
// ncol - column number of p-values if input wave w is a 2D matrix
//*********************************************************************
function statsBC(w,[ncol])

    //-----------------------------------------------------------------
    // parameters
    //-----------------------------------------------------------------
    wave w;
    variable ncol;
   
    //-----------------------------------------------------------------
    // variables
    //-----------------------------------------------------------------
    variable npnts,ndim;
   
    //-----------------------------------------------------------------
    // init variables
    //-----------------------------------------------------------------
   
    if(paramIsDefault(ncol))
        ncol = 0;
    endif
   
    npnts = dimsize(w,0);
    ndim = dimsize(w,1);

    //-----------------------------------------------------------------
    // main function
    //-----------------------------------------------------------------
   
    if(ndim)
        w[][ncol] *= npnts;
    else
        w *= npnts;
    endif

    // set p-values greater than 1 to 1
    w[] = w[p] > 1 ? 1 : w[p];

end