在R中創建事件變量的時間

Question

親愛的所有人：我一直在嘗試創建一個時間到事件變量。事實上，前一段時間，我要求here尋求幫助。但是，我發現它並未完全達到我的目的。

以下是我的數據和我想創建「Time-to-event」的變量。

df2 = structure(list(Country = c("USA", "USA", "USA", "USA", "USA", 
"USA", "USA", "USA", "USA", "USA", "USA", "USA", "USA"), year = 2000:2012, 
    Event = c(0L, 0L, 1L, 0L, 0L, 0L, 1L, 0L, 1L, 0L, 0L, 0L, 
    0L), `**Time-to-event**` = c(0L, 1L, 2L, 0L, 1L, 2L, 3L, 
    0L, 1L, 0L, 1L, 2L, 3L)), .Names = c("Country", "year", "Event", 
"**Time-to-event**"), row.names = c(NA, -13L), class = "data.frame") 

Country year    Event  **Time-to-event** 
USA  2000    0   0 
USA  2001    0   1 
USA  2002    1   2 
USA  2003    0   0 
USA  2004    0   1 
USA  2005    0   2 
USA  2006    1   3 
USA  2007    0   0 
USA  2008    1   1 
USA  2009    0   0 
USA  2010    0   1 
USA  2011    0   2 
USA  2012    0   3 

，我建議使用下面的代碼，以創建時間 - 事件varaible

i1 <- with(df2, ave(Event, Country, FUN= 
     function(x) cumsum(c(TRUE, diff(x)<0)))) 
df2$Time_to_event <- with(df2, ave(i1, i1, Country, FUN= seq_along)-1) 

它運作良好，但與此代碼的問題是，它計算的情況下事件=連續多年。請參閱下面的例子：

Country year    Event  **Time-to-event** 
USA  2000    0   0 
USA  2001    0   1 
USA  2002    1   2 
USA  2003    0   0 
USA  2004    1   **1** 
USA  2005    1   **2** 
USA  2006    1   **3** 
USA  2007    0   0 
USA  2008    1   1 

相反，我想它對於該事件變量是在隨後的幾年中1例（0）給出的值爲零，而不是計數1的 需要明確的是，這是我希望看到「時間到事件」變量的方式。

Country year    Event  **Time-to-event** 
USA  2000    0   0 
USA  2001    0   1 
USA  2002    1   2 
USA  2003    0   0 
USA  2004    0   1 
USA  2005    1   2 
USA  2006    1   0 
USA  2007    1   0 
USA  2008    1   0 
USA  2009    0   0 
USA  2010    0   1 

Answer 1

您可以使用data.table如下：

require(data.table) 
setDT(dat)[,tte := seq.int(0,.N-1L), by = cumsum(Event)-Event] 

所以你最終獲得：

> dat 
    Country year Event **Time-to-event** tte 
1:  USA 2000  0     0 0 
2:  USA 2001  0     1 1 
3:  USA 2002  1     2 2 
4:  USA 2003  0     0 0 
5:  USA 2004  0     1 1 
6:  USA 2005  1     2 2 
7:  USA 2006  1     0 0 
8:  USA 2007  1     0 0 
9:  USA 2008  1     0 0 
10:  USA 2009  0     0 0 
11:  USA 2010  0     1 1 

爲什麼？

讓我們看看組件：

> dat[,.(Event, cumsum = cumsum(Event), run = cumsum(Event)-Event)] 
    Event cumsum run 
1:  0  0 0 
2:  0  0 0 
3:  1  1 0 
4:  0  1 1 
5:  0  1 1 
6:  1  2 1 
7:  1  3 2 
8:  1  4 3 
9:  1  5 4 
10:  0  5 5 
11:  0  5 5 

事件+ cumsum加起來建設運行的數量。按照這個順序進行分組使其工作。