Update est_pow.R

ttriche · Mar 20, 2023 · f2a2570 · f2a2570
1 parent dfcf190
commit f2a2570
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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.