# rmeanDiff: rmeanDiff In mbrueckner/AHR: Estimation and Testing of Average Hazard Ratios

## Description

Estimate difference of restricted mean survival based on (weighted) Kaplan-Meier estimates of the survival functions in each group.

## Usage

 `1` ```rmeanDiff(L, formula, data, rr.subset = rep(TRUE, nrow(data))) ```

## Arguments

 `L` time-limit specifying up to which time restricted mean will be calculated `formula` an object of class '"formula"' specifying the conditional survival model `data` data frame containing the variables in formula `rr.subset` logical vector defining subset of observations to use for response rate estimation (default: use all observations)

## Value

An object of class '"rmd"', i.e. a list containing:

 `L` time limit, i.e. restricted mean up to time L is calculated `rmean1` restricted mean in group 1 `rmean2` restricted mean in group 2 `rmean.diff` estimated restricted mean difference `var.rmean1` an estimate of the asymptotic variance of the restricted mean in group 1 `var.rmean2` an estimate of the asymptotic variance of the restricted mean in group 2 `var.rmean.diff` an estimate of the asymptotic variance of the restricted mean difference `Z.rmean` the standardized test statistic for testing rmean.diff=0 `p.value` p-value corresponding to Z.rmean

## Examples

 ```1 2 3 4 5 6``` ```T <- c(rexp(100, 1), rexp(100, 2)) C <- c(rexp(100, 1), rexp(100, 2)) Y <- pmin(T, C) D <- T <= C Z <- rep(c(0,1), c(100, 100)) fit <- rmeanDiff(2, formula=Surv(Y, D) ~ Z, data.frame(Y=Y, D=D, Z=Z)) ```

mbrueckner/AHR documentation built on May 22, 2017, 12:16 a.m.