# Generate sequences of the Poisson events and plot histogram to check predictive distribution

include("./metadata.jl")
using Metadata

function run_sim()
    n_sim = 500
    n = 100
    w = 30
    seqs = cumsum(rand(1 : w, n))
    for i in 1 : n_sim
        new_dat = cumsum(rand(1 : w, n))
        seqs = vcat(seqs, new_dat)
    end  
    println("sim_poisson_hist.jl > Save generated data ino ", outdir*"sim_poisson_hist.txt")
    writedlm(outdir*"sim_poisson_hist.txt", seqs)
end 
run_sim() 

