时间戳

首页 > 解决方案 > 在函数调用内部构造的 glms 上的 step/stepAIC

r - 在函数调用内部构造的 glms 上的 step/stepAIC

问题描述

我遇到了与可见性/环境相关的问题。简而言之，在函数内部构造的 glms 不能使用 step/stepAIC 来简化：

foo = function(model) {
  m = glm(y~x, family=model$family, data = dframe)
  return(m)
}

y = rbinom(100, 1, 0.5)
x = y*rnorm(100) + rnorm(100)
dframe = data.frame(y, x)
m = glm(y~x, family='binomial', data = dframe)
m2 = foo(m)
library(MASS)
summary(m2)
print(m2$family)
m3 = stepAIC(m2, k = 2)

这会导致以下错误：

Error in glm(formula = y ~ 1, family = model$family, data = dframe) : 
  object 'model' not found

尽管 m2 看起来很合身并且家庭已定义，但这仍然存在。对不起，如果这个例子有点做作。

标签： rfunctionscope

解决方案

找到了解决办法——原来的glm需要用do.call来构造。

foo = function(model) {
  form.1<-as.formula(y ~ x)
  dat = model$data
  fam = model$family
  m <- do.call("glm", list(form.1, data=dat, family=fam))
  ##m = glm(y~x, family='binomial', data = model$dframe)
  return(m)
}

y = rbinom(100, 1, 0.5)
x = y*rnorm(100) + rnorm(100)
dframe = data.frame(y, x)
m = glm(y~x, family='binomial', data = dframe)
m2 = foo(m)
library(MASS)
summary(m2)
print(m2$family)
m3 = stepAIC(m2, k = 2)

推荐阅读