Exemplo de Correlação

Introdução

Este é um exemplo de análise de correlação entre severidade de ferrugem (FER) e produtividade (PROD) em ensaio de campo.

Código

library(gsheet)
library(tidyverse)

── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
✔ dplyr     1.1.4     ✔ readr     2.1.5
✔ forcats   1.0.0     ✔ stringr   1.5.1
✔ ggplot2   3.5.2     ✔ tibble    3.2.1
✔ lubridate 1.9.4     ✔ tidyr     1.3.1
✔ purrr     1.0.4     
── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
✖ dplyr::filter() masks stats::filter()
✖ dplyr::lag()    masks stats::lag()
ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

url_campo <- "https://docs.google.com/spreadsheets/d/1bq2N19DcZdtax2fQW9OHSGMR0X2__Z9T/edit?gid=866852711#gid=866852711"
campo <- gsheet2tbl(url_campo)

Correlação

cor(campo$FER, campo$PROD)

[1] -0.6258321

cor.test(campo$FER, campo$PROD)

    Pearson's product-moment correlation

data:  campo$FER and campo$PROD
t = -4.3949, df = 30, p-value = 0.0001277
alternative hypothesis: true correlation is not equal to 0
95 percent confidence interval:
 -0.7999565 -0.3544981
sample estimates:
       cor 
-0.6258321 

Gráfico de dispersão com linha de tendência

ggplot(campo, aes(x = FER, y = PROD)) +
  geom_point(color = "tomato", size = 3) +
  geom_smooth(method = "lm", color = "black") +
  theme_minimal() +
  labs(
    x = "Severidade de Ferrugem (%)",
    y = "Produtividade (kg/ha)",
    title = "Correlação entre severidade e produtividade"
  )

`geom_smooth()` using formula = 'y ~ x'