# NOT RUN {
SG <- CGGPcreate(d=3, batchsize=100)
y <- apply(SG$design, 1, function(x){x[1]+x[2]^2})
SG <- CGGPfit(SG, Y=y)
MSE_MAP <- outer(1:SG$d, 1:8,
Vectorize(function(dimlcv, lcv1) {
CGGP_internal_calcMSE(SG$xb[1:SG$sizest[dimlcv]],
theta=SG$thetaMAP[(dimlcv-1)*SG$numpara+1:SG$numpara],
CorrMat=SG$CorrMat)
}))
CGGP_internal_calcMSEde(SG$po[1:SG$poCOUNT, ], MSE_MAP)
# }
Run the code above in your browser using DataLab