Hello, I have a cohort with approx 1,200 patients at the ages of 30-65 that had their first myocardial infarction during 1992: · They were in a follow up until 2005. · About 400 of them died during this period of time (right censored) · Each one of them had up to 4 mi recurrent events. I am using the semi-parametric model in order to assess the relationship of predictors to rate of occurrence allowing for multiple events per subject. The exposure variable ("ESHKOL") is the socio-economic of the neighborhood while adjusting for additional variables such as age, sex, etc. In addition each patient has geo- statics -code for his neighborhood ("stat_unit"). I ran an Andersen Gill model: Ag<-coxph(Surv(TStart,TStop,Status)~ ESHKOL+VAR1+ VAR2+cluster(id),data=recmi) My questions are: · I would like to check the unique effect of the "stat_unit" variable. How can I write it in code? · In your opinion does the AG model is the best model for this type of data? Or should I use models that consider the order of the events such as PWP or WLW? Thank you, Avi. [[alternative HTML version deleted]]
Hello, I have a cohort with approx 1,200 patients at the ages of 30-65 that had their first myocardial infarction during 1992: ? They were in a follow up until 2005. ? About 400 of them died during this period of time (right censored) ? Each one of them had up to 4 mi recurrent events. I am using the semi-parametric model in order to assess the relationship of predictors to rate of occurrence allowing for multiple events per subject. The exposure variable ("ESHKOL") is the socio-economic of the neighborhood while adjusting for additional variables such as age, sex, etc. In addition each patient has geo- statics -code for his neighborhood ("stat_unit"). I ran an Andersen Gill model: Ag<-coxph(Surv(TStart,TStop,Status)~ ESHKOL+VAR1+ VAR2+cluster(id),data=recmi) My questions are: ? I would like to check the unique effect of the "stat_unit" variable. How can I write it in code? ? In your opinion does the AG model is the best model for this type of data? Or should I use models that consider the order of the events such as PWP or WLW? Thank you, Avi. -- View this message in context: http://r.789695.n4.nabble.com/recurrent-events-tp2318296p2318296.html Sent from the R help mailing list archive at Nabble.com.
Hello everyone, I have a question about using "frailtypack" to run recurrent events data. I have a file with 1,580 events. The file structure fits the andersen gill model. When I run the file with all records, the R stops and shuts down. Once I run only part of the file, about 400 events, R runs without a problem. Any idea why it happens? In addition when I run a "Nested frailty model" with the below code, R stops and shuts down as well: cgd<-frailtyPenal(Surv(TStart, TStop, Status)~cluster(Center)+subcluster(ID)+ Treatment,data=cgd.ag, n.knots=8,kappa1=50000) I would greatly appreciate it if anyone can help me out. Thank you very much, Avi -- View this message in context: http://r.789695.n4.nabble.com/recurrent-events-tp2329539p2329539.html Sent from the R help mailing list archive at Nabble.com.