在眼科学领域的随机临床试验中，我想找到由线性混合效应模型给出的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
谢谢