在眼科领域的随机临床试验中，我想找到由
给出的线性混合效应模型的R实现。log(y_ijk) = b0 + b1_j + b2_k + b3_jk + w_i + e_ijk
其中y_ijk是残留抗生素的数量，b1_j是结构的影响（j=1如果是角膜，j=2是水样的，b2_k是治疗效果（k=1,2,3），b3是互动效应，w_i是随机效应，e_ijk是随机错误。我们对每个个体都有成对的观察结果（每个受试者的两种结构都测量残留的抗生素）。
我的主要困难是模型应该包含错误的异方差结构，因此每个治疗组的方差不同。这样，
e_ijk ~ Normal(0,sigma_k), k = 1,2,3
由于