Title: | Parameter Estimation in Conditional GEE for Recurrent Event Gap Times |
---|---|
Description: | Solves for the mean parameters, the variance parameter, and their asymptotic variance in a conditional GEE for recurrent event gap times, as described by Clement and Strawderman (2009) in the journal Biostatistics. Makes a parametric assumption for the length of the censored gap time. |
Authors: | David Clement |
Maintainer: | David Clement <[email protected]> |
License: | GPL (>= 2) |
Version: | 0.2.0 |
Built: | 2024-12-14 06:32:00 UTC |
Source: | CRAN |
This data set gives the start and stop times of recurrent asthma events in children. It also provides a subject ID, treatment indicator, censoring indicator, number of events per subject and a first event indicator.
A data frame with 1037 rows and 7 columns. See asthma.txt header for details.
http://www.blackwellpublishing.com/rss/
Duchateau et al. JRSSC 2003. Volume 52, 355–363.
Solves for the mean parameters (), the
variance parameter (
), and their asymptotic variance
in a conditional GEE for recurrent event gap times, as described by
Clement, D. Y. and Strawderman, R. L. (2009)
condGEE( data, start, mu.fn = MU, mu.d = MU.d, var.fn = V, k1 = K1.norm, k2 = K2.norm, robust = TRUE, asymp.var = TRUE, maxiter = 100, rtol = 1e-06, atol = 1e-08, ctol = 1e-08, useFortran = TRUE )
condGEE( data, start, mu.fn = MU, mu.d = MU.d, var.fn = V, k1 = K1.norm, k2 = K2.norm, robust = TRUE, asymp.var = TRUE, maxiter = 100, rtol = 1e-06, atol = 1e-08, ctol = 1e-08, useFortran = TRUE )
data |
matrix of data with one row for each gap time; the first column should be a subject ID, the second column the gap time, the third column a completeness indicator equal to 1 if the gap time is complete and 0 if the gap time is censored, and the remaining columns the covariates for use in the mean and variance functions |
start |
vector containing initial guesses for the unknown parameter vector |
mu.fn |
the specification for the mean of the gap time; the default is
a linear combination of the covariates; the function should take two arguments
( |
mu.d |
the derivative of |
var.fn |
the specification for |
k1 |
the function to solve for the conditional mean length of the censored
gap times; its sole argument should be the vector of standardized (i.e.\
|
k2 |
the function to solve for the conditional mean length of the square
of the censored gap times; its sole argument should be the vector of
standardized (i.e.\ |
robust |
logical, if |
asymp.var |
logical, if |
maxiter |
see |
rtol |
see |
atol |
see |
ctol |
see |
useFortran |
see |
conditional expectation
David Clement
E(Y|Y>w)
where Y is exponential dist with mean 0
and variance 1
K1.exp(w)
K1.exp(w)
w |
real value |
conditional expectation
David Clement
E(Y|Y>w)
where Y is normal
K1.norm(w)
K1.norm(w)
w |
real value |
conditional expectation
David Clement
E(Y|Y>w)
where Y is t dist with 3 df
K1.t3(w)
K1.t3(w)
w |
real value |
conditional expectation
David Clement
E(Y^2|Y>w)
where Y is exponential dist with mean 0
and variance 1
K2.exp(w)
K2.exp(w)
w |
real value |
conditional expectation
David Clement
E(Y^2|Y>w)
where Y is normal
K2.norm(w)
K2.norm(w)
w |
real value |
conditional expectation
David Clement
E(Y^2|Y>w)
where Y is t dist with 3 df
K2.t3(w)
K2.t3(w)
w |
real value |
conditional expectation
David Clement