dati0<- dati[dati$varB==”0”,]
dati1<- dati[dati$varB==”1”,]
p<- diffPerc<- vector()
for (i in 1:1000){
  datiboot<- rbind(dati0[sample(nrow(dati0),nrow(dati1)),],dati1)
  glmm1boot<- lmer(varA~ varB+varC+varD+(1|Year),data=dati,REML = F)
  glmm1bootVarB<- lmer(varA~ varC+varD+(1|Year),data=dati,REML = F)
  p[i]<- (anova(glmm1boot,glmm1boot))[2,8]
  diffPerc[i]<-((effect("varB",glmm1boot)$fit[1] effect("varB",glmm1boot)$fit[2])/effect("varB",glmm1boot)$fit[1])*100
}