Uso do Pacote ExpDes.pt

Experimento fatorial (A x B) em DBC

• carregar pacote

require(gsheet)

require(ExpDes.pt) # Pacote utilizado para análise de comparação de médias

• Dados para os exemplos, baixados de documento online compartilhado.

url <- 'https://docs.google.com/spreadsheets/d/19yFydp3pdUo48T59ERA6UDmPNaqyOekEkD-SEMfk6-Q/edit?usp=sharing'
dados <- gsheet2tbl(url)

• Exemplo de chamada de dados do computador

  dados=read.table("C:/Rafael/Documentos/Experimentos/AprendizadoR/exemplos/mortalidade.txt", header = T, dec=".", sep="\t")

• fatorial duplo dbc - Teste tukey (a linha de comando é praticamente igual para todos os tipos de delineamento)

  fat2.dbc(dados$tempo,dados$temperatura,dados$bloco,dados$mortas, quali = c(TRUE, TRUE), mcomp = "tukey", fac.names = c("Tempo", "Temperatura"), sigT = 0.05, sigF = 0.05)

## ------------------------------------------------------------------------
## Legenda:
## FATOR 1:  Tempo 
## FATOR 2:  Temperatura 
## ------------------------------------------------------------------------
## 
## 
## Quadro da analise de variancia
## ------------------------------------------------------------------------
##                   GL     SQ     QM     Fc   Pr>Fc
## Bloco              3  31.84 10.615 0.6882 0.56933
## Tempo              3   8.09  2.698 0.1749 0.91217
## Temperatura        1  87.78 87.781 5.6912 0.02655
## Tempo*Temperatura  3   9.84  3.281 0.2127 0.88644
## Residuo           21 323.91 15.424               
## Total             31 461.47                      
## ------------------------------------------------------------------------
## CV = 103.86 %
## 
## ------------------------------------------------------------------------
## Teste de normalidade dos residuos (Shapiro-Wilk)
## p-valor:  0.00195293 
## ATENCAO: a 5% de significancia, os residuos nao podem ser considerados normais!
## ------------------------------------------------------------------------
## 
## Interacao nao significativa: analisando os efeitos simples
## ------------------------------------------------------------------------
## Tempo
## De acordo com o teste F, as medias desse fator sao estatisticamente iguais.
##   Niveis Medias
## 1      0  3.000
## 2      1  4.375
## 3      2  3.750
## 4      3  4.000
## ------------------------------------------------------------------------
## Temperatura
## Teste de Tukey
## ------------------------------------------------------------------------
## Grupos Tratamentos Medias
## a     2   5.4375 
##  b    1   2.125 
## ------------------------------------------------------------------------

• dados não normais - Transformação Box-Cox (variável com +1 pra eliminar 0)

require(MASS)

## Loading required package: MASS

## 
## Attaching package: 'MASS'

## The following object is masked from 'package:ExpDes.pt':
## 
##     ginv

bxc <- boxcox(dados$mortas+1 ~ dados$tempo+dados$temperatura+dados$tempo:temperatura, data = dados, plotit = T)

• Encontrar o valor de lambda que corresponde ao maior y.

bxc$x[which.max(bxc$y)]

## [1] -0.02020202

• ou pode-se aproximar gráfico para ver todos os valores dentro da área de 95%

boxcox(dados$mortas+1 ~ dados$tempo+dados$temperatura+dados$tempo:temperatura, data = dados, lambda = seq(-0.5,0.5, 1/10))

lambda =~ -0.02

fat2.dbc(dados$tempo,dados$temperatura,dados$bloco,dados$mortas ^ (0.1), quali = c(TRUE, TRUE), mcomp = "tukey", fac.names = c("Tempo", "Temperatura"), sigT = 0.05, sigF = 0.05)

## ------------------------------------------------------------------------
## Legenda:
## FATOR 1:  Tempo 
## FATOR 2:  Temperatura 
## ------------------------------------------------------------------------
## 
## 
## Quadro da analise de variancia
## ------------------------------------------------------------------------
##                   GL     SQ      QM     Fc   Pr>Fc
## Bloco              3 0.5707 0.19025 2.2195 0.11580
## Tempo              3 0.4108 0.13692 1.5973 0.21996
## Temperatura        1 0.5683 0.56831 6.6301 0.01766
## Tempo*Temperatura  3 0.2785 0.09283 1.0829 0.37790
## Residuo           21 1.8001 0.08572               
## Total             31 3.6283                       
## ------------------------------------------------------------------------
## CV = 28.73 %
## 
## ------------------------------------------------------------------------
## Teste de normalidade dos residuos (Shapiro-Wilk)
## p-valor:  0.351395 
## De acordo com o teste de Shapiro-Wilk a 5% de significancia, os residuos podem ser considerados normais.
## ------------------------------------------------------------------------
## 
## Interacao nao significativa: analisando os efeitos simples
## ------------------------------------------------------------------------
## Tempo
## De acordo com o teste F, as medias desse fator sao estatisticamente iguais.
##   Niveis    Medias
## 1      0 0.9819125
## 2      1 1.1232554
## 3      2 0.8497683
## 4      3 1.1212765
## ------------------------------------------------------------------------
## Temperatura
## Teste de Tukey
## ------------------------------------------------------------------------
## Grupos Tratamentos Medias
## a     2   1.152319 
##  b    1   0.8857874 
## ------------------------------------------------------------------------

Rafael Henrique Pertille