gamllik: Log likelihood from the gamma distribution

In RPEXE.RPEXT: Reduced Piecewise Exponential Estimate/Test Software

Description

A function computing the log likelihood from the gamma distribution under an order restriction reduction

Usage

gamllik(structtime,structttot,structdeaths,time_die,ttot,deaths)

Arguments

structtime	change-point times to be used to compute the likelihood value
structttot	total time on test (ttot) between each time point and the previous time point (or 0) corresponding to structtime
structdeaths	number of deaths corresponding to structttot
time_die	all event and censoring times from small to large
ttot	total time on test corresponding to time_die
deaths	the number of deaths corresponding to "ttot"

Value

log of the likelihood

Examples

time_die <- c(0.05,0.08,0.38,0.41,0.64)
ttot <- c(9.2,5.8,52.1,5.8,40.0)
deaths <- c(1,1,1,1,1)
structtime <- c(0.05,0.64)
structttot <- c(9.2, 40.0)
structdeaths = c(1, 5)
gamllik(structtime,structttot,structdeaths,time_die,ttot,deaths)

RPEXE.RPEXT documentation built on July 1, 2020, 6:02 p.m.