search

drizopoulos / JM

Joint Models for Longitudinal & Survival Data under Maximum Likelihood
34 stars 7 forks source link

Error in aeqSurv(Y) : aeqSurv exception, an interval has effective length 0 #30

Closed ytripodis closed 2 years ago

ytripodis commented 3 years ago

I'm running the following in a dataset with 2641 subjects with follow-up in years and time to event between 0.01 and 4.9 years. The lme and the coxfit run with no problems but the jointModel returns

Error in aeqSurv(Y) : aeqSurv exception, an interval has effective length 0

Time-to-event Min. 1st Qu. Median Mean 3rd Qu. Max. 0.01369 1.14716 2.00685 2.05442 2.71869 4.86790

lme.NACC<-lme(Y ~ yrs + AGE + SEX + EDUC, random = ~ yrs| NACCID, data=longdata_complete3)

coxfit.NACC<-coxph(Surv(time=timetoevent, event=event) ~ AGE + SEX + EDUC, data=survdata_complete2, x=TRUE)

jointfit.NACC<-jointModel(lme.NACC,coxfit.NACC, timeVar="yrs", method = "piecewise-PH-aGH")

drizopoulos commented 3 years ago

I cannot replicate this error.

From: ytripodis @.> Sent: Thursday, April 15, 2021 2:03 AM To: drizopoulos/JM @.> Cc: Subscribed @.***> Subject: [drizopoulos/JM] Error in aeqSurv(Y) : aeqSurv exception, an interval has effective length 0 (#30)

I'm running the following in a dataset with 2641 subjects with follow-up in years and time to event between 0.01 and 4.9 years. The lme and the coxfit run with no problems but the jointModel returns

Error in aeqSurv(Y) : aeqSurv exception, an interval has effective length 0

Time-to-event Min. 1st Qu. Median Mean 3rd Qu. Max. 0.01369 1.14716 2.00685 2.05442 2.71869 4.86790

lme.NACC<-lme(Y ~ yrs + AGE + SEX + EDUC, random = ~ yrs| NACCID, data=longdata_complete3)

coxfit.NACC<-coxph(Surv(time=timetoevent, event=event) ~ AGE + SEX + EDUC, data=survdata_complete2, x=TRUE)

jointfit.NACC<-jointModel(lme.NACC,coxfit.NACC, timeVar="yrs", method = "piecewise-PH-aGH")

- You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHubhttps://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fdrizopoulos%2FJM%2Fissues%2F30&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C9ef15840cfda471294e208d8ffa1e6d7%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C0%7C637540418115309331%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=eID4%2BEvc3H6AwQPuQYD8M6jd%2Fh4Rr%2BE%2FgZL3wlzKZ9k%3D&reserved=0, or unsubscribehttps://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FADE7TT7YHRKGIWA6WVTB63LTIYUNBANCNFSM426L3MEQ&data=04%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C9ef15840cfda471294e208d8ffa1e6d7%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C0%7C637540418115319294%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C1000&sdata=9wPgl8gFmZAv%2Fq8EyE2wtq1%2Bi0chqmviG4UEWc610EA%3D&reserved=0.