4

以下で見つかった実行中model0に、「eval（expr、envir、enclos）のエラー：関数「Lag」が見つかりませんでした」というエラーメッセージが表示されました。このメッセージを投稿する前に、このフォーラムとWebをスキャンしましたが、関連する解決策が見つかりませんでした。Lag以下に示す他のモデル（1と2）は問題なく実行できる ため、私のモデルのエラーはおそらく機能によるものではないと思います。

私の主な動機は、説明変数とそのラグのリストをループしてGAMモデルを実行することでした。

library(quantmod)
library(gamair) 
library(mgcv) 
data(chicago) 

names(chicago)

varlist0 <- c("pm10median", "pm25median", "o3median", "so2median")

model0<- lapply(varlist0, function(x) {
  gam(substitute(death ~ s(time,bs="cr",k=200)+ s(tmpd,bs="cr") + Lag(i,0:4) , list(i = as.name(x))),family=quasipoisson,na.action=na.omit, data=chicago)

})

エラーメッセージのないGam+Lag：

model1<- gam(death ~ s(time,bs="cr",k=200)+ s(tmpd,bs="cr") + Lag(pm10median, 0:4),family=quasipoisson,na.action=na.omit, data=chicago)

ラグがあり、エラーメッセージがないLm：

hsb2 <- read.csv("http://www.ats.ucla.edu/stat/data/hsb2.csv")
varlist <- names(hsb2)[8:11]   
models <- lapply(varlist, function(x) {
    lm(substitute(read ~ Lag(i,0:4) , list(i = as.name(x))), data = hsb2)
})

このエラーの原因を特定できません。最初のモデルで何を間違えましたか？

4

1 に答える 1

3

正確な原因はわかりませんがmgcv、式を評価している環境と関係があります。substitute次のようなことを行うよりも安全なアプローチ：

varlist0 <- c("pm10median", "pm25median", "o3median", "so2median")
model0 <- lapply(varlist0,function(v) {
    f <- sprintf("death ~ s(time,bs='cr',k=200)+s(tmpd,bs='cr') + Lag(%s,0:4)",v)
    gam(as.formula(f),family=quasipoisson,na.action=na.omit,data=chicago)
})
于 2012-11-23T15:33:42.250 に答える