diff --git a/R/summarize_BSM_tables_v1.R b/R/summarize_BSM_tables_v1.R
index 7c0f32a..8a87518 100644
--- a/R/summarize_BSM_tables_v1.R
+++ b/R/summarize_BSM_tables_v1.R
@@ -583,7 +583,7 @@ count_ana_dispersal_events <- function(ana_events_table, areanames, actual_names
 			} # END for (j in 1:length(areanames))
 
 		# Extinction/extirpation vector
-		TF = events_df2$extirpation_from == areanames[i]
+		TF = extirp_df2$extirpation_from == areanames[i]
 		if (sum(TF) > 0)
 			{
 			e_counts_list[,i] = sum(TF)
@@ -1318,6 +1318,7 @@ count_ana_clado_events <- function(clado_events_tables, ana_events_tables, arean
 				{
 				e_counts_rectangle[i,] = as.matrix(e_counts_df)
 				e_totals_list[i] = sum(e_counts_df)
+				ana_totals_list[i] = ana_totals_list[i] + e_totals_list[i]
 				} else {
 				e_counts_rectangle[i,] = 0
 				e_totals_list[i] = 0