Parcelas Subdivididas em Delineamento de Blocos ao Acaso

# para realizar uma análise de variância em esquema de parcelas subdivididas vamos requerer o pacote "easyanova" fazendo

require(easyanova)

 

# observação: este pacote pode ser utilizado para análise de parcelas subdivididas nos delineamentos inteiramente ao acaso, blocos ao acaso e quadrado latino, mesmo em caso de perda de unidades experimentais

 

# vamos carregar exemplo do pacote

data(data8)

 

# para visualizar os dados 

data8

 

# observação: estes dados foram obtidos de Kaps and Lamberson (2009) página 386

 

# vamos usar a função "ea2()" e gravar o resultado em um objeto chamado "resultado"

# o primeiro argumento da função define o conjunto de dados, "design" define o delinemanto e "cov" define o tipo de matriz de variância e covariância (veja documentação da função digitando ?ea2) 

resultado=ea2(data8, design=5, cov=4)

 

# observação: na entrada de dados (veja em data8) a primeira coluna deve conter as parcelas, a segunda coluna as repetições das parcelas (neste caso os blocos), a terceira coluna as subparcelas e a quarta coluna a variável resposta

 

# para visualizar o resultado

resultado

 

# na análise de variância observa-se interação não significativa e portanto não é necessário avaliar o desdobramento da mesma

# o fator isolado que foi designado na parcela (plot) foi significativo e deve ser avaliado mais detalhadamente

resultado[6]

resultado[8]

 

# no teste no teste de normalidade e de homogeneidade de variâncias observa-se (p-value>0,05) que os resíduos podem ser considerados aproximadamente normais e com variâncias homogêneas (homocedasticidade). 

 resultado[10]

 

 

# para salvar o resultado em arquivo doc (word)

capture.output(resultado, file="Resultado.doc") 

 

# para salvar o resultado em arquivo txt (bloco de notas)

capture.output(resultado, file="Resultado.txt") 

 

# outros formatos de arquivo para salvar são odt, ods, docx, xls, xlsx

 

# para rodar a análise com várias variáveis resposta altere o argumento "list" para TRUE

resultado=ea2(data8, design=5, cov=4, list=TRUE)


 

# observe a documentação da função fazendo 

?ea2

 
Referências
KAPS, M. and LAMBERSON, W. R. Biostatistics for Animal Science: an introductory text. 2nd Edition. CABI Publishing, Wallingford, Oxfordshire, UK, 2009. 504p.