binregG: G-estimator for binomial regression model (Standardized...

In kkholst/mets: Analysis of Multivariate Event Times

G-estimator for binomial regression model (Standardized estimates)

Description

Computes G-estimator

\hat F(t,A=a) = n^{-1} \sum_i \hat F(t,A=a,Z_i)

. Assumes that the first covariate is $A$. Gives influence functions of these risk estimates and SE's are based on these. If first covariate is a factor then all contrast are computed, and if continuous then considered covariate values are given by Avalues.

Usage

binregG(x, data, Avalues = c(0, 1), varname = NULL)

Arguments

x	phreg or cifreg object
data	data frame for risk averaging
Avalues	values to compare for first covariate A
varname	if given then averages for this variable, default is first variable

Author(s)

Thomas Scheike

Examples

data(bmt); bmt$time <- bmt$time+runif(408)*0.001
bmt$event <- (bmt$cause!=0)*1

b1 <- binreg(Event(time,cause)~age+tcell+platelet,bmt,cause=1,time=50)
sb1 <- binregG(b1,bmt,Avalues=c(0,1,2))
summary(sb1)

kkholst/mets documentation built on April 24, 2024, 11:33 a.m.