def generate_dataset_rep(theta):
hp = HillRepressilator(alpha=theta[0],alpha0=theta[1], beta=theta[2], n=theta[3])
ds = hp.run(T=50)
dataset = np.zeros([data_points, 3])
for ind, time in enumerate(times):
dataset[ind, :] = ds[time, :]
return dataset
def generate_dataset_hp(theta):
hp = HillRepressilator(*theta, n=3.)
return hp.run(20)[:, 3:]