some update to get mkulazi to work

This commit is contained in:
Martin Folkerts 2024-11-18 13:41:58 +01:00
parent 576250982c
commit 3bfa7e1750
2 changed files with 6 additions and 2 deletions

View file

@ -275,7 +275,11 @@ The below table shows estimates of the biomass if you would harvest them now.
```{r message=FALSE, warning=FALSE, include=FALSE}
CI_quadrant <- readRDS(here(cumulative_CI_vals_dir,"All_pivots_Cumulative_CI_quadrant_year_v2.rds"))
#CI_quadrant <- readRDS(here(cumulative_CI_vals_dir,"All_pivots_Cumulative_CI_quadrant_year_v2.rds"))
CI_quadrant <- readRDS(here(cumulative_CI_vals_dir,"All_pivots_Cumulative_CI_quadrant_year_v2.rds")) %>%
group_by(model) %>% # Group by model or other relevant columns
fill(field, sub_field, .direction = "downup") %>% # Fill down then up within each group
ungroup()
harvesting_data <- harvesting_data %>% rename(season = year)
@ -348,7 +352,7 @@ ggplot(pred_rf_current_season, aes(x = Age_days, y = predicted_Tcha)) +
labs(title = "Predicted Yields for Fields Over 300 Days \nOld Yet to Be Harvested",
x = "Age (days)",
y = "Predicted tonnage/ha (Tcha)") +
facet_wrap(~sub_area) +
# facet_wrap(~sub_area) +
scale_y_continuous(limits = c(0, 200)) + # Optional: Set limits for y-axis
theme_minimal()

Binary file not shown.