Re: [R] Joint modelling of survival data

2012-01-05 Thread ryusuke
Hi Dr Terry, Thank you for your reply. Step(1) - Lets assume Generalized Poisson model (GLM) as basic model where constant hazards ratio as time goes by. Below are two correlated GLM. X_ij = Poisson( lambda_1 = \gamma * \alpha_i * \delta_j ) Y_ij = Poisson( lambda_2 =

Re: [R] Joint modelling of survival data

2012-01-03 Thread Terry Therneau
--- begin included messge --- Assume that we collect below data : - subjects = 20 males + 20 females, every single individual is independence, and difference events = 1, 2, 3... n covariates = 4 blood types A, B, AB, O http://r.789695.n4.nabble.com/file/n4245397/CodeCogsEqn.jpeg ?m = hazards ra

[R] Joint modelling of survival data

2011-12-30 Thread ryusuke
Assume that we collect below data : - subjects = 20 males + 20 females, every single individual is independence, and difference events = 1, 2, 3... n covariates = 4 blood types A, B, AB, O http://r.789695.n4.nabble.com/file/n4245397/CodeCogsEqn.jpeg λm = hazards rates for male λn = hazards rate