# aalenfrailty: Aalen frailty model In mets: Analysis of Multivariate Event Times

 aalenfrailty R Documentation

## Aalen frailty model

### Description

Additive hazards model with (gamma) frailty

### Usage

```aalenfrailty(time, status, X, id, theta, B = NULL, ...)
```

### Arguments

 `time` Time variable `status` Status variable (0,1) `X` Covariate design matrix `id` cluster variable `theta` list of thetas (returns score evaluated here), or starting point for optimization (defaults to magic number 0.1) `B` (optional) Cumulative coefficients (update theta by fixing B) `...` Additional arguments to lower level functions

### Details

Aalen frailty model

### Value

Parameter estimates

Klaus K. Holst

### Examples

```library("timereg")
dd <- simAalenFrailty(5000)
f <- ~1##+x
X <- model.matrix(f,dd) ## design matrix for non-parametric terms
system.time(out<-timereg::aalen(update(f,Surv(time,status)~.),dd,n.sim=0,robust=0))
dix <- which(dd\$status==1)
t1 <- system.time(bb <- .Call("Bhat",as.integer(dd\$status),
X,0.2,as.integer(dd\$id),NULL,NULL,
PACKAGE="mets"))
spec <- 1
##plot(out,spec=spec)
## plot(dd\$time[dix],bb\$B2[,spec],col="red",type="s",
##      ylim=c(0,max(dd\$time)*c(beta0,beta)[spec]))
## abline(a=0,b=c(beta0,beta)[spec])
##'

## Not run:
thetas <- seq(0.1,2,length.out=10)
Us <- unlist(aalenfrailty(dd\$time,dd\$status,X,dd\$id,as.list(thetas)))
##plot(thetas,Us,type="l",ylim=c(-.5,1)); abline(h=0,lty=2); abline(v=theta,lty=2)
op <- aalenfrailty(dd\$time,dd\$status,X,dd\$id)
op

## End(Not run)
```

mets documentation built on Jan. 17, 2023, 5:12 p.m.