簡體   English   中英

在 glmulti 中定義一個包裝器 function

[英]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.

 
粵ICP備18138465號  © 2020-2024 STACKOOM.COM