Minicurso: Estatística Experimental com Software R

O planejamento, execução e análise de experimentos é uma competência do profissional das Ciências Agrárias. Para agilizar a análise dos dados obtidos nos experimentos, a utilização de recursos computacionais pode ser de grande utilidade. O curso objetiva analisar dados vinculados a estatística experimental, utilizando o software livre R.

• Prof. Paulo Justiniano Ribeiro Junior, LEG (Laboratório de Estatística e Geoinformação) da UFPR (Universidade Federal do Paraná) CV
  
• Acad. Éder David Borges da Silva CV email: ederdbs@gmail.com
  
• Acad. Renato Gonçalves Oliveira CV
  
• Prof. Jomar Antonio Camarinha Filho CV
  
• Prof. Volnei Pauletti CV
  
• Prof. Ricardo Augusto de Oliveira CV
  

* Viagem sai do agrárias Quinta-feira às 23:30 e retornando sábado ás 7:00.

1. Materiais aulas R Prof. Paulo Justiniano
   • Dia 1:
     • arquivo de dados fictícios
     • alguns comandos do R utilizados na aula
     • Analise Descritiva
     • R Embrapa ler Sessões 4 a 9
   • Dia 2:
     • alguns comandos do R utilizados na aula
     • R Embrapa ler Sessões 4 a 9
2. Apresentação Prof. Volnei Pauletti
   • Apresentação
3. Apresentações
   • Delineamentos
   • Exemplos Analisados
   • Regressão
4. Exemplos de análises de diferentes tipos de experimentos
   1. Delineamento Inteiramente Casualizado
      • Script DICDados DIC
   2. Script Geral do curso
      • Script Geral
   3. Delineamento Blocos ao Acaso
      • Dados DBA
   4. Transformação de dados
      • Dados Trans
   5. Esquema Fatorial
      • Dados Tripla
   6. Delineamento em Parcelas sub-divididas
      • Dados Sub-Dividida
   7. Regressão Linear
      • Dados no Script
   8. Regressão Não Linear
      • Dados no Script
   9. Algumas formas de salvar análises e figuras

## Salvando arquivo com resultados de análises
sink("AnovaDIC.txt")# Com este comando é criado um arquivo de texto no Workspace do R
dados.av<-aov(resp~trat)
summary(dados.av)#Todas as saidas que apareceriam na tela agora vão ser salvas no arquivo
sink()#Com este comando é fechado o arquivo texto com os resultados das analises
 
## Salvando graficos
## em jpg
jpeg("grafico1.jpg")
plot(x,y)
dev.off()
 
## em pdf
pdf("grafico2.pdf")
plot(x,y)
dev.off()

Com o conjunto de dados "exemplo01E.txt", dados retirados do livro "Estatística Aplicada a pesquisa agícola" de Zimmermann (2006), o exemplo encontra-se na pagina 62 e foi modificado, os tratamentos são:

1. Solo sem cultivo previo de arroz, com furadam
   
2. Solo sem cultivo previo de arroz, sem furadam
   
3. Solo sem cultivo previo de arroz, com herbicida
   
4. Solo sem cultivo previo de arroz, sem herbicida
   
5. Solo com cultivo previo de arroz, com herbicida
   

A variavel é materia seca por planta, em gramas, o experimento foi conduzido em casa de vegetação em vasos. Dado O exemplo, obtenha: a análise exploratoria dos dados, análise de variância, e teste de comparação multipla caso for necessario, e obtenha a comparação por contraste:

1. Com versus sem cultivo previo
   
2. Furadam versus herbicida
   
3. Sem furadam versus sem herbida
   

• DIC Exe Dados

Com o conjunto de dados "exemplo02E.txt" que contem dados de produtividade de trigo sobre diferentes doses de esterco liquido bovino, foi conduzido com 4 blocos e 4 doses (0,60,120,180). Obtenha: Análise exploratoria, Análise de variância, transformação via Box-Cox, teste de comparação multipla, regressão linear, e recomendações.

• DBA Exe Dados

A análise de variância de um experimento conduzido no Delineamento em Quadrado Latino. O experimento é sobre produção de grãos de feijão (Kg/parcela) de diferentes variedades de feijão. (Zimmermann) Obtenha: Análise exploratoria, Análise de variância, teste de comparação multipla, e recomendações.

• DQL Exe Dados

O experimento refere-se a um estudo de mudas de Eucalipto onde foi estudado o crescimento (cm) de duas espécies (E1 e E2) em três tipos de recipiente (R1, R2 e R3). O delineamento foi completamente casualizado com 4 repetições. extraido "Experimentação Agricola" Kronha & Banzatto pg 106. Obtenha: Análise exploratoria, Análise de variância, Desdodramento da Interação, teste de comparação multipla, e recomendações.

• Fatorial DIC Exe Dados

O exemplo a seguir é de um delineamento em blocos completos casualizados com 6 repetições no esquema fatorial 2x2x2, no qual foi estudado o efeito da adubação mineral no cafeiro, com Nitrogênio (N), Fósforo (P) e Potássio(K), sendo cada contendo 2 níveis: 0 (ausência) e 1 (presença). Cada parcela era constituída por 12 covas no espaçamento 3,5 x 2,5, com uma área de 105 m², sendo que os dados obtidos foram transformados para kg/ha. Extraido "Experimentação Agricola" Kronha & Banzatto pg 112.

• Fatorial DBA Exe Dados

No arquivo "exemplo07E.txt", estão os dados retirados do livro "Experimentação Agrícola" Kronha & Banzatto, este experimento foi conduzido em bloco ao acaso em parcela subdivididas, o tratamento principal foi adubos:

1. A1-Salitre do Chile
2. A2-Sulfato de amônio
3. A3-Uréia
4. A4-Calnitro

Como tratamento secundario foi adotado 3 Doses:

1. D1- 10 kg de N/ha
2. D2- 40 kg de N/ha
3. D3- 70 kg de N/ha

Estes tratamentos foram avaliados na cultura do milho, onde como variavel resposta foi mensurado a produtividade em kg/ha. Obtenha: ANOVA, desdobramento da ANOVA, testes de comparação multiplas e conclusões sobre o experimento.

• Sub DBA Exe Dados

O exemplo a seguir e retirados do livro "Experimentação Agrícola" Kronha & Banzatto pg 161, o experimento avaliou t/acre de beteraba submetida a 0,80,160 e 320 libras/acre de N e 5 épocas de plantio conduzidos em um experimento em faixas.

• Faixa DBA Exe Dados

No arquivo "exemplo08E.txt" existem dados oriundos da produção de milho em kg/parcela em diferentes espaçamentos e diferentes densidades, sendo o experimento conduzido em faixas. Obtenha: ANOVA e as comparações multiplas caso for necessario.

• Faixa2 DBA Exe Dados

• Solução

• Software R
• agricolae Implementa uma serie de funções interessantes em estatistíca experimental no R.
• ScottKnott Executa o teste ScottKnott no R
• nlme Pacote para análise de modelos mistos no R

• Box-Cox
• Uso do teste de Scott-Knott e da análise de agrupamentos
• PODER E TAXAS DE ERRO TIPO I DOS TESTES SCOTT-KNOTT, TUKEY E STUDENT-NEWMANKEULS, SOB DISTRIBUIÇÕES NORMAL E NÃO NORMAIS DOS RESÍDUOS
• Avaliação de procedimentos de comparações múltiplas em trabalhos publicados na revista Horticultura Brasileira de 1.983 a 2.000.
• AVALIAÇÃO DO PODER E TAXAS DE ERRO TIPO I DO TESTE DE SCOTT-KNOTT POR MEIO DO MÉTODO DE MONTE CARLO
• A precisão experimental relacionada ao uso de bordaduras nas extremidades das fileiras em ensaios de milho
• APLICAÇÃO DA METODOLOGIA DE SUPERFÍCIE DE RESPOSTA NO ESTUDO DA PRODUÇÃO DE PECTINASE POR FERMENTAÇÃO EM ESTADO SÓLIDO DO PEDÚNCULO DE CAJU
• Análise de Superfície de Resposta na Produção de Biodiesel a partir dos Blends do Óleo de Mamona e Babaçu
• EMPREGO DE SUPERFÍCIE DE RESPOSTA NO ESTUDO DE ESPAÇAMENTO E DENSIDADE DE PLANTAS EM FILEIRAS DUPLAS E SIMPLES, EM CULTIVARES DE ALGODOEIRO

Centro Academico de Agronomia Lycio Vellozo - UFPR

• LEG: Laboratório de Estatística e Geoinformação
  
• LABEST: Laboratório de Estatística
  
• C3SL: Centro de Computação Científica e Software Livre
  

• Wiki do R
• A galeria de gráficos do R é uma página que ilustra a diversidade de recursos gráficos do R
• Material sobre o uso do R de Prof. Paulo Justiniano
• Apostila Estatistíca Coputacional Utilizando R, Prof. Daniel Furtado Ferreira - UFLA
• Apostila Recursos Computacionais Utilizando R, Prof. Daniel Furtado Ferreira - UFLA
• Tinn-R GUI/Editor para o ambiente R sob Windows.
• O Xemacs é uma outra opção de editor que facilita a edição de arquivos do e R e disponível para plataformas Linux e Windows e deve ser usado junto com o pacote ESS (Emacs Speaks Statisitcs)
• R_STAT é uma lista de discussão em português sobre o uso do R no yahoo groups
• Reference cards:
  • Cartão de referência rápida de Jonathan Baron
  • Cartão de referência rápida de Tom Short
  • 3-D Graphics in R

• Local: LABEST - Laboratório de Estatística - Centro Politécnico - UFPR

Localizado no centro do mapa abaixo