R/CalculateGapK.R

In AssocTests: Genetic Association Studies

## Gap Statistics
CalculateGapK <- function(W, WStar, kG, n.monteCarlo)
{
    WLog <- log(W)
    WStarLog <- log(WStar)
    Gap <- apply(WStarLog, 2, mean) - WLog
    s <- apply( WStarLog,2,sd) * sqrt( (n.monteCarlo-1)/n.monteCarlo ) * sqrt(1+1/n.monteCarlo)

    I <- 2;
    gap_kHat <- 1
    for (k in I:(kG-1))
    {
        if ( Gap[k] >= Gap[k+1]-s[k+1] ) 
    	{
    	    gap_kHat <- k
    	    break
    	}
    }

    gap_kHat
}