lmec-package: Linear Mixed-Effects Models with Censored Responses

In lmec: Linear Mixed-Effects Models with Censored Responses

Description

This package includes a function to fit a linear mixed-effects model in the formulation described in Laird and Ware (1982) but allowing for censored normal responses. In this version, the with-in group errors are assumed independent and identically distributed.

Details

Package:	lmec
Type:	Package
Version:	1.0
Date:	2009-01-28

Author(s)

Florin Vaida and Lin Liu

Maintainer: Florin Vaida <[e-mail]) and Lin Liu <[e-mail]>

References

Vaida, Florin and Liu, Lin, Fast Implementation For Normal Mixed Effects Models with Censored Response (submitted).

Vaida, Florin and Fitzgerald, Anthony and DeGruttola, Victor (2007), Efficient Hybrid EM for nonlinear mixed effects models with censored response, Computational Statistics and Data Analysis, 51, 5718-5730.

Examples

data(UTIdata)
UTIdata <- subset(UTIdata, !is.na(RNA))
o <- order(UTIdata$Patid, UTIdata$Fup)
UTIdata <- UTIdata[o,]
cens <- (UTIdata$RNAcens==1)+0
yL<-  log10(UTIdata$RNA)
X<- cbind((UTIdata$Fup==0)+0, (UTIdata$Fup==1)+0, (UTIdata$Fup==3)+0, (UTIdata$Fup==6)+0, (UTIdata$Fup==9)+0, (UTIdata$Fup==12)+0, (UTIdata$Fup==18)+0, (UTIdata$Fup==24)+0)
Z <- matrix(rep(1, length(yL)), ncol=1)
cluster<-  as.numeric(UTIdata$Patid)
fit <- lmec(yL,cens, X, Z, cluster, method='ML', maxstep=40)

Example output

Loading required package: mvtnorm
[1] "1"                 "0"                 "-466.568101849883"
[1] "2"                 "0"                 "48.7180426043605" 
[4] "-417.850059245523"
[1] "3"                 "0"                 "4.80013344709465" 
[4] "-413.049925798428"
[1] "4"                 "0"                 "0.83934928265063" 
[4] "-412.210576515777"
[1] "5"                 "0"                 "0.139392293211529"
[4] "-412.071184222566"
[1] "6"                  "0"                  "0.0290105840074375"
[4] "-412.042173638558" 
[1] "7"                   "0"                   "-0.0028765939925961"
[4] "-412.045050232551"  
[1] "8"                   "0"                   "0.00196282729200448"
[4] "-412.043087405259"  
[1] "9"                   "0"                   "0.00164700355088598"
[4] "-412.041440401708"  
[1] "10"                   "0"                    "-0.00142453318204616"
[4] "-412.04286493489"    
[1] "11"                  "0"                   "0.00373711926431497"
[4] "-412.039127815626"  
[1] "12"                  "0"                   "0.00316452889950369"
[4] "-412.035963286726"  
[1] "13"                   "1"                    "-0.00165990116903458"
[4] "-412.037623187895"   
[1] "14"                  "2"                   "-0.0051415821660612"
[4] "-412.042764770062"  

lmec documentation built on May 1, 2019, 9:15 p.m.