Learn R Programming

lmeNB (version 1.3)

Compute the Personalized Activity Index Based on a Negative Binomial Model

Description

The functions in this package implement the safety monitoring procedures proposed in the paper titled "Detection of unusual increases in MRI lesion counts in individual multiple sclerosis patients" by Zhao, Y., Li, D.K.B., Petkau, A.J., Riddehough, A., Traboulsee, A., published in Journal of the American Statistical Association in 2013. The procedure first models longitudinally collected count variables with a negative binomial mixed-effect regression model. To account for the correlation among repeated measures from the same patient, the model has subject-specific random intercept, which can be modelled with a gamma or log-normal distributions. One can also choose the semi-parametric option which does not assume any distribution for the random effect. These mixed-effect models could be useful beyond the application of the safety monitoring. The maximum likelihood methods are used to estimate the unknown fixed effect parameters of the model. Based on the fitted model, the personalized activity index is computed for each patient. Lastly, this package is companion to R package lmeNBBayes, which contains the functions to compute the Personalized Activity Index in Bayesian framework.

Copy Link

Version

Install

install.packages('lmeNB')

Monthly Downloads

11

Version

1.3

License

GPL (>= 2)

Maintainer

Yumi Kondo

Last Published

February 2nd, 2015

Functions in lmeNB (1.3)

fitParaAR1

Performs the maximum likelihood estimation for the negative binomial mixed-effect AR(1) model
index.batch

The main function to compute the point estimates and 95% confidence intervals (for a parametric model) of the conditional probabilities $Pr(q(\boldsymbol{Y}_{i,new}) \ge q(\boldsymbol{y}_{i,new})| \boldsymbol{Y}_{i,pre}=\boldsymbol{y}_{i,pre})$ for multiple subjects.
fitSemiAR1

Fit the semi-parametric negative binomial mixed-effect AR(1) model.
CP.se

Compute a conditional probability of observing a set of counts as extreme as the new observations of a subject given the previous observations from the same subject based on the negative binomial mixed effect independent model.
lmeNB-internal

Internal lmeNB functions
fitParaIND

Performs the maximum likelihood estimation for the negative binomial mixed-effect independent model
CP.ar1.se

Compute a conditional probability of observing a set of counts as extreme as the new observations of a subjectvisit given the previous observations of the same subject based on the negative binomial mixed-effect AR(1) model.
rNBME.R

Simulate a dataset from the negative binomial mixed-effect independent/AR(1) model
fitSemiIND

Fit the semi-parametric negative binomial mixed-effect independent model.
lmeNB

RElmeNB

Calculate predicted values of E(Gi|Yi) given the estimates of parameters