diff --git a/blockmodels_correct_sampling_effort.R b/blockmodels_correct_sampling_effort.R
index 9802ee6..0f6dc76 100644
--- a/blockmodels_correct_sampling_effort.R
+++ b/blockmodels_correct_sampling_effort.R
@@ -26,6 +26,10 @@ message("At level: ", level, " starting at: ", epoch)
 
 save_path <- here("results", "seq-effort", paste0(author, "_", level, "_", epoch, ".Rds"))
 
+if (!dir.exists(save_path)) {
+    dir.create(save_path, recursive = TRUE)
+}
+
 per_taxa_networks <- collapse_otu_at_taxo(the_data)
 
 covar_mat <- t(colSums(per_taxa_networks[[level]]) %*% t(rep(1, nrow(per_taxa_networks[[level]]))))