data(colon)
normal.group <- colon[,1:40]
colon.group <- colon[,41:62]
Sigma.normal <- shrinkcovmat.identity(normal.group)
Sigma.normal
Sigma.colon <- shrinkcovmat.identity(colon.group)
Sigma.colonRun the code above in your browser using DataLab