r语言 - 与dplyr和broom模型比较::glance:如果产生错误如何继续

Question

我想在数据集中运行每个变量作为使用r中的lme4包的单变量glmer模型。我想用dplyr/tidyr包准备数据，并使用扫把包(即do(glance(glmer…))组织每个模型的结果。我将非常感谢在这个框架内提供的帮助。我在R中不是那么好，但能够产生一个抛出错误的数据集，并且与我使用的数据具有相同的结构:

library(lme4)
library(dplyr)
library(tidyr)
library(broom)
Bird<-c(rep(c(0,0,0,0,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,0),10))
Stop<-c(rep(seq(1,10), 20))
Count<-c(rep(c(rep(c(1,2), each=10)), each=10))
Route<-c(rep(seq(1,10), each=20))
X1<-rnorm(200, 50, 10)
X2<-rnorm(200, 10, 1)
X3<-c(rep(c(0),200))#trouble maker variable
Data<-data.frame(cbind(Bird, Stop, Count, Route, X1, X2, X3))
Data%>%
  gather(Variable, Value, 5:7)%>%
  group_by(Variable)%>%
  do(glance(glmer(Bird~Value+Stop+(1+Stop|Route/Count), data=.,     family=binomial)))

最后一个变量产生错误，因此没有输出。如果发生这种情况，我想要的是在输出中产生NA值，或者直接跳过这个变量。我试过用'try'来吹过麻烦制造者变量:

do(try(glance(glmer(Bird~Value+Stop+(1+Stop|Route/Count), data=.,       family=binomial))))

，它做了，但仍然没有产生输出，因为它不能强制'try-error'到data.frame。不幸的是，没有tryharder函数。我试过一些对我有意义的if语句，但对计算机没有意义。我肯定我写得不对，但是如果我使用:

try(glance(glmer(Bird~Value+Stop+(1+Stop|Route/Count), data=., family=binomial)))->mod
if(is.data.frame(mod)){do(mod)}

我得到下标越界错误。非常感谢您提供的任何意见!

Answer 1

在调用glance之前使用tryCatch:

zz = Data %>%
  gather(Variable, Value, 5:7) %>%
  group_by(Variable) %>%
  do(aa = tryCatch(glmer(Bird~Value+Stop+(1+Stop|Route/Count), data=., 
                  family=binomial), error = function(e) data.frame(NA))) 

zz %>% 
  glance(aa)