用FIML和R中的bootstrap拟合中介路径模型中的负二项式

我一直在分析R中中介路径模型的纵向数据。但我的结果是高度正偏态的，所以正则回归不起作用（偏态太大，无法拟合泊松分布）

Lavan不支持负二项式（我与Lavan小组进行了核对，Lavan的开发人员回答了问题）

我已经尝试过将结果二分法并使用拉万。但是Lavan不能将FIML用于二进制选项（仅成对），因此我以这种方式失去了1/3的参与者，而且我也不喜欢因为二分法而失去信息

我想继续使用引导，而不是手动检查a、b、c和c'

我的合作者和我正在努力避免多重插补

我没有MPLU，我想继续使用R

任何意见都将不胜感激

以下是我使用Lavan的模型：

p16a18c24 <- ' 
# direct effect
myOutcome ~ c*myPredictor
# mediator
myMediator1 ~ a1*myPredictor
myOutcome ~ b1*myMediator1
myMediator2 ~ a2*myPredictor
myOutcome ~ b2*myMediator2
myMediator3 ~ a3*myPredictor
myOutcome ~ b3*myMediator3
# variances: any variable that was listed as DV
myOutcome ~~ myOutcome 
myMediator1 ~~ myMediator1
myMediator2 ~~ myMediator2
myMediator3 ~~ myMediator3
# covariances: mediators who are (positively) correlated with each other 
myMediator1 ~~ myMediator2
myMediator2 ~~ myMediator3
myMediator3 ~~ myMediator1
# intercepts: *1, only for DVs 
myOutcome ~ 1 
myMediator1 ~ 1
myMediator2 ~ 1
myMediator3 ~ 1
# indirect effects (a*b)
ab1 := a1*b1
ab2 := a2*b2
ab3 := a3*b3
# total effect
total := c + (a1*b1) + (a2*b2) + (a3*b3)
'
fit.p16a18c24 <- sem(p16a18c24,data=m,missing="fiml",se="bootstrap")
summary(fit.p16a18c24,fit.measures=TRUE,standardized=TRUE,rsquare=TRUE)

p16a18c24