Compare survival distributions
Description
Computes a weighted version of the logrank and stratified logrank tests
for comparing two or more survival distributions. The generalization to
complex samples is based on the characterization of the logrank test as
the score test in a Cox model, Under simple random sampling with
replacement, this function with rho=0
and gamma=0
is almost identical to the robust score test
in the survival package.
Usage
1  svylogrank(formula, design, rho=0,gamma=0,method=c("small","large","score"), ...)

Arguments
formula 
Model formula with a single predictor and optionally a 
design 
A survey design object 
rho,gamma 
Coefficients for the Harrington/Fleming Grhogamma
tests. The default is the logrank test, 
method 

... 
for future expansion. 
Value
A vector containing the zstatistic for comparing each level of the variable to the lowest, the chisquared statistic for the logrank test, and the pvalue.
References
Rader K (2014)
See Also
svykm
, svycoxph
.
Examples
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19  library("survival")
data(nwtco)
## stratified on case status
dcchs<twophase(id=list(~seqno,~seqno), strata=list(NULL,~rel),
subset=~I(in.subcohort  rel), data=nwtco, method="simple")
svylogrank(Surv(edrel,rel)~factor(stage),design=dcchs)
data(pbc, package="survival")
pbc$randomized < with(pbc, !is.na(trt) & trt>0)
biasmodel<glm(randomized~age*edema,data=pbc)
pbc$randprob<fitted(biasmodel)
dpbc<svydesign(id=~1, prob=~randprob, strata=~edema, data=subset(pbc,randomized))
svylogrank(Surv(time,status==2)~trt,design=dpbc)
svylogrank(Surv(time,status==2)~trt,design=dpbc,rho=1)
rpbc<as.svrepdesign(dpbc)
svylogrank(Surv(time,status==2)~trt,design=rpbc)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker. Vote for new features on Trello.