search for: workerc

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...