Displaying 1 result from an estimated 1 matches for "workerc".
Did you mean:
worker
2011 May 16
4
Problem on glmer
...#############################################################
## Code to fit the model
#############################################################
library(MASS)
library(lme4)
address <- "http://instruction.bus.wisc.edu/jfrees/jfreesbooks/Longitudinal%20and%20Panel%20Data/Book/Data/TXTData/WorkerC.txt"
wc <- read.table(address,header=T)
# clean data
cl <- unique(wc$CL[wc$LOSS==0])
wc2 <- subset(wc, !(CL %in% cl) )
# fit model
fit1 <- glmmPQL(LOSS~ YR+offset(log(PR)), random=~1|CL,family=Gamma(link="log"),
data=wc2)
fit2 <- glmer(LOSS~ (1|CL) + YR+offse...