Comments (7)
drizopoulos
commented on August 16, 2024
Have you tried to use coxph instead of survreg?
——
Dimitris Rizopoulos
Professor of Biostatistics
Erasmus University Medical Center
The Netherlands
…________________________________
From: agstone90 <[email protected]>
Sent: Friday, August 28, 2020 3:15:20 AM
To: drizopoulos/JM <[email protected]>
Cc: Subscribed <[email protected]>
Subject: [drizopoulos/JM] “aeqSurv exception, an interval has effective length 0” (#19)
Hi Dr. Rizopoulos,
Like a few of the commentors below, I am having an issue with getting my joint model to run, receiving an error
"Error in aeqSurv(Y) :
aeqSurv exception, an interval has effective length 0"
My survival object uses the survreg function with weeks as my unit of time.
lmefit<-lme(marker~ns(week,3),random=~week|ID, data=Data,control=lmeControl(opt="optim")) survfit<-survreg(Surv(time,status)~covariate,data=Data.id,x=TRUE) fitjoint<-jointModel(lmefit,survfit,timeVar="week")
I haven't been able to get any of the fixes listed here<https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstackoverflow.com%2Fquestions%2F46988426%2Ferror-in-aeqsurvy-aeqsurv-exception-an-interval-has-effective-length-0&data=02%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C53c45118febc4b99f6e708d84aefd569%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C0%7C637341741220342748&sdata=l%2BiPngHbcuCSzGPV9Fu%2BHfbR6PRHu7DWxj7BfelRtG4%3D&reserved=0> to work. The issue doesn't seem to be one of 2 weeks vs 2.0 weeks. Entering "timefix=FALSE" into the survreg call gets an error for an unused argument, and adding
aeqSurv(survfit, tolerance = sqrt(.Machine$double.eps))
returns an error that survfit isn't a Surv object.
Any advice for how to get my model to run would be very appreciated!
Thanks!
—
You are receiving this because you are subscribed to this thread.
Reply to this email directly, view it on GitHub<https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fdrizopoulos%2FJM%2Fissues%2F19&data=02%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C53c45118febc4b99f6e708d84aefd569%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C0%7C637341741220352741&sdata=U5fr%2FmoKfsFgm2JCfXRIsanO13MjwKjRASSU2qxb0R4%3D&reserved=0>, or unsubscribe<https://eur01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FADE7TT5NNC6KUFZOCU44WYLSC4AKRANCNFSM4QNS2HDQ&data=02%7C01%7Cd.rizopoulos%40erasmusmc.nl%7C53c45118febc4b99f6e708d84aefd569%7C526638ba6af34b0fa532a1a511f4ac80%7C0%7C0%7C637341741220352741&sdata=wUfqs2mEp6T8JbBCooupub7BS9W6uZxETjGq1%2Bnawv4%3D&reserved=0>.
from jm.
agstone90
commented on August 16, 2024
Hi Dr Rizopoulos,
Thank you so much for responding! It does seem to be running now that I've switched to coxph. Is there a reason that one should work and not the other? Given a time-independent covariate, would both produce the same model?
from jm.
jinjinzhou
commented on August 16, 2024
Even tried coxph model, I still have this error
from jm.
rberch
commented on August 16, 2024
Hi @jinjinzhou, please, did you find a solution to this? I have run into a similar situation.
Thanks
from jm.
jinjinzhou
commented on August 16, 2024
Hi @jinjinzhou, please, did you find a solution to this? I have run into a similar situation. Thanks
For my situation, it is related to how I bin my longitudinal markers, i.e., some of them were measured on the same day or even more frequently. I filtered those when Tstart-stop < threshold
from jm.
rberch
commented on August 16, 2024
Hi @jinjinzhou, please, did you find a solution to this? I have run into a similar situation. Thanks
For my situation, it is related to how I bin my longitudinal markers, i.e., some of them were measured on the same day or even more frequently. I filtered those when Tstart-stop < threshold
Thank you @jinjinzhou ,
In my case, I am doing a simulation study. I simulate the longitudinal data and survival data (which has the ID, event status, and follow-up time) using simsurv.
And then when I fit the linear mixed model and coxph, they run without errors.
However, when I run the jointModel, that is when I get the error:
Error in aeqSurv(Y) :
aeqSurv exception, an interval has effective length 0
I have searched online for possible solutions but I have not found any.
from jm.
mariellabogoni
commented on August 16, 2024
@rberch Did you find a solution? I have the same problem. Thank you!
from jm.
Related Issues (20)
- Confusing point estimates for association parameters in current values + current slopes model HOT 5
- Potential bug in survfitJM.jointModel.R HOT 2
- error when applying weight function to ins()
- Error in aeqSurv(Y) : aeqSurv exception, an interval has effective length 0 HOT 1
- Query on getting back association parameter from data simulation in weibull-ph-agh
- Error in exp(coefs[1]) : non-numeric argument to mathematical function- fail to run CompRisk HOT 5
- Error in if (w_dif < specs$tol || iter == specs$maxiter) break : missing value where TRUE/FALSE needed
- Error in if (w_dif < specs$tol || iter == specs$maxiter) break : missing value where TRUE/FALSE needed
- Simulating data and GLMMPQL using JM
- issues in running survfitJM
- Error in aeqSurv(Y) : aeqSurv exception, an interval has effective length 0
- Error in d * log.hazard + log.survival : non-conformable arrays In addition: Warning message: In wk * rep(x$P, each = nk) : longer object length is not a multiple of shorter object length HOT 3
- Error in optim(thetas, opt.survPC, gr.survPC, method = "BFGS", control = list(maxit = if (it < : non-finite value supplied by optim HOT 3
- `survfitJM` expects newdata to be ordered by subject, and by time within each subject, else things go wrong
- Error in in optim(thetas, LogLik.splineGH, Score.splineGH, method = "BFGS", HOT 1
- ***caught segfault*** caused by Error in names(fit$coef) <- cname : 'names' attribute [4] must be the same length as the vector [2] HOT 5
- problem running fuction jointModel() in JM package HOT 1
- problem running fuction jointModel() in JM package with method = "Cox-PH-aGH" option HOT 1
- Question about the calculation of standard deviation in the "random effects parameterization"
- Extract fitted lme object from the jointModel HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from jm.