如何使用Stata&x27对标准错误进行聚类；s stcrprep

我正在使用Stata对存在竞争风险的累积发病率函数进行建模。每一个咒语都以一个事件（“失败”）的发生而结束，但是受试者有可能在之后进入一个新的咒语。因此，咒语不是独立的，而是嵌套在主体中

标准Stata命令stcrreg可以通过对主题级别聚集的标准错误进行建模来处理此结构。然而，我的数据集非常庞大（超过300万次观测），计算时间也非常长。其他用户建议使用用户编写的程序stcrprep，该程序还具有其他功能

---

有人知道如何使用stcrprep对标准错误进行集群吗

我从

stcrprep

的帮助文件中找到了解决方案，并在命令编写者的演示文稿中进行了解释（）

集群标准错误可以很容易地实现如下

使用

stcrprep

为生存分析准备数据

使用

stset

创建正确的格式，将权重指定为pweights。将它们指定为iweights将使稳健和集群稳健标准误差的估计变得不可能

使用

stcox

执行竞争风险回归，并指定所需的

，vce（）

选项。请注意，通过这种数据准备方式，将产生与

stcrreg

完全相同的结果

对于感兴趣的读者，以下是一个在存在竞争性心血管疾病和其他死亡原因的情况下，通过审查后的观察，以及聚集在治疗中的患者，估计癌症风险的示例：

use "http://www.stata-journal.com/software/sj4-2/st0059/prostatecancer", clear
stset time, failure(status=1,2,3) id(id) //specify all possible events
stcrprep, events(status) keep(treatment age therapy) trans(1) //trans creates failcode (event of interest)
gen event = failcode == status
stset tstop [pw=weight_c], failure(event==1) enter(tstart)
stcox treatment age, vce(cluster therapy)