single.impute.dyads() should not mess with networks when the constrai…

…nts are dyad-dependent.

DESCRIPTION

@@ -1,6 +1,6 @@
 Package: ergm
-Version: 4.7-7367
-Date: 2024-06-05
+Version: 4.7-7368
+Date: 2024-06-11
 Title: Fit, Simulate and Diagnose Exponential-Family Models for Networks
 Authors@R: c(
              person(c("Mark", "S."), "Handcock", role=c("aut"), email="[email protected]"),

R/ergm.utility.R

@@ -144,6 +144,12 @@ single.impute.dyads <- function(nw, constraints=NULL, constraints.obs=NULL, min_
   response <- nw %ergmlhs% "response"
   stopifnot(!is.null(constraints)||is.null(constraints.obs))
 
+  if(!is.dyad.independent(constraints) || !is.dyad.independent(constraints.obs)){
+    message("Model and/or observational constraints are not dyad-independent. Dyad imputation cannot be used. Please ensure your LHS network satisfies all constraints.")
+    if(output=="network") return(nw)
+    else return(ergm_state(nw))
+  }
+
   if(!is.null(constraints)){
     imputable <- as.rlebdm(constraints, constraints.obs, "missing")
     nae <- NVL3(imputable, sum(.), 0)