有没有类似Stata的metobit命令的R函数？

Question

我正在尝试在 Rstudio 中运行多级混合效果 Tobit 回归，但我找不到方法。

我的数据集包括 1989 年至 2018 年向意大利地区的转移（以欧元为单位）以及其他几个经济和政治变量。因变量是某一年某个地区收到的转移支付金额。该分布在 0 处有一个峰值，968 个观测值中有 312 个观测值假设值为 0。因此，我想通过运行 Tobit 模型来解释这个在 0 处的峰值。此外，我想在区域层面应用随机效应。 我找不到在 Rstudio 中适合这种多级混合效果 Tobit 模型的方法。

我发现在 Stata 中有一个

metobit

命令可以执行此操作，我想在 Rstudio 中复制它。

我尝试了几种包和功能的组合，但似乎都不起作用。 基本上我想在 Rstudio 中复制以下 Stata 命令。

metobit Total_fund GDP_capita || region:, ll(0)

我尝试使用

gamlss

包。下面是我运行的代码示例。

model1 <- gamlss(Total_fund ~ GDP_capita,
                 sigma.formula = ~ 1,  
                 data = data_regression_2_ventile,
                 family = TF)  

model1$method$REML <- TRUE 
model1$method$REMLscope <- list(bsc = "none")  
model1$method$REMLconstraint <- NULL  

model1$method$na.action <- na.exclude  
model1$method$REMLscope$region <- "~ 1"  

发生的情况是我收到错误（

Warning message: In RS() : Algorithm RS has not yet converged

）并且代码的第二部分不会更改回归输出。

我还看到使用

glmmTMB

包可能会有帮助，但我收到了几个错误。 警告消息：

1: In checkMatrixPackageVersion() :
  Package version inconsistency detected.
TMB was built with Matrix version 1.5.4
Current Matrix version is 1.5.1
Please re-install 'TMB' from source using install.packages('TMB', 
        type = 'source') or ask CRAN for a binary version of 'TMB' matching 
        CRAN's 'Matrix' package
2: In checkDepPackageVersion(dep_pkg = "TMB") :
  Package version inconsistency detected.
  glmmTMB was built with TMB version 1.9.3
  Current TMB version is 1.9.4
  Please re-install glmmTMB from source or restore original ‘TMB’ package   
  (see '?reinstalling' for more information)
> 

我尝试重新安装，但没有帮助。

关于如何在 Rstudio 中拟合多级混合效果 Tobit 模型有什么建议吗？

Answer 1

我也有同样的问题。 有 R 相当于 STATA 的 metobit 吗？