Estrutura, croqui de campo com FielDHub, ANOVA e comparação de médias
Que cada gole desperte uma nova ideia.
Que cada script abra uma nova conversa.
Que o Café com R se torne um ponto de encontro nosso.
SITE OFICIAL > https://jenniferlopes.github.io/meu_site/
FielDHubmetanemmeans e apresentar o resultado com letrasNo Delineamento Inteiramente Casualizado, DIC, os tratamentos são atribuídos às unidades experimentais de forma completamente aleatória, sem nenhuma restrição ou agrupamento.
Contexto: 5 cultivares de milho são avaliados quanto à produtividade, em gramas por planta, em casa de vegetação, ambiente controlado e homogêneo, sem gradiente ambiental relevante a controlar.
É o delineamento apropriado quando as unidades experimentais são, de fato, homogêneas entre si, condição mais comum em ambientes controlados do que em campo aberto.
\[Y_{ij} = \mu + \tau_i + \epsilon_{ij}\]
| Termo | Significado |
|---|---|
| \(Y_{ij}\) | Observação da repetição \(j\) no tratamento \(i\) |
| \(\mu\) | Média geral |
| \(\tau_i\) | Efeito do tratamento \(i\) |
| \(\epsilon_{ij}\) | Erro experimental, assumido \(\sim N(0, \sigma^2)\) |
Sortear tratamentos “à mão” ou por critérios informais compromete a validade da aleatorização, que é a base teórica que sustenta a inferência estatística de qualquer delineamento.
O pacote FielDHub gera a aleatorização formalmente, documenta a semente utilizada para reprodutibilidade, e produz o croqui de campo, o mapa de onde cada tratamento deve ser instalado.
ID LOCATION PLOT REP TREATMENT
1 1 1 101 3 Cultivar 3
2 2 1 102 6 Cultivar 2
3 3 1 103 6 Cultivar 4
4 4 1 104 2 Cultivar 3
5 5 1 105 4 Cultivar 5
6 6 1 106 4 Cultivar 4
7 7 1 107 2 Cultivar 1
8 8 1 108 4 Cultivar 2
9 9 1 109 3 Cultivar 2
10 10 1 110 5 Cultivar 2
fieldBook documenta exatamente qual cultivar ocupa cada parcela, com a aleatorização registrada e reprodutível pela semente utilizada. O croqui visual facilita a instalação correta do experimento em campo, evitando erro de troca de parcelas.# A tibble: 5 × 11
cultivar variable cv max mean median min sd.amo se ci.t n.valid
<chr> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
1 Cultivar 1 produtiv… 9.33 193. 176. 179. 147. 16.4 6.69 17.2 6
2 Cultivar 2 produtiv… 2.88 212. 202. 201. 195. 5.82 2.38 6.11 6
3 Cultivar 3 produtiv… 11.5 195. 173. 175. 137. 19.8 8.09 20.8 6
4 Cultivar 4 produtiv… 6.75 252. 229. 226. 209. 15.4 6.30 16.2 6
5 Cultivar 5 produtiv… 6.13 227. 213. 216. 199. 13.1 5.34 13.7 6
# A tibble: 2 × 6
fonte_variacao Df `Sum Sq` `Mean Sq` `F value` `Pr(>F)`
<chr> <dbl> <dbl> <dbl> <dbl> <dbl>
1 "cultivar " 4 13848. 3462. 15.7 0
2 "Residuals " 25 5522. 221. NA NA
O quadrado médio de cultivar, dividido pelo quadrado médio do resíduo, produz a estatística F, comparada à distribuição F com os graus de liberdade do modelo.
Um valor de p pequeno indica que ao menos um cultivar difere dos demais quanto à produtividade média, mas não indica, sozinho, quais pares diferem entre si, questão respondida no próximo bloco.
# A tibble: 5 × 7
cultivar emmean SE df lower.CL upper.CL .group
<fct> <dbl> <dbl> <dbl> <dbl> <dbl> <chr>
1 Cultivar 3 173. 6.07 25 156. 190. " a "
2 Cultivar 1 176. 6.07 25 159. 193. " a "
3 Cultivar 2 202. 6.07 25 185 219. " b "
4 Cultivar 5 213. 6.07 25 196. 230. " bc"
5 Cultivar 4 229. 6.07 25 212. 246. " c"
| Etapa | Função no R |
|---|---|
| Croqui de campo | FielDHub::CRD() |
| Estatística descritiva | metan::desc_stat() |
| ANOVA | aov() |
| Comparação de médias | emmeans::cld() |
Continue praticando e explorando.
Esta aula é parte do projeto Café com R. É open source. Use, compartilhe e adapte.
Que cada gole desperte uma nova ideia.
Que cada script abra uma nova conversa.
Que o Café com R se torne um ponto de encontro nosso.
SITE OFICIAL > https://jenniferlopes.github.io/meu_site/

Jennifer Lopes | Café com R