struc_eval_impact: Evaluate impact equations numerically

In jandraor/structuR: Structural analysis on System Dynamics models

Description

Evaluate impact equations numerically

Usage

struc_eval_impact(impact_df, sim_df)

Arguments

impact_df	Data frame returned by the function struc_impact_on
sim_df	Data frame with the simulation results

Examples

I_z_x <- "-(R * x) * (a * y - z) / (-R * x * z)"
I_x_x <- "-(R * z) * (-R * x * z) / (-R * x * z)"

impact_df <- data.frame(from    = c("z", "x"),
                        to      = c("x", "x"),
                        through = c("f1", "f1"),
                        impact  = c(I_z_x, I_x_x))

x0  <- 0.974182771
x10 <- 0.466656479
y0  <- 0.009831036
y10 <- 0.078797261
z0  <- 0.007933189
z10 <- 0.078490487

sim_df <- data.frame(
  time = c(0L, 10L),
  x = c(x0, x10),
  y = c(y0, y10),
  z = c(z0, z10),
  R = c(2L, 2L),
  a = c(1L, 1L))
struc_eval_impact(impact_df, sim_df)

jandraor/structuR documentation built on Feb. 11, 2022, 10:24 p.m.