lmm1 <- lmer((total) ~ DISTANCE+POS + (1|NO_UNIT),data=TURN)
TURN$POS<-as.factor(TURN$POS)#Change position and distance to factors
TURN$DISTANCE<-as.factor(TURN$DISTANCE)
TURN$NO_UNIT <- as.factor(TURN$NO_UNIT)

mySumm <- function(.) { s <- sigma(.)
    c(beta =getME(., "beta"), sigma = s, sig01 = unname(s * getME(., "theta"))) }

# run bootstrap analysis for calculation of confidence intervals of parameter estimates
mod_lmm1_boot <- bootMer(lmm1,mySumm, nsim=300)
boot.ci(mod_lmm1_boot,type="perc",conf=.95,index=1)
boot.ci(mod_lmm1_boot,type="perc",conf=.95,index=2)
boot.ci(mod_lmm1_boot,type="perc",conf=.95,index=3)#I want to plot these 3 indexes with the mean!