changed default initializtion for mixSQP

stephens999 · Feb 18, 2020 · 082b660 · 082b660
1 parent fe9680c
commit 082b660
Showing 1 changed file with 4 additions and 0 deletions.
diff --git a/R/ash.R b/R/ash.R
@@ -369,6 +369,10 @@ ash.workhorse <-
     k = length(mixsd)
     prior = setprior(prior, k, nullweight, null.comp)
     pi = initpi(k, length(data$x), null.comp)
+    if(optmethod=="mixSQP"){ # we found that using a constant initialization for mixSQP works better than
+      #initpi, which was aimed at initializing EM algorithm
+      pi = rep(1,k)
+    }
 
     if (mixcompdist == "normal")
       g = normalmix(pi, rep(mode, k), mixsd)