DIC e DBS

Author

Isabela Oliveira

DIC- delineamento inteiramente casualizado

Atividade 13-6 (p. 323) - analisar a quantidade (%) de radônio liberado (variável resposta) em chuveiros com diferentes diâmetros no orifício de saida (tratamentos).

A=0,37 (testemunha?) B=0,51 C=0,71 D= 1,02 E= 1,40 F= 1,99

Analise Descritiva

  Repeticoes Media Variancia DesvioPadrao CoefVariacao
A          4 82.75  4.250000     2.061553          2.5
B          4 77.00  5.333333     2.309401          3.0
C          4 75.00  3.333333     1.825742          2.4
D          4 71.75 10.916667     3.304038          4.6
E          4 65.00 12.666667     3.559026          5.5
F          4 62.75  7.583333     2.753785          4.4

Teste de Normalidade da variável resposta

[1] 60
[1] 85
[1] 20.0 27.5 35.0 42.5 50.0 57.5

QQplot (variável resposta)

Teste de Shapiro para normalidade (resposta?)

         Teste Estatística_W Valor_p
W Shapiro-Wilk        0.9544  0.3369

Teste da homogeneidade de variâncias

        Teste Estatistica Valor_p
       Levene      0.8605  0.5261
F Max Hartley      3.8000  0.8953

Teste de hipótese

                  Fonte GL  Soma SQ Quadrado Médio       F Valor-p
tratamento  tratamento   5 1133.375       226.6750 30.8518 < 0.001
Residuals   Residuals   18  132.250         7.3472      NA    <NA>

Resíduos ordinários

Residuos x valores ajustados

Residuos padronizados x Valores ajustados

Teste de Homogeneidade de variâncias para os residuos

        Teste Estatistica Valor_p
       Levene      0.8605  0.5261
F Max Hartley      3.8000  0.8953

Gráfico para Normalidade dos residuos

QQplot (residuos)

Teste de hipótese (normalidade de resíduos)

         Teste Estatística_W Valor_p
W Shapiro-Wilk        0.9348  0.1247

Comparações múltiplas

Tukey

  Tukey multiple comparisons of means
    95% family-wise confidence level
    factor levels have been ordered

Fit: aov(formula = radônio_liberado ~ tratamento, data = dados)

$tratamento
     diff        lwr       upr     p adj
E-F  2.25 -3.8412336  8.341234 0.8432736
D-F  9.00  2.9087664 15.091234 0.0021152
C-F 12.25  6.1587664 18.341234 0.0000650
B-F 14.25  8.1587664 20.341234 0.0000089
A-F 20.00 13.9087664 26.091234 0.0000001
D-E  6.75  0.6587664 12.841234 0.0249971
C-E 10.00  3.9087664 16.091234 0.0007059
B-E 12.00  5.9087664 18.091234 0.0000841
A-E 17.75 11.6587664 23.841234 0.0000004
C-D  3.25 -2.8412336  9.341234 0.5513482
B-D  5.25 -0.8412336 11.341234 0.1153360
A-D 11.00  4.9087664 17.091234 0.0002404
B-C  2.00 -4.0912336  8.091234 0.8968057
A-C  7.75  1.6587664 13.841234 0.0084181
A-B  5.75 -0.3412336 11.841234 0.0707511

Ducan


Study: mod_aov ~ "tratamento"

Duncan's new multiple range test
for radônio_liberado 

Mean Square Error:  7.347222 

tratamento,  means

  radônio_liberado      std r       se Min Max   Q25  Q50   Q75
A            82.75 2.061553 4 1.355288  80  85 82.25 83.0 83.50
B            77.00 2.309401 4 1.355288  75  79 75.00 77.0 79.00
C            75.00 1.825742 4 1.355288  73  77 73.75 75.0 76.25
D            71.75 3.304038 4 1.355288  67  74 70.75 73.0 74.00
E            65.00 3.559026 4 1.355288  62  69 62.00 64.5 67.50
F            62.75 2.753785 4 1.355288  60  66 60.75 62.5 64.50

Alpha: 0.05 ; DF Error: 18 

Critical Range
       2        3        4        5        6 
4.026767 4.224954 4.350007 4.436661 4.499999 

Means with the same letter are not significantly different.

  radônio_liberado groups
A            82.75      a
B            77.00      b
C            75.00     bc
D            71.75      c
E            65.00      d
F            62.75      d

Scheffé


Study: mod_aov ~ "tratamento"

Scheffe Test for radônio_liberado 

Mean Square Error  : 7.347222 

tratamento,  means

  radônio_liberado      std r       se Min Max   Q25  Q50   Q75
A            82.75 2.061553 4 1.355288  80  85 82.25 83.0 83.50
B            77.00 2.309401 4 1.355288  75  79 75.00 77.0 79.00
C            75.00 1.825742 4 1.355288  73  77 73.75 75.0 76.25
D            71.75 3.304038 4 1.355288  67  74 70.75 73.0 74.00
E            65.00 3.559026 4 1.355288  62  69 62.00 64.5 67.50
F            62.75 2.753785 4 1.355288  60  66 60.75 62.5 64.50

Alpha: 0.05 ; DF Error: 18 
Critical Value of F: 2.772853 

Minimum Significant Difference: 7.13666 

Means with the same letter are not significantly different.

  radônio_liberado groups
A            82.75      a
B            77.00     ab
C            75.00      b
D            71.75     bc
E            65.00     cd
F            62.75      d

Scott-Knott

List of 2
 $ new.fact.lev:List of 1
  ..$ tratamento: chr [1:6] "A" "B" "C" "D" ...
 $ grid.data   :'data.frame':   6 obs. of  1 variable:
  ..$ tratamento: chr [1:6] "A" "B" "C" "D" ...
Goups of means at sig.level = 0.05 
  Means G1 G2 G3 G4
A 82.75  a         
B 77.00     b      
C 75.00     b      
D 71.75        c   
E 65.00           d
F 62.75           d

Apresentação final dos resultados

DBC-delineamento em blocos casualizados

Experimento para minimizar variações nos níveis da glicose no sangue. Tratamento=10 ou 20 min de exercício; Variável bloco=tarde,manhã. Duas repetições por tratamento

T1-10 min/manhã

T2-10 min/Tarde

T3-20 min/manhã

T4-20 min/Tarde

  bloco tratamento glicose_média
1 tarde         10          71.5
2 manhã         10         103.0
3 tarde         10          83.5
4 manhã         10         126.0
5 tarde         20         125.5
6 manhã         20         129.5
7 tarde         20          95.0
8 manhã         20          93.0

Análise descritiva da variável resposta (glicose_média)

Analise do glicose por tratamento

$`10`
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
  71.50   80.50   93.25   96.00  108.75  126.00 

$`20`
   Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
   93.0    94.5   110.2   110.8   126.5   129.5 
# A tibble: 2 × 7
  tratamento `N de repeticoes` Media Mediana `Desvio Padrao` Variancia
  <fct>                  <int> <dbl>   <dbl>           <dbl>     <dbl>
1 10                         4   96     93.2            23.8      568.
2 20                         4  111.   110.             19.4      377.
# ℹ 1 more variable: `Coef. de variacao` <dbl>
# A tibble: 2 × 14
  tratamento variable         n   min   max median    q1    q3   iqr   mad  mean
  <fct>      <fct>        <dbl> <dbl> <dbl>  <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
1 10         glicose_méd…     4  71.5  126    93.2  80.5  109.  28.2  23.4   96 
2 20         glicose_méd…     4  93    130.  110.   94.5  126.  32    24.1  111.
# ℹ 3 more variables: sd <dbl>, se <dbl>, ci <dbl>

Grafico em pontos com o traco da mediana

Grafico em pontos com o traco da media

Analise no tratamento e no bloco

Ajuste do modelo

                  Fonte GL  Soma SQ Quadrado Médio      F Valor-p
bloco       bloco        1  722.000        722.000 1.7063  0.2483
tratamento  tratamento   1  435.125        435.125 1.0283  0.3571
Residuals   Residuals    5 2115.750        423.150     NA    <NA>

Oq pode ser extraído

                     Componente
1         Nomes dos Componentes
2 Graus de Liberdade (Resíduos)
3                      Resíduos
                                                                                                                   Valor
1 coefficients, residuals, effects, rank, fitted.values, assign, qr, df.residual, contrasts, xlevels, call, terms, model
2                                                                                                                      5
3                                                                        -15, -2.5, -3, 20.5, 24.25, 9.25, -6.25, -27.25

Resíduos padronizados

         1          2          3          4          5          6          7 
-0.9223677 -0.1537279 -0.1844735  1.2605692  1.4911611  0.5687934 -0.3843199 
         8 
-1.6756346 

Valores ajustados pelo modelo

Grafico Residuos Padronizados x Valores ajustados

Gráfico Q-Q plot: para verificação da suposição de Normalidade

[1] 8 5

Variância dos residuos em cada tratamento

      10       20 
220.1667 485.0833 

teste de Homogeneidade de variancias dos residuos

        Teste Estatistica Valor_p
       Levene      1.0169  0.3522
F Max Hartley      2.2033  0.5332

Teste para normalidade dos residuos

         Teste Estatística_W Valor_p
W Shapiro-Wilk        0.9627  0.8356

Teste de tukey

  Tukey multiple comparisons of means
    95% family-wise confidence level
    factor levels have been ordered

Fit: aov(formula = glicose_média ~ bloco + tratamento)

$tratamento
       diff      lwr     upr     p adj
20-10 14.75 -22.6407 52.1407 0.3570873