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] 重复

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

有人可以帮助我吗？

4

0 回答 0