Dear R users,

 

I'd like to fit different longevities for males and females (痠 vs 番 below) with a Weibull interval-censored survival model.

Unfortunately more than 50% of the individuals in the experiment died before their sex could be assigned.

Hence I was wondering about fitting the proportion of females in the model ( "prop" below) with a likelihood function such as:

L(prop, 番, 痠, shape)=prop*L(番, shape)*(1-prop)L(痠, shape).

 

I tried to follow the Cauchy example, but I realized that I don't have a proper density function to be passed on to the survreg.fit function.

Does anyone has an advice on how to implement that in R?

 

Thank you very much!

 

Best regards,

 

N. Rode  

 

CEFE (Centre d'Ecologie Fonctionnelle et Evolutive)
UMR 5175 CNRS, 1919 route de Mende
34293 Montpellier, France


	[[alternative HTML version deleted]]

