0

我收到了一个奇怪的错误，例如 Dredge: MuMIn在某些情况下在 MuMIn 包中使用 dredge() 或 pdredge() 的奇怪错误。 

Error in while ((iComb <- iComb + 1L) < ncomb) { : 
  missing value where TRUE/FALSE needed
In addition: Warning message:
In iComb + 1L : NAs produced by integer overflow

当全局模型包含 31 个不同的项（环境的主效应、5 个基因的主效应、5 个基因-环境相互作用、10 个基因-基因相互作用和 10 个基因-基因-环境相互作用）时，就会出现此错误。31 是 dredge() 允许的最大项数。

该错误会在大约 30 小时后终止该过程，这可能已接近完成。当全局效应减少到 15 项（5 个基因的主效应，10 个基因-基因相互作用）时，没有错误。我的数据集没有任何缺失数据。

有谁知道这个错误是什么意思或如何解决它？

这是我的脚本：

# Full model has 31 terms
Full <- lm(DTH ~ Environment + PpdD1 + PpdB1 + PpdA1 + RhtB1 + RhtD1 + 
               PpdD1:PpdB1 + PpdD1:PpdA1 + PpdD1:RhtB1 + PpdD1:RhtD1 + PpdB1:PpdA1 + 
               PpdB1:RhtB1 + PpdB1:RhtD1 + PpdA1:RhtB1 + PpdA1:RhtD1 + RhtB1:RhtD1 + 
               PpdD1*PpdB1*Environment + PpdD1*PpdA1*Environment + PpdD1*RhtB1*Environment + 
               PpdD1*RhtD1*Environment + PpdB1*PpdA1*Environment + PpdB1*RhtB1*Environment + 
               PpdB1*RhtD1*Environment + PpdA1*RhtB1*Environment + PpdA1*RhtD1*Environment + 
               RhtB1*RhtD1*Environment, data=data)    
# Model selection using MuMIN package
options(na.action = "na.fail")
AllModels <- dredge(Full, rank="AIC") # Not run in parallel
AllModels
4

1 回答 1

0

dredge使用整数的位来指定变量的组合。事实证明，R 的整数被限制为 2^31 - 1，所以错误只发生在最后一轮有 31 个变量。目前，我已将当前版本（在 R-forge 上）中的允许限制更改为 30。在某个时候，我会寻找一种允许更多变量的解决方案，但这不是现在的优先事项。

于 2016-01-21T23:50:24.740 回答