simData <- simulatedData(p = 10, n = 20, mu = 100, sigma = 0.25,
ppower = 0.73, noise = FALSE)
counts <- simData$counts
adjMat <- simData$adjMat
miML <- parMIEstimate(counts, method = "ML", unit = "nat", nchips = 2)
miML <- aracne.a(miML, eps = 0.5)
valML <- performanceNET(miML, adjMat)
YI <- YoudenIndex(valML)
Run the code above in your browser using DataLab