asreml 4r 報錯

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)$varcomp
               component 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