nM <- 30
nP <- 5
nO <- 3
y <- c(7.650117e+01, 7.650499e+01, 7.585679e+01, 7.585256e+01, 7.453181e+01, 7.451245e+01, 8.103738e+01, 8.104761e+01, 8.033791e+01, 8.030458e+01, 7.906916e+01, 7.905566e+01, 1.443217e+02, 1.443228e+02, 1.436584e+02, 1.436182e+02, 1.423920e+02, 1.423718e+02, 4.158867e+01, 4.160192e+01, 4.089084e+01, 4.089339e+01, 3.965006e+01, 3.962394e+01, 1.609838e+02, 1.609655e+02, 1.602413e+02, 1.602854e+02, 1.589873e+02, 1.589909e+02)
part <- c(1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 4, 4, 4, 4, 4, 4, 5, 5, 5, 5, 5, 5)
oper <- c(1, 1, 2, 2, 3, 3, 1, 1, 2, 2, 3, 3, 1, 1, 2, 2, 3, 3, 1, 1, 2, 2, 3, 3, 1, 1, 2, 2, 3, 3)
sigmaPest <- 5.000175e+01
sigmaOest <- 9.987473e-01
sigmaRepEst <- 1.197054e-02