[Developers] Random effects without penalties work: should they?

[Hans J. Skaug]

Hi,

I have noted that if you define random effects in your model,
but do not use them for anything [u(M+1) and u(M+2)) in
the attach model] the program runs fine. This
used to cause problems, because the Laplace approximation
is not well defined.

My guess is that the reason why this now works is that
the effort to avoid the "infinite loop" in the inner problem
that was carried out some years ago has introduced 
some penalty on the random effects.

Maybe the new behaviour is OK, but it means that you are no longer
told if you introduce a random effect that is not used for anything.
(Probably you did not do this unintentionally.

Mostly I was curious about finding out what has caused the change.

Hans
-------------- next part --------------
A non-text attachment was scrubbed...
Name: orange.dat
Type: application/octet-stream
Size: 541 bytes
Desc: orange.dat
URL: <http://lists.admb-project.org/pipermail/developers/attachments/20121214/6a611fb0/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: orange.pin
Type: application/octet-stream
Size: 59 bytes
Desc: orange.pin
URL: <http://lists.admb-project.org/pipermail/developers/attachments/20121214/6a611fb0/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: orange.tpl
Type: application/octet-stream
Size: 2046 bytes
Desc: orange.tpl
URL: <http://lists.admb-project.org/pipermail/developers/attachments/20121214/6a611fb0/attachment-0002.obj>