This commit is contained in:
Martin Folkerts 2024-10-31 22:35:59 +01:00
parent c1a32bf6ad
commit 3093c86e3a
2 changed files with 5 additions and 3 deletions

View file

@ -272,8 +272,10 @@ print(" PRINT")
# Yield prediction # Yield prediction
The below table shows estimates of the biomass if you would harvest them now. The below table shows estimates of the biomass if you would harvest them now.
```{r eval=FALSE, message=FALSE, warning=FALSE, include=TRUE}
CI_quadrant <- readRDS(here(cumulative_CI_vals_dir,"All_pivots_Cumulative_CI_quadrant_year_v2.rds")) ```{r message=FALSE, warning=FALSE, include=FALSE}
CI_quadrant <- readRDS(here(cumulative_CI_vals_dir,"All_pivots_Cumulative_CI_quadrant_year_v2.rds"))
harvesting_data <- harvesting_data %>% rename(season = year) harvesting_data <- harvesting_data %>% rename(season = year)
@ -330,7 +332,7 @@ pred_rf_current_season <- prepare_predictions(predict(model_ffs_rf, newdata = pr
mutate(CI_per_day = round(total_CI / Age_days, 1)) mutate(CI_per_day = round(total_CI / Age_days, 1))
``` ```
```{r yield_plaatjes, eval=FALSE, include=TRUE} ```{r yield_plaatjes, echo=FALSE }
ggplot(pred_ffs_rf, aes(y = predicted_Tcha, x = tonnage_ha, col = sub_area)) + ggplot(pred_ffs_rf, aes(y = predicted_Tcha, x = tonnage_ha, col = sub_area)) +
geom_point(size = 2, alpha = 0.6) + # Adjust point size and transparency geom_point(size = 2, alpha = 0.6) + # Adjust point size and transparency
geom_abline(intercept = 0, slope = 1, linetype = "dashed", color = "red") + # Reference line geom_abline(intercept = 0, slope = 1, linetype = "dashed", color = "red") + # Reference line

Binary file not shown.