N <- 200 y_obs <- c(-3.057663e+00, 3.934837e-01, 9.644205e-01, -7.624387e-01, 3.226636e+00, -3.626843e-01, 5.364191e+00, 1.296407e+01, 1.959366e+01, 1.825506e+01, 2.594398e+01, 2.786073e+01, 1.006841e+01, 3.458381e+01, 1.310577e+01, 4.101316e+01, 2.270777e+01, 2.288842e+01, 2.403384e+01, 2.342435e+01, 4.616151e+01, 4.558888e+01, 4.358490e+01, 2.226437e+01, 1.757128e+01, 3.206889e+01, 2.977210e+01, -6.044165e-01, 3.295767e+01, 2.333601e+01, 1.721693e+01, -6.163232e+00, 2.139006e+01, 2.177696e+01, 3.365025e+00, 1.095881e+01, 1.146749e+01, 1.859508e+01, 1.079218e+01, 9.567876e+00, -4.242011e+00, 1.252658e+01, -4.643001e+00, 1.623040e+00, -2.013934e+00, -4.834483e+00, 1.594055e+01, 1.415615e+01, -7.618235e+00, 1.794551e+01, 8.338405e-01, 1.814164e+01, -6.761922e+00, 1.428407e+01, 1.608445e-01, 1.123914e+01, 4.902329e+00, 1.430029e+00, 6.134329e+00, 1.759869e+00, 9.261259e-01, 5.058128e+00, 2.180423e+01, 1.690808e+01, 1.083946e+01, 2.506721e+01, 2.277659e+01, 2.262838e+01, 2.861009e+01, 1.508836e+01, 1.239550e+01, 1.772265e+01, 1.358432e+01, 5.960574e+00, 2.070645e+00, 7.709353e+00, 2.654013e+00, -6.467063e+00, 2.423809e+00, -1.510695e+01, -1.632846e+01, 4.994124e-01, -1.103995e+00, -1.534157e+00, -9.516422e+00, -6.937091e+00, -2.582503e+00, 2.220652e+00, -2.147436e+01, -2.233727e+01, -3.201241e+00, -2.239910e+01, -1.537922e+01, 1.642056e+00, -2.323601e+01, -6.469554e+00, -3.616788e+01, -2.294741e+01, -2.573157e+01, -3.433442e+01, -1.673511e+01, -1.771076e+01, -2.374033e+01, -3.305220e+01, -1.602951e+01, -1.693395e+01, -3.500640e+01, -4.140945e+00, -8.679148e+00, -1.829626e+01, -1.236252e+01, -1.739295e+01, -1.586489e+01, -1.134302e+01, -8.159659e+00, -1.003693e+01, -9.856867e+00, -1.244714e+00, -9.699456e+00, -1.205075e+00, 8.746308e+00, 7.569233e+00, 7.351682e+00, -1.930489e+00, -5.074708e+00, -2.565565e+00, -2.209159e+00, -8.271702e-01, -6.087701e+00, -2.840574e+00, -1.155985e+01, -3.932338e+00, 1.524788e+00, -1.930625e+01, -2.393919e+01, -1.814665e+01, 9.623492e+00, 1.206280e+01, 1.257934e+01, 1.309632e+01, 9.464990e+00, 9.455126e+00, 5.664491e+00, -2.709795e+01, -2.879625e+01, -2.435599e+01, 5.156196e+00, 1.243802e+01, 9.730232e+00, 4.831957e+00, -1.056241e+00, -2.719436e+01, -2.336785e+01, -2.164253e+01, 9.348593e+00, -6.714123e+00, -8.206383e+00, -5.475411e-01, 1.555114e+01, 1.340368e+01, -1.049673e+01, -1.106348e+01, -1.597448e+01, -1.874201e+01, -8.722378e+00, 1.818545e+01, -6.343910e+00, 1.840087e+01, 1.779397e+01, 2.046048e+01, 9.470498e+00, 3.716778e+01, 4.443878e+01, 4.277739e+01, 4.747872e+01, 2.129356e+01, 4.839371e+01, 5.517677e+01, 5.337970e+01, 5.999150e+01, 1.820266e+01, 2.622632e+01, 3.013788e+01, 2.096102e+01, 5.193941e+01, 3.092789e+01, 3.201391e+01, 6.016909e+01, 6.451852e+01, 5.359153e+01, 5.348913e+01, 5.189238e+01, 4.574828e+01, 4.513777e+01, 4.709225e+01, 4.765625e+01, 4.528252e+01, 4.693325e+01, 4.198988e+01, 3.749127e+01) obs <- c(1, 2, 1, 1, 2, 2, 2, 2, 1, 1, 2, 2, 1, 2, 1, 2, 1, 1, 1, 1, 2, 2, 2, 1, 1, 2, 2, 1, 2, 2, 2, 1, 2, 2, 1, 2, 2, 2, 2, 2, 1, 2, 1, 1, 1, 1, 2, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 1, 1, 1, 1, 1, 2, 1, 1, 2, 2, 2, 2, 1, 1, 2, 2, 2, 1, 1, 2, 2, 1, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 2, 2, 1, 2, 1, 2, 1, 1, 2, 1, 1, 2, 2, 1, 1, 2, 1, 1, 2, 2, 2, 2, 2, 1, 2, 2, 2, 1, 2, 2, 2, 1, 2, 1, 1, 2, 2, 1, 2, 1, 2, 2, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 1, 1, 2, 2, 1, 1, 1, 1, 1, 2, 1, 2, 2, 2, 1, 2, 2, 2, 2, 1, 2, 2, 2, 2, 1, 1, 1, 1, 2, 1, 1, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 2, 1, 1)