Fits a multiplicative relative risk model to interval censored data.

Description

Utility function.

The model fitted assumes a piecewise constant baseline rate in intervals specified by the argument breaks, and a multiplicative relative risk function.

Usage

1
2
  fit.mult( y, rates.frame, cov.frame, start )
  

Arguments

y

Binary vector of outcomes

rates.frame

Dataframe expanded from the original data by expand.data, cooresponding to covariates for the rate parameters.

cov.frame

do., but covariates corresponding to the formula argument of Icens

start

Starting values for the rate parameters. If not supplied, then starting values are generated.

Details

The model is fitted by alternating between two generalized linear models where one estimates the underlying rates in the intervals, and the other estimates the log-relative risks.

Value

A list with three components:

rates

A glm object from a binomial model with log-link, estimating the baseline rates.

cov

A glm object from a binomial model with complementary log-log link, estimating the log-rate-ratios

niter

Nuber of iterations, a scalar

Author(s)

Martyn Plummer, plummer@iarc.fr, Bendix Carstensen, bxc@steno.dk

References

B Carstensen: Regression models for interval censored survival data: application to HIV infection in Danish homosexual men. Statistics in Medicine, 15(20):2177-2189, 1996.

CP Farrington: Interval censored survival data: a generalized linear modelling approach. Statistics in Medicine, 15(3):283-292, 1996.

See Also

Icens fit.add

Examples

1
2
  data( HIV.dk ) 
  

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