1. 报错信息

示例代码：

library(asreml)
data("harvey")
head(harvey)
ainv = ainverse(harvey)
head(ainv)
m2 = asreml(y3 ~ Line, random = ~ vm(Calf = ainv), data=harvey)

报错信息：

Error in asr_vm(Calf = ainv, data = data) : unused argument (Calf = ainv)
In addition: Warning messages:
1: In `[.data.table`(data, , `:=`((vv), lapply(vv, function(x) eval(parse(text = x),  :3 column matrix RHS of := will be treated as one vector
2: In `[.data.table`(data, , `:=`((vv), lapply(vv, function(x) eval(parse(text = x),  :Supplied 417 items to be assigned to 65 items of column 'ainv' (352 unused)

2. 问题解析

这典型的是asreml3r转到asreml4r的不适合综合征，vm中ID和ainv之间应该用逗号，而不是等于号：
Right:
vm(Calf,ainv)
Wrong:
vm(Calf=ainv)

3. 正确代码示例

library(asreml)
data("harvey")
head(harvey)
ainv = ainverse(harvey)
head(ainv)
m2 = asreml(y3 ~ Line, random = ~ vm(Calf = ainv), data=harvey)

> summary(m2)$varcompcomponent std.error  z.ratio bound %ch
vm(Calf, ainv)  500.3245  498.1786 1.004307     P 0.3
units!R         273.5984  409.7263 0.667759     P 0.0