R TIME-TO-EVENT

Question

Я новичок в R (начал учить себя на прошлой неделе) и новичок в этом форуме.

Я работаю над набором данных, где мне нужно определить время-событие. У меня есть три переменные (лекарство, пациент и дата), и мне нужно разработать разницу во времени для каждого пациента, когда они переключают наркотики. Следует также отметить, что пациент может перейти от препарата А к препарату В и обратно к препарату А, поэтому последовательность даты имеет важное значение.

я в настоящее время эти данные являются dataframe называемые «данные»

data <-data.frame(Drug=sample(c("Drug A","Drug B", "Drug C", "Drug D", "Drug E"),100, 
replace=TRUE),PAT=rep(1:10,each=10),Date=sample(seq(as.Date('2003/01/01'), 
as.Date('2016/01/01'), by="day"), 100)) 
data<-data[order(data$PAT,data$Date),] 

Любая помощь будет большим, и я сожалею, что не было вещи ясно или отклонился от протокола данного форума ...

Требуемый результат - это доля пациентов, которые перешли от своего первого препарата к своему второму препарату (а также от 2-го, 3-го, 4-го и 4-го до 5-го) в каждый момент времени ... время 0 (когда существует первая запись каждого пациента), 100% пациентов не переключились, следующая строка будет днем ​​2 (где 99% пациентов могут не иметь Witched), затем 3-й день, то день 4 и т.д.

Приветствия

Answer 1

.... ответ на «часть Deux», если я правильно понял намерение. Примечание. Я уменьшил количество дат от 100 до 10, чтобы получить последовательности одновременно.

library(TraMineR) 

data <- data.frame(Drug=sample(c("Drug A","Drug B", "Drug C", "Drug D", "Drug E"),100,replace=TRUE), 
       PAT=rep(1:10,each=10), 
       Date=sample(seq(as.Date('2003/01/01'),as.Date('2016/01/01'), by="day"), 10)) 

data <- data[order(data$PAT,data$Date),] 

data2 <- reshape(data, idvar="PAT", timevar="Date", direction="wide") 

data.seq <- seqdef(data2,2:11) 

seqiplot(data.seq) 
seqdplot(data.seq) 

Последняя функция сюжет дает ....