Avoid potential bug with Matrix library, by ensuring that `design$sub…

…set` is a logical vector.
bschneidr · Jul 9, 2024 · 76beceb · 76beceb
1 parent 6b4f361
commit 76beceb
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/R/quadratic_forms_of_survey_design_objects.R b/R/quadratic_forms_of_survey_design_objects.R
@@ -318,7 +318,10 @@ get_design_quad_form.twophase2 <- function(design, variance_estimator,
     variance_estimator = variance_estimator[[1]],
     aux_var_names = aux_var_names
   )
-  Sigma_phase1 <- Sigma_phase1[design$subset, design$subset]
+  Sigma_phase1 <- Sigma_phase1[
+    as.logical(design$subset), 
+    as.logical(design$subset)
+  ]
 
   # Extract the quadratic form for second phase
   # (conditional on first phase sample)