list(theta0 = -1, theta1 = -3, theta2 = 0,
gamma = c(0.00406, 0.00518, 0.00916, 0.0117, 0.0101, 0.00249, 0.0178, 0.00112, 0.000703,
0.00427, 0.0112, 0.00599, 0.00131, 0.00599, 0.000738, 0.00259))