Thanks a lot. Both solutions seem to work, although the BW option with LAPLACE method needs a nested design of random effects -birdidentity(region)- to correct the d.f. As we are not especially interested in obtaining very accurate estimates or AIC values the DDFM=KR (or SAT) with residual ML or PL is also a good option.
... View more