Nothing
This R package contains multivariate two-sample survival permutation tests, based on the logrank and Gehan statistics. The tests are described in Persson et al. (2019).
To install the development version from GitHub:
library(devtools)
install_github("lukketotte/MultSurvTests")
Example usage, comparing the bivariate survival times of the two treatment groups in the diabetes
data (included in the package):
library(MultSurvTests)
# Diabetes data:
?diabetes
# Survival times for the two groups:
x <- as.matrix(subset(diabetes, LASER==1)[c(6,8)])
y <- as.matrix(subset(diabetes, LASER==2)[c(6,8)])
# Censoring status for the two groups:
delta.x <- as.matrix(subset(diabetes, LASER==1)[c(7,9)])
delta.y <- as.matrix(subset(diabetes, LASER==2)[c(7,9)])
# Create the input for the test:
z <- rbind(x, y)
delta.z <- rbind(delta.x, delta.y)
# Run the tests with 99 permutations:
perm_gehan(B = 99, z, delta.z, n1 = nrow(x))
perm_mvlogrank(B = 99, z, delta.z, n1 = nrow(x))
# In most cases, it is preferable to use more than 99
# permutations for computing p-values. choose_B() can
# be used to determine how many permutations are needed.
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.