growth_and_yield_evaluation: Extensive Growth and Yield Evaluation

View source: R/growth_and_yield_evaluation.R

growth_and_yield_evaluationR Documentation

Extensive Growth and Yield Evaluation

Description

Provide an extensive compilation of growth and yield related results.

Usage

growth_and_yield_evaluation(sim_agg, concept_def, detailed_out = FALSE)

Arguments

sim_agg

Aggregated simulated area dynamics, more precisely, the output of aggregate_raw_sim_rslt

concept_def

Concept definition matching sim_agg, a c4c_concept object

detailed_out

Boolean, if TRUE, also pre-evaluation output (calculated by the internal function growth_and_yield_pre_eval) will be part of the result list (default = FALSE)

Details

The result object, a list, contains the following elements:

  • gyield_summary: A tibble that contains phase overarching growth and yield information. In this tibble, each row is a point in time. The columns (in addition to time) are: vol_standing, the standing volume on the total area of interest; vol_rmv_cont,the continuous removals that take place, as long as an area is in a given phase; vol_rmv_trans, the volume removals occurring at phase transitions from phases with more to such with less standing volume; vol_rmv_damage, the volume losses due to damage events; vol_rmv_harvest, regular harvest volume, the sum of vol_rmv_cont and vol_rmv_trans; vol_rmv_total, all removed volume, the sum of vol_rmv_harvest and vol_rmv_damage; vol_mort, the mortality volume (normal, not event-triggered); vol_inc, the volume increment on the whole area resulting from the vol_standing, vol_rmv_total, and vol_mort.

  • gyield_phases: A list of tibbles, one for each variable as contained in gyield_summary (except the volume increment) which makes no sense in a phase-wise context), but in phase-wise resolution. In each tibble, each row is a point in time, the columns represent the stand development phases.

  • gyield_pre (in case the user has chosen detailed_out = TRUE): A tibble, where each row is a point in time. The columns (in addition to time) are: vol_standing, the standing volume on the total area of interest; vol_rmv_cont,the continuous removals that take place, as long as an area is in a given phase; vol_rmv_trans, the volume removals occurring at phase transitions from phases with more to such with less standing volume;vol_rmv_damage, the volume losses due to damage events; vol_rmv_harvest, regular harvest volume, the sum of vol_rmv_cont and vol_rmv_trans; vol_rmv_total, all removed volume, the sum of vol_rmv_harvest and vol_rmv_damage; vol_mort, the mortality volume (normal, not event-triggered); vol_inc, the volume increment on the whole area resulting from the vol_standing, vol_rmv_total, and vol_mort.

Value

A list with two elements (see also details),

  • gyield_summary: A tibble that contains phase overarching growth and yield information.

  • gyield_phases: A list of tibbles, each one representing one of the growth and yield variables also found in gyield_summary, but here on the level of the single stand development phases.

If the user has selected detailed_out = TRUE, there will also be another list element, gyield_pre, which contains the interim information from which gyield_summary and gyield_phases are calculated.

Examples

  # Run a simulation and store the aggregated outcome
  state_vars <- setup_statevars(
    pine_thinning_from_above_1, c(1000, 0, 0, 0, 0, 0)
  )
  time_span <- 200
  parms     <- setup_parms(pine_thinning_from_above_1)
  parms$risk_mat <- setup_risk_events(
    time_span, avg_event_strength = 1, parms$risk
  )

  sim_areas_agg <- sim_area_single_concept_with_risk(
    state_vars,
    parms       = parms,
    event_times = c(0:time_span),
    time_span   = time_span
  ) |>
    aggregate_raw_sim_rslt(pine_thinning_from_above_1)

  # Growth and yield evaluation
  growth_and_yield_evaluation(sim_areas_agg, pine_thinning_from_above_1)


care4cmodel documentation built on April 4, 2025, 1:39 a.m.