Tuesday, 3 September 2013

frailtyPenal for two recurrent risks and memory not mapped

frailtyPenal for two recurrent risks and memory not mapped

I am a beginner in R and I am using frailtyPenal() to model a joint
frailty model with two competing risks. Firstly, my two competing risk
events may be understood as recurrent, but in frailtyPenal it is supposed
that one competing risk is recurrent and the other one is terminal. So, I
do not know if it might be a problem. Secondly, after running 3 hours (on
mac) this function it has appeared a message " caught segfault" cause
memory not mapped" I do not know if it is either because my database is
large, or because there is something wrong in the function... This is the
function I am running:
fit.joint.05_males_services<-frailtyPenal(Surv(.t0,.t,delta.unstable)~cluster(id)+u.rate+as.factor(h.skill)+as.factor(m.skill)+as.factor(non.manual)
+as.factor(municipio)+as.factor(spanish.speakers)+
as.factor(no.spanish.speaker)+as.factor(Aged.16.19)+as.factor(Aged.20.24)+as.factor(Aged.25.29)+as.factor(Aged.30.34)+as.factor(Aged.35.39)+
as.factor(Aged.40.44)+as.factor(Aged.45.51)+as.factor(older61)+
as.factor(responsabilities)+terminal(delta.stable),formula.terminalEvent=~
u.rate+as.factor(h.skill)+as.factor(m.skill)+as.factor(municipio)+as.factor(spanish.speakers)+
as.factor(no.spanish.speaker)+as.factor(Aged.16.19)+as.factor(Aged.20.24)+as.factor(Aged.25.29)+as.factor(Aged.30.34)+as.factor(Aged.35.39)+as.factor(Aged.40.44)+as.factor(Aged.45.51)+as.factor(older61)+
as.factor(responsabilities),
data=cr_05_males_services,n.knots=8,kappa1=1000,kappa2=1000,maxit=1000,
Frailty=TRUE,joint=TRUE, recurrentAG=TRUE)
Thank you very much

No comments:

Post a Comment