Skip to contents

Generate and propagate Monte Carlo Simulations based on a template input file.

Source: R/fct_combine_mcs_E.R
fct_combine_mcs_E.Rd

TBD

Usage

fct_combine_mcs_E(.ad, .cs, .usr, .time)

Arguments

.ad

Activity Data input table for the shiny app (AD_lu_transitions)

.cs

Carbon Stock input table for the shiny app (c_stocks)

.usr

User inputs' table for the shiny app (user_inputs). Contains the number of iterations of the MCS, carbon fraction if needed and if truncated PDFs should be used when necessary.

.time

the time period table from the input file (see template)

Value

A data frame with Monte Carlo simulations of CO2 emissions for each land use transition, REDD+ activity or emission reductions level.

Examples

library(mocaredd)
library(readxl)
library(dplyr)

path <- system.file("extdata/example2-with-sims.xlsx", package = "mocaredd")

cs   <- read_xlsx(path = path, sheet = "c_stocks", na = "NA")
ad   <- read_xlsx(path = path, sheet = "AD_lu_transitions", na = "NA")
usr  <- read_xlsx(path = path, sheet = "user_inputs", na = "NA")
time <- read_xlsx(path = path, sheet = "time_periods", na = "NA")

cs_clean <- cs |> filter(!is.na(c_value) | !is.na(c_pdf_a))
time_p   <- time |> dplyr::mutate(nb_years = year_end - year_start + 1)

res      <- fct_combine_mcs_E(.ad = ad, .cs = cs_clean, .usr = usr, .time = time_p)

get_trans <- sample(res$trans_id, 1)
res_sub <- res |> filter(trans_id == get_trans)

hist(res_sub$E)