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)
valML <- performanceIndex(miML, adjMat)
plotPR(valML, col = "red")
Run the code above in your browser using DataLab