这是我的问题：我需要在 MCMCglmm 包中使用 hupossion。这是我之前的：

> prior = list(R = list(V = diag(2), nu = 0.002),
 +                G = list(G1 = list(V = diag(10), nu = 0.002),
 +                           G2 = list(V = 1, nu = 0.002)))

这是我的训练：

> myTrain <- MCMCglmm(pipenumyearnosameday ~ trait - 1 +
+                         at.level(trait,1):lastyearnumnosameday + 
+                         at.level(trait,1):pipeInsideSizemm,
+                     random = ~idh(pipematerialgroup1):installyeargroup1,
+                     data = train,
+                     prior = prior,
+                     verbose = F,
+                     pr=T,
+                     family = "hupoisson",
+                     rcov = ~idh(trait):units)

属性 pipematerialroup1 中有 5 个级别。我知道之前的事情一定是错的。我得到以下信息：

Error in MCMCglmm(pipenumyearnosameday ~ trait - 1 + at.level(trait, 1):lastyearnumnosameday + : prior$G has the wrong number of structures

我还需要说一件事：

myPredict <- predict(myTrain,newdata = myTest,type = "response",marginal = ~idh(trait):pipematerialgroup1,interval = "confidence") 

当我只是将 pipematerialgroup1 作为随机效应时，这是我的预测部分。如果结合 pipematerialgroup1 和 installyeargroup1，这个预测会是什么？

有人可以帮我获得优先权吗？非常感谢。