PROC PHREG聚类

Question

我有一个眼睛的数据集，每位患者有2只眼睛。数据是聚类的，因为对于那一个患者，双眼的年龄，性别和种族等协变量是相同的。

到目前为止，这是我的代码：

PROC PHREG data = eyes covs(aggregate) plots(overlay)=(survival);
id patientid;
class Var1 (ref = "0") Gender (ref = "M") Ethnic Agegroup (ref = "0") / param = ref;
model TimeToTherapy*therapy(0) = Var1 Agegroup Gender Ethnic NumA1c/ ties=discrete rl;
strata patientid;
run;

然而，当我运行它时，我会得到错误和风险比，p值等。

然而，当我拿出strata patientid;声明时，一切正常。

是否有必要将此声明保留在那里？