list(theta=c(3700, 3600, 3550, 3504, 3586, 3800, 3529, 3525, 3500, 3500, 3402, 3542, 3492, 3618, 3148, 3638))