Dear
R users,



I
created a matrix that tells me the first day of use of a category by
id.



#Calculate
time difference
test$tdiff<-as.numeric(difftime(as.Date("2002-09-01"), test$ftime, units = 
"days"))



#
obtain the index date per person and dcategory
index.date.test<-tapply(test$tdiff,
list(test$id, test$rcat), max)



Nonetheless,
at the moment I think will be more useful to create a column in my
data that tells me which row is the index date. 




Something
like:



ti<-function(x){
        ifelse(x==max(x),
"i", "n") # x = test$tdiff
}



tapply(test$tdiff,
list(test$rcat, test$id), FUN=ti)



I
have been testing different things for few days but I am in a loop 
and I do not see my mistake.



It
should be simple but I can't get it



Bellow
a test data



Thanks in advance for your time,
Jose
Back ground info:  I want to use the index.date to obtain information from 
other df for every id six month previous the index date. 
Then I should normalize the ftime to a common time frame and look form patterns 
in that time frame.
(Do not know yet how I will do it. )






###
test data ####



structure(list(id = c(1L, 1L, 1L, 46L, 80L, 80L, 80L, 80L, 88L, 
160L, 179L, 179L, 179L, 179L, 179L, 179L, 192L, 192L, 192L, 204L, 
204L, 204L, 204L, 205L, 211L, 233L, 233L, 272L, 272L, 272L, 272L, 
309L, 309L, 309L, 310L, 310L, 314L, 314L, 315L, 316L, 320L, 320L, 
320L, 320L, 324L, 324L, 324L, 329L, 329L, 339L, 354L, 354L, 354L, 
357L, 358L, 359L, 364L, 366L, 377L, 377L, 377L, 377L, 377L, 377L, 
377L, 377L, 377L, 377L, 377L, 377L, 379L, 383L, 383L, 387L, 387L, 
391L, 395L, 398L, 401L, 401L, 401L, 401L, 401L, 407L, 407L, 407L, 
409L, 414L, 414L, 414L, 434L, 434L, 434L, 437L, 437L, 437L, 437L, 
437L, 439L, 439L, 439L, 439L, 442L, 443L, 450L, 452L, 452L, 459L, 
459L, 468L, 472L, 472L, 472L, 478L, 478L, 484L, 484L, 484L, 484L, 
484L, 486L, 486L, 486L, 487L, 487L, 487L, 487L, 487L), ftime = 
structure(c(11761, 
11824, 11925, 11852, 11814, 11814, 11929, 11929, 11902, 11857, 
11779, 11779, 11807, 11841, 11871, 11899, 11831, 11894, 11925, 
11761, 11801, 11843, 11905, 11832, 11877, 11838, 11901, 11783, 
11783, 11818, 11850, 11750, 11782, 11905, 11852, 11877, 11852, 
11922, 11855, 11838, 11845, 11878, 11901, 11927, 11795, 11817, 
11837, 11901, 11928, 11853, 11751, 11751, 11877, 11922, 11760, 
11914, 11857, 11912, 11752, 11752, 11785, 11785, 11825, 11825, 
11862, 11862, 11891, 11891, 11926, 11926, 11919, 11907, 11907, 
11842, 11873, 11842, 11922, 11865, 11782, 11829, 11858, 11888, 
11912, 11750, 11803, 11897, 11871, 11787, 11787, 11787, 11764, 
11817, 11882, 11778, 11808, 11863, 11894, 11918, 11771, 11817, 
11851, 11907, 11799, 11766, 11794, 11765, 11828, 11788, 11884, 
11897, 11810, 11852, 11922, 11810, 11846, 11801, 11835, 11859, 
11891, 11922, 11771, 11884, 11925, 11765, 11765, 11801, 11843, 
11892), class = "Date"), rcat = structure(c(1L, 1L, 1L, 1L, 1L, 
2L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 
1L, 1L, 1L, 1L, 1L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 2L, 3L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 
1L, 1L, 1L, 2L, 2L, 2L, 1L, 2L, 1L, 1L, 1L), .Label = c("ICS", 
"LABA", "MCSs"), class = "factor"), tdiff = c(170, 107, 6, 79, 
117, 117, 2, 2, 29, 74, 152, 152, 124, 90, 60, 32, 100, 37, 6, 
170, 130, 88, 26, 99, 54, 93, 30, 148, 148, 113, 81, 181, 149, 
26, 79, 54, 79, 9, 76, 93, 86, 53, 30, 4, 136, 114, 94, 30, 3, 
78, 180, 180, 54, 9, 171, 17, 74, 19, 179, 179, 146, 146, 106, 
106, 69, 69, 40, 40, 5, 5, 12, 24, 24, 89, 58, 89, 9, 66, 149, 
102, 73, 43, 19, 181, 128, 34, 60, 144, 144, 144, 167, 114, 49, 
153, 123, 68, 37, 13, 160, 114, 80, 24, 132, 165, 137, 166, 103, 
143, 47, 34, 121, 79, 9, 121, 85, 130, 96, 72, 40, 9, 160, 47, 
6, 166, 166, 130, 88, 39)), .Names = c("id", "ftime", "rcat", 
"tdiff"), row.names = c(11L, 4L, 13L, 25L, 39L, 41L, 35L, 44L, 
54L, 57L, 96L, 98L, 88L, 107L, 80L, 77L, 118L, 136L, 124L, 146L, 
150L, 157L, 153L, 169L, 196L, 210L, 214L, 225L, 230L, 221L, 222L, 
258L, 266L, 281L, 311L, 324L, 333L, 334L, 358L, 372L, 400L, 419L, 
423L, 434L, 439L, 437L, 443L, 479L, 465L, 496L, 517L, 516L, 519L, 
525L, 539L, 598L, 606L, 634L, 658L, 649L, 637L, 655L, 640L, 644L, 
645L, 636L, 647L, 646L, 639L, 654L, 665L, 673L, 680L, 701L, 688L, 
712L, 737L, 738L, 784L, 766L, 785L, 753L, 773L, 799L, 791L, 808L, 
818L, 826L, 821L, 820L, 838L, 830L, 837L, 841L, 840L, 844L, 850L, 
845L, 886L, 875L, 887L, 868L, 899L, 912L, 915L, 931L, 929L, 934L, 
939L, 957L, 988L, 975L, 981L, 1015L, 1003L, 1043L, 1051L, 1056L, 
1034L, 1031L, 1073L, 1068L, 1065L, 1079L, 1101L, 1092L, 1089L, 
1096L), class = "data.frame")
 

                                          
        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to