对于我的博士学位,我目前正在使用以下命令在 AS Reml R 中运行双变量方差分量估计:
Test <- asreml(cbind(Reflex,lakt_1) ~ trait + Betriebsnummer +
at(trait,1):gestLength + at(trait,1):halfyear,
random= ~us(trait):vm(SampleID,Gs),
residual= ~units:diag(trait),
data= phenW,
workspace= 128e06)
但每次我尝试运行它时,我都会收到以下错误消息:
错误
levels <- (*tmp*, value = as.character(levels)):因子级别 [4] 重复
我的固定因子(Betriebsnummer和gestLength)halfyear都被定义为因子,我还将它们的属性设置为“唯一”。我也用我的随机术语做到了这一点,SampleID是一个因素和Gs我的遗传关系矩阵。phenW是我的数据集,它是一个数据框。Reflex并且lakt_1是我的响应变量。
有人可以帮助我吗?