用於統計生存時間中的關係的高效計數器（R）

Question

生存時間是在大小爲1 x n的有序向量x中給出的數值時間測量。對於x中的兩個時間點，聯繫被定義爲x[i] == x[j]，其中如果所有中間時間都是關係，則排序強制i和j相鄰，i+1 == j或k指數分開。

鑑於x我需要一個有效的計數器，返回（1）唯一時間向量（無雙），（2）每次關係計數。

下面是我的判斷代碼。有沒有更有效的方法來做到這一點？

我們有一些時間對指數5,6和9和10

x = 1:10 
x[5] = 6 
x[10] = 9 

領帶indentifier關係：

x_new <- ties <- numeric() 
n <- length(x) 
i = 1; k=1 

while(i <= n){ 
    d = NA 
    if(i+1<=n) {d = x[i] - x[i+1]} 
    j = 0 

    while(d==0 & (i+j+1)<=n){ 
    j = j + 1 
    d = x[i] - x[i+j+1] 
    } 

    x_new[k] = x[i] 
    ties[k] = j + 1 
    i = i + j + 1 
    k = k + 1 
} 

給予期望的結果：

cbind(x_new,ties) 
    x_new ties 
[1,]  1 1 
[2,]  2 1 
[3,]  3 1 
[4,]  4 1 
[5,]  6 2 
[6,]  7 1 
[7,]  8 1 
[8,]  9 2 

Answer 1

只需使用rle

sapply(rle(sort(x)), cbind) 
#  lengths values 
#[1,]  1  1 
#[2,]  1  2 
#[3,]  1  3 
#[4,]  1  4 
#[5,]  2  6 
#[6,]  1  7 
#[7,]  1  8 
#[8,]  2  9