Experimentos em Parcelas Subdivididas©cnicas experiment… · Vantagens e Desvantagens Modelo...
Transcript of Experimentos em Parcelas Subdivididas©cnicas experiment… · Vantagens e Desvantagens Modelo...
Parcelas SubdivididasDesdobramentos
Experimentos em Parcelas Subdivididas
Lucas Santana da Cunhahttp://www.uel.br/pessoal/lscunha
08 de novembro de 2018Londrina
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Introducao
Tal como no caso de fatorial, o termo parcelas subdivididasnao se refere a um tipo de delineamento e sim ao esquema doexperimento, ou seja, a maneira pela qual os tratamentos saoorganizados.
Nos experimentos em parcelas subdivididas, em geral, estuda-sesimultaneamente dois tipos de fatores os quais sao geralmentedenominados de fatores primarios e fatores secundarios.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Na instalacao os nıveis do fator primario (A) sao distribuıdos asparcelas segundo um tipo de delineamento experimental: DIC,DBC, DQL.
Posteriormente os nıveis do fator secundario (B) sao distribuıdosao acaso as subparcerlas de cada parcela, quando possıvel.
Tal disposicao permite obter uma estimativa geral de maiorprecisao para os efeitos dos nıveis do segundo fator.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Nos experimentos em parcelas subdivididas tem-se dois resıduosdistintos: um correspondente as parcelas e outro as subparcelasdentro das parcelas.
Em casos mais complexos, as subparcelas podem, tambem, serrepartidas em subsubparcelas. Tem-se, neste caso, tres resıduosdistintos:
Resıduo (a), referente as parcelas;
Resıduo (b), a subparcelas e
Resıduo (c), correspondendo as subsubparcelas.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Eem parcela subdividida com a nıveis primarios, b nıveis se-cundarios e r repeticoes, temos a seguinte decomposicao dosgraus de liberdade:
Tabela 1: Parcela subdividida no delineamento inteiramentecasualizado
CV GL
Tratamento A a− 1
Resıduo(a) a(r − 1)
Parcelas ar − 1
Tratamento B b − 1
A× B (a− 1)(b − 1)
Resıduo(b) a(b − 1)(r − 1)
Total abr − 1
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Tabela 2: Parcela subdividida no delineamento em blocoscasualizados.
CV GL
Blocos r − 1
Tratamento A a− 1
Resıduo(a) (a− 1)(r − 1)
Parcelas ar − 1
Tratamento B b − 1
A× B (a− 1)(b − 1)
Resıduo(b) a(r − 1)(b − 1)
Total abr − 1
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Tabela 3: Parcela subdividida no delineamento em quadrado latino.
CV GL
Linhas a− 1
Colunas a− 1
Tratamento A a− 1
Resıduo(a) (a− 1)(a− 2)
Parcelas a2 − 1
Tratamento B b − 1
A× B (a− 1)(b − 1)
Resıduo(b) a(a− 1)(b − 1)
Total a2b − 1
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Vantagens
1 Em comparacao com experimentos fatoriais, experimentos emparcelas subdivididas sao mais faceis de instalar;
2 Quando os tratamentos associados aos nıveis de um dos fatoresexigem maior quantidade de material na unidade experimentaldo que os tratamentos do outro fator.
3 O esquema pode ser utilizado quando um fator adicional e in-corporado num experimento, para ampliar seu objetivo.
4 Atraves da previa informacao, sabe-se que maiores diferencaspodem ser esperadas entre os nıveis de um certo fator do queentre os nıveis do outro fator.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Desvantagens
1 Do ponto de vista estatıstico, os fatoriais sao, em geral, maiseficientes que os em parcelas subdivididas;
2 Enquanto nos fatoriais temos um so resıduo para todos os F ecomparacoes de medias, no ”split-plot” ha dois resıduos, umpara comparacoes de parcelas e outro para subparcelas;
3 Para parcela, o numero de GL geralmente e pequeno, levandoa pouca sensibilidade na analise;
4 Sempre que possıvel, e preferıvel utilizar experimentos fatoriaisem lugar dos experimentos em parcelas subdivididas.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Modelo estatıstico
O modelo linear para o experimento em parcelas subdivididas no delineamentoem blocos ao acaso e dado por:
yijk = µ+ τi + γk + eik + βj + (τβ)ij + εijk ,
i = 1, 2, . . . , aj = 1, 2, . . . , bk = 1, 2, . . . , r
(1)
em que:
yijk e o valor observado no i-esimo tratamento, k-esimo bloco e j-esimasubparcela;
µ e uma constante;
τi e o efeito do i-esimo fator A;
γk e o efeito do k-esimo bloco;
eik e o resıduo (a) da parcela;
βj e o efeito do j-esimo fator B;
(τβ)ij e a interacao entre o i-esimo fator A e o j-esimo fator B;
εijk e o resıduo (b) da subparcela;
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
No experimento em parcelas subdivididas, em geral, deseja-se testarprimeiramente a significancia da interacao entre os fatores. No casode dois fatores, tem-se:
H0 : (τβ)ij = 0 para todo i , j
H1 : Pelo menos um (τβ)ij 6= 0
Caso a interacao nao seja significativa, testa-se os efeitos principais:
H0 : τ1 = τ2 = . . . τa = 0
H1 : Pelo menos um τi 6= 0
H0 : β1 = β2 = . . . βb = 0
H1 : Pelo menos um βj 6= 0
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Analise de Variancia
Tabela 4: Quadro da Analise de Variancia em um delineamento emblocos.
C.V. G.L. S.Q. Q.M. Fcalc
Blocos r − 1 SQBlocosSQBlocos
r−1QMBlocosQMRes(a)
A a − 1 SQASQAa−1
QMAQMRes(a)
Resıduo(a) (a − 1)(r − 1) SQRes(a)
SQRes(a)(a−1)(r−1)
(Parcelas) (ar − 1) (SQParcelas )
B b − 1 SQBSQBb−1
QMBQMRes(b)
A× B (a − 1)(b − 1) SQAxBSQAxB
(a−1)(b−1)
QMAxBQMRes(b)
Resıduo(b) a(r − 1)(b − 1) SQRes(b)
SQRes(b)a(r−1)(b−1)
Total abr − 1 SQTotal
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Em que as somas de quadrados sao dadas por:
SQTotal =a∑
i=1
b∑j=1
r∑k=1
y 2ijk − C C =
(∑ai=1
∑bj=1
∑rk=1 yijk
)2
abr
SQA =1
r × b
a∑i=1
T 2Ai
− C
SQBlocos =1
a× b
r∑k=1
T 2Blocok − C
SQParcelas =1
b
a∑i=1
r∑k=1
T 2Parcela − C
SQRes(a) = SQParcelas − SQA − SQBlocos
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
SQB =1
a× r
b∑j=1
T 2B − C
SQA,B =1
r
a∑i=1
b∑j=1
T 2Ai ,Bj− C
SQAxB = SQA,B − SQA − SQB
SQRes(b) = SQTotal − SQParcelas − SQA − SQAxB
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Exemplo 1
Suponha o caso de um experimento com tres racoes (A,B, e C ), em seis
blocos casualizados, cada parcela constituıda por dois animais. Em uma
determinada fase do ensaio, os bovinos, dentro de cada parcela, passaram
a receber, por sorteio, um dos tipos de suplementos minerais (M ou P).
Os ganhos de pesos individuais, ao final do experimento, sao apresentados
na tabela abaixo.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Tabela 5: Ganhos de pesos, em quilos, ao final do experimento.Tipos de Racao
Blocos A B C TotaisM P M P M P
I 107 89 116 101 90 96 599II 117 101 136 110 112 89 665III 122 98 130 104 99 92 645IV 111 101 122 91 105 78 608V 90 95 117 100 110 90 602VI 116 90 114 94 114 93 621
Totais 663 574 735 600 630 538 3.740
A um nıvel de significancia de 5%, faca a analise de variancia e con-siderando um experimento em parcela subdividida no delineamentoem blocos ao acaso, em que o tipo de suplemento mineral esta nasubparcela.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Tabela 6: Tabela auxiliar para calculo das somas de quadrados dasparcelas.
Blocos (2)Tipos de Racao
TotaisA B C
I 196 217 186 599
II 218 246 201 665
III 220 234 191 645
IV 212 213 183 608
V 185 217 200 602
VI 206 208 207 621
Totais 1.237 (12) 1.335 (12) 1.168 (12) 3.740
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
C =(107 + 117 + · · · + 90 + 93)2
3 × 2 × 6= 388.544,4
SQTotal = (1072 + 1172 + · · · + 902 + 932) − 388.544, 4 = 6.061,556
SQRac =1
2 × 6× (1.2372 + 1.3352 + 1.1682) − 388.544, 4 = 1.173,722
SQBlocos =1
2 × 3× (5992 + · · · + 6212) − 388.544, 4 = 582,2222
SQParcelas =1
2× (1962 + · · · + 2002 + 2072) − 388.544, 4 = 2.377,556
SQRes(a) = SQParcelas − SQTrat − SQBlocos
= 2.377, 556 − 1.173, 722 − 582, 2222 = 621,6111
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Tabela 7: Tabela auxiliar para calculo das somas de quadrados dasSubparcelas.
Suplementos (6)Tipos de Racao
TotaisA B C
M 663 735 630 2.028
P 574 600 538 1.712
Totais 1.237 (12) 1.335 (12) 1.168 (12) 3.740
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
SQSup =1
3 × 6× (2.0282 + 1.7122) − 388.544, 4 = 2.773,778
SQRac,Sup =1
3 × 2(6632 + 5742 + · · · + 6302 + 5382) − 388.544, 4
= 4.057.889
SQInter = 4.057, 889 − 1.173, 722 − 2.773, 778
= 110,3889
SQRes(b) = SQTotal − SQParcelas − SQSup − SQInter
= 6.061, 556 − 2.377, 556 − 2.773, 778 − 110, 3889
= 799,8333
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Tabela 8: Quadro da analise de variancia do experimento em parcelassubdivididas no delineamento em blocos ao acaso.
Causa da Variacao S.Q. g.l. Q.M. Fcalc Pr(> F )
Blocos 582, 22 5 116, 44
Racao 1.173, 72 2 586, 86 9, 441 0, 004976∗∗
Resıduo(a) 621, 61 10 62, 16
(Parcelas) 2.377, 556 17
Suplementos 2.773, 78 1 2.773, 78 52, 0192 3, 011 × 10(−6)∗∗∗
Racao × Suplementos 110, 39 2 55, 19 1, 0351 0, 3792
Resıduo(b) 799, 83 15 53, 32
Total 6.061, 556 35
Os efeitos das Racoes e dos Blocos sao testados usando o Resıduo(a).
Os efeitos dos Suplementos e da Interacao sao testados usando o Resıduo(b).
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Efeito da Interacao
Verifica-se da Tabela 8 que a interacao entre os tipos de Racaoe Suplementos nao foi significativa, havendo efeito dos fatoresprincipais: Racao e Suplemento.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Efeito Racao
No caso de racao, verifica-se que o efeito e significativo, eassim, pelo teste de tukey, temos:
∆ = q(3,10,5%)
√QMRes(a)
br
= 3, 88×√
62, 16
12∆ = 8,8 kg
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Construindo-se a tabela das medias ordenadas em ordem de-crescente, tem-se:
Racao Medias (kg)
B 111,25 a
A 103,0833 ab
C 97,3333 b
em que letras iguais indicam medias semelhantes.
Portanto, a racao B difere da racao C, a um nıvel de signi-ficancia de 5%, no ganho de peso dos bovinos, em que a Bproporcionou um ganho maior, em kg.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoVantagens e DesvantagensModelo estatısticoAnalise de VarianciaExemplo
Efeito Suplementos
No caso dos suplementos, basta observar que a media de ga-nho de peso dos animais que foram alimentados com o supri-mento M foi de
yM = 112, 7 kg
e com o suprimento F foi de
yF = 95, 1 kg
Assim, o suprimento M foi mais eficiente no ganho de peso, emKg.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoDesdobramento A/BDesdobramento B/A
Introducao
Quando a hipotese H0 para a interacao entre os fatores e re-jeitada, entao dizemos que a interacao e significativa.
Este resultado implica que os efeitos dos fatores atuam de formadependente, ou seja, o efeito de um fator depende do nıvel dooutro fator.
Assim, nao e recomendado realizar o teste F para cada fatorisoladamente tal como foi apresentado para o caso da interacaonao significativa.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoDesdobramento A/BDesdobramento B/A
O procedimento recomendado e realizar o desdobramento doefeito da interacao.
Para realizar este desdobramento deve-se fazer uma nova analisede variancia em que os nıveis de um fator sao comparados den-tro de cada nıvel do outro fator.
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoDesdobramento A/BDesdobramento B/A
Desdobramento A/B
Tabela 9: Analise de variancia para o desdobramento do fator A dentrode cada nıvel de B.
C.V. G.L. S.Q Q.M. Fcal
B b − 1 SQB QMB =SQBb−1
Fcalc =QMB
QMRes(b)
A|B1 a − 1 SQA|B1QMA|B1
=SQA|B1a−1
Fcal =QMA|B1
QMResComb
A|B2 a − 1 SQA|B2QMA|B2
=SQA|B2a−1
Fcal =QMA|B2
QMResComb
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
A|Bj a − 1 SQA|Bj QMA|Bj =SQA|Bja−1
Fcal =QMA|Bj
QMResComb
ResComb n∗ SQResComb QMResComb =SQResComb
n∗ -
Total SQTotal abn − 1 - -
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoDesdobramento A/BDesdobramento B/A
Para comparar os nıveis de um fator principal em cada nıveldo fator secundario, e necessario fazer uma combinacao dasduas estimativas obtidas para o erro experimental bem comodo numero de graus de liberdade associado as mesmas.
Esta combinacao e denominada de resıduo combinado (ResComb).
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoDesdobramento A/BDesdobramento B/A
A estimativa do quadrado medio deste resıduo combinado eobtida por
QMResComb =QMRes(a) + (b − 1)QMRes(b)
b
O numero de graus de liberdade associado a esta estimativa eobtido pela formula dos graus de liberdade de Satterhwaitte(n∗) dada por
n∗ =
[QMRes(a) + (b − 1)QMRes(b)
]2[QMRes(a)]
2
GLRes(a)+
[(b−1)QMRes(b)]2
GLRes(b)
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais
Parcelas SubdivididasDesdobramentos
IntroducaoDesdobramento A/BDesdobramento B/A
Desdobramento B/A
Tabela 10: Analise de variancia para o desdobramento do fator B dentrode cada nıvel de A.
C.V. G.L. S.Q Q.M. Fcal
A a − 1 SQA QMA =SQAa−1
Fcalc =QMA
QMRes(a)
B|A1 b − 1 SQB|A1QMB|A1
=SQB|A1b−1
Fcal =QMB|A1QMRes(b)
B|A2 b − 1 SQB|A2QMB|A2
=SQB|A2b−1
Fcal =QMB|A2QMRes(b)
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
B|Aj b − 1 SQB|Aj QMB|Aj =SQB|Ajb−1
Fcal =QMB|AjQMRes(b)
Res(b) ab(n − 1) SQRes(b) QMRes(b) =SQRes(b)ab(n−1)
-
Total SQTotal abn − 1 - -
Lucas Santana da Cunha http://www.uel.br/pessoal/lscunha Tecnicas Experimentais com Animais