0

我想使用 R 中的 lme4 包将数据集中的每个变量作为单变量 glmer 模型运行。我想用 dplyr/tidyr 包准备数据，并用 broom 包组织每个模型的结果（即（一目了然（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 值，或者只是跳过该变量。我尝试使用“尝试”来克服麻烦制造者变量：

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

它确实如此，但仍然没有产生输出，因为它不能将“尝试错误”强制到 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)}

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

4

1 回答 1

1

tryCatch在调用之前使用glance：

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)
于 2016-08-19T17:17:29.983 回答