simulate_sir.jl

using MarkovProcesses

load_model("SIR")

σ = simulate(SIR)

d1 = Dict(:S => 1, :I=> 2, :R => 3)
d2 = Dict(:I => 1)

bool_test = SIR.g == [:I] && SIR._g_idx == [2] && 
            SIR.map_var_idx == d1 && 
            SIR._map_obs_var_idx == d2

return bool_test