powered by
compInd(S,group=1:ncol(S),alpha=1)
dcovustat
energy::dcov
library(steadyICA) library(combinat) set.seed(100) S = matrix(rnorm(40),ncol=4) group = c(1,2,3,3) compInd(S,group,1)
Run the code above in your browser using DataLab