data(leukemia)
group = leukemia$Y # 38 patients status with a value of 1 or 2
x = leukemia$X # 38 x 3051 genes
# apply the function
dist_mah(x, group) # 26.8
# default alpha
d = 3051; n = 38
alpha = sqrt(log(d)/n)
dist_mah(x, group, alpha) # 26.8
Run the code above in your browser using DataLab