[英]Defining a wrapper function in glmulti
我是 R 腳本的新手,我想使用包裝器 function 運行 glmulti,以獲得帶有偏移的負二項式系列。 由於收到錯誤消息,因此我想就以下代碼提出建議。 變量 a、b 和 c 是因子和協變量。 變量 e 用作偏移量。
glmnb.glmulti <- function(formula, data, offset, ...)
glm.nb(formula, offset=offset, data=data, ...)
selection<- glmulti(y~ a + b + c + d,level=2,
fitfunction=glmnb.glmulti, offset=log(e), crit="aic")
Initialization...
Error in model.frame.default(formula = formula, data = data, offset = offset, :
invalid type (closure) for variable '(offset)'
任何建議將不勝感激。
十分感謝! 我也嘗試了以下方法,它也有效
eqt <- as.formula(y~ a + b + c + d)
glmnb <- 函數(公式,數據,權重,...) glm.nb(公式,數據=數據,權重=偏移(log(e)),...)
選擇 <- glmulti(eqt, data=data, fitfunction=glmnb, crit = "aic")
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.