我想知道是否有办法从R中的robu()
包中为robumeta
模型获得AIC或BIC(请参见下面的示例)?
library(robumeta)
fit <- robu(effectsize ~ binge + followup + sreport
+ age,data = get(data(hierdat)),study = studyid,var = var,model = "HIER",small = TRUE)
AIC(fit) # fails
BIC(fit) # fails
logLik(fit)# fails