Make it possible to disable the MPLE existence check.

references #570

DESCRIPTION

@@ -1,5 +1,5 @@
 Package: ergm
-Version: 4.7-7393
+Version: 4.7-7394
 Date: 2024-09-29
 Title: Fit, Simulate and Diagnose Exponential-Family Models for Networks
 Authors@R: c(

R/control.ergm.R

@@ -165,6 +165,9 @@
 #'   simulate, the MCMC burn-in, and the MCMC interval for `Godambe`
 #'   and `bootstrap` methods.
 #'
+#' @param MPLE.check If `TRUE` (the default), perform the MPLE
+#'   existence check described by \insertCite{ScHu23c;textual}{ergm}.
+#'
 #' @param MPLE.constraints.ignore If `TRUE`, MPLE will ignore all
 #'   dyad-independent constraints except for those due to attributes
 #'   missingness. This can be used to avert evaluating and storing the
@@ -523,7 +526,8 @@ control.ergm<-function(drop=TRUE,
                        MPLE.covariance.method ="invHess",
                        MPLE.covariance.sim.burnin = 1024,
                        MPLE.covariance.sim.interval = 1024,
-                       MPLE.constraints.ignore=FALSE,
+                       MPLE.check = TRUE,
+                       MPLE.constraints.ignore = FALSE,
 
                        MCMC.prop=trim_env(~sparse + .triadic),
                        MCMC.prop.weights="default", MCMC.prop.args=list(),

R/ergm.mple.R

@@ -72,7 +72,7 @@ ergm.mple<-function(s, s.obs, init=NULL,
 
   # test whether the MPLE actually exists
   # FIXME: Figure out how to test for MPLE's existence in penalised and curved MPLEs.
-  if(! control$MPLE.type%in%c("penalized","logitreg"))  mple.existence(pl)
+  if(control$MPLE.check && ! control$MPLE.type%in%c("penalized","logitreg"))  mple.existence(pl)
 
   message("Maximizing the pseudolikelihood.")
   if(control$MPLE.type=="penalized"){