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.