Python simulate_age_group_data примеры использования

Пример #1

Файл: age_group_integrate.py Проект: aflaxman/gbd

import scipy.interpolate

import graphics
import book_graphics

reload(book_graphics)


import age_group_models as agm

reload(agm)

m = {}

mc.np.random.seed(12345)
f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [0, 0.1, 0.9, 0.2, 0])
model = agm.simulate_age_group_data(N=25, delta_true=5e2, pi_true=f)
agm.fit_age_integrating_model(model)

m[0] = model

mc.np.random.seed(12345)
f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [0, 0.1, 0.9, 0.2, 0])
model = agm.simulate_age_group_data(N=25, delta_true=5, pi_true=f)
agm.fit_age_integrating_model(model)
m[1] = model

agm.plot_fits(m)
pl.savefig("age_group_integrate.pdf")
pl.show()

Пример #2

import graphics
import book_graphics

reload(book_graphics)

import age_group_models as agm

reload(agm)

mc.np.random.seed(123456)

m = {}

f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [0, .1, .9, .2, 0])
model = agm.simulate_age_group_data(N=25, delta_true=5e2, pi_true=f)
agm.fit_midpoint_covariate_model(model)
m[0] = model

mc.np.random.seed(123456)

f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [0, .1, .9, .2, 0])
model = agm.simulate_age_group_data(N=25, delta_true=5, pi_true=f)
agm.fit_midpoint_covariate_model(model)
m[1] = model

agm.plot_fits(m)
pl.savefig('book/graphics/age_group_midpoint_covariate.pdf')

pl.show()

Пример #3

Файл: age_group_midpoint.py Проект: aflaxman/gbd

import scipy.interpolate

import graphics
import book_graphics

reload(book_graphics)


import age_group_models as agm

reload(agm)

mc.np.random.seed(123567)

m = {}

f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [0.4, 0.425, 0.6, 0.5, 0.4])
model = agm.simulate_age_group_data(N=20, delta_true=50, pi_true=f)
agm.fit_midpoint_covariate_model(model)
m[0] = model

f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [0, 0.1, 0.9, 0.2, 0])
model = agm.simulate_age_group_data(N=20, delta_true=50, pi_true=f)
agm.fit_midpoint_covariate_model(model)
m[1] = model

agm.plot_fits(m)
pl.savefig("book/graphics/age_group_midpoint.pdf")

pl.show()

Пример #4

Файл: age_group_midpoint.py Проект: studentmicky/gbd

import pylab as pl
import pymc as mc
import scipy.interpolate

import graphics
import book_graphics
reload(book_graphics)


import age_group_models as agm
reload(agm)

mc.np.random.seed(123567)

m = {}

f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [.4, .425, .6, .5, .4])
model = agm.simulate_age_group_data(N=20, delta_true=50, pi_true=f)
agm.fit_midpoint_covariate_model(model)
m[0] = model

f = scipy.interpolate.interp1d([0, 20, 40, 60, 100], [0, .1, .9, .2, 0])
model = agm.simulate_age_group_data(N=20, delta_true=50, pi_true=f)
agm.fit_midpoint_covariate_model(model)
m[1] = model

agm.plot_fits(m)
pl.savefig('book/graphics/age_group_midpoint.pdf')

pl.show()