测量干预效果的长度

机器算法验证 时间序列 生存 干预分析
2022-03-28 18:06:26

我进行了一项研究,参与者被随机分配到对照组或干预组,结果以事件发生时间数据的形式出现。虽然干预组的总体事件发生时间比对照组短(时序检验p<0.05),我有兴趣检测研究开始后是否有一段时间干预没有效果。更正式地说,我有兴趣检测最短时间t超过此干预组的风险函数值没有显着高于对照组的风险函数值。

我的一个想法是重新计算对数秩检验 p 值,忽略一段时间之前发生的所有事件t在干预组和对照组(左删失事件发生时间数据)。对一系列重复此过程t值,我可以确定最小值t除此以外没有显着的干预效果。

我想知道是否有更多标准方法来确定时间t超出此范围,干预对事件发生的时间没有可察觉的影响。

2个回答

鉴于您似乎有一组您随时间关注的人,其中一些人得到治疗,而另一些人则没有,您可以在差异分析中进行差异分析。你可以像这样运行回归

yit=β1(treati)+β2(interventiont)+β3(treatiinterventiont)+ϵit
在哪里treati是一个假人是否个人i在治疗组,interventiont是治疗后时期的虚拟变量,两者之间的相互作用捕捉到治疗效果β3.

如果您现在要估计淡出时间,请改为估计

yit=γ=0mβγ(treatmentit)+ηit
在哪里treatmentit是一个虚拟变量,如果单独,则等于 1i在治疗组和时间t在治疗日期或之后。这估计了第一个方程,但m处理的滞后,您可以选择从处理开始到样本期结束的时间段数。

然后β0是干预日期的治疗效果,β1是干预日期后第一个时期的干预效果,以此类推。这种方法的好处是

  • 它很容易在任何统计软件中实现(您只需要创建假人并运行回归)
  • β0,...,βm系数将具有标准误差和置信区间,您可以使用它们来查看干预停止产生效果的时间(治疗滞后)
  • β0,...,βm系数将为您提供对后续时期干预幅度的估计

如果您还有其他控制变量,例如研究参与者的特征Xit您可以轻松地将它们包含在回归中,

yit=γ=0mβγ(treatmentit)+Xitρ+ηit
这不会影响干预效果的估计(因为识别来自治疗组和对照组之间的组差异),但它有助于减少残余方差,从而提高精度。

干预检测http://www.unc.edu/~jbhill/tsay.pdf和其他地方可以在有或没有用户建议的干预变量的情况下使用。在任何一种情况下,都需要处理可能存在的任何自投影过程,即 ARIMA 结构。识别 ARIMA 结构和对任何用户建议变量和“新干预系列”的响应都需要一些试验和错误,因为搜索过程/方法是确定是否以及何时出现这种等待被发现的干预。此搜索过程可以通过软件解决/辅助,但需要确认软件包含任何所需的影响,即用户指定的可能预测变量系列/变量,包括它们的滞后结构,并且 ARIMA 过程识别阶段没有受到干预的损坏/影响那是有待发现的,这当然会产生有害影响。

长度将是用户指定的干预变量的点之间的差异......如果它存在,或者是发现/发现/未屏蔽的两个识别的 Level Shift/Step Shift 变量之间的差异。