def sim_counts(freq_param_ts, freq_rv):
    """
    Simulate count process with non-stationary parameter (currently only 1d)
    """
    counts = [freq_rv(l_t).rvs() for l_t in freq_param_ts]
    return(counts)