Modelo Linear Generalizado Misto - GLIMMIX

data tatiana;
input tamanho modelo etapa rep y;
cards;
1 1 1 1 0
.
.
.
.
3 2 4 3 1
;
proc sort data=tatiana;
by tamanho;
proc npar1way wilcoxon dscf;
class modelo;
var y;
by tamanho;
run;

proc glimmix data=tatiana plots=(residualpanel);
class modelo etapa rep tamanho;
model y=modelo etapa /dist=lognormal;
random tamanho rep;
lsmeans modelo etapa /adjust=tukey pdiff=all;
run;