From f2a2570adadd2a0ba360fed8387a92a6ebe2f349 Mon Sep 17 00:00:00 2001
From: madajz <31075974+madajz@users.noreply.github.com>
Date: Mon, 20 Mar 2023 15:21:03 -0400
Subject: [PATCH] Update est_pow.R

---
 R/est_pow.R | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)

diff --git a/R/est_pow.R b/R/est_pow.R
index 50e953b..561b6cb 100644
--- a/R/est_pow.R
+++ b/R/est_pow.R
@@ -59,8 +59,8 @@ est_pow = function(n,alpha,nsim,dist,params,tests){
   
   if("mclust" %in% tests){
     pwr.df = rbind(pwr.df,data.frame(N = n, Test = "Mclust",
-                                     power = sum(sapply(n.dfs, function(s) I(mclust::mclustBootstrapLRT(as.data.frame(s),modelName="V",verbose=F,maxG=1)$p.value<alpha)))/nsim,
-                                     FP = sum(sapply(a.dfs, function(s) I(mclust::mclustBootstrapLRT(as.data.frame(s),modelName="V",verbose=F,maxG=1)$p.value<alpha)))/nsim ))
+                                     power = sum(sapply(n.dfs, function(s) I(mclust::mclustBootstrapLRT(as.data.frame(s),modelName="E",verbose=F,maxG=1)$p.value<alpha)))/nsim,
+                                     FP = sum(sapply(a.dfs, function(s) I(mclust::mclustBootstrapLRT(as.data.frame(s),modelName="E",verbose=F,maxG=1)$p.value<alpha)))/nsim ))
   }
   
   ## sigclust is giving me issues so I've dropped it for now.