* Delineamento e Análise de Experimentos 1
* Prof. George von Borries
* Exemplos de Modelos Mistos;

options nocenter nodate;

* Exemplo 1: Turbinas - Ott e Longnecker, pg. 996;

Title1 'Exemplo das Turbinas';

data turbinas;
    do temperaturas = 15 to 35 by 10;
        do medidores = 1 to 4;
            do rep = 1 to 6;
                input pressao @@;
                output;
            end;
        end;
    end;
    datalines;
    40 40 37 47 42 41
    43 34 38 42 39 35
    42 35 35 41 43 36
    47 47 40 36 41 47
    57 57 65 67 63 59
    49 43 51 49 45 43
    44 45 49 45 46 43
    36 49 38 45 38 42
    35 35 35 46 41 42
    41 43 44 36 42 41
    42 41 34 35 39 36
    41 44 35 46 44 46
    ;

ods graphics on;

proc glm data=turbinas;
  class temperaturas medidores;
  model pressao = medidores|temperaturas;
  random medidores temperaturas*medidores / test;
run;

proc mixed data=turbinas covtest;
  class temperaturas medidores;
  model pressao = temperaturas;
  random medidores temperaturas*medidores;
* lsmeans temperaturas;
run;


* Exemplo 2: Kuehl, Exemplo 7.2, pág. 237;

title 'Métodos Químicos';

data eval;
    input method day serum @@;
    datalines;
1 1 142.3 1 1 144.0
1 2 134.9 1 2 146.3
1 3 148.6 1 3 156.5
1 4 152.0 1 4 151.4
2 1 142.9 2 1 147.4
2 2 125.9 2 2 127.6
2 3 135.5 2 3 138.9
2 4 142.9 2 4 142.3
;
run;

proc glm data=eval;
    class method day;
    model serum=method day method*day;
    test h=method e=method*day;
    test h=day e=method*day;
    random day method*day / test;
run;

proc mixed cl covtest;
    class method day;
    model serum = method;
    random day method*day;
run;

quit;