Description Usage Arguments Details Value Source References See Also Examples

`fit.nb.gf`

fits data using the pseudo maximum likelihood of a Gamma frailty model

1 2 3 4 5 6 7 8 9 10 11 12 13 |

`dataC` |
a matrix containing count data from the control group, which is to be fitted. Columns correspond to time points, rows to observations. |

`dataE` |
a matrix containing count data from the experiment group, which is to be fitted. Columns correspond to time points, rows to observations. |

`trend` |
the trend which assumed to underlying in the data. |

`lower` |
vector of lower bounds for estimated parameters |

`upper` |
vector of upper bounds for estimated parameters |

`method` |
algorithm used for minimization of the likelihood, see |

`start` |
vector of starting values for estimated parameters |

`approx` |
numer of iterations in numerical calculation of the sandwich estimator, see 'Details'. |

`rho` |
indicates whether or not to calculate the correlation coefficient of Gamma frailties. Must be TRUE or FALSE. |

`H0` |
indicates whether or not to calculate the hessian and outer gradient matrix under the null hypothesis, see 'Details'. |

`h0` |
the value against which is tested under the null |

the function `fit.nb.gf`

fits a Gamma frailty model as found in Fiocco (2009). The fitting function allows for incomplete follow up,
but not for intermittent missingness.

When calculating the expected sandwich estimator required for the sample size, certain terms can not be computed analytically and have
to be approximated numerically. The value `approx`

defines how close the approximation is to the true expected sandwich estimator.
High values of `approx`

provide better approximations but are compuationally more expensive.

If parameter H0 is set to TRUE, the hessian and outer gradient are calculated under the assumption that `lambda[2]`

*≥q* `h0`

if
`trend = "constant"`

or `lambda[3]`

*≥q* `h0`

if `trend = "exponential"`

.

`fit.nb.gf`

returns estimates of the trend parameters `lambda`

, dispersion parameter `size`

,
Hessian matrix `hessian`

, outer gradient product matrix `ogradient`

and, if inquired, correlation coefficient `rho`

.

`fit.nb.gf`

uses code contributed by Thomas Asendorf.

Fiocco M, Putter H, Van Houwelingen JC, (2009), A new serially correlated gamma-frailty process for longitudinal count data *Biostatistics* Vol. 10, No. 2, pp. 245-257.

`rnbinom.gf`

for information on the Gamma frailty model, `n.nb.gf`

for calculating
initial sample size required when performing inference, `bssr.nb.gf`

for blinded
sample size reestimation within a running trial, `optim`

for more information on the used minimization algorithms.

1 2 3 4 |

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.