[R] multistate data w/out individual ID

Emilio A. Laca ealaca at ucdavis.edu
Wed Nov 6 22:06:25 CET 2013


Dear R help community,

I have a version of the data below generated by counting the number of individuals in each state every 2-3 time units for a series of cohorts (sets).
States are A, B, C, E and T.
It is known that the following transitions are possible A->B, B->C, C->E, B->T, C->T and E->T. 
The identity of the individuals in each set was not recorded.
I am interested in estimating the transition rates or probabilities and sojourn times as a function of covariates.
We assume that rates (probabilities) are independent of time.

Is there a package to estimate those quantities?
Is there a method to estimate those quantities?
Do you think a Bayesian approach integrating over all possible paths might work?
Everything I read about multi-state modeling packages seems to indicate that it is necessary to know the identity of individuals.

time,N,A,B,C,E,T,set,cov1,cov2
199,24,22,2,0,0,0,s1,1,4.65816667
201,24,22,2,0,0,0,s1,3,0.29922222
203,24,22,1,1,0,0,s1,10,9.570125
205,24,21,2,1,0,0,s1,6,9.94370139
207,24,20,3,1,0,0,s1,4,1.34693056
209,24,15,8,1,0,0,s1,2,1.20429167
212,24,14,9,1,0,0,s1,9,8.3008125
214,24,12,10,2,0,0,s1,15,9.14613194
216,24,10,11,3,0,0,s1,8,8.94250694
219,24,4,12,8,0,0,s1,5,0.97334722
221,24,3,12,9,0,0,s1,7,7.29372917
223,24,3,6,12,2,1,s1,17,1.28570139
225,24,1,7,8,5,3,s1,16,1.41599306
227,24,1,4,5,11,3,s1,18,7.07947222
229,24,1,3,2,15,3,s1,14,6.09359028
231,24,0,4,0,17,3,s1,12,4.70666667
233,24,0,4,0,17,3,s1,13,5.31577083
235,24,0,2,0,19,3,s1,11,4.62228472
199,24,19,5,0,0,0,s2,1,4.65816667
201,24,23,1,0,0,0,s2,3,0.29922222
203,24,22,2,0,0,0,s2,10,9.570125
205,24,22,2,0,0,0,s2,6,9.94370139
207,23,21,2,0,0,0,s2,4,1.34693056
209,23,19,3,1,0,0,s2,2,1.20429167
212,23,12,10,1,0,0,s2,9,8.3008125
214,23,10,12,1,0,0,s2,15,9.14613194
216,23,3,18,2,0,0,s2,8,8.94250694
219,23,0,16,7,0,0,s2,5,0.97334722
221,23,0,14,8,0,1,s2,7,7.29372917
223,23,0,5,14,2,2,s2,17,1.28570139
225,23,0,0,18,2,3,s2,16,1.41599306
227,23,0,2,13,2,6,s2,18,7.07947222
229,23,0,0,6,7,10,s2,14,6.09359028
231,23,0,0,3,8,12,s2,12,4.70666667
233,26,0,0,1,13,12,s2,13,5.31577083
235,28,0,0,0,16,12,s2,11,4.62228472
199,24,20,4,0,0,0,s3,1,4.65816667
201,24,22,2,0,0,0,s3,3,0.29922222
203,24,18,5,1,0,0,s3,10,9.570125
205,24,16,7,1,0,0,s3,6,9.94370139
207,24,13,10,1,0,0,s3,4,1.34693056
209,24,8,12,4,0,0,s3,2,1.20429167
212,24,5,14,5,0,0,s3,9,8.3008125
214,24,4,10,10,0,0,s3,15,9.14613194
216,24,2,9,12,1,0,s3,8,8.94250694
219,24,0,8,13,3,0,s3,5,0.97334722
221,24,2,6,5,11,0,s3,7,7.29372917
223,22,0,5,4,12,1,s3,17,1.28570139
225,24,1,3,1,17,2,s3,16,1.41599306
227,24,0,2,1,18,3,s3,18,7.07947222
229,25,0,2,0,20,3,s3,14,6.09359028
231,24,0,2,0,19,3,s3,12,4.70666667
233,24,0,1,1,19,3,s3,13,5.31577083
235,24,0,1,0,20,3,s3,11,4.62228472

Any comment will be appreciated.
Regards,

Emilio A. Laca, Professor                                                
One Shields Avenue, 2306 PES Bldg.            
Plant Sciences  Mail Stop 1                        ealaca at ucdavis.edu 
University of California                                 fax: (530) 752-4361     
Davis, California  95616                              voice: (530) 754-4083
						 		          mobile: (530) 220-5315
Include the Mail stop in the address or mail will be delayed.



More information about the R-help mailing list