[R-sig-ME] coding a certain model

Lin, Shin Shin.Lin at uphs.upenn.edu
Wed Nov 12 03:31:25 CET 2008

Does anyone know how to code the following?
As on page 162 of Mixed Effects Models in S and S-Plus by Pinheiro and Bates with the same b_i** and PSI** but with W_i = 
1 1 0 0                         1
1 0 1 0    Kronecker product    1 
1 1 1 1                         1
The difference is in the last row of the 3x4 matrix. 
Granted, it's a weird model, but this situation actually mimics a biological scenario.  I'm guessing that the solution, if there is one, would be in the manipulation of the groupedData object (Oats in this case).

The information contained in this e-mail message is intended only for the personal and confidential use of the recipient(s) named above. If the reader of this message is not the intended recipient or an agent responsible for delivering it to the intended recipient, you are hereby notified that you have received this document in error and that any review, dissemination, distribution, or copying of this message is strictly prohibited. If you have received this communication in error, please notify us immediately by e-mail, and delete the original message.

More information about the R-sig-mixed-models mailing list