nM <- 30 nP <- 5 nO <- 3 y <- c(8.352189e+01, 8.346839e+01, 8.443409e+01, 8.431714e+01, 8.563262e+01, 8.558676e+01, 1.084811e+02, 1.084628e+02, 1.094307e+02, 1.093453e+02, 1.104854e+02, 1.103567e+02, 1.027943e+02, 1.028578e+02, 1.036348e+02, 1.036713e+02, 1.048034e+02, 1.048445e+02, 9.505933e+01, 9.491565e+01, 9.568071e+01, 9.610842e+01, 9.704310e+01, 9.693448e+01, 1.054701e+02, 1.054411e+02, 1.062438e+02, 1.064447e+02, 1.071204e+02, 1.074093e+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 <- 1.000309e+01 sigmaOest <- 9.890137e-01 sigmaRepEst <- 8.431680e-02