Add density summary for later-burned interior

youngdjn · Sep 29, 2024 · e4fcdcd · e4fcdcd
1 parent a6f98e1
commit e4fcdcd
Showing 1 changed file with 10 additions and 0 deletions.
diff --git a/scripts/analysis/summary-stats-for-paper.R b/scripts/analysis/summary-stats-for-paper.R
@@ -42,6 +42,16 @@ d_sp |>
             dens_median = median(seedl_dens_sp),
             n = n())
 
+## Specifically later-burned interior plots
+d_sp |>
+  filter(grn_vol_abs_sp == 0,
+         ((is.na(dist_grn_sp) | dist_grn_sp > minimum_dist_green) & sight_line > minimum_dist_green),
+         plot_type %in% c("core", "delayed"),
+         day_of_burning > 210) |>
+  summarize(dens_mean = mean(seedl_dens_sp),
+            dens_median = median(seedl_dens_sp),
+            n = n())
+
 
 ## Seed wall plots burning BEfore 1 Aug
 d_sp |>