0

我正在处理一个包含 184,113 名参与者赌博行为的纵向数据的大型数据集。该数据基于对赌博运营商内电子赌博行为的完整跟踪。赌博行为数据按月汇总，共 70 个月。我有一个用于分隔参与者的 ID 变量、一个时间变量（月份），以及许多赌博行为变量，例如给定月份的活跃天数、给定月份的投注、给定月份的总损失等。参与者的时间不同一直在积极赌博。一位参与者可能在第 2、3、4 和 7 个月赌博，另一位参与者在 3、5 和 7 个月赌博，第三位参与者在 23、24、48、65 个月等。

我试图在 glmmTMB 中运行一个负二项式 2 截断模型，我想知道包如何处理缺少 0。我有关于赌博行为的纵向数据，每个月玩的天数（总共 70 个月）。该变量可以取 1-31 之间的值（取决于月份），没有 0。数据集中不存在参与者的月份为 0。仅由两个参与者构成数据的示例：

# Example variables and data frame in long form 

  # Includes id variable, time variable and example variable 

id <- c(1, 1, 1, 1, 2, 2, 2) 

time <- c(2, 3, 4, 7, 3, 5, 7) 

daysPlayed <- c(2, 2, 3, 3, 2, 2, 2) 

dfLong <- data.frame(id = id, time = time, daysPlayed = daysPlayed)

我的问题：如何指定 glmmTMB 中截断的位置？它默认为0吗？我想截断 0 并运行以下代码（我将比较模型，第一个是简单的无条件模型）：

DaysPlayedUnconditional <- glmmTMB(daysPlayed ~ 1 + (1 | id), dfLong, family = truncated_nbinom2)

它会成功吗？

4

1 回答 1

0

从 Ben Bolker 到 r-sig-mixed-models@r-project.org：“我不是 100% 清楚你的问题，但是：glmmTMB 只进行零截断，而不是 k>0 的 k-截断，即你只能指定模型 Prob(x==0) = 0 Prob(x>0) = Prob(NBinom(x))/Prob(NBinom(x>0)) （可怕的符号，但希望你明白）“

于 2021-01-20T08:25:10.843 回答