1

对于我的博士学位,我目前正在使用以下命令在 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] 重复

我的固定因子(BetriebsnummergestLengthhalfyear都被定义为因子,我还将它们的属性设置为“唯一”。我也用我的随机术语做到了这一点,SampleID是一个因素和Gs我的遗传关系矩阵。phenW是我的数据集,它是一个数据框。Reflex并且lakt_1是我的响应变量。

有人可以帮助我吗?

4

0 回答 0