兽医流行病学，数据为4000个奶牛场10年产量和疾病发生数据，已知有些没有发生疾病，有些场发病，想要分析发病后1年、2年、3年...对奶产量的影响。现在用混合线性模型分析，奶牛场ID作为随机变量，发病后年份year.bd作为固定变量，奶产量作为因变量。但因奶产量随时间变动，需剔除正常年份year对奶产量影响，不知道应如何在R里建模。求指导！