最近在学习Mixed model 的理论和在R中的运用,但是碰到了一些问题,我的数据是个体*环境*区组(200*2*3)有关基因与环境互作的。如下表所示

,我的分析思路是这样的,
1.首先分单个环境逐一分析,每个环境下的数据就做简单的区组、个体两项方差分析,得出两个环境各自的方差,作为权重因子进入下一步的计算;
2.接下来进行个体与环境的mixed model的分析,此时个体用的就是平均值了。
具体思路是这样的,但是代码可能有些问题,具体代码如下

,我想问是不是权重的加入搞错了,因为翻阅了一些英文的教程书,也没详细写怎么加入权重,只是很粗略的写weights = varFixed(~age)这种根本不知道它在说什么玩意的东西。。

很是头大!麻烦高手们有空帮忙看看!!!谢谢谢谢谢谢!!困扰我好久好久了!!