def nburst(n=100, rt=15.0, dur=100, prec=15, sd=2000, latsd=5000, slp=2, mpt=3, sprob=.9):
dur = int(dur * 1000)
prec = int(prec * 1000)
m = T.BURST(dur, rt, prec, 1, sd=sd, latsd=latsd, sprob=sprob, slp=slp, mpt=mpt)
d = gd.Doc()
d['condburst'] = T.draw_single(m, 0, n)
rt = d['condburst.evts']
rt = len(flat(rt)) / float(len(rt))
rt = rt * 1e6 / dur
print(rt)
m = T.AC(dur, rt, prec, 3)
d['condhpois'] = T.draw_single(m, 1, n)
return d
def rburst(n=200, prate=10):
dur = 100000
d = gd.Doc()
mp = T.AC(dur, prate, 1, 3)
d['condpois'] = T.draw_single(mp, 0, n)
m = T.BURST(dur, 1, 15000, 1, sd=1000, latsd=2000, mfst=10000, sprob=.8, slp=2, mpt=3)
b = T.draw_single(m, 0, n)
rt = 1e6 * len(flat(b['evts'])) / float(len(b['evts'])) / dur
print('burst1 rate', rt)
p = T.draw_single(mp, 0, n)
rt = 1e6 * len(flat(p['evts'])) / float(len(p['evts'])) / dur
print('pois rate', rt)
d['condb1'] = combinedraws(b, p)
m = T.BURST(dur, 1, 10000, 1, sd=2000, latsd=2000, mfst=15000, sprob=.9, slp=3, mpt=2)
b = T.draw_single(m, 0, n)
rt = 1e6 * len(flat(b['evts'])) / float(len(b['evts'])) / dur
print('burst2 rate', rt)
p = T.draw_single(mp, 0, n)
d['condb2'] = combinedraws(b, p)
return d
