# <markdowncell> # ### Rerun the model # <codecell> sim.run(0, 12) # <markdowncell> # ### Plot the posterior histogram # <codecell> model.plot_distributions() # <markdowncell> # ### Plot the posterior histogram with the normal approximation # <codecell> model.plot_distributions(show_normal=True) # <markdowncell> # ### Fit the model parameter, $a$, and the initial value of the variable, $h$ # <codecell>
# <markdowncell> # ### Rerun the model # <codecell> sim.run(0,12) # <markdowncell> # ### Plot the posterior histogram # <codecell> model.plot_distributions() # <markdowncell> # ### Plot the posterior histogram with the normal approximation # <codecell> model.plot_distributions(show_normal=True) # <markdowncell> # ### Fit the model parameter, $a$, and the initial value of the variable, $h$ # <codecell>
# <codecell> model=MCMCModel(sim,{'b':[0,10]}) model.fit(iter=25000) # <codecell> model.b # <codecell> sim.run(5) # <codecell> model.plot_distributions() # <codecell> model.sigma # <codecell> print "The parameter b has best estimate",model.mu['b'],' +/- ',model.sigma['b'] # <markdowncell> # ## A linear growth example # # Data from [http://www.seattlecentral.edu/qelp/sets/009/009.html](http://www.seattlecentral.edu/qelp/sets/009/009.html)
# <codecell> model = MCMCModel(sim, {'b': [0, 10]}) model.fit(iter=25000) # <codecell> model.b # <codecell> sim.run(5) # <codecell> model.plot_distributions() # <codecell> model.sigma # <codecell> print "The parameter b has best estimate", model.mu['b'], ' +/- ', model.sigma[ 'b'] # <markdowncell> # ## A linear growth example # # Data from [http://www.seattlecentral.edu/qelp/sets/009/009.html](http://www.seattlecentral.edu/qelp/sets/009/009.html)
sim.run(0, 350) # <codecell> model = MCMCModel( sim, { 'beta0': [0, 1], 'invGamma': [3.5, 10.7], 'beta1': [0, 1], 'q': [0, 100], 'tau': [100, 150], 'invk': [5, 22] }) #model = MCMCModel(sim,{'invGamma':[3.5,10.7],'q':[0,10]}) model.fit(iter=500) model.plot_distributions() # <codecell> Beta0 = model.beta0 # <codecell> Beta1 = model.beta1 # <codecell> #beta = (Beta0 + Beta1)/2 # <codecell>