Analisis de datos bovinos sometidos a diferentes tratamientos.

Este experimento se baso en la alimentacion suplementaria de bovinos usando extractos acuosos de Oregano, tomillo y ajo con el objetivo de verificar carga parasitaria antes y despues de los tratamientos.

library(readr)
library(tidyverse)
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ purrr     1.0.2
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.2     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.1
## ── 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
coprologico <- read_csv("data/coprologico.csv")
## Rows: 30 Columns: 6
## ── Column specification ────────────────────────────────────────────────────────
## Delimiter: ","
## chr (3): Tratamiento, pareja, moo
## dbl (3): Codigo Animal, Coprologico 1, Coprologico 2
## 
## ℹ Use `spec()` to retrieve the full column specification for this data.
## ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message.
coprologico$`Codigo Animal` <- as.factor(coprologico$`Codigo Animal`)
coprologico$pareja <- as.factor(coprologico$pareja)

Visualizaciones de la data de coprologicos

Coprologico 1

coprologico %>% ggplot(aes(x = Tratamiento, y = `Coprologico 1`, fill = `Codigo Animal`)) +
  geom_col(col = "black") +
  facet_grid(pareja~moo, scales = "free_y")

Coprologico 2

coprologico %>% ggplot(aes(x = Tratamiento, y = `Coprologico 2`, fill = `Codigo Animal`)) +
 geom_col(col = "black") +
  facet_grid(pareja~moo, scales = "free_y")
## Warning: Removed 3 rows containing missing values or values outside the scale range
## (`geom_col()`).

Coprologico 1 y 2 juntos HPG (Trichostrongylidae sp.)

coprologico %>% pivot_longer(cols = c("Coprologico 1", "Coprologico 2"), names_to = "test_copro",
                             values_to = "conteo") %>% 
  # group_by(Tratamiento, moo, `Codigo Animal`, test_copro) %>% 
  # summarise(conteo = mean(conteo)) %>% 
  filter(moo == "HPG (Trichostrongylidae sp.)") %>% 
  ggplot(aes(x = Tratamiento, y = conteo)) +
  facet_grid(test_copro ~pareja) +
  geom_col(aes(fill = `Codigo Animal`), col = "black") +
  labs(title = "HPG (Trichostrongylidae sp.)") +
  theme(axis.text.x = element_text(face = "bold", color = "red", angle = 90, hjust = 1))

Coprologico 1 y 2 juntos Oopg (Eimeria sp.)

coprologico %>% pivot_longer(cols = c("Coprologico 1", "Coprologico 2"), names_to = "test_copro",
                             values_to = "conteo") %>% 
  # group_by(Tratamiento, moo, `Codigo Animal`, test_copro) %>% 
  # summarise(conteo = mean(conteo)) %>% 
  filter(moo == "Oopg (Eimeria sp.)") %>% 
  ggplot(aes(x = Tratamiento, y = conteo)) +
  facet_grid(test_copro ~pareja) +
  geom_col(aes(fill = `Codigo Animal`), col = "black") +
  labs(title = "Oopg (Eimeria sp.)") +
  theme(axis.text.x = element_text(face = "bold", color = "red", angle = 45, hjust = 1))

Coprologico 1 y 2 juntos HPG (Moniezia sp.)

coprologico %>% pivot_longer(cols = c("Coprologico 1", "Coprologico 2"), names_to = "test_copro",
                             values_to = "conteo") %>% 
  # group_by(Tratamiento, moo, `Codigo Animal`, test_copro) %>% 
  # summarise(conteo = mean(conteo), na.rm = T) %>% 
  filter(moo == "HPG (Moniezia sp.)") %>% 
  ggplot(aes(x = Tratamiento, y = conteo)) +
  facet_grid(test_copro ~pareja) +
  geom_col(aes(fill = `Codigo Animal`), col = "black") +
  labs(title = "HPG (Moniezia sp.)") +
  theme(axis.text.x = element_text(face = "bold", color = "red", angle = 45, hjust = 1))