survival 1.1.33

A high-performance survival analysis library written in Rust with Python bindings
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86

use pyo3::prelude::*;
#[pyfunction]
pub fn tmerge(
    id: Vec<i32>,
    time1: Vec<f64>,
    newx: Vec<f64>,
    nid: Vec<i32>,
    ntime: Vec<f64>,
    x: Vec<f64>,
) -> Vec<f64> {
    let n1 = id.len();
    let n2 = nid.len();
    let mut result = newx;
    let mut k = 0;
    let mut current_id = -1;
    let mut csum = 0.0;
    let mut has_one = false;
    for i in 0..n1 {
        if id[i] != current_id {
            current_id = id[i];
            csum = 0.0;
            has_one = false;
            while k < n2 && nid[k] < current_id {
                k += 1;
            }
        }
        let start_time = time1[i];
        let mut local_k = k;
        while local_k < n2 && nid[local_k] == current_id && ntime[local_k] <= start_time {
            csum += x[local_k];
            has_one = true;
            local_k += 1;
        }
        if has_one {
            result[i] = if result[i].is_nan() {
                csum
            } else {
                result[i] + csum
            };
        }
    }
    result
}
#[pyfunction]
pub fn tmerge2(id: Vec<i32>, time1: Vec<f64>, nid: Vec<i32>, ntime: Vec<f64>) -> Vec<usize> {
    let n1 = id.len();
    let n2 = nid.len();
    let mut result = vec![0; n1];
    let mut k = 0;
    for i in 0..n1 {
        let current_id = id[i];
        let start_time = time1[i];
        result[i] = 0;
        while k < n2 && nid[k] < current_id {
            k += 1;
        }
        let mut last_valid = 0;
        let mut local_k = k;
        while local_k < n2 && nid[local_k] == current_id && ntime[local_k] <= start_time {
            last_valid = local_k + 1;
            local_k += 1;
        }
        result[i] = last_valid;
    }
    result
}
#[pyfunction]
pub fn tmerge3(id: Vec<i32>, miss: Vec<bool>) -> Vec<usize> {
    let n = id.len();
    let mut result = vec![0; n];
    let mut last_good = 0;
    let mut current_id = -1;
    for (i, (&current, is_missing)) in id.iter().zip(miss).enumerate() {
        if current != current_id {
            current_id = current;
            last_good = 0;
        }
        if is_missing {
            result[i] = last_good;
        } else {
            result[i] = i + 1;
            last_good = i + 1;
        }
    }
    result
}