Вейбулл с зависящим от времени исходом

Question

Я пытаюсь запустить модель Вейбулла для продольных измерений. То есть Вейбулла с зависящим от времени исходом. Я нашел следующие функции в библиотеке ("Эха") R. Я не уверен, если можно добавить cluster(ID) к этим функциям? Если ответ "нет", то как я могу соответствовать модели? В чем разница между этими функциями?

fit1=aftreg(Surv(start, stop, event) ~ x1+x2+x3+cluster(ID),dist="weibull"
    ,param = "lifeExp", data = df)

fit2=phreg(Surv(start,stop, event) ~ x1+x2+x3+cluster(ID),dist = "weibull" 
    ,data =df)


fit3=weibreg(Surv(start,stop, event) ~ x1+x2+x3+cluster(ID) ,data =nadf)