library(dplyr)
library(tidyr)
library(broom)
library(agricolae)
library(lme4)
# library(lmerTest)
library(emmeans)

# Open and structure Data
# Change path
datos <- read.csv("C:/Users/Data/Cultivo_To_Sy.csv", sep=";", dec=",")
datos$Tratamiento <- factor(datos$Treatment, levels= c("Control", "EsSy", "EsTo"))
datos_cat <- split(datos, list(datos$Target))

# Variance per category
aov <- lapply(datos_cat, function(x) lmer(Area ~ Treatment + (1|Dish),  data = x))
tidy_aov_To_Sy_cross <- lapply(aov, function(x) summary(x)$coefficients)
aov_matrix <- do.call(rbind, tidy_aov_To_Sy_cross)
names(aov_matrix) <- names(datos_cat)

hsd_To_Sy_cross  <- lapply(aov, function(aov_tmp) {emmeans(aov_tmp, list(pairwise ~ Treatment), adjust = "tukey")})

hsd_To_Sy_cross