Extensões de Modelos de Regressão

leg.ufpr.br/~walmes/ensino/extensoes

1 Delineamento de blocos completos casualizados

# Lista de pacotes.
pkg <- c("lattice",
         "latticeExtra",
         "gridExtra",
         "reshape",
         "car",
         "lme4",
         "doBy",
         "multcomp",
         "contrast",
         "wzRfun",
         "labestData")
sapply(pkg,
       FUN = library,
       character.only = TRUE,
       logical.return = TRUE)

1.1 Modelo de efeitos fixos

# da <- read.table("http://www.leg.ufpr.br/~walmes/data/pimentel_batatinha.txt",
#                  header = TRUE, sep = "\t")
# labestDataView()

# help(PimentelEg5.2, help_type = "html")

xtabs(~variedade + bloco, data = PimentelEg5.2)
##              bloco
## variedade     I II III IV
##   B 116-51    1  1   1  1
##   B 1-52      1  1   1  1
##   B 25-50 E   1  1   1  1
##   B 72-53 A   1  1   1  1
##   Buena Vista 1  1   1  1
##   Huinkul     1  1   1  1
##   Kennebec    1  1   1  1
##   S. Rafalela 1  1   1  1
xyplot(producao ~ variedade,
       groups = bloco,
       data = PimentelEg5.2,
       type = "o",
       ylab = expression(Produção ~ (t ~ ha^{-1})),
       xlab = "Variedades de batatinha")

da <- PimentelEg5.2
str(da)
## 'data.frame':    32 obs. of  3 variables:
##  $ bloco    : Factor w/ 4 levels "I","II","III",..: 1 1 1 1 1 1 1 1 2 2 ...
##  $ variedade: Factor w/ 8 levels "B 116-51","B 1-52",..: 1 2 3 4 5 6 7 8 1 2 ...
##  $ producao : num  23.1 20 12.7 18 15.4 21.1 9.2 22.6 24.2 21.1 ...
#-----------------------------------------------------------------------

# Ajuste o modelo linear de efeitos fixos.
m0 <- lm(producao ~ bloco + variedade, data = da)

# Diagnóstico.
par(mfrow = c(2, 2))
plot(m0)