Summary of Aspen Data & Variables

Example subset of stand-level Aspen data from TreeMap:
MgmtID StandID Year Species DBHClass LiveTpa MrtTPA HrvTPA RsdTPA LiveBA MrtBA HrvBA RsdBA LiveTCuFt MrtTCuFt HrvTCuFt RsdTCuFt LiveMCuFt MrtMCuFt HrvMCuFt RsdMCuFt LiveBdFt MrtBdFt HrvBdFt RsdBdFt RunTitle CaseID StandBA pBA
CRbs 10354 2025 POTR5 All 17.035972 0.5629876 0 17.035972 17.366452 0.4733497 0 17.366452 456.22493 11.969565 0 456.22493 436.88563 11.429029 0 436.88563 1996.8072 46.10289 0 1996.8072 eastern-sugarpine-base 07ab727d-e545-43ed-a55a-6e2b7a12e271 103.51123 0.1677736
CRbs 10516 2025 POTR5 All 23.120510 0.5350774 0 23.120510 20.069320 0.4496264 0 20.069320 633.13865 14.332618 0 633.13865 601.31755 13.600376 0 601.31755 2816.5587 62.39936 0 2816.5587 eastern-sugarpine-base 02c6a721-00d7-40d5-a8ab-769996a4eb92 185.17822 0.1083784
CRbs 11423 2025 POTR5 All 5.482844 0.3715488 0 5.482844 2.486656 0.1685052 0 2.486656 64.50596 4.371164 0 64.50596 58.50206 3.964305 0 58.50206 163.5948 11.08559 0 163.5948 eastern-sugarpine-base 2ec14edf-3bc8-4e55-a178-eb8ac91d3a0d 263.38471 0.0094412
CRbs 11433 2025 POTR5 All 57.147223 1.6777010 0 57.147223 21.452504 0.6046028 0 21.452504 563.22390 15.918083 0 563.22390 497.38292 13.950890 0 497.38292 786.9044 18.32762 0 786.9044 eastern-sugarpine-base 5d5a51d6-12fc-45f2-b89f-25446c944db3 138.91260 0.1544317
CRbs 11434 2025 POTR5 All 5.440652 0.3211918 0 5.440652 8.295720 0.4894587 0 8.295720 308.51935 18.203379 0 308.51935 298.11681 17.589016 0 298.11681 1561.8632 92.14618 0 1561.8632 eastern-sugarpine-base 30bfc306-690a-42af-830a-3982b1a029cc 66.68517 0.1244013
CRbs 11437 2025 POTR5 All 124.777616 4.4138805 0 124.777616 107.294122 3.4441384 0 107.294122 3652.11449 115.086686 0 3652.11449 3476.95339 109.122781 0 3476.95339 15752.5105 475.13556 0 15752.5105 eastern-sugarpine-base 407afc2c-f34a-40bb-b3f6-258530026c70 202.08057 0.5309473

Explore unfiltered Live BA for all tm_id pixels that contain Aspen

fig <- plot_ly(x = ~Asp.25$LiveBA, type = "histogram")%>%
  layout(title='Aspen BA (unfiltered)')

fig
fig <- plot_ly(x = ~Asp.25$pBA, type = "histogram")%>%
  layout(title='Aspen Proportional BA (unfiltered)')

fig
asp.plot=plot_ly()%>%
  add_trace(data = Asp.25, y = ~ LiveTpa, x = ~Year, type = 'box', color=I("forestgreen"))
asp.plot
asp.tpa=plot_ly()%>%
  add_trace(data = Asp.25, y = ~ LiveBA, x = ~Year, type = 'box', color=I("goldenrod"))
asp.tpa
ID.top = ID.All %>% 
  filter(Species %in% c("PICO","POTR5", "PIPO", "PSME"))%>%
  arrange(desc(pBA))#, 
p1 = ID.top %>% 
  ggplot(aes(x = pBA, y = Species, fill = Species)) +
  
  # add half-violin from {ggdist} package
  stat_halfeye(
    # adjust bandwidth
    adjust = 0.5,
    # move to the right
    justification = -0.2,
    # remove the slub interval
    .width = 0,
    point_colour = NA
  ) +
  geom_boxplot(
    width = 0.12,
    # removing outliers
    outlier.color = NA,
    alpha = 0.5
  ) +
  # Themes and Labels
  scale_fill_tq() +
  theme_tq() +
  labs(
    title = "Proportional BA by Species",
    x = "pBA",
    y = "Species",
    fill = "Species"
  )

p1