将 Stata xtmelogit 转换为 R glmer (lme4) 包

编程问答 2022-05-27

问题描述

我尝试在 R 中重新计算准确发布的 Stata 代码。在第一步中，我从 Stata 导出了相同的数据集并将其导入到 R 中。尽管如此，我还是在与代码中的错误作斗争。 我做错了什么？

原始状态代码：

xtmelogit redpref1 c.incomedif c.incomedif#c.forpop forpop i.year if (brncntr==1) || country:

我在 R 中的方法是：

dataset <- dataset%>%
  filter(brncntr==1) %>%
  mutate(c.incomedif = factor(incomedif))%>%
  mutate(c.forpop = factor(forpop)) %>%
  mutate(i.year = as.integer(year))) 

library(lme4) 
    logit <- glmer(redpref1~ c.incomedif+ i.year|country,family=binomial,rueda4,nAGQ=0L)
    summary(logit)

这种方法会产生如下错误：

Error in if (ctrl$npt > (2 * n + 1)) warning("Setting npt > 2 * length(par) + 1 is not recommended.") : 
  missing value where TRUE/FALSE needed
In addition: Warning message:
In (n + 1L) * (n + 2L) : NAs produced by integer overflow
>

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

lme4mixed-modelsrrstata