Niwot Ridge Dendrobands

2024 Update

# Growth for all by year
ggplot(data = dendro_reg) +
  geom_boxplot(mapping = aes(x = as.factor(Year), y = log(Growth), fill = as.factor(Year)), alpha = .5) +
  geom_jitter(mapping = aes(x = as.factor(Year), y = log(Growth), col = as.factor(Year)), alpha = .2) +
  theme_bw() +
  labs(x = "Year", fill = "", col = "", title = "Pooled Growth by Year")

# Growth for all by plot
ggplot(data = dendro_reg) +
  geom_boxplot(mapping = aes(x = as.factor(Plot), y = log(Growth), fill = as.factor(Plot)), alpha = .5) +
  geom_jitter(mapping = aes(x = as.factor(Plot), y = log(Growth), col = as.factor(Plot)), alpha = .2) +
  theme_bw() +
  labs(x = "Plot", fill = "", col = "", title = "Pooled Growth by Plot")

# Growth for all by plot
ggplot(data = dendro_reg) +
  geom_boxplot(mapping = aes(x = as.factor(Plot), y = log(Growth), fill = as.factor(Plot)), alpha = .6) +
  geom_jitter(mapping = aes(x = as.factor(Plot), y = log(Growth), col = as.factor(Plot)), alpha = .2) +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 50, vjust =1, hjust=1)) +
  labs(x = "Plot", col = "", fill = "", title = "Plot Growth by Year") +
  facet_wrap(~Year)

# Growth for Species by Year
ggplot(data = dendro_reg) +
  geom_boxplot(mapping = aes(x = as.factor(Year), y = log(Growth), fill = as.factor(Species)), alpha = .6) +
  geom_jitter(mapping = aes(x = as.factor(Year), y = log(Growth), col = as.factor(Species)), alpha = .2) +
  theme_bw() +
  labs(x = "Year", col = "", fill = "", title = "Species Growth by Year")

# Growth for Species by Year/Plot
ggplot(data = dendro_reg) +
  geom_boxplot(mapping = aes(x = as.factor(Year), y = log(Growth), fill = as.factor(Species)), alpha = .6) +
  geom_jitter(mapping = aes(x = as.factor(Year), y = log(Growth), col = as.factor(Species)), alpha = .2) +
  facet_wrap(~Plot) +
  theme_bw() +
  theme(axis.text.x = element_text(angle = 50, vjust =1, hjust=1)) +
  labs(x = "Year", col = "", fill = "", title = "Species Growth by Plot + Year")

# Growth Density by Year/Species
ggplot(data = dendro_reg) +
  geom_density(mapping = aes(x = log(Growth), fill = as.factor(Year)), alpha = .75) +
  theme_bw() +
  facet_wrap(~Species) +
  labs(fill = "Year", title = "Species Growth Density by Year")