test new interface

[hyunjimoon]

list of test:

• data/lookup structure with subscript
• different combination of driving data, estimated and assumed parameters
• parameter's lower and upper bound

tested, not pass

1. how to estiamte adj_frac1? the following didn't work

v2s_code = """
ss_obs[timesteps] ~ normal(ss[timesteps], sigma);
adj_frac1<lower=0> ~ normal(0.2, 0.02);
sigma<lower=0> ~ normal(0.1, 0.01);
"""

2. how to use two data variables? passing ["ss_obs", "s_obs"] to Vensim2Stan() doesn't work well

3. when data structure has exo_demand[reg]

[Dashadower]

1. how to estiamte adj_frac1? the following didn't work
   v2s_code = """
   ss_obs[timesteps] ~ normal(ss[timesteps], sigma);
   adj_frac1<lower=0> ~ normal(0.2, 0.02);
   sigma<lower=0> ~ normal(0.1, 0.01);
   """

Running the above model with show_console=True for sample() yields the following error:

Chain [1] Gradient evaluated at the initial value is not finite.
Chain [1] Stan can't start sampling from this initial value.
Chain [1] 
Chain [1] Initialization between (-2, 2) failed after 100 attempts.
Chain [1] Try specifying initial values, reducing ranges of constrained values, or reparameterizing the model.
Chain [1] Initialization failed.