Pour plus de précisions, voilà ce que je viens d'essayer : [global.model.nmds1 = glmer.nb(nmds1~area+age+nnd+F100+ppt+pptwar+ (1|region/window_ty+(1|patch),interval=log(th)+c(-3,3),verbose=F,data=data3)] Ce qui me donne le message d'erreur suivant : Error in eval(expr, envir, enclos) : negative valu...