我与写调用使用get()
为什么无法在gamm中找到对象,但它在lm中工作?
假设我有这个数据帧的变量的R功能挣扎:
mydat = data.frame(y = rnorm(100),
x = rnorm(100),
day = sample(90:260, 100, replace = T),
r1 = sample(seq(2008,2015,1), 100, replace = T),
r2 = sample(letters, 100, replace = T),
r3 = sample(letters, 100, replace = T))
,我想编写一个返回GAMM模型的摘要函数,例如像这样
gamm_summary = function(data, response = "y"){
require(mgcv)
gamm_model = gamm(get(response) ~ s(day),
random = list(r1=~1, r2=~1, r3 =~1), data = data, method = "REML")
summary(gamm_model$gam)
}
gamm_summary(mydat)
为什么这给我的错误:
项Error in get(response) : object 'response' not found
但以下几项工作:
lm_summary = function(data, response = "y") {
lm_model = lm(get(response) ~ x, data = data)
summary(lm_model)
}
lm_summary(mydat)
问:为什么无法得到我GAMM功能工作,我怎么可以重写,以便它能发挥的作用?
尝试使用此功能'?reconulate' –