R Graphical Manual

Browse All

Last data update: 2014.03.03

R: Fits a multiplicative relative risk model to interval...

fit.mult

R Documentation

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

  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.

Examples

  data( HIV.dk )

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(Epi)

Attaching package: 'Epi'

The following object is masked from 'package:base':

    merge.data.frame

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Epi/fit.mult.Rd_%03d_medium.png", width=480, height=480)
> ### Name: fit.mult
> ### Title: Fits a multiplicative relative risk model to interval censored
> ###   data.
> ### Aliases: fit.mult
> ### Keywords: models regression survival
> 
> ### ** Examples
> 
>   data( HIV.dk ) 
Warning message:
In data(HIV.dk) : data set 'HIV.dk' not found
>   
> 
> 
> 
> 
> dev.off()
null device 
          1 
>