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