setwd("/Volumes/GoogleDrive/Mi unidad/Agrosavia/Colaboraciones/Laura/BD Estadística SGR Santander/data")
datos5<-read.table("veinteclosep.csv", header=T, sep=',')
datos5$gen<-as.factor(datos5$gen)
datos5$bloque<-as.factor(datos5$bloque)
datos5$semana<-as.factor(datos5$semana)
attach(datos5)
library(ggplot2)
library(emmeans)
#Gráfica altura
ggplot(datos5, aes(semana, alt, group = gen, colour = gen)) +
  geom_smooth(method="lm", se=F) +
  theme_classic() +
  xlab ("Semana") +
  ylab ("Altura") +
  labs(colour = "Genotipo") +
  theme_linedraw() +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    strip.text = element_text(size = 16),
    axis.title.y = element_text(size = 16),
    axis.title.x = element_text(size = 16),
    axis.text.x = element_text(size = 14),
    axis.text.y = element_text(size = 14)
  ) 
## `geom_smooth()` using formula 'y ~ x'
## Warning: Removed 54 rows containing non-finite values (stat_smooth).

# Gráfica diámetro patrón
ggplot(datos5, aes(semana, patrodia, group = gen, colour = gen)) +
  geom_smooth(method="lm", se=F) +
  theme_classic() +
  xlab ("Semana") +
  ylab ("Diámetro del patrón") +
  labs(colour = "Genotipo") +
  theme_linedraw() +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    strip.text = element_text(size = 16),
    axis.title.y = element_text(size = 16),
    axis.title.x = element_text(size = 16),
    axis.text.x = element_text(size = 14),
    axis.text.y = element_text(size = 14)
  )  
## `geom_smooth()` using formula 'y ~ x'
## Warning: Removed 58 rows containing non-finite values (stat_smooth).

#Gráfica diámetro injerto
ggplot(datos5, aes(semana, injedia, group = gen, colour = gen)) +
  geom_smooth(method="lm", se=F) +
  theme_classic() +
  xlab ("Semana") +
  ylab ("Diámetro del injerto") +
  labs(colour = "Genotipo") +
  theme_linedraw() +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    strip.text = element_text(size = 16),
    axis.title.y = element_text(size = 16),
    axis.title.x = element_text(size = 16),
    axis.text.x = element_text(size = 14),
    axis.text.y = element_text(size = 14)
  ) 
## `geom_smooth()` using formula 'y ~ x'
## Warning: Removed 58 rows containing non-finite values (stat_smooth).

# Gráfica área copa
ggplot(datos5, aes(semana, coparea, group = gen, colour = gen)) +
  geom_smooth(method="lm", se=F) +
  theme_classic() +
  xlab ("Semana") +
  ylab ("Área de la copa") +
  labs(colour = "Genotipo") +
  theme_linedraw() +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    strip.text = element_text(size = 16),
    axis.title.y = element_text(size = 16),
    axis.title.x = element_text(size = 16),
    axis.text.x = element_text(size = 14),
    axis.text.y = element_text(size = 14)
  )
## `geom_smooth()` using formula 'y ~ x'
## Warning: Removed 54 rows containing non-finite values (stat_smooth).

# Gráfica índice de área foliar
ggplot(datos5, aes(semana, iaf, group = gen, colour = gen)) +
  geom_smooth(method="lm", se=F) +
  theme_classic() +
  xlab ("Semana") +
  ylab ("Índice de área foliar") +
  labs(colour = "Genotipo") +
  theme_linedraw() +
  theme(
    plot.title = element_text(hjust = 0.5, size = 16),
    strip.text = element_text(size = 16),
    axis.title.y = element_text(size = 16),
    axis.title.x = element_text(size = 16),
    axis.text.x = element_text(size = 14),
    axis.text.y = element_text(size = 14)
  )
## `geom_smooth()` using formula 'y ~ x'
## Warning: Removed 54 rows containing non-finite values (stat_smooth).

# Anova general
aov.alt<-aov(alt~semana+gen+bloque,na.action=na.exclude)
aov.patrodia<-aov(patrodia~semana+gen+bloque, na.action=na.exclude)
aov.injedia<-aov(injedia~semana+gen+bloque, na.action=na.exclude)
aov.coparea<-aov(coparea~semana+gen+bloque, na.action=na.exclude)
aov.iaf<-aov(iaf~semana+gen+bloque, na.action=na.exclude)
#Análisis para altura
library(nlme)
fit.compsym.alt <- gls(alt ~ semana+gen+bloque, data=datos5, corr=corCompSymm(, form= ~ 1 | gen),na.action=na.exclude)
fit.ar1.alt <- gls(alt ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), na.action=na.exclude)
fit.ar1het.alt <- gls(alt ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), weights=varIdent(form = ~ 1 | semana), na.action=na.exclude)
anova(fit.compsym.alt, fit.ar1.alt, fit.ar1het.alt) #compares the models
##                 Model df      AIC      BIC    logLik   Test  L.Ratio p-value
## fit.compsym.alt     1 26 15780.11 15918.96 -7864.056                        
## fit.ar1.alt         2 26 15742.79 15881.64 -7845.397                        
## fit.ar1het.alt      3 28 15736.65 15886.18 -7840.325 2 vs 3 10.14399  0.0063
anova(fit.ar1.alt)
## Denom. DF: 1541 
##             numDF   F-value p-value
## (Intercept)     1 26681.673  <.0001
## semana          2   237.152  <.0001
## gen            19    28.832  <.0001
## bloque          2   145.935  <.0001
anova(fit.ar1het.alt)
## Denom. DF: 1541 
##             numDF   F-value p-value
## (Intercept)     1 26226.565  <.0001
## semana          2   246.629  <.0001
## gen            19    29.650  <.0001
## bloque          2   146.673  <.0001
# Análisis para diámetro del patrón
fit.compsym.patrodia <- gls(patrodia ~ semana+gen+bloque, data=datos5, corr=corCompSymm(, form= ~ 1 | gen),na.action=na.exclude)
fit.ar1.patrodia <- gls(patrodia ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), na.action=na.exclude)
fit.ar1het.patrodia <- gls(patrodia ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), weights=varIdent(form = ~ 1 | semana), na.action=na.exclude)
anova(fit.compsym.patrodia, fit.ar1.patrodia, fit.ar1het.patrodia) #compares the models
##                      Model df      AIC      BIC    logLik   Test  L.Ratio
## fit.compsym.patrodia     1 26 11869.22 12008.00 -5908.611                
## fit.ar1.patrodia         2 26 11814.36 11953.13 -5881.178                
## fit.ar1het.patrodia      3 28 11787.00 11936.45 -5865.500 2 vs 3 31.35614
##                      p-value
## fit.compsym.patrodia        
## fit.ar1.patrodia            
## fit.ar1het.patrodia   <.0001
anova(fit.ar1.patrodia)
## Denom. DF: 1537 
##             numDF   F-value p-value
## (Intercept)     1 21687.142  <.0001
## semana          2   214.630  <.0001
## gen            19    15.558  <.0001
## bloque          2   194.097  <.0001
anova(fit.ar1het.patrodia)
## Denom. DF: 1537 
##             numDF   F-value p-value
## (Intercept)     1 21449.078  <.0001
## semana          2   234.225  <.0001
## gen            19    15.763  <.0001
## bloque          2   193.952  <.0001
# Análisis para diámetro del injerto
fit.compsym.injedia <- gls(injedia ~ semana+gen+bloque, data=datos5, corr=corCompSymm(, form= ~ 1 | gen),na.action=na.exclude)
fit.ar1.injedia <- gls(injedia ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), na.action=na.exclude)
fit.ar1het.injedia <- gls(injedia ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), weights=varIdent(form = ~ 1 | semana), na.action=na.exclude)
anova(fit.compsym.injedia, fit.ar1.injedia, fit.ar1het.injedia) #compares the models
##                     Model df      AIC      BIC    logLik   Test  L.Ratio
## fit.compsym.injedia     1 26 11159.49 11298.27 -5553.746                
## fit.ar1.injedia         2 26 11072.18 11210.95 -5510.089                
## fit.ar1het.injedia      3 28 11037.89 11187.34 -5490.945 2 vs 3 38.28768
##                     p-value
## fit.compsym.injedia        
## fit.ar1.injedia            
## fit.ar1het.injedia   <.0001
anova(fit.ar1.injedia)
## Denom. DF: 1537 
##             numDF   F-value p-value
## (Intercept)     1 15301.997  <.0001
## semana          2   162.277  <.0001
## gen            19    17.822  <.0001
## bloque          2   228.427  <.0001
anova(fit.ar1het.injedia)
## Denom. DF: 1537 
##             numDF   F-value p-value
## (Intercept)     1 15170.675  <.0001
## semana          2   179.653  <.0001
## gen            19    19.843  <.0001
## bloque          2   240.340  <.0001
# Análisis para área de la copa
fit.compsym.coparea <- gls(coparea ~ semana+gen+bloque, data=datos5, corr=corCompSymm(, form= ~ 1 | gen),na.action=na.exclude)
fit.ar1.coparea <- gls(coparea ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), na.action=na.exclude)
fit.ar1het.coparea <- gls(coparea ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), weights=varIdent(form = ~ 1 | semana), na.action=na.exclude)
anova(fit.compsym.coparea, fit.ar1.coparea, fit.ar1het.coparea) #compares the models
##                     Model df      AIC      BIC    logLik   Test  L.Ratio
## fit.compsym.coparea     1 26 1207.267 1346.112 -577.6334                
## fit.ar1.coparea         2 26 1147.563 1286.408 -547.7814                
## fit.ar1het.coparea      3 28 1128.961 1278.486 -536.4806 2 vs 3 22.60163
##                     p-value
## fit.compsym.coparea        
## fit.ar1.coparea            
## fit.ar1het.coparea   <.0001
anova(fit.ar1.coparea)
## Denom. DF: 1541 
##             numDF   F-value p-value
## (Intercept)     1 14565.860  <.0001
## semana          2   150.721  <.0001
## gen            19    22.740  <.0001
## bloque          2   134.157  <.0001
anova(fit.ar1het.coparea)
## Denom. DF: 1541 
##             numDF   F-value p-value
## (Intercept)     1 14319.360  <.0001
## semana          2   167.581  <.0001
## gen            19    23.436  <.0001
## bloque          2   140.124  <.0001
# Análisis para índice de área foliar
fit.compsym.iaf <- gls(iaf ~ semana+gen+bloque, data=datos5, corr=corCompSymm(, form= ~ 1 | gen),na.action=na.exclude)
fit.ar1.iaf <- gls(iaf ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), na.action=na.exclude)
fit.ar1het.iaf <- gls(iaf ~ semana+gen+bloque, data=datos5, corr=corAR1(, form= ~ 1 | gen), weights=varIdent(form = ~ 1 | semana), na.action=na.exclude)
anova(fit.compsym.iaf, fit.ar1.iaf, fit.ar1het.iaf) #compares the models
##                 Model df      AIC      BIC    logLik   Test  L.Ratio p-value
## fit.compsym.iaf     1 26 11261.69 11400.54 -5604.848                        
## fit.ar1.iaf         2 26 11144.21 11283.06 -5546.105                        
## fit.ar1het.iaf      3 28 10614.21 10763.74 -5279.105 2 vs 3 533.9997  <.0001
anova(fit.ar1.iaf)
## Denom. DF: 1541 
##             numDF   F-value p-value
## (Intercept)     1 1079.6228  <.0001
## semana          2  123.8303  <.0001
## gen            19    5.6458  <.0001
## bloque          2  137.4081  <.0001
anova(fit.ar1het.iaf)
## Denom. DF: 1541 
##             numDF   F-value p-value
## (Intercept)     1 1034.3624  <.0001
## semana          2  104.1960  <.0001
## gen            19    7.3376  <.0001
## bloque          2  130.8218  <.0001
#Tukey altura
library(multcompView)
gen.tuk.alt<-TukeyHSD(aov.alt, "gen", ordered = TRUE)
#Tukey diámetro del patrón
gen.tuk.patrodia<-TukeyHSD(aov.patrodia, "gen", ordered = TRUE)
#Tukey diámetro del injerto
gen.tuk.injedia<-TukeyHSD(aov.injedia, "gen", ordered = TRUE)
#Tukey área de la copa
gen.tuk.coparea<-TukeyHSD(aov.coparea, "gen", ordered = TRUE)
#Tukey índice de área foliar
gen.tuk.iaf<-TukeyHSD(aov.iaf, "gen", ordered = TRUE)
#Etiquetas Tukey altura
#Genotipos
generate_label_df_gen_alt <- function(gen.tuk.alt, variable){
  Tukey.levels <- gen.tuk.alt[[variable]][,4]
  Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])
  Tukey.labels$treatment=rownames(Tukey.labels)
  Tukey.labels=Tukey.labels[order(Tukey.labels$treatment) , ]
  return(Tukey.labels)
}
labels.gen.alt <- generate_label_df_gen_alt(gen.tuk.alt, "gen")
labels.gen.alt
##       Letters treatment
## TCS01      ab     TCS01
## TCS02      bc     TCS02
## TCS03      de     TCS03
## TCS04      de     TCS04
## TCS05      de     TCS05
## TCS08     bcd     TCS08
## TCS10     def     TCS10
## TCS11     efg     TCS11
## TCS12      de     TCS12
## TCS20       h     TCS20
## TCS40      bc     TCS40
## TCS41       a     TCS41
## TCS42       a     TCS42
## TCS43      de     TCS43
## TCS44      de     TCS44
## TCS45      cd     TCS45
## TCS46      fg     TCS46
## TCS47      de     TCS47
## TCS48       g     TCS48
## TCS49      de     TCS49
#Etiquetas Tukey diámetro del patrón
#Genotipos
generate_label_df_gen_patrodia <- function(gen.tuk.patrodia, variable){
  Tukey.levels <- gen.tuk.patrodia[[variable]][,4]
  Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])
  Tukey.labels$treatment=rownames(Tukey.labels)
  Tukey.labels=Tukey.labels[order(Tukey.labels$treatment) , ]
  return(Tukey.labels)
}
labels.gen.patrodia <- generate_label_df_gen_patrodia(gen.tuk.patrodia, "gen")
labels.gen.patrodia
##       Letters treatment
## TCS01      ab     TCS01
## TCS02     abc     TCS02
## TCS03     efg     TCS03
## TCS04     abc     TCS04
## TCS05     efg     TCS05
## TCS08      ab     TCS08
## TCS10     cde     TCS10
## TCS11      fg     TCS11
## TCS12     efg     TCS12
## TCS20       g     TCS20
## TCS40      ab     TCS40
## TCS41       a     TCS41
## TCS42       a     TCS42
## TCS43    cdef     TCS43
## TCS44    defg     TCS44
## TCS45     abc     TCS45
## TCS46     cde     TCS46
## TCS47     bcd     TCS47
## TCS48      fg     TCS48
## TCS49     def     TCS49
#Etiquetas Tukey diámetro del injerto
#Genotipos
generate_label_df_gen_injedia <- function(gen.tuk.injedia, variable){
  Tukey.levels <- gen.tuk.injedia[[variable]][,4]
  Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])
  Tukey.labels$treatment=rownames(Tukey.labels)
  Tukey.labels=Tukey.labels[order(Tukey.labels$treatment) , ]
  return(Tukey.labels)
}
labels.gen.injedia <- generate_label_df_gen_injedia(gen.tuk.injedia, "gen")
labels.gen.injedia
##       Letters treatment
## TCS01      ab     TCS01
## TCS02    bcde     TCS02
## TCS03     fgh     TCS03
## TCS04    cdef     TCS04
## TCS05     fgh     TCS05
## TCS08      ab     TCS08
## TCS10    defg     TCS10
## TCS11      gh     TCS11
## TCS12    defg     TCS12
## TCS20       h     TCS20
## TCS40      ab     TCS40
## TCS41       a     TCS41
## TCS42       a     TCS42
## TCS43     fgh     TCS43
## TCS44     efg     TCS44
## TCS45      bc     TCS45
## TCS46     efg     TCS46
## TCS47     bcd     TCS47
## TCS48     fgh     TCS48
## TCS49     efg     TCS49
#Etiquetas Tukey área de la copa
#Genotipos
generate_label_df_gen_coparea <- function(gen.tuk.coparea, variable){
  Tukey.levels <- gen.tuk.coparea[[variable]][,4]
  Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])
  Tukey.labels$treatment=rownames(Tukey.labels)
  Tukey.labels=Tukey.labels[order(Tukey.labels$treatment) , ]
  return(Tukey.labels)
}
labels.gen.coparea <- generate_label_df_gen_coparea(gen.tuk.coparea, "gen")
labels.gen.coparea
##       Letters treatment
## TCS01     bcd     TCS01
## TCS02     abc     TCS02
## TCS03    efgh     TCS03
## TCS04     bcd     TCS04
## TCS05     fgh     TCS05
## TCS08      ab     TCS08
## TCS10     efg     TCS10
## TCS11     fgh     TCS11
## TCS12     fgh     TCS12
## TCS20       h     TCS20
## TCS40       a     TCS40
## TCS41       a     TCS41
## TCS42       a     TCS42
## TCS43     cde     TCS43
## TCS44      gh     TCS44
## TCS45    efgh     TCS45
## TCS46      gh     TCS46
## TCS47     efg     TCS47
## TCS48     fgh     TCS48
## TCS49     def     TCS49
#Etiquetas Tukey índice de área foliar
#Genotipos
generate_label_df_gen_iaf <- function(gen.tuk.iaf, variable){
  Tukey.levels <- gen.tuk.iaf[[variable]][,4]
  Tukey.labels <- data.frame(multcompLetters(Tukey.levels)['Letters'])
  Tukey.labels$treatment=rownames(Tukey.labels)
  Tukey.labels=Tukey.labels[order(Tukey.labels$treatment) , ]
  return(Tukey.labels)
}
labels.gen.iaf <- generate_label_df_gen_iaf(gen.tuk.iaf, "gen")
labels.gen.iaf
##       Letters treatment
## TCS01   abcde     TCS01
## TCS02   acdef     TCS02
## TCS03    defg     TCS03
## TCS04    abcd     TCS04
## TCS05   abcde     TCS05
## TCS08      ab     TCS08
## TCS10     efg     TCS10
## TCS11     efg     TCS11
## TCS12   cdefg     TCS12
## TCS20      fg     TCS20
## TCS40     abc     TCS40
## TCS41       b     TCS41
## TCS42      ab     TCS42
## TCS43   cdefg     TCS43
## TCS44   cdefg     TCS44
## TCS45   abcde     TCS45
## TCS46       g     TCS46
## TCS47   abcde     TCS47
## TCS48    defg     TCS48
## TCS49       g     TCS49
## Gráficas contrastes de medias Altura
#Gen
contrast <- emmeans(aov.alt, ~gen)
plot(contrast, comparisons = TRUE, xlab ="Altura")

medias.gen <- emmeans(aov.alt, pairwise ~ gen)
medias.gen
## $emmeans
##  gen   emmean   SE   df lower.CL upper.CL
##  TCS01    150 4.50 1541      142      159
##  TCS02    162 4.27 1541      154      171
##  TCS03    191 4.30 1541      182      199
##  TCS04    184 4.27 1541      175      192
##  TCS05    195 4.27 1541      187      203
##  TCS08    170 5.33 1541      160      181
##  TCS10    196 4.27 1541      188      205
##  TCS11    205 4.27 1541      196      213
##  TCS12    193 4.27 1541      184      201
##  TCS20    246 4.27 1541      238      255
##  TCS40    161 4.35 1541      153      170
##  TCS41    140 4.27 1541      131      148
##  TCS42    132 4.30 1541      124      141
##  TCS43    187 4.67 1541      178      196
##  TCS44    189 4.27 1541      180      197
##  TCS45    181 4.27 1541      172      189
##  TCS46    218 4.27 1541      209      226
##  TCS47    189 4.27 1541      181      197
##  TCS48    220 4.27 1541      212      229
##  TCS49    189 4.27 1541      180      197
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast      estimate   SE   df t.ratio p.value
##  TCS01 - TCS02  -11.813 6.20 1541  -1.905  0.9362
##  TCS01 - TCS03  -40.245 6.22 1541  -6.469  <.0001
##  TCS01 - TCS04  -33.344 6.20 1541  -5.376  <.0001
##  TCS01 - TCS05  -44.616 6.20 1541  -7.193  <.0001
##  TCS01 - TCS08  -19.923 6.98 1541  -2.856  0.3188
##  TCS01 - TCS10  -45.986 6.20 1541  -7.414  <.0001
##  TCS01 - TCS11  -54.369 6.20 1541  -8.766  <.0001
##  TCS01 - TCS12  -42.208 6.20 1541  -6.805  <.0001
##  TCS01 - TCS20  -95.826 6.20 1541 -15.450  <.0001
##  TCS01 - TCS40  -10.827 6.26 1541  -1.730  0.9749
##  TCS01 - TCS41   10.779 6.20 1541   1.738  0.9737
##  TCS01 - TCS42   18.150 6.22 1541   2.918  0.2801
##  TCS01 - TCS43  -36.748 6.48 1541  -5.674  <.0001
##  TCS01 - TCS44  -38.221 6.20 1541  -6.162  <.0001
##  TCS01 - TCS45  -30.196 6.20 1541  -4.868  0.0002
##  TCS01 - TCS46  -67.356 6.20 1541 -10.860  <.0001
##  TCS01 - TCS47  -38.554 6.20 1541  -6.216  <.0001
##  TCS01 - TCS48  -69.850 6.20 1541 -11.262  <.0001
##  TCS01 - TCS49  -38.443 6.20 1541  -6.198  <.0001
##  TCS02 - TCS03  -28.432 6.06 1541  -4.695  0.0005
##  TCS02 - TCS04  -21.531 6.04 1541  -3.566  0.0474
##  TCS02 - TCS05  -32.802 6.04 1541  -5.433  <.0001
##  TCS02 - TCS08   -8.109 6.83 1541  -1.187  0.9998
##  TCS02 - TCS10  -34.173 6.04 1541  -5.660  <.0001
##  TCS02 - TCS11  -42.556 6.04 1541  -7.049  <.0001
##  TCS02 - TCS12  -30.395 6.04 1541  -5.035  0.0001
##  TCS02 - TCS20  -84.012 6.04 1541 -13.916  <.0001
##  TCS02 - TCS40    0.986 6.10 1541   0.162  1.0000
##  TCS02 - TCS41   22.593 6.04 1541   3.742  0.0261
##  TCS02 - TCS42   29.963 6.06 1541   4.948  0.0002
##  TCS02 - TCS43  -24.934 6.32 1541  -3.942  0.0126
##  TCS02 - TCS44  -26.407 6.04 1541  -4.374  0.0022
##  TCS02 - TCS45  -18.383 6.04 1541  -3.045  0.2094
##  TCS02 - TCS46  -55.543 6.04 1541  -9.200  <.0001
##  TCS02 - TCS47  -26.741 6.04 1541  -4.429  0.0017
##  TCS02 - TCS48  -58.037 6.04 1541  -9.613  <.0001
##  TCS02 - TCS49  -26.630 6.04 1541  -4.411  0.0018
##  TCS03 - TCS04    6.901 6.06 1541   1.139  0.9999
##  TCS03 - TCS05   -4.371 6.06 1541  -0.722  1.0000
##  TCS03 - TCS08   20.322 6.85 1541   2.968  0.2508
##  TCS03 - TCS10   -5.741 6.06 1541  -0.948  1.0000
##  TCS03 - TCS11  -14.124 6.06 1541  -2.332  0.7117
##  TCS03 - TCS12   -1.963 6.06 1541  -0.324  1.0000
##  TCS03 - TCS20  -55.580 6.06 1541  -9.177  <.0001
##  TCS03 - TCS40   29.418 6.11 1541   4.812  0.0003
##  TCS03 - TCS41   51.024 6.06 1541   8.425  <.0001
##  TCS03 - TCS42   58.395 6.07 1541   9.612  <.0001
##  TCS03 - TCS43    3.498 6.34 1541   0.551  1.0000
##  TCS03 - TCS44    2.024 6.06 1541   0.334  1.0000
##  TCS03 - TCS45   10.049 6.06 1541   1.659  0.9839
##  TCS03 - TCS46  -27.111 6.06 1541  -4.477  0.0014
##  TCS03 - TCS47    1.691 6.06 1541   0.279  1.0000
##  TCS03 - TCS48  -29.605 6.06 1541  -4.888  0.0002
##  TCS03 - TCS49    1.802 6.06 1541   0.298  1.0000
##  TCS04 - TCS05  -11.272 6.04 1541  -1.867  0.9469
##  TCS04 - TCS08   13.421 6.83 1541   1.965  0.9162
##  TCS04 - TCS10  -12.642 6.04 1541  -2.094  0.8599
##  TCS04 - TCS11  -21.025 6.04 1541  -3.482  0.0620
##  TCS04 - TCS12   -8.864 6.04 1541  -1.468  0.9962
##  TCS04 - TCS20  -62.481 6.04 1541 -10.349  <.0001
##  TCS04 - TCS40   22.517 6.10 1541   3.694  0.0309
##  TCS04 - TCS41   44.123 6.04 1541   7.308  <.0001
##  TCS04 - TCS42   51.494 6.06 1541   8.503  <.0001
##  TCS04 - TCS43   -3.403 6.32 1541  -0.538  1.0000
##  TCS04 - TCS44   -4.877 6.04 1541  -0.808  1.0000
##  TCS04 - TCS45    3.148 6.04 1541   0.521  1.0000
##  TCS04 - TCS46  -34.012 6.04 1541  -5.634  <.0001
##  TCS04 - TCS47   -5.210 6.04 1541  -0.863  1.0000
##  TCS04 - TCS48  -36.506 6.04 1541  -6.047  <.0001
##  TCS04 - TCS49   -5.099 6.04 1541  -0.845  1.0000
##  TCS05 - TCS08   24.693 6.83 1541   3.615  0.0404
##  TCS05 - TCS10   -1.370 6.04 1541  -0.227  1.0000
##  TCS05 - TCS11   -9.753 6.04 1541  -1.615  0.9881
##  TCS05 - TCS12    2.407 6.04 1541   0.399  1.0000
##  TCS05 - TCS20  -51.210 6.04 1541  -8.482  <.0001
##  TCS05 - TCS40   33.788 6.10 1541   5.543  <.0001
##  TCS05 - TCS41   55.395 6.04 1541   9.175  <.0001
##  TCS05 - TCS42   62.766 6.06 1541  10.364  <.0001
##  TCS05 - TCS43    7.868 6.32 1541   1.244  0.9996
##  TCS05 - TCS44    6.395 6.04 1541   1.059  1.0000
##  TCS05 - TCS45   14.420 6.04 1541   2.388  0.6702
##  TCS05 - TCS46  -22.741 6.04 1541  -3.767  0.0239
##  TCS05 - TCS47    6.062 6.04 1541   1.004  1.0000
##  TCS05 - TCS48  -25.235 6.04 1541  -4.180  0.0049
##  TCS05 - TCS49    6.173 6.04 1541   1.022  1.0000
##  TCS08 - TCS10  -26.063 6.83 1541  -3.815  0.0201
##  TCS08 - TCS11  -34.446 6.83 1541  -5.043  0.0001
##  TCS08 - TCS12  -22.286 6.83 1541  -3.262  0.1191
##  TCS08 - TCS20  -75.903 6.83 1541 -11.111  <.0001
##  TCS08 - TCS40    9.095 6.88 1541   1.322  0.9990
##  TCS08 - TCS41   30.702 6.83 1541   4.494  0.0013
##  TCS08 - TCS42   38.073 6.85 1541   5.560  <.0001
##  TCS08 - TCS43  -16.825 7.08 1541  -2.375  0.6799
##  TCS08 - TCS44  -18.298 6.83 1541  -2.679  0.4449
##  TCS08 - TCS45  -10.273 6.83 1541  -1.504  0.9949
##  TCS08 - TCS46  -47.434 6.83 1541  -6.944  <.0001
##  TCS08 - TCS47  -18.631 6.83 1541  -2.727  0.4085
##  TCS08 - TCS48  -49.928 6.83 1541  -7.309  <.0001
##  TCS08 - TCS49  -18.520 6.83 1541  -2.711  0.4205
##  TCS10 - TCS11   -8.383 6.04 1541  -1.388  0.9981
##  TCS10 - TCS12    3.778 6.04 1541   0.626  1.0000
##  TCS10 - TCS20  -49.840 6.04 1541  -8.255  <.0001
##  TCS10 - TCS40   35.159 6.10 1541   5.768  <.0001
##  TCS10 - TCS41   56.765 6.04 1541   9.402  <.0001
##  TCS10 - TCS42   64.136 6.06 1541  10.590  <.0001
##  TCS10 - TCS43    9.239 6.32 1541   1.461  0.9964
##  TCS10 - TCS44    7.765 6.04 1541   1.286  0.9993
##  TCS10 - TCS45   15.790 6.04 1541   2.615  0.4935
##  TCS10 - TCS46  -21.370 6.04 1541  -3.540  0.0517
##  TCS10 - TCS47    7.432 6.04 1541   1.231  0.9996
##  TCS10 - TCS48  -23.864 6.04 1541  -3.953  0.0121
##  TCS10 - TCS49    7.543 6.04 1541   1.249  0.9996
##  TCS11 - TCS12   12.160 6.04 1541   2.014  0.8968
##  TCS11 - TCS20  -41.457 6.04 1541  -6.867  <.0001
##  TCS11 - TCS40   43.542 6.10 1541   7.143  <.0001
##  TCS11 - TCS41   65.148 6.04 1541  10.791  <.0001
##  TCS11 - TCS42   72.519 6.06 1541  11.974  <.0001
##  TCS11 - TCS43   17.621 6.32 1541   2.786  0.3664
##  TCS11 - TCS44   16.148 6.04 1541   2.675  0.4479
##  TCS11 - TCS45   24.173 6.04 1541   4.004  0.0099
##  TCS11 - TCS46  -12.988 6.04 1541  -2.151  0.8292
##  TCS11 - TCS47   15.815 6.04 1541   2.620  0.4903
##  TCS11 - TCS48  -15.481 6.04 1541  -2.564  0.5335
##  TCS11 - TCS49   15.926 6.04 1541   2.638  0.4761
##  TCS12 - TCS20  -53.617 6.04 1541  -8.881  <.0001
##  TCS12 - TCS40   31.381 6.10 1541   5.148  0.0001
##  TCS12 - TCS41   52.988 6.04 1541   8.777  <.0001
##  TCS12 - TCS42   60.359 6.06 1541   9.966  <.0001
##  TCS12 - TCS43    5.461 6.32 1541   0.863  1.0000
##  TCS12 - TCS44    3.988 6.04 1541   0.661  1.0000
##  TCS12 - TCS45   12.012 6.04 1541   1.990  0.9067
##  TCS12 - TCS46  -25.148 6.04 1541  -4.165  0.0052
##  TCS12 - TCS47    3.654 6.04 1541   0.605  1.0000
##  TCS12 - TCS48  -27.642 6.04 1541  -4.579  0.0009
##  TCS12 - TCS49    3.765 6.04 1541   0.624  1.0000
##  TCS20 - TCS40   84.998 6.10 1541  13.945  <.0001
##  TCS20 - TCS41  106.605 6.04 1541  17.658  <.0001
##  TCS20 - TCS42  113.976 6.06 1541  18.820  <.0001
##  TCS20 - TCS43   59.078 6.32 1541   9.341  <.0001
##  TCS20 - TCS44   57.605 6.04 1541   9.541  <.0001
##  TCS20 - TCS45   65.630 6.04 1541  10.871  <.0001
##  TCS20 - TCS46   28.469 6.04 1541   4.716  0.0005
##  TCS20 - TCS47   57.272 6.04 1541   9.486  <.0001
##  TCS20 - TCS48   25.975 6.04 1541   4.302  0.0029
##  TCS20 - TCS49   57.383 6.04 1541   9.505  <.0001
##  TCS40 - TCS41   21.607 6.10 1541   3.545  0.0509
##  TCS40 - TCS42   28.977 6.11 1541   4.740  0.0004
##  TCS40 - TCS43  -25.920 6.38 1541  -4.063  0.0079
##  TCS40 - TCS44  -27.393 6.10 1541  -4.494  0.0013
##  TCS40 - TCS45  -19.369 6.10 1541  -3.178  0.1498
##  TCS40 - TCS46  -56.529 6.10 1541  -9.274  <.0001
##  TCS40 - TCS47  -27.727 6.10 1541  -4.549  0.0010
##  TCS40 - TCS48  -59.023 6.10 1541  -9.683  <.0001
##  TCS40 - TCS49  -27.616 6.10 1541  -4.531  0.0011
##  TCS41 - TCS42    7.371 6.06 1541   1.217  0.9997
##  TCS41 - TCS43  -47.527 6.32 1541  -7.514  <.0001
##  TCS41 - TCS44  -49.000 6.04 1541  -8.116  <.0001
##  TCS41 - TCS45  -40.975 6.04 1541  -6.787  <.0001
##  TCS41 - TCS46  -78.136 6.04 1541 -12.942  <.0001
##  TCS41 - TCS47  -49.333 6.04 1541  -8.171  <.0001
##  TCS41 - TCS48  -80.630 6.04 1541 -13.355  <.0001
##  TCS41 - TCS49  -49.222 6.04 1541  -8.153  <.0001
##  TCS42 - TCS43  -54.898 6.34 1541  -8.655  <.0001
##  TCS42 - TCS44  -56.371 6.06 1541  -9.308  <.0001
##  TCS42 - TCS45  -48.346 6.06 1541  -7.983  <.0001
##  TCS42 - TCS46  -85.507 6.06 1541 -14.119  <.0001
##  TCS42 - TCS47  -56.704 6.06 1541  -9.363  <.0001
##  TCS42 - TCS48  -88.000 6.06 1541 -14.531  <.0001
##  TCS42 - TCS49  -56.593 6.06 1541  -9.345  <.0001
##  TCS43 - TCS44   -1.473 6.32 1541  -0.233  1.0000
##  TCS43 - TCS45    6.552 6.32 1541   1.036  1.0000
##  TCS43 - TCS46  -30.609 6.32 1541  -4.839  0.0003
##  TCS43 - TCS47   -1.807 6.32 1541  -0.286  1.0000
##  TCS43 - TCS48  -33.103 6.32 1541  -5.234  <.0001
##  TCS43 - TCS49   -1.695 6.32 1541  -0.268  1.0000
##  TCS44 - TCS45    8.025 6.04 1541   1.329  0.9989
##  TCS44 - TCS46  -29.136 6.04 1541  -4.826  0.0003
##  TCS44 - TCS47   -0.333 6.04 1541  -0.055  1.0000
##  TCS44 - TCS48  -31.630 6.04 1541  -5.239  <.0001
##  TCS44 - TCS49   -0.222 6.04 1541  -0.037  1.0000
##  TCS45 - TCS46  -37.160 6.04 1541  -6.155  <.0001
##  TCS45 - TCS47   -8.358 6.04 1541  -1.384  0.9982
##  TCS45 - TCS48  -39.654 6.04 1541  -6.568  <.0001
##  TCS45 - TCS49   -8.247 6.04 1541  -1.366  0.9985
##  TCS46 - TCS47   28.802 6.04 1541   4.771  0.0004
##  TCS46 - TCS48   -2.494 6.04 1541  -0.413  1.0000
##  TCS46 - TCS49   28.914 6.04 1541   4.789  0.0003
##  TCS47 - TCS48  -31.296 6.04 1541  -5.184  <.0001
##  TCS47 - TCS49    0.111 6.04 1541   0.018  1.0000
##  TCS48 - TCS49   31.407 6.04 1541   5.202  <.0001
## 
## Results are averaged over the levels of: semana, bloque 
## P value adjustment: tukey method for comparing a family of 20 estimates
cld_gen <-multcomp::cld(contrast, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_gen
##  gen   emmean   SE   df lower.CL upper.CL .group     
##  TCS20    246 4.27 1541      238      255  A         
##  TCS48    220 4.27 1541      212      229   B        
##  TCS46    218 4.27 1541      209      226   BC       
##  TCS11    205 4.27 1541      196      213   BCD      
##  TCS10    196 4.27 1541      188      205    CDE     
##  TCS05    195 4.27 1541      187      203     DE     
##  TCS12    193 4.27 1541      184      201     DEF    
##  TCS03    191 4.30 1541      182      199     DEF    
##  TCS47    189 4.27 1541      181      197     DEF    
##  TCS49    189 4.27 1541      180      197     DEF    
##  TCS44    189 4.27 1541      180      197     DEF    
##  TCS43    187 4.67 1541      178      196     DEF    
##  TCS04    184 4.27 1541      175      192     DEF    
##  TCS45    181 4.27 1541      172      189      EFG   
##  TCS08    170 5.33 1541      160      181       FGH  
##  TCS02    162 4.27 1541      154      171        GH  
##  TCS40    161 4.35 1541      153      170        GHI 
##  TCS01    150 4.50 1541      142      159         HIJ
##  TCS41    140 4.27 1541      131      148          IJ
##  TCS42    132 4.30 1541      124      141           J
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 20 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.
## Gráficas contrastes de medias diámetro del patrón
#Gen
contrast <- emmeans(aov.patrodia, ~gen)
plot(contrast, comparisons = TRUE, xlab ="Diámetro del patrón")

medias.gen <- emmeans(aov.patrodia, pairwise ~ gen)
medias.gen
## $emmeans
##  gen   emmean   SE   df lower.CL upper.CL
##  TCS01   41.6 1.28 1537     39.1     44.1
##  TCS02   45.8 1.24 1537     43.4     48.2
##  TCS03   54.8 1.22 1537     52.4     57.2
##  TCS04   44.7 1.21 1537     42.4     47.1
##  TCS05   54.8 1.21 1537     52.4     57.2
##  TCS08   40.1 1.51 1537     37.2     43.1
##  TCS10   50.3 1.21 1537     48.0     52.7
##  TCS11   56.8 1.21 1537     54.4     59.2
##  TCS12   54.2 1.21 1537     51.8     56.6
##  TCS20   58.4 1.21 1537     56.1     60.8
##  TCS40   42.8 1.24 1537     40.4     45.3
##  TCS41   41.3 1.21 1537     38.9     43.6
##  TCS42   39.8 1.22 1537     37.4     42.2
##  TCS43   49.8 1.33 1537     47.2     52.4
##  TCS44   52.6 1.21 1537     50.2     55.0
##  TCS45   44.7 1.21 1537     42.4     47.1
##  TCS46   50.4 1.21 1537     48.0     52.8
##  TCS47   47.6 1.21 1537     45.2     50.0
##  TCS48   57.3 1.21 1537     54.9     59.7
##  TCS49   51.5 1.21 1537     49.1     53.9
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast       estimate   SE   df t.ratio p.value
##  TCS01 - TCS02  -4.23322 1.78 1537  -2.374  0.6809
##  TCS01 - TCS03 -13.26505 1.77 1537  -7.510  <.0001
##  TCS01 - TCS04  -3.16283 1.76 1537  -1.796  0.9634
##  TCS01 - TCS05 -13.23591 1.76 1537  -7.516  <.0001
##  TCS01 - TCS08   1.44661 1.98 1537   0.730  1.0000
##  TCS01 - TCS10  -8.76431 1.76 1537  -4.977  0.0001
##  TCS01 - TCS11 -15.20603 1.76 1537  -8.635  <.0001
##  TCS01 - TCS12 -12.63999 1.76 1537  -7.178  <.0001
##  TCS01 - TCS20 -16.87554 1.76 1537  -9.583  <.0001
##  TCS01 - TCS40  -1.27588 1.78 1537  -0.718  1.0000
##  TCS01 - TCS41   0.31359 1.76 1537   0.178  1.0000
##  TCS01 - TCS42   1.79151 1.77 1537   1.014  1.0000
##  TCS01 - TCS43  -8.27095 1.84 1537  -4.498  0.0013
##  TCS01 - TCS44 -11.01406 1.76 1537  -6.254  <.0001
##  TCS01 - TCS45  -3.16665 1.76 1537  -1.798  0.9629
##  TCS01 - TCS46  -8.83505 1.76 1537  -5.017  0.0001
##  TCS01 - TCS47  -6.04863 1.76 1537  -3.435  0.0719
##  TCS01 - TCS48 -15.71591 1.76 1537  -8.924  <.0001
##  TCS01 - TCS49  -9.94739 1.76 1537  -5.649  <.0001
##  TCS02 - TCS03  -9.03183 1.74 1537  -5.185  <.0001
##  TCS02 - TCS04   1.07039 1.74 1537   0.616  1.0000
##  TCS02 - TCS05  -9.00270 1.74 1537  -5.184  <.0001
##  TCS02 - TCS08   5.67982 1.96 1537   2.898  0.2922
##  TCS02 - TCS10  -4.53109 1.74 1537  -2.609  0.4982
##  TCS02 - TCS11 -10.97282 1.74 1537  -6.319  <.0001
##  TCS02 - TCS12  -8.40677 1.74 1537  -4.841  0.0003
##  TCS02 - TCS20 -12.64233 1.74 1537  -7.280  <.0001
##  TCS02 - TCS40   2.95733 1.75 1537   1.687  0.9808
##  TCS02 - TCS41   4.54681 1.74 1537   2.618  0.4912
##  TCS02 - TCS42   6.02472 1.74 1537   3.459  0.0668
##  TCS02 - TCS43  -4.03774 1.82 1537  -2.221  0.7870
##  TCS02 - TCS44  -6.78085 1.74 1537  -3.905  0.0145
##  TCS02 - TCS45   1.06656 1.74 1537   0.614  1.0000
##  TCS02 - TCS46  -4.60183 1.74 1537  -2.650  0.4667
##  TCS02 - TCS47  -1.81541 1.74 1537  -1.045  1.0000
##  TCS02 - TCS48 -11.48270 1.74 1537  -6.613  <.0001
##  TCS02 - TCS49  -5.71418 1.74 1537  -3.291  0.1100
##  TCS03 - TCS04  10.10223 1.72 1537   5.875  <.0001
##  TCS03 - TCS05   0.02914 1.72 1537   0.017  1.0000
##  TCS03 - TCS08  14.71166 1.94 1537   7.566  <.0001
##  TCS03 - TCS10   4.50074 1.72 1537   2.617  0.4919
##  TCS03 - TCS11  -1.94098 1.72 1537  -1.129  0.9999
##  TCS03 - TCS12   0.62506 1.72 1537   0.364  1.0000
##  TCS03 - TCS20  -3.61049 1.72 1537  -2.100  0.8570
##  TCS03 - TCS40  11.98917 1.74 1537   6.906  <.0001
##  TCS03 - TCS41  13.57864 1.72 1537   7.897  <.0001
##  TCS03 - TCS42  15.05656 1.72 1537   8.729  <.0001
##  TCS03 - TCS43   4.99410 1.80 1537   2.773  0.3756
##  TCS03 - TCS44   2.25099 1.72 1537   1.309  0.9991
##  TCS03 - TCS45  10.09840 1.72 1537   5.873  <.0001
##  TCS03 - TCS46   4.43000 1.72 1537   2.576  0.5240
##  TCS03 - TCS47   7.21642 1.72 1537   4.197  0.0046
##  TCS03 - TCS48  -2.45086 1.72 1537  -1.425  0.9974
##  TCS03 - TCS49   3.31766 1.72 1537   1.929  0.9284
##  TCS04 - TCS05 -10.07309 1.71 1537  -5.876  <.0001
##  TCS04 - TCS08   4.60943 1.94 1537   2.377  0.6791
##  TCS04 - TCS10  -5.60148 1.71 1537  -3.268  0.1173
##  TCS04 - TCS11 -12.04321 1.71 1537  -7.026  <.0001
##  TCS04 - TCS12  -9.47716 1.71 1537  -5.529  <.0001
##  TCS04 - TCS20 -13.71272 1.71 1537  -8.000  <.0001
##  TCS04 - TCS40   1.88695 1.73 1537   1.090  0.9999
##  TCS04 - TCS41   3.47642 1.71 1537   2.028  0.8909
##  TCS04 - TCS42   4.95433 1.72 1537   2.881  0.3027
##  TCS04 - TCS43  -5.10813 1.80 1537  -2.844  0.3265
##  TCS04 - TCS44  -7.85123 1.71 1537  -4.580  0.0009
##  TCS04 - TCS45  -0.00383 1.71 1537  -0.002  1.0000
##  TCS04 - TCS46  -5.67222 1.71 1537  -3.309  0.1044
##  TCS04 - TCS47  -2.88580 1.71 1537  -1.684  0.9812
##  TCS04 - TCS48 -12.55309 1.71 1537  -7.323  <.0001
##  TCS04 - TCS49  -6.78457 1.71 1537  -3.958  0.0118
##  TCS05 - TCS08  14.68252 1.94 1537   7.570  <.0001
##  TCS05 - TCS10   4.47161 1.71 1537   2.609  0.4988
##  TCS05 - TCS11  -1.97012 1.71 1537  -1.149  0.9999
##  TCS05 - TCS12   0.59593 1.71 1537   0.348  1.0000
##  TCS05 - TCS20  -3.63963 1.71 1537  -2.123  0.8446
##  TCS05 - TCS40  11.96003 1.73 1537   6.911  <.0001
##  TCS05 - TCS41  13.54951 1.71 1537   7.904  <.0001
##  TCS05 - TCS42  15.02742 1.72 1537   8.739  <.0001
##  TCS05 - TCS43   4.96496 1.80 1537   2.765  0.3814
##  TCS05 - TCS44   2.22185 1.71 1537   1.296  0.9993
##  TCS05 - TCS45  10.06926 1.71 1537   5.874  <.0001
##  TCS05 - TCS46   4.40086 1.71 1537   2.567  0.5311
##  TCS05 - TCS47   7.18728 1.71 1537   4.193  0.0047
##  TCS05 - TCS48  -2.48000 1.71 1537  -1.447  0.9968
##  TCS05 - TCS49   3.28852 1.71 1537   1.918  0.9319
##  TCS08 - TCS10 -10.21091 1.94 1537  -5.265  <.0001
##  TCS08 - TCS11 -16.65264 1.94 1537  -8.586  <.0001
##  TCS08 - TCS12 -14.08659 1.94 1537  -7.263  <.0001
##  TCS08 - TCS20 -18.32215 1.94 1537  -9.447  <.0001
##  TCS08 - TCS40  -2.72249 1.95 1537  -1.393  0.9980
##  TCS08 - TCS41  -1.13301 1.94 1537  -0.584  1.0000
##  TCS08 - TCS42   0.34490 1.94 1537   0.177  1.0000
##  TCS08 - TCS43  -9.71756 2.01 1537  -4.832  0.0003
##  TCS08 - TCS44 -12.46067 1.94 1537  -6.425  <.0001
##  TCS08 - TCS45  -4.61326 1.94 1537  -2.379  0.6776
##  TCS08 - TCS46 -10.28166 1.94 1537  -5.301  <.0001
##  TCS08 - TCS47  -7.49524 1.94 1537  -3.864  0.0168
##  TCS08 - TCS48 -17.16252 1.94 1537  -8.849  <.0001
##  TCS08 - TCS49 -11.39400 1.94 1537  -5.875  <.0001
##  TCS10 - TCS11  -6.44173 1.71 1537  -3.758  0.0247
##  TCS10 - TCS12  -3.87568 1.71 1537  -2.261  0.7612
##  TCS10 - TCS20  -8.11124 1.71 1537  -4.732  0.0004
##  TCS10 - TCS40   7.48843 1.73 1537   4.327  0.0027
##  TCS10 - TCS41   9.07790 1.71 1537   5.296  <.0001
##  TCS10 - TCS42  10.55581 1.72 1537   6.139  <.0001
##  TCS10 - TCS43   0.49335 1.80 1537   0.275  1.0000
##  TCS10 - TCS44  -2.24975 1.71 1537  -1.312  0.9991
##  TCS10 - TCS45   5.59765 1.71 1537   3.266  0.1180
##  TCS10 - TCS46  -0.07074 1.71 1537  -0.041  1.0000
##  TCS10 - TCS47   2.71568 1.71 1537   1.584  0.9905
##  TCS10 - TCS48  -6.95160 1.71 1537  -4.055  0.0081
##  TCS10 - TCS49  -1.18309 1.71 1537  -0.690  1.0000
##  TCS11 - TCS12   2.56605 1.71 1537   1.497  0.9951
##  TCS11 - TCS20  -1.66951 1.71 1537  -0.974  1.0000
##  TCS11 - TCS40  13.93015 1.73 1537   8.049  <.0001
##  TCS11 - TCS41  15.51963 1.71 1537   9.054  <.0001
##  TCS11 - TCS42  16.99754 1.72 1537   9.885  <.0001
##  TCS11 - TCS43   6.93508 1.80 1537   3.862  0.0170
##  TCS11 - TCS44   4.19198 1.71 1537   2.446  0.6266
##  TCS11 - TCS45  12.03938 1.71 1537   7.024  <.0001
##  TCS11 - TCS46   6.37099 1.71 1537   3.717  0.0285
##  TCS11 - TCS47   9.15741 1.71 1537   5.342  <.0001
##  TCS11 - TCS48  -0.50988 1.71 1537  -0.297  1.0000
##  TCS11 - TCS49   5.25864 1.71 1537   3.068  0.1981
##  TCS12 - TCS20  -4.23556 1.71 1537  -2.471  0.6068
##  TCS12 - TCS40  11.36411 1.73 1537   6.566  <.0001
##  TCS12 - TCS41  12.95358 1.71 1537   7.557  <.0001
##  TCS12 - TCS42  14.43149 1.72 1537   8.393  <.0001
##  TCS12 - TCS43   4.36903 1.80 1537   2.433  0.6363
##  TCS12 - TCS44   1.62593 1.71 1537   0.949  1.0000
##  TCS12 - TCS45   9.47333 1.71 1537   5.527  <.0001
##  TCS12 - TCS46   3.80494 1.71 1537   2.220  0.7881
##  TCS12 - TCS47   6.59136 1.71 1537   3.845  0.0180
##  TCS12 - TCS48  -3.07593 1.71 1537  -1.794  0.9637
##  TCS12 - TCS49   2.69259 1.71 1537   1.571  0.9914
##  TCS20 - TCS40  15.59966 1.73 1537   9.014  <.0001
##  TCS20 - TCS41  17.18914 1.71 1537  10.028  <.0001
##  TCS20 - TCS42  18.66705 1.72 1537  10.856  <.0001
##  TCS20 - TCS43   8.60459 1.80 1537   4.792  0.0003
##  TCS20 - TCS44   5.86148 1.71 1537   3.419  0.0754
##  TCS20 - TCS45  13.70889 1.71 1537   7.997  <.0001
##  TCS20 - TCS46   8.04049 1.71 1537   4.691  0.0005
##  TCS20 - TCS47  10.82691 1.71 1537   6.316  <.0001
##  TCS20 - TCS48   1.15963 1.71 1537   0.677  1.0000
##  TCS20 - TCS49   6.92815 1.71 1537   4.042  0.0086
##  TCS40 - TCS41   1.58947 1.73 1537   0.918  1.0000
##  TCS40 - TCS42   3.06739 1.74 1537   1.767  0.9689
##  TCS40 - TCS43  -6.99507 1.81 1537  -3.861  0.0170
##  TCS40 - TCS44  -9.73818 1.73 1537  -5.627  <.0001
##  TCS40 - TCS45  -1.89077 1.73 1537  -1.093  0.9999
##  TCS40 - TCS46  -7.55917 1.73 1537  -4.368  0.0022
##  TCS40 - TCS47  -4.77275 1.73 1537  -2.758  0.3864
##  TCS40 - TCS48 -14.44003 1.73 1537  -8.344  <.0001
##  TCS40 - TCS49  -8.67151 1.73 1537  -5.011  0.0001
##  TCS41 - TCS42   1.47791 1.72 1537   0.859  1.0000
##  TCS41 - TCS43  -8.58455 1.80 1537  -4.780  0.0003
##  TCS41 - TCS44 -11.32765 1.71 1537  -6.608  <.0001
##  TCS41 - TCS45  -3.48025 1.71 1537  -2.030  0.8899
##  TCS41 - TCS46  -9.14864 1.71 1537  -5.337  <.0001
##  TCS41 - TCS47  -6.36222 1.71 1537  -3.712  0.0291
##  TCS41 - TCS48 -16.02951 1.71 1537  -9.351  <.0001
##  TCS41 - TCS49 -10.26099 1.71 1537  -5.986  <.0001
##  TCS42 - TCS43 -10.06246 1.80 1537  -5.587  <.0001
##  TCS42 - TCS44 -12.80557 1.72 1537  -7.447  <.0001
##  TCS42 - TCS45  -4.95816 1.72 1537  -2.883  0.3013
##  TCS42 - TCS46 -10.62655 1.72 1537  -6.180  <.0001
##  TCS42 - TCS47  -7.84014 1.72 1537  -4.560  0.0010
##  TCS42 - TCS48 -17.50742 1.72 1537 -10.182  <.0001
##  TCS42 - TCS49 -11.73890 1.72 1537  -6.827  <.0001
##  TCS43 - TCS44  -2.74310 1.80 1537  -1.528  0.9938
##  TCS43 - TCS45   5.10430 1.80 1537   2.842  0.3279
##  TCS43 - TCS46  -0.56409 1.80 1537  -0.314  1.0000
##  TCS43 - TCS47   2.22233 1.80 1537   1.238  0.9996
##  TCS43 - TCS48  -7.44496 1.80 1537  -4.146  0.0056
##  TCS43 - TCS49  -1.67644 1.80 1537  -0.934  1.0000
##  TCS44 - TCS45   7.84741 1.71 1537   4.578  0.0009
##  TCS44 - TCS46   2.17901 1.71 1537   1.271  0.9994
##  TCS44 - TCS47   4.96543 1.71 1537   2.897  0.2929
##  TCS44 - TCS48  -4.70185 1.71 1537  -2.743  0.3972
##  TCS44 - TCS49   1.06667 1.71 1537   0.622  1.0000
##  TCS45 - TCS46  -5.66840 1.71 1537  -3.307  0.1050
##  TCS45 - TCS47  -2.88198 1.71 1537  -1.681  0.9815
##  TCS45 - TCS48 -12.54926 1.71 1537  -7.321  <.0001
##  TCS45 - TCS49  -6.78074 1.71 1537  -3.956  0.0119
##  TCS46 - TCS47   2.78642 1.71 1537   1.626  0.9872
##  TCS46 - TCS48  -6.88086 1.71 1537  -4.014  0.0095
##  TCS46 - TCS49  -1.11235 1.71 1537  -0.649  1.0000
##  TCS47 - TCS48  -9.66728 1.71 1537  -5.640  <.0001
##  TCS47 - TCS49  -3.89877 1.71 1537  -2.274  0.7521
##  TCS48 - TCS49   5.76852 1.71 1537   3.365  0.0886
## 
## Results are averaged over the levels of: semana, bloque 
## P value adjustment: tukey method for comparing a family of 20 estimates
cld_gen <-multcomp::cld(contrast, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_gen
##  gen   emmean   SE   df lower.CL upper.CL .group   
##  TCS20   58.4 1.21 1537     56.1     60.8  A       
##  TCS48   57.3 1.21 1537     54.9     59.7  AB      
##  TCS11   56.8 1.21 1537     54.4     59.2  AB      
##  TCS03   54.8 1.22 1537     52.4     57.2  ABC     
##  TCS05   54.8 1.21 1537     52.4     57.2  ABC     
##  TCS12   54.2 1.21 1537     51.8     56.6  ABC     
##  TCS44   52.6 1.21 1537     50.2     55.0  ABCD    
##  TCS49   51.5 1.21 1537     49.1     53.9   BCDE   
##  TCS46   50.4 1.21 1537     48.0     52.8    CDEF  
##  TCS10   50.3 1.21 1537     48.0     52.7    CDEF  
##  TCS43   49.8 1.33 1537     47.2     52.4    CDEF  
##  TCS47   47.6 1.21 1537     45.2     50.0     DEFG 
##  TCS02   45.8 1.24 1537     43.4     48.2      EFGH
##  TCS45   44.7 1.21 1537     42.4     47.1       FGH
##  TCS04   44.7 1.21 1537     42.4     47.1       FGH
##  TCS40   42.8 1.24 1537     40.4     45.3        GH
##  TCS01   41.6 1.28 1537     39.1     44.1        GH
##  TCS41   41.3 1.21 1537     38.9     43.6         H
##  TCS08   40.1 1.51 1537     37.2     43.1         H
##  TCS42   39.8 1.22 1537     37.4     42.2         H
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 20 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.
## Gráficas contrastes de medias diámetro del injerto
#Gen
contrast <- emmeans(aov.injedia, ~gen)
plot(contrast, comparisons = TRUE, xlab ="Diámetro del injerto")

medias.gen <- emmeans(aov.injedia, pairwise ~ gen)
medias.gen
## $emmeans
##  gen   emmean    SE   df lower.CL upper.CL
##  TCS01   28.8 1.014 1537     26.8     30.8
##  TCS02   33.2 0.987 1537     31.3     35.2
##  TCS03   38.3 0.968 1537     36.4     40.2
##  TCS04   35.0 0.962 1537     33.1     36.9
##  TCS05   39.1 0.962 1537     37.2     41.0
##  TCS08   28.1 1.202 1537     25.8     30.5
##  TCS10   36.2 0.962 1537     34.3     38.0
##  TCS11   40.6 0.962 1537     38.7     42.5
##  TCS12   36.2 0.962 1537     34.3     38.0
##  TCS20   43.0 0.962 1537     41.1     44.9
##  TCS40   28.2 0.981 1537     26.2     30.1
##  TCS41   24.5 0.962 1537     22.6     26.4
##  TCS42   24.7 0.968 1537     22.8     26.6
##  TCS43   38.8 1.052 1537     36.7     40.8
##  TCS44   37.3 0.962 1537     35.4     39.2
##  TCS45   30.3 0.962 1537     28.4     32.2
##  TCS46   37.5 0.962 1537     35.6     39.4
##  TCS47   31.7 0.962 1537     29.8     33.6
##  TCS48   38.9 0.962 1537     37.0     40.8
##  TCS49   36.8 0.962 1537     35.0     38.7
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast       estimate   SE   df t.ratio p.value
##  TCS01 - TCS02 -4.44e+00 1.42 1537  -3.138  0.1661
##  TCS01 - TCS03 -9.55e+00 1.40 1537  -6.810  <.0001
##  TCS01 - TCS04 -6.19e+00 1.40 1537  -4.429  0.0017
##  TCS01 - TCS05 -1.03e+01 1.40 1537  -7.387  <.0001
##  TCS01 - TCS08  6.58e-01 1.57 1537   0.419  1.0000
##  TCS01 - TCS10 -7.37e+00 1.40 1537  -5.272  <.0001
##  TCS01 - TCS11 -1.18e+01 1.40 1537  -8.441  <.0001
##  TCS01 - TCS12 -7.37e+00 1.40 1537  -5.272  <.0001
##  TCS01 - TCS20 -1.42e+01 1.40 1537 -10.168  <.0001
##  TCS01 - TCS40  6.29e-01 1.41 1537   0.446  1.0000
##  TCS01 - TCS41  4.33e+00 1.40 1537   3.094  0.1856
##  TCS01 - TCS42  4.06e+00 1.40 1537   2.894  0.2947
##  TCS01 - TCS43 -9.97e+00 1.46 1537  -6.833  <.0001
##  TCS01 - TCS44 -8.52e+00 1.40 1537  -6.096  <.0001
##  TCS01 - TCS45 -1.54e+00 1.40 1537  -1.104  0.9999
##  TCS01 - TCS46 -8.68e+00 1.40 1537  -6.208  <.0001
##  TCS01 - TCS47 -2.93e+00 1.40 1537  -2.096  0.8587
##  TCS01 - TCS48 -1.01e+01 1.40 1537  -7.213  <.0001
##  TCS01 - TCS49 -8.05e+00 1.40 1537  -5.756  <.0001
##  TCS02 - TCS03 -5.11e+00 1.38 1537  -3.694  0.0309
##  TCS02 - TCS04 -1.75e+00 1.38 1537  -1.270  0.9994
##  TCS02 - TCS05 -5.89e+00 1.38 1537  -4.269  0.0034
##  TCS02 - TCS08  5.10e+00 1.56 1537   3.278  0.1140
##  TCS02 - TCS10 -2.93e+00 1.38 1537  -2.124  0.8443
##  TCS02 - TCS11 -7.36e+00 1.38 1537  -5.338  <.0001
##  TCS02 - TCS12 -2.93e+00 1.38 1537  -2.124  0.8440
##  TCS02 - TCS20 -9.77e+00 1.38 1537  -7.090  <.0001
##  TCS02 - TCS40  5.07e+00 1.39 1537   3.644  0.0366
##  TCS02 - TCS41  8.77e+00 1.38 1537   6.360  <.0001
##  TCS02 - TCS42  8.50e+00 1.38 1537   6.147  <.0001
##  TCS02 - TCS43 -5.53e+00 1.44 1537  -3.834  0.0188
##  TCS02 - TCS44 -4.08e+00 1.38 1537  -2.960  0.2549
##  TCS02 - TCS45  2.90e+00 1.38 1537   2.103  0.8554
##  TCS02 - TCS46 -4.24e+00 1.38 1537  -3.073  0.1955
##  TCS02 - TCS47  1.51e+00 1.38 1537   1.096  0.9999
##  TCS02 - TCS48 -5.64e+00 1.38 1537  -4.093  0.0070
##  TCS02 - TCS49 -3.60e+00 1.38 1537  -2.615  0.4937
##  TCS03 - TCS04  3.36e+00 1.37 1537   2.459  0.6158
##  TCS03 - TCS05 -7.78e-01 1.37 1537  -0.570  1.0000
##  TCS03 - TCS08  1.02e+01 1.54 1537   6.613  <.0001
##  TCS03 - TCS10  2.18e+00 1.37 1537   1.597  0.9896
##  TCS03 - TCS11 -2.25e+00 1.37 1537  -1.649  0.9850
##  TCS03 - TCS12  2.18e+00 1.37 1537   1.596  0.9896
##  TCS03 - TCS20 -4.67e+00 1.37 1537  -3.418  0.0756
##  TCS03 - TCS40  1.02e+01 1.38 1537   7.386  <.0001
##  TCS03 - TCS41  1.39e+01 1.37 1537  10.165  <.0001
##  TCS03 - TCS42  1.36e+01 1.37 1537   9.938  <.0001
##  TCS03 - TCS43 -4.25e-01 1.43 1537  -0.297  1.0000
##  TCS03 - TCS44  1.03e+00 1.37 1537   0.752  1.0000
##  TCS03 - TCS45  8.01e+00 1.37 1537   5.865  <.0001
##  TCS03 - TCS46  8.71e-01 1.37 1537   0.638  1.0000
##  TCS03 - TCS47  6.62e+00 1.37 1537   4.849  0.0002
##  TCS03 - TCS48 -5.35e-01 1.37 1537  -0.392  1.0000
##  TCS03 - TCS49  1.50e+00 1.37 1537   1.101  0.9999
##  TCS04 - TCS05 -4.14e+00 1.36 1537  -3.039  0.2124
##  TCS04 - TCS08  6.85e+00 1.54 1537   4.449  0.0016
##  TCS04 - TCS10 -1.18e+00 1.36 1537  -0.865  1.0000
##  TCS04 - TCS11 -5.61e+00 1.36 1537  -4.121  0.0062
##  TCS04 - TCS12 -1.18e+00 1.36 1537  -0.866  1.0000
##  TCS04 - TCS20 -8.02e+00 1.36 1537  -5.896  <.0001
##  TCS04 - TCS40  6.82e+00 1.37 1537   4.965  0.0001
##  TCS04 - TCS41  1.05e+01 1.36 1537   7.729  <.0001
##  TCS04 - TCS42  1.02e+01 1.37 1537   7.509  <.0001
##  TCS04 - TCS43 -3.78e+00 1.43 1537  -2.653  0.4645
##  TCS04 - TCS44 -2.33e+00 1.36 1537  -1.713  0.9774
##  TCS04 - TCS45  4.65e+00 1.36 1537   3.416  0.0761
##  TCS04 - TCS46 -2.49e+00 1.36 1537  -1.827  0.9567
##  TCS04 - TCS47  3.26e+00 1.36 1537   2.397  0.6639
##  TCS04 - TCS48 -3.89e+00 1.36 1537  -2.860  0.3162
##  TCS04 - TCS49 -1.85e+00 1.36 1537  -1.363  0.9985
##  TCS05 - TCS08  1.10e+01 1.54 1537   7.135  <.0001
##  TCS05 - TCS10  2.96e+00 1.36 1537   2.174  0.8162
##  TCS05 - TCS11 -1.47e+00 1.36 1537  -1.082  0.9999
##  TCS05 - TCS12  2.96e+00 1.36 1537   2.173  0.8166
##  TCS05 - TCS20 -3.89e+00 1.36 1537  -2.857  0.3182
##  TCS05 - TCS40  1.10e+01 1.37 1537   7.975  <.0001
##  TCS05 - TCS41  1.47e+01 1.36 1537  10.768  <.0001
##  TCS05 - TCS42  1.44e+01 1.37 1537  10.538  <.0001
##  TCS05 - TCS43  3.53e-01 1.43 1537   0.248  1.0000
##  TCS05 - TCS44  1.80e+00 1.36 1537   1.326  0.9990
##  TCS05 - TCS45  8.78e+00 1.36 1537   6.455  <.0001
##  TCS05 - TCS46  1.65e+00 1.36 1537   1.212  0.9997
##  TCS05 - TCS47  7.40e+00 1.36 1537   5.436  <.0001
##  TCS05 - TCS48  2.43e-01 1.36 1537   0.179  1.0000
##  TCS05 - TCS49  2.28e+00 1.36 1537   1.676  0.9821
##  TCS08 - TCS10 -8.03e+00 1.54 1537  -5.214  <.0001
##  TCS08 - TCS11 -1.25e+01 1.54 1537  -8.091  <.0001
##  TCS08 - TCS12 -8.03e+00 1.54 1537  -5.214  <.0001
##  TCS08 - TCS20 -1.49e+01 1.54 1537  -9.660  <.0001
##  TCS08 - TCS40 -2.93e-02 1.55 1537  -0.019  1.0000
##  TCS08 - TCS41  3.67e+00 1.54 1537   2.382  0.6751
##  TCS08 - TCS42  3.40e+00 1.54 1537   2.203  0.7987
##  TCS08 - TCS43 -1.06e+01 1.60 1537  -6.660  <.0001
##  TCS08 - TCS44 -9.18e+00 1.54 1537  -5.963  <.0001
##  TCS08 - TCS45 -2.20e+00 1.54 1537  -1.430  0.9973
##  TCS08 - TCS46 -9.34e+00 1.54 1537  -6.064  <.0001
##  TCS08 - TCS47 -3.59e+00 1.54 1537  -2.331  0.7125
##  TCS08 - TCS48 -1.07e+01 1.54 1537  -6.977  <.0001
##  TCS08 - TCS49 -8.70e+00 1.54 1537  -5.654  <.0001
##  TCS10 - TCS11 -4.43e+00 1.36 1537  -3.256  0.1212
##  TCS10 - TCS12 -7.41e-04 1.36 1537  -0.001  1.0000
##  TCS10 - TCS20 -6.85e+00 1.36 1537  -5.031  0.0001
##  TCS10 - TCS40  8.00e+00 1.37 1537   5.822  <.0001
##  TCS10 - TCS41  1.17e+01 1.36 1537   8.595  <.0001
##  TCS10 - TCS42  1.14e+01 1.37 1537   8.372  <.0001
##  TCS10 - TCS43 -2.60e+00 1.43 1537  -1.827  0.9567
##  TCS10 - TCS44 -1.15e+00 1.36 1537  -0.847  1.0000
##  TCS10 - TCS45  5.83e+00 1.36 1537   4.282  0.0032
##  TCS10 - TCS46 -1.31e+00 1.36 1537  -0.962  1.0000
##  TCS10 - TCS47  4.44e+00 1.36 1537   3.262  0.1192
##  TCS10 - TCS48 -2.71e+00 1.36 1537  -1.995  0.9047
##  TCS10 - TCS49 -6.77e-01 1.36 1537  -0.498  1.0000
##  TCS11 - TCS12  4.43e+00 1.36 1537   3.255  0.1214
##  TCS11 - TCS20 -2.42e+00 1.36 1537  -1.775  0.9674
##  TCS11 - TCS40  1.24e+01 1.37 1537   9.047  <.0001
##  TCS11 - TCS41  1.61e+01 1.36 1537  11.850  <.0001
##  TCS11 - TCS42  1.59e+01 1.37 1537  11.617  <.0001
##  TCS11 - TCS43  1.83e+00 1.43 1537   1.281  0.9994
##  TCS11 - TCS44  3.28e+00 1.36 1537   2.408  0.6550
##  TCS11 - TCS45  1.03e+01 1.36 1537   7.538  <.0001
##  TCS11 - TCS46  3.12e+00 1.36 1537   2.294  0.7386
##  TCS11 - TCS47  8.87e+00 1.36 1537   6.518  <.0001
##  TCS11 - TCS48  1.72e+00 1.36 1537   1.261  0.9995
##  TCS11 - TCS49  3.75e+00 1.36 1537   2.758  0.3862
##  TCS12 - TCS20 -6.84e+00 1.36 1537  -5.030  0.0001
##  TCS12 - TCS40  8.00e+00 1.37 1537   5.822  <.0001
##  TCS12 - TCS41  1.17e+01 1.36 1537   8.595  <.0001
##  TCS12 - TCS42  1.14e+01 1.37 1537   8.372  <.0001
##  TCS12 - TCS43 -2.60e+00 1.43 1537  -1.826  0.9568
##  TCS12 - TCS44 -1.15e+00 1.36 1537  -0.847  1.0000
##  TCS12 - TCS45  5.83e+00 1.36 1537   4.282  0.0032
##  TCS12 - TCS46 -1.31e+00 1.36 1537  -0.961  1.0000
##  TCS12 - TCS47  4.44e+00 1.36 1537   3.263  0.1190
##  TCS12 - TCS48 -2.71e+00 1.36 1537  -1.994  0.9049
##  TCS12 - TCS49 -6.77e-01 1.36 1537  -0.497  1.0000
##  TCS20 - TCS40  1.48e+01 1.37 1537  10.805  <.0001
##  TCS20 - TCS41  1.85e+01 1.36 1537  13.625  <.0001
##  TCS20 - TCS42  1.83e+01 1.37 1537  13.387  <.0001
##  TCS20 - TCS43  4.24e+00 1.43 1537   2.975  0.2465
##  TCS20 - TCS44  5.69e+00 1.36 1537   4.183  0.0048
##  TCS20 - TCS45  1.27e+01 1.36 1537   9.312  <.0001
##  TCS20 - TCS46  5.54e+00 1.36 1537   4.069  0.0077
##  TCS20 - TCS47  1.13e+01 1.36 1537   8.293  <.0001
##  TCS20 - TCS48  4.13e+00 1.36 1537   3.036  0.2139
##  TCS20 - TCS49  6.17e+00 1.36 1537   4.533  0.0011
##  TCS40 - TCS41  3.70e+00 1.37 1537   2.691  0.4358
##  TCS40 - TCS42  3.43e+00 1.38 1537   2.488  0.5933
##  TCS40 - TCS43 -1.06e+01 1.44 1537  -7.373  <.0001
##  TCS40 - TCS44 -9.15e+00 1.37 1537  -6.661  <.0001
##  TCS40 - TCS45 -2.17e+00 1.37 1537  -1.581  0.9907
##  TCS40 - TCS46 -9.31e+00 1.37 1537  -6.775  <.0001
##  TCS40 - TCS47 -3.56e+00 1.37 1537  -2.591  0.5125
##  TCS40 - TCS48 -1.07e+01 1.37 1537  -7.798  <.0001
##  TCS40 - TCS49 -8.68e+00 1.37 1537  -6.315  <.0001
##  TCS41 - TCS42 -2.68e-01 1.37 1537  -0.196  1.0000
##  TCS41 - TCS43 -1.43e+01 1.43 1537 -10.031  <.0001
##  TCS41 - TCS44 -1.28e+01 1.36 1537  -9.442  <.0001
##  TCS41 - TCS45 -5.87e+00 1.36 1537  -4.313  0.0028
##  TCS41 - TCS46 -1.30e+01 1.36 1537  -9.556  <.0001
##  TCS41 - TCS47 -7.26e+00 1.36 1537  -5.333  <.0001
##  TCS41 - TCS48 -1.44e+01 1.36 1537 -10.589  <.0001
##  TCS41 - TCS49 -1.24e+01 1.36 1537  -9.092  <.0001
##  TCS42 - TCS43 -1.40e+01 1.43 1537  -9.814  <.0001
##  TCS42 - TCS44 -1.26e+01 1.37 1537  -9.216  <.0001
##  TCS42 - TCS45 -5.60e+00 1.37 1537  -4.103  0.0067
##  TCS42 - TCS46 -1.27e+01 1.37 1537  -9.330  <.0001
##  TCS42 - TCS47 -6.99e+00 1.37 1537  -5.120  0.0001
##  TCS42 - TCS48 -1.41e+01 1.37 1537 -10.360  <.0001
##  TCS42 - TCS49 -1.21e+01 1.37 1537  -8.868  <.0001
##  TCS43 - TCS44  1.45e+00 1.43 1537   1.018  1.0000
##  TCS43 - TCS45  8.43e+00 1.43 1537   5.914  <.0001
##  TCS43 - TCS46  1.30e+00 1.43 1537   0.909  1.0000
##  TCS43 - TCS47  7.04e+00 1.43 1537   4.941  0.0002
##  TCS43 - TCS48 -1.10e-01 1.43 1537  -0.077  1.0000
##  TCS43 - TCS49  1.93e+00 1.43 1537   1.352  0.9987
##  TCS44 - TCS45  6.98e+00 1.36 1537   5.129  0.0001
##  TCS44 - TCS46 -1.56e-01 1.36 1537  -0.114  1.0000
##  TCS44 - TCS47  5.59e+00 1.36 1537   4.109  0.0065
##  TCS44 - TCS48 -1.56e+00 1.36 1537  -1.147  0.9999
##  TCS44 - TCS49  4.76e-01 1.36 1537   0.350  1.0000
##  TCS45 - TCS46 -7.14e+00 1.36 1537  -5.243  <.0001
##  TCS45 - TCS47 -1.39e+00 1.36 1537  -1.020  1.0000
##  TCS45 - TCS48 -8.54e+00 1.36 1537  -6.277  <.0001
##  TCS45 - TCS49 -6.50e+00 1.36 1537  -4.779  0.0003
##  TCS46 - TCS47  5.75e+00 1.36 1537   4.224  0.0041
##  TCS46 - TCS48 -1.41e+00 1.36 1537  -1.033  1.0000
##  TCS46 - TCS49  6.31e-01 1.36 1537   0.464  1.0000
##  TCS47 - TCS48 -7.15e+00 1.36 1537  -5.257  <.0001
##  TCS47 - TCS49 -5.12e+00 1.36 1537  -3.760  0.0245
##  TCS48 - TCS49  2.04e+00 1.36 1537   1.497  0.9951
## 
## Results are averaged over the levels of: semana, bloque 
## P value adjustment: tukey method for comparing a family of 20 estimates
cld_gen <-multcomp::cld(contrast, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_gen
##  gen   emmean    SE   df lower.CL upper.CL .group    
##  TCS20   43.0 0.962 1537     41.1     44.9  A        
##  TCS11   40.6 0.962 1537     38.7     42.5  AB       
##  TCS05   39.1 0.962 1537     37.2     41.0  ABC      
##  TCS48   38.9 0.962 1537     37.0     40.8  ABC      
##  TCS43   38.8 1.052 1537     36.7     40.8  ABC      
##  TCS03   38.3 0.968 1537     36.4     40.2  ABC      
##  TCS46   37.5 0.962 1537     35.6     39.4   BCD     
##  TCS44   37.3 0.962 1537     35.4     39.2   BCD     
##  TCS49   36.8 0.962 1537     35.0     38.7   BCD     
##  TCS12   36.2 0.962 1537     34.3     38.0   BCDE    
##  TCS10   36.2 0.962 1537     34.3     38.0   BCDE    
##  TCS04   35.0 0.962 1537     33.1     36.9    CDEF   
##  TCS02   33.2 0.987 1537     31.3     35.2     DEFG  
##  TCS47   31.7 0.962 1537     29.8     33.6      EFGH 
##  TCS45   30.3 0.962 1537     28.4     32.2       FGH 
##  TCS01   28.8 1.014 1537     26.8     30.8        GHI
##  TCS40   28.2 0.981 1537     26.2     30.1         HI
##  TCS08   28.1 1.202 1537     25.8     30.5        GHI
##  TCS42   24.7 0.968 1537     22.8     26.6          I
##  TCS41   24.5 0.962 1537     22.6     26.4          I
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 20 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.
## Gráficas contrastes de medias área de la copa
#Gen
contrast <- emmeans(aov.coparea, ~gen)
plot(contrast, comparisons = TRUE, xlab ="Área de la copa")

medias.gen <- emmeans(aov.coparea, pairwise ~ gen)
medias.gen
## $emmeans
##  gen   emmean     SE   df lower.CL upper.CL
##  TCS01  1.121 0.0398 1541    1.043    1.199
##  TCS02  1.010 0.0377 1541    0.936    1.084
##  TCS03  1.394 0.0380 1541    1.319    1.468
##  TCS04  1.140 0.0377 1541    1.066    1.215
##  TCS05  1.427 0.0377 1541    1.353    1.501
##  TCS08  0.894 0.0471 1541    0.801    0.986
##  TCS10  1.341 0.0377 1541    1.267    1.415
##  TCS11  1.431 0.0377 1541    1.357    1.505
##  TCS12  1.448 0.0377 1541    1.374    1.522
##  TCS20  1.532 0.0377 1541    1.458    1.606
##  TCS40  0.928 0.0385 1541    0.853    1.003
##  TCS41  0.884 0.0377 1541    0.810    0.958
##  TCS42  0.874 0.0380 1541    0.800    0.949
##  TCS43  1.168 0.0413 1541    1.087    1.249
##  TCS44  1.505 0.0377 1541    1.431    1.579
##  TCS45  1.365 0.0377 1541    1.291    1.439
##  TCS46  1.494 0.0377 1541    1.420    1.568
##  TCS47  1.337 0.0377 1541    1.263    1.411
##  TCS48  1.460 0.0377 1541    1.386    1.534
##  TCS49  1.304 0.0377 1541    1.230    1.378
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast      estimate     SE   df t.ratio p.value
##  TCS01 - TCS02  0.11101 0.0548 1541   2.024  0.8925
##  TCS01 - TCS03 -0.27294 0.0550 1541  -4.963  0.0001
##  TCS01 - TCS04 -0.01974 0.0548 1541  -0.360  1.0000
##  TCS01 - TCS05 -0.30665 0.0548 1541  -5.592  <.0001
##  TCS01 - TCS08  0.22683 0.0617 1541   3.678  0.0326
##  TCS01 - TCS10 -0.21998 0.0548 1541  -4.012  0.0096
##  TCS01 - TCS11 -0.30986 0.0548 1541  -5.651  <.0001
##  TCS01 - TCS12 -0.32751 0.0548 1541  -5.973  <.0001
##  TCS01 - TCS20 -0.41159 0.0548 1541  -7.506  <.0001
##  TCS01 - TCS40  0.19278 0.0553 1541   3.484  0.0617
##  TCS01 - TCS41  0.23631 0.0548 1541   4.310  0.0029
##  TCS01 - TCS42  0.24630 0.0550 1541   4.478  0.0014
##  TCS01 - TCS43 -0.04682 0.0573 1541  -0.818  1.0000
##  TCS01 - TCS44 -0.38430 0.0548 1541  -7.008  <.0001
##  TCS01 - TCS45 -0.24455 0.0548 1541  -4.460  0.0015
##  TCS01 - TCS46 -0.37307 0.0548 1541  -6.803  <.0001
##  TCS01 - TCS47 -0.21665 0.0548 1541  -3.951  0.0122
##  TCS01 - TCS48 -0.33961 0.0548 1541  -6.193  <.0001
##  TCS01 - TCS49 -0.18344 0.0548 1541  -3.345  0.0940
##  TCS02 - TCS03 -0.38394 0.0535 1541  -7.171  <.0001
##  TCS02 - TCS04 -0.13074 0.0534 1541  -2.449  0.6235
##  TCS02 - TCS05 -0.41765 0.0534 1541  -7.825  <.0001
##  TCS02 - TCS08  0.11582 0.0604 1541   1.918  0.9321
##  TCS02 - TCS10 -0.33099 0.0534 1541  -6.201  <.0001
##  TCS02 - TCS11 -0.42086 0.0534 1541  -7.885  <.0001
##  TCS02 - TCS12 -0.43852 0.0534 1541  -8.216  <.0001
##  TCS02 - TCS20 -0.52259 0.0534 1541  -9.791  <.0001
##  TCS02 - TCS40  0.08178 0.0539 1541   1.518  0.9943
##  TCS02 - TCS41  0.12531 0.0534 1541   2.348  0.7004
##  TCS02 - TCS42  0.13529 0.0535 1541   2.527  0.5629
##  TCS02 - TCS43 -0.15783 0.0559 1541  -2.823  0.3412
##  TCS02 - TCS44 -0.49531 0.0534 1541  -9.280  <.0001
##  TCS02 - TCS45 -0.35556 0.0534 1541  -6.661  <.0001
##  TCS02 - TCS46 -0.48407 0.0534 1541  -9.069  <.0001
##  TCS02 - TCS47 -0.32765 0.0534 1541  -6.139  <.0001
##  TCS02 - TCS48 -0.45062 0.0534 1541  -8.442  <.0001
##  TCS02 - TCS49 -0.29444 0.0534 1541  -5.517  <.0001
##  TCS03 - TCS04  0.25320 0.0535 1541   4.729  0.0004
##  TCS03 - TCS05 -0.03371 0.0535 1541  -0.630  1.0000
##  TCS03 - TCS08  0.49977 0.0605 1541   8.254  <.0001
##  TCS03 - TCS10  0.05295 0.0535 1541   0.989  1.0000
##  TCS03 - TCS11 -0.03692 0.0535 1541  -0.690  1.0000
##  TCS03 - TCS12 -0.05458 0.0535 1541  -1.019  1.0000
##  TCS03 - TCS20 -0.13865 0.0535 1541  -2.590  0.5137
##  TCS03 - TCS40  0.46572 0.0541 1541   8.616  <.0001
##  TCS03 - TCS41  0.50925 0.0535 1541   9.511  <.0001
##  TCS03 - TCS42  0.51923 0.0537 1541   9.668  <.0001
##  TCS03 - TCS43  0.22611 0.0561 1541   4.032  0.0089
##  TCS03 - TCS44 -0.11137 0.0535 1541  -2.080  0.8669
##  TCS03 - TCS45  0.02839 0.0535 1541   0.530  1.0000
##  TCS03 - TCS46 -0.10013 0.0535 1541  -1.870  0.9460
##  TCS03 - TCS47  0.05629 0.0535 1541   1.051  1.0000
##  TCS03 - TCS48 -0.06667 0.0535 1541  -1.245  0.9996
##  TCS03 - TCS49  0.08950 0.0535 1541   1.672  0.9826
##  TCS04 - TCS05 -0.28691 0.0534 1541  -5.375  <.0001
##  TCS04 - TCS08  0.24656 0.0604 1541   4.083  0.0073
##  TCS04 - TCS10 -0.20025 0.0534 1541  -3.752  0.0252
##  TCS04 - TCS11 -0.29012 0.0534 1541  -5.436  <.0001
##  TCS04 - TCS12 -0.30778 0.0534 1541  -5.766  <.0001
##  TCS04 - TCS20 -0.39185 0.0534 1541  -7.341  <.0001
##  TCS04 - TCS40  0.21252 0.0539 1541   3.944  0.0125
##  TCS04 - TCS41  0.25605 0.0534 1541   4.797  0.0003
##  TCS04 - TCS42  0.26603 0.0535 1541   4.969  0.0001
##  TCS04 - TCS43 -0.02709 0.0559 1541  -0.484  1.0000
##  TCS04 - TCS44 -0.36457 0.0534 1541  -6.830  <.0001
##  TCS04 - TCS45 -0.22481 0.0534 1541  -4.212  0.0043
##  TCS04 - TCS46 -0.35333 0.0534 1541  -6.620  <.0001
##  TCS04 - TCS47 -0.19691 0.0534 1541  -3.689  0.0314
##  TCS04 - TCS48 -0.31988 0.0534 1541  -5.993  <.0001
##  TCS04 - TCS49 -0.16370 0.0534 1541  -3.067  0.1984
##  TCS05 - TCS08  0.53348 0.0604 1541   8.833  <.0001
##  TCS05 - TCS10  0.08667 0.0534 1541   1.624  0.9874
##  TCS05 - TCS11 -0.00321 0.0534 1541  -0.060  1.0000
##  TCS05 - TCS12 -0.02086 0.0534 1541  -0.391  1.0000
##  TCS05 - TCS20 -0.10494 0.0534 1541  -1.966  0.9157
##  TCS05 - TCS40  0.49943 0.0539 1541   9.268  <.0001
##  TCS05 - TCS41  0.54296 0.0534 1541  10.173  <.0001
##  TCS05 - TCS42  0.55295 0.0535 1541  10.327  <.0001
##  TCS05 - TCS43  0.25983 0.0559 1541   4.647  0.0006
##  TCS05 - TCS44 -0.07765 0.0534 1541  -1.455  0.9966
##  TCS05 - TCS45  0.06210 0.0534 1541   1.163  0.9998
##  TCS05 - TCS46 -0.06642 0.0534 1541  -1.244  0.9996
##  TCS05 - TCS47  0.09000 0.0534 1541   1.686  0.9809
##  TCS05 - TCS48 -0.03296 0.0534 1541  -0.618  1.0000
##  TCS05 - TCS49  0.12321 0.0534 1541   2.308  0.7286
##  TCS08 - TCS10 -0.44681 0.0604 1541  -7.398  <.0001
##  TCS08 - TCS11 -0.53669 0.0604 1541  -8.887  <.0001
##  TCS08 - TCS12 -0.55434 0.0604 1541  -9.179  <.0001
##  TCS08 - TCS20 -0.63842 0.0604 1541 -10.571  <.0001
##  TCS08 - TCS40 -0.03405 0.0608 1541  -0.560  1.0000
##  TCS08 - TCS41  0.00949 0.0604 1541   0.157  1.0000
##  TCS08 - TCS42  0.01947 0.0605 1541   0.322  1.0000
##  TCS08 - TCS43 -0.27365 0.0626 1541  -4.370  0.0022
##  TCS08 - TCS44 -0.61113 0.0604 1541 -10.119  <.0001
##  TCS08 - TCS45 -0.47138 0.0604 1541  -7.805  <.0001
##  TCS08 - TCS46 -0.59990 0.0604 1541  -9.933  <.0001
##  TCS08 - TCS47 -0.44348 0.0604 1541  -7.343  <.0001
##  TCS08 - TCS48 -0.56644 0.0604 1541  -9.379  <.0001
##  TCS08 - TCS49 -0.41027 0.0604 1541  -6.793  <.0001
##  TCS10 - TCS11 -0.08988 0.0534 1541  -1.684  0.9812
##  TCS10 - TCS12 -0.10753 0.0534 1541  -2.015  0.8966
##  TCS10 - TCS20 -0.19160 0.0534 1541  -3.590  0.0439
##  TCS10 - TCS40  0.41276 0.0539 1541   7.660  <.0001
##  TCS10 - TCS41  0.45630 0.0534 1541   8.549  <.0001
##  TCS10 - TCS42  0.46628 0.0535 1541   8.709  <.0001
##  TCS10 - TCS43  0.17316 0.0559 1541   3.097  0.1844
##  TCS10 - TCS44 -0.16432 0.0534 1541  -3.079  0.1929
##  TCS10 - TCS45 -0.02457 0.0534 1541  -0.460  1.0000
##  TCS10 - TCS46 -0.15309 0.0534 1541  -2.868  0.3111
##  TCS10 - TCS47  0.00333 0.0534 1541   0.062  1.0000
##  TCS10 - TCS48 -0.11963 0.0534 1541  -2.241  0.7742
##  TCS10 - TCS49  0.03654 0.0534 1541   0.685  1.0000
##  TCS11 - TCS12 -0.01765 0.0534 1541  -0.331  1.0000
##  TCS11 - TCS20 -0.10173 0.0534 1541  -1.906  0.9358
##  TCS11 - TCS40  0.50264 0.0539 1541   9.327  <.0001
##  TCS11 - TCS41  0.54617 0.0534 1541  10.233  <.0001
##  TCS11 - TCS42  0.55616 0.0535 1541  10.387  <.0001
##  TCS11 - TCS43  0.26304 0.0559 1541   4.704  0.0005
##  TCS11 - TCS44 -0.07444 0.0534 1541  -1.395  0.9980
##  TCS11 - TCS45  0.06531 0.0534 1541   1.224  0.9997
##  TCS11 - TCS46 -0.06321 0.0534 1541  -1.184  0.9998
##  TCS11 - TCS47  0.09321 0.0534 1541   1.746  0.9724
##  TCS11 - TCS48 -0.02975 0.0534 1541  -0.557  1.0000
##  TCS11 - TCS49  0.12642 0.0534 1541   2.369  0.6850
##  TCS12 - TCS20 -0.08407 0.0534 1541  -1.575  0.9911
##  TCS12 - TCS40  0.52030 0.0539 1541   9.655  <.0001
##  TCS12 - TCS41  0.56383 0.0534 1541  10.564  <.0001
##  TCS12 - TCS42  0.57381 0.0535 1541  10.717  <.0001
##  TCS12 - TCS43  0.28069 0.0559 1541   5.020  0.0001
##  TCS12 - TCS44 -0.05679 0.0534 1541  -1.064  1.0000
##  TCS12 - TCS45  0.08296 0.0534 1541   1.554  0.9924
##  TCS12 - TCS46 -0.04556 0.0534 1541  -0.854  1.0000
##  TCS12 - TCS47  0.11086 0.0534 1541   2.077  0.8683
##  TCS12 - TCS48 -0.01210 0.0534 1541  -0.227  1.0000
##  TCS12 - TCS49  0.14407 0.0534 1541   2.699  0.4294
##  TCS20 - TCS40  0.60437 0.0539 1541  11.215  <.0001
##  TCS20 - TCS41  0.64790 0.0534 1541  12.139  <.0001
##  TCS20 - TCS42  0.65788 0.0535 1541  12.287  <.0001
##  TCS20 - TCS43  0.36476 0.0559 1541   6.523  <.0001
##  TCS20 - TCS44  0.02728 0.0534 1541   0.511  1.0000
##  TCS20 - TCS45  0.16704 0.0534 1541   3.130  0.1698
##  TCS20 - TCS46  0.03852 0.0534 1541   0.722  1.0000
##  TCS20 - TCS47  0.19494 0.0534 1541   3.652  0.0356
##  TCS20 - TCS48  0.07198 0.0534 1541   1.348  0.9987
##  TCS20 - TCS49  0.22815 0.0534 1541   4.274  0.0033
##  TCS40 - TCS41  0.04353 0.0539 1541   0.808  1.0000
##  TCS40 - TCS42  0.05352 0.0541 1541   0.990  1.0000
##  TCS40 - TCS43 -0.23961 0.0564 1541  -4.248  0.0037
##  TCS40 - TCS44 -0.57709 0.0539 1541 -10.709  <.0001
##  TCS40 - TCS45 -0.43733 0.0539 1541  -8.116  <.0001
##  TCS40 - TCS46 -0.56585 0.0539 1541 -10.500  <.0001
##  TCS40 - TCS47 -0.40943 0.0539 1541  -7.598  <.0001
##  TCS40 - TCS48 -0.53239 0.0539 1541  -9.880  <.0001
##  TCS40 - TCS49 -0.37622 0.0539 1541  -6.981  <.0001
##  TCS41 - TCS42  0.00998 0.0535 1541   0.186  1.0000
##  TCS41 - TCS43 -0.28314 0.0559 1541  -5.064  0.0001
##  TCS41 - TCS44 -0.62062 0.0534 1541 -11.628  <.0001
##  TCS41 - TCS45 -0.48086 0.0534 1541  -9.009  <.0001
##  TCS41 - TCS46 -0.60938 0.0534 1541 -11.417  <.0001
##  TCS41 - TCS47 -0.45296 0.0534 1541  -8.486  <.0001
##  TCS41 - TCS48 -0.57593 0.0534 1541 -10.790  <.0001
##  TCS41 - TCS49 -0.41975 0.0534 1541  -7.864  <.0001
##  TCS42 - TCS43 -0.29312 0.0561 1541  -5.227  <.0001
##  TCS42 - TCS44 -0.63060 0.0535 1541 -11.778  <.0001
##  TCS42 - TCS45 -0.49085 0.0535 1541  -9.168  <.0001
##  TCS42 - TCS46 -0.61937 0.0535 1541 -11.568  <.0001
##  TCS42 - TCS47 -0.46295 0.0535 1541  -8.646  <.0001
##  TCS42 - TCS48 -0.58591 0.0535 1541 -10.943  <.0001
##  TCS42 - TCS49 -0.42974 0.0535 1541  -8.026  <.0001
##  TCS43 - TCS44 -0.33748 0.0559 1541  -6.035  <.0001
##  TCS43 - TCS45 -0.19773 0.0559 1541  -3.536  0.0523
##  TCS43 - TCS46 -0.32625 0.0559 1541  -5.834  <.0001
##  TCS43 - TCS47 -0.16983 0.0559 1541  -3.037  0.2133
##  TCS43 - TCS48 -0.29279 0.0559 1541  -5.236  <.0001
##  TCS43 - TCS49 -0.13662 0.0559 1541  -2.443  0.6284
##  TCS44 - TCS45  0.13975 0.0534 1541   2.618  0.4912
##  TCS44 - TCS46  0.01123 0.0534 1541   0.210  1.0000
##  TCS44 - TCS47  0.16765 0.0534 1541   3.141  0.1648
##  TCS44 - TCS48  0.04469 0.0534 1541   0.837  1.0000
##  TCS44 - TCS49  0.20086 0.0534 1541   3.763  0.0242
##  TCS45 - TCS46 -0.12852 0.0534 1541  -2.408  0.6555
##  TCS45 - TCS47  0.02790 0.0534 1541   0.523  1.0000
##  TCS45 - TCS48 -0.09506 0.0534 1541  -1.781  0.9663
##  TCS45 - TCS49  0.06111 0.0534 1541   1.145  0.9999
##  TCS46 - TCS47  0.15642 0.0534 1541   2.931  0.2723
##  TCS46 - TCS48  0.03346 0.0534 1541   0.627  1.0000
##  TCS46 - TCS49  0.18963 0.0534 1541   3.553  0.0496
##  TCS47 - TCS48 -0.12296 0.0534 1541  -2.304  0.7319
##  TCS47 - TCS49  0.03321 0.0534 1541   0.622  1.0000
##  TCS48 - TCS49  0.15617 0.0534 1541   2.926  0.2750
## 
## Results are averaged over the levels of: semana, bloque 
## P value adjustment: tukey method for comparing a family of 20 estimates
cld_gen <-multcomp::cld(contrast, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_gen
##  gen   emmean     SE   df lower.CL upper.CL .group   
##  TCS20  1.532 0.0377 1541    1.458    1.606  A       
##  TCS44  1.505 0.0377 1541    1.431    1.579  AB      
##  TCS46  1.494 0.0377 1541    1.420    1.568  AB      
##  TCS48  1.460 0.0377 1541    1.386    1.534  ABC     
##  TCS12  1.448 0.0377 1541    1.374    1.522  ABC     
##  TCS11  1.431 0.0377 1541    1.357    1.505  ABC     
##  TCS05  1.427 0.0377 1541    1.353    1.501  ABC     
##  TCS03  1.394 0.0380 1541    1.319    1.468  ABC     
##  TCS45  1.365 0.0377 1541    1.291    1.439  ABCD    
##  TCS10  1.341 0.0377 1541    1.267    1.415   BCD    
##  TCS47  1.337 0.0377 1541    1.263    1.411   BCD    
##  TCS49  1.304 0.0377 1541    1.230    1.378    CDE   
##  TCS43  1.168 0.0413 1541    1.087    1.249     DEF  
##  TCS04  1.140 0.0377 1541    1.066    1.215      EF  
##  TCS01  1.121 0.0398 1541    1.043    1.199      EFG 
##  TCS02  1.010 0.0377 1541    0.936    1.084       FGH
##  TCS40  0.928 0.0385 1541    0.853    1.003        GH
##  TCS08  0.894 0.0471 1541    0.801    0.986         H
##  TCS41  0.884 0.0377 1541    0.810    0.958         H
##  TCS42  0.874 0.0380 1541    0.800    0.949         H
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 20 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.
## Gráficas contrastes de medias índice de área foliar
#Gen
contrast <- emmeans(aov.iaf, ~gen)
plot(contrast, comparisons = TRUE, xlab ="Índice de área foliar")

medias.gen <- emmeans(aov.iaf, pairwise ~ gen)
medias.gen
## $emmeans
##  gen   emmean    SE   df lower.CL upper.CL
##  TCS01   7.48 1.039 1541     5.45     9.52
##  TCS02   9.66 0.985 1541     7.73    11.59
##  TCS03  11.42 0.992 1541     9.48    13.37
##  TCS04   7.05 0.985 1541     5.11     8.98
##  TCS05   8.36 0.985 1541     6.42    10.29
##  TCS08   3.69 1.231 1541     1.28     6.11
##  TCS10  12.19 0.985 1541    10.26    14.12
##  TCS11  12.24 0.985 1541    10.30    14.17
##  TCS12  10.56 0.985 1541     8.62    12.49
##  TCS20  13.51 0.985 1541    11.58    15.44
##  TCS40   6.03 1.004 1541     4.06     8.00
##  TCS41   4.69 0.985 1541     2.76     6.62
##  TCS42   4.90 0.992 1541     2.95     6.84
##  TCS43  10.39 1.077 1541     8.28    12.51
##  TCS44  10.61 0.985 1541     8.68    12.54
##  TCS45   8.26 0.985 1541     6.33    10.19
##  TCS46  14.84 0.985 1541    12.90    16.77
##  TCS47   8.24 0.985 1541     6.31    10.18
##  TCS48  11.89 0.985 1541     9.95    13.82
##  TCS49  15.31 0.985 1541    13.38    17.24
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## 
## $contrasts
##  contrast      estimate   SE   df t.ratio p.value
##  TCS01 - TCS02  -2.1753 1.43 1541  -1.519  0.9942
##  TCS01 - TCS03  -3.9391 1.44 1541  -2.743  0.3970
##  TCS01 - TCS04   0.4370 1.43 1541   0.305  1.0000
##  TCS01 - TCS05  -0.8733 1.43 1541  -0.610  1.0000
##  TCS01 - TCS08   3.7919 1.61 1541   2.355  0.6950
##  TCS01 - TCS10  -4.7072 1.43 1541  -3.288  0.1109
##  TCS01 - TCS11  -4.7525 1.43 1541  -3.319  0.1013
##  TCS01 - TCS12  -3.0744 1.43 1541  -2.147  0.8314
##  TCS01 - TCS20  -6.0276 1.43 1541  -4.210  0.0043
##  TCS01 - TCS40   1.4562 1.44 1541   1.008  1.0000
##  TCS01 - TCS41   2.7941 1.43 1541   1.952  0.9209
##  TCS01 - TCS42   2.5853 1.44 1541   1.800  0.9625
##  TCS01 - TCS43  -2.9097 1.49 1541  -1.946  0.9227
##  TCS01 - TCS44  -3.1252 1.43 1541  -2.183  0.8108
##  TCS01 - TCS45  -0.7747 1.43 1541  -0.541  1.0000
##  TCS01 - TCS46  -7.3542 1.43 1541  -5.137  0.0001
##  TCS01 - TCS47  -0.7610 1.43 1541  -0.532  1.0000
##  TCS01 - TCS48  -4.4026 1.43 1541  -3.075  0.1946
##  TCS01 - TCS49  -7.8265 1.43 1541  -5.466  <.0001
##  TCS02 - TCS03  -1.7638 1.40 1541  -1.262  0.9995
##  TCS02 - TCS04   2.6123 1.39 1541   1.875  0.9448
##  TCS02 - TCS05   1.3020 1.39 1541   0.934  1.0000
##  TCS02 - TCS08   5.9672 1.58 1541   3.784  0.0225
##  TCS02 - TCS10  -2.5319 1.39 1541  -1.817  0.9590
##  TCS02 - TCS11  -2.5772 1.39 1541  -1.849  0.9514
##  TCS02 - TCS12  -0.8991 1.39 1541  -0.645  1.0000
##  TCS02 - TCS20  -3.8523 1.39 1541  -2.764  0.3818
##  TCS02 - TCS40   3.6315 1.41 1541   2.581  0.5203
##  TCS02 - TCS41   4.9694 1.39 1541   3.566  0.0475
##  TCS02 - TCS42   4.7606 1.40 1541   3.405  0.0786
##  TCS02 - TCS43  -0.7344 1.46 1541  -0.503  1.0000
##  TCS02 - TCS44  -0.9499 1.39 1541  -0.682  1.0000
##  TCS02 - TCS45   1.4006 1.39 1541   1.005  1.0000
##  TCS02 - TCS46  -5.1789 1.39 1541  -3.716  0.0286
##  TCS02 - TCS47   1.4143 1.39 1541   1.015  1.0000
##  TCS02 - TCS48  -2.2273 1.39 1541  -1.598  0.9895
##  TCS02 - TCS49  -5.6512 1.39 1541  -4.055  0.0081
##  TCS03 - TCS04   4.3762 1.40 1541   3.130  0.1694
##  TCS03 - TCS05   3.0658 1.40 1541   2.193  0.8046
##  TCS03 - TCS08   7.7311 1.58 1541   4.891  0.0002
##  TCS03 - TCS10  -0.7680 1.40 1541  -0.549  1.0000
##  TCS03 - TCS11  -0.8133 1.40 1541  -0.582  1.0000
##  TCS03 - TCS12   0.8647 1.40 1541   0.619  1.0000
##  TCS03 - TCS20  -2.0885 1.40 1541  -1.494  0.9953
##  TCS03 - TCS40   5.3954 1.41 1541   3.823  0.0196
##  TCS03 - TCS41   6.7332 1.40 1541   4.816  0.0003
##  TCS03 - TCS42   6.5245 1.40 1541   4.653  0.0006
##  TCS03 - TCS43   1.0294 1.46 1541   0.703  1.0000
##  TCS03 - TCS44   0.8140 1.40 1541   0.582  1.0000
##  TCS03 - TCS45   3.1645 1.40 1541   2.264  0.7594
##  TCS03 - TCS46  -3.4150 1.40 1541  -2.443  0.6286
##  TCS03 - TCS47   3.1782 1.40 1541   2.273  0.7528
##  TCS03 - TCS48  -0.4634 1.40 1541  -0.332  1.0000
##  TCS03 - TCS49  -3.8874 1.40 1541  -2.781  0.3701
##  TCS04 - TCS05  -1.3104 1.39 1541  -0.940  1.0000
##  TCS04 - TCS08   3.3549 1.58 1541   2.128  0.8423
##  TCS04 - TCS10  -5.1442 1.39 1541  -3.691  0.0312
##  TCS04 - TCS11  -5.1895 1.39 1541  -3.724  0.0278
##  TCS04 - TCS12  -3.5115 1.39 1541  -2.520  0.5686
##  TCS04 - TCS20  -6.4647 1.39 1541  -4.639  0.0007
##  TCS04 - TCS40   1.0192 1.41 1541   0.724  1.0000
##  TCS04 - TCS41   2.3570 1.39 1541   1.691  0.9803
##  TCS04 - TCS42   2.1483 1.40 1541   1.537  0.9933
##  TCS04 - TCS43  -3.3467 1.46 1541  -2.292  0.7398
##  TCS04 - TCS44  -3.5622 1.39 1541  -2.556  0.5399
##  TCS04 - TCS45  -1.2117 1.39 1541  -0.869  1.0000
##  TCS04 - TCS46  -7.7912 1.39 1541  -5.591  <.0001
##  TCS04 - TCS47  -1.1980 1.39 1541  -0.860  1.0000
##  TCS04 - TCS48  -4.8396 1.39 1541  -3.473  0.0639
##  TCS04 - TCS49  -8.2636 1.39 1541  -5.930  <.0001
##  TCS05 - TCS08   4.6653 1.58 1541   2.959  0.2559
##  TCS05 - TCS10  -3.8338 1.39 1541  -2.751  0.3913
##  TCS05 - TCS11  -3.8791 1.39 1541  -2.784  0.3681
##  TCS05 - TCS12  -2.2011 1.39 1541  -1.579  0.9908
##  TCS05 - TCS20  -5.1543 1.39 1541  -3.699  0.0304
##  TCS05 - TCS40   2.3296 1.41 1541   1.656  0.9843
##  TCS05 - TCS41   3.6674 1.39 1541   2.632  0.4809
##  TCS05 - TCS42   3.4587 1.40 1541   2.474  0.6043
##  TCS05 - TCS43  -2.0364 1.46 1541  -1.395  0.9980
##  TCS05 - TCS44  -2.2519 1.39 1541  -1.616  0.9881
##  TCS05 - TCS45   0.0986 1.39 1541   0.071  1.0000
##  TCS05 - TCS46  -6.4809 1.39 1541  -4.650  0.0006
##  TCS05 - TCS47   0.1123 1.39 1541   0.081  1.0000
##  TCS05 - TCS48  -3.5293 1.39 1541  -2.532  0.5585
##  TCS05 - TCS49  -6.9532 1.39 1541  -4.989  0.0001
##  TCS08 - TCS10  -8.4991 1.58 1541  -5.390  <.0001
##  TCS08 - TCS11  -8.5444 1.58 1541  -5.419  <.0001
##  TCS08 - TCS12  -6.8664 1.58 1541  -4.355  0.0024
##  TCS08 - TCS20  -9.8196 1.58 1541  -6.227  <.0001
##  TCS08 - TCS40  -2.3357 1.59 1541  -1.470  0.9961
##  TCS08 - TCS41  -0.9978 1.58 1541  -0.633  1.0000
##  TCS08 - TCS42  -1.2066 1.58 1541  -0.763  1.0000
##  TCS08 - TCS43  -6.7016 1.63 1541  -4.099  0.0068
##  TCS08 - TCS44  -6.9171 1.58 1541  -4.387  0.0021
##  TCS08 - TCS45  -4.5666 1.58 1541  -2.896  0.2933
##  TCS08 - TCS46 -11.1461 1.58 1541  -7.069  <.0001
##  TCS08 - TCS47  -4.5529 1.58 1541  -2.887  0.2988
##  TCS08 - TCS48  -8.1945 1.58 1541  -5.197  <.0001
##  TCS08 - TCS49 -11.6185 1.58 1541  -7.368  <.0001
##  TCS10 - TCS11  -0.0453 1.39 1541  -0.033  1.0000
##  TCS10 - TCS12   1.6327 1.39 1541   1.172  0.9998
##  TCS10 - TCS20  -1.3205 1.39 1541  -0.948  1.0000
##  TCS10 - TCS40   6.1634 1.41 1541   4.380  0.0021
##  TCS10 - TCS41   7.5012 1.39 1541   5.383  <.0001
##  TCS10 - TCS42   7.2925 1.40 1541   5.217  <.0001
##  TCS10 - TCS43   1.7975 1.46 1541   1.231  0.9996
##  TCS10 - TCS44   1.5820 1.39 1541   1.135  0.9999
##  TCS10 - TCS45   3.9325 1.39 1541   2.822  0.3417
##  TCS10 - TCS46  -2.6470 1.39 1541  -1.899  0.9378
##  TCS10 - TCS47   3.9462 1.39 1541   2.832  0.3351
##  TCS10 - TCS48   0.3046 1.39 1541   0.219  1.0000
##  TCS10 - TCS49  -3.1194 1.39 1541  -2.238  0.7761
##  TCS11 - TCS12   1.6780 1.39 1541   1.204  0.9997
##  TCS11 - TCS20  -1.2752 1.39 1541  -0.915  1.0000
##  TCS11 - TCS40   6.2087 1.41 1541   4.413  0.0018
##  TCS11 - TCS41   7.5465 1.39 1541   5.415  <.0001
##  TCS11 - TCS42   7.3378 1.40 1541   5.249  <.0001
##  TCS11 - TCS43   1.8428 1.46 1541   1.262  0.9995
##  TCS11 - TCS44   1.6273 1.39 1541   1.168  0.9998
##  TCS11 - TCS45   3.9778 1.39 1541   2.854  0.3200
##  TCS11 - TCS46  -2.6017 1.39 1541  -1.867  0.9469
##  TCS11 - TCS47   3.9915 1.39 1541   2.864  0.3136
##  TCS11 - TCS48   0.3499 1.39 1541   0.251  1.0000
##  TCS11 - TCS49  -3.0741 1.39 1541  -2.206  0.7968
##  TCS12 - TCS20  -2.9532 1.39 1541  -2.119  0.8468
##  TCS12 - TCS40   4.5307 1.41 1541   3.220  0.1337
##  TCS12 - TCS41   5.8685 1.39 1541   4.211  0.0043
##  TCS12 - TCS42   5.6598 1.40 1541   4.049  0.0083
##  TCS12 - TCS43   0.1647 1.46 1541   0.113  1.0000
##  TCS12 - TCS44  -0.0507 1.39 1541  -0.036  1.0000
##  TCS12 - TCS45   2.2998 1.39 1541   1.650  0.9849
##  TCS12 - TCS46  -4.2798 1.39 1541  -3.071  0.1965
##  TCS12 - TCS47   2.3135 1.39 1541   1.660  0.9839
##  TCS12 - TCS48  -1.3281 1.39 1541  -0.953  1.0000
##  TCS12 - TCS49  -4.7521 1.39 1541  -3.410  0.0776
##  TCS20 - TCS40   7.4839 1.41 1541   5.319  <.0001
##  TCS20 - TCS41   8.8217 1.39 1541   6.330  <.0001
##  TCS20 - TCS42   8.6130 1.40 1541   6.161  <.0001
##  TCS20 - TCS43   3.1179 1.46 1541   2.136  0.8379
##  TCS20 - TCS44   2.9025 1.39 1541   2.083  0.8655
##  TCS20 - TCS45   5.2530 1.39 1541   3.769  0.0237
##  TCS20 - TCS46  -1.3265 1.39 1541  -0.952  1.0000
##  TCS20 - TCS47   5.2667 1.39 1541   3.779  0.0229
##  TCS20 - TCS48   1.6251 1.39 1541   1.166  0.9998
##  TCS20 - TCS49  -1.7989 1.39 1541  -1.291  0.9993
##  TCS40 - TCS41   1.3378 1.41 1541   0.951  1.0000
##  TCS40 - TCS42   1.1291 1.41 1541   0.800  1.0000
##  TCS40 - TCS43  -4.3659 1.47 1541  -2.964  0.2525
##  TCS40 - TCS44  -4.5814 1.41 1541  -3.256  0.1211
##  TCS40 - TCS45  -2.2309 1.41 1541  -1.586  0.9904
##  TCS40 - TCS46  -8.8104 1.41 1541  -6.262  <.0001
##  TCS40 - TCS47  -2.2172 1.41 1541  -1.576  0.9910
##  TCS40 - TCS48  -5.8588 1.41 1541  -4.164  0.0052
##  TCS40 - TCS49  -9.2828 1.41 1541  -6.598  <.0001
##  TCS41 - TCS42  -0.2087 1.40 1541  -0.149  1.0000
##  TCS41 - TCS43  -5.7038 1.46 1541  -3.907  0.0144
##  TCS41 - TCS44  -5.9193 1.39 1541  -4.247  0.0037
##  TCS41 - TCS45  -3.5688 1.39 1541  -2.561  0.5362
##  TCS41 - TCS46 -10.1483 1.39 1541  -7.282  <.0001
##  TCS41 - TCS47  -3.5551 1.39 1541  -2.551  0.5439
##  TCS41 - TCS48  -7.1967 1.39 1541  -5.164  <.0001
##  TCS41 - TCS49 -10.6206 1.39 1541  -7.621  <.0001
##  TCS42 - TCS43  -5.4950 1.46 1541  -3.753  0.0251
##  TCS42 - TCS44  -5.7105 1.40 1541  -4.085  0.0072
##  TCS42 - TCS45  -3.3600 1.40 1541  -2.404  0.6588
##  TCS42 - TCS46  -9.9395 1.40 1541  -7.110  <.0001
##  TCS42 - TCS47  -3.3463 1.40 1541  -2.394  0.6662
##  TCS42 - TCS48  -6.9879 1.40 1541  -4.999  0.0001
##  TCS42 - TCS49 -10.4119 1.40 1541  -7.448  <.0001
##  TCS43 - TCS44  -0.2155 1.46 1541  -0.148  1.0000
##  TCS43 - TCS45   2.1350 1.46 1541   1.462  0.9964
##  TCS43 - TCS46  -4.4445 1.46 1541  -3.044  0.2097
##  TCS43 - TCS47   2.1487 1.46 1541   1.472  0.9961
##  TCS43 - TCS48  -1.4929 1.46 1541  -1.023  1.0000
##  TCS43 - TCS49  -4.9168 1.46 1541  -3.368  0.0880
##  TCS44 - TCS45   2.3505 1.39 1541   1.687  0.9808
##  TCS44 - TCS46  -4.2290 1.39 1541  -3.035  0.2146
##  TCS44 - TCS47   2.3642 1.39 1541   1.696  0.9796
##  TCS44 - TCS48  -1.2774 1.39 1541  -0.917  1.0000
##  TCS44 - TCS49  -4.7014 1.39 1541  -3.374  0.0865
##  TCS45 - TCS46  -6.5795 1.39 1541  -4.721  0.0004
##  TCS45 - TCS47   0.0137 1.39 1541   0.010  1.0000
##  TCS45 - TCS48  -3.6279 1.39 1541  -2.603  0.5030
##  TCS45 - TCS49  -7.0519 1.39 1541  -5.060  0.0001
##  TCS46 - TCS47   6.5932 1.39 1541   4.731  0.0004
##  TCS46 - TCS48   2.9516 1.39 1541   2.118  0.8475
##  TCS46 - TCS49  -0.4723 1.39 1541  -0.339  1.0000
##  TCS47 - TCS48  -3.6416 1.39 1541  -2.613  0.4953
##  TCS47 - TCS49  -7.0656 1.39 1541  -5.070  0.0001
##  TCS48 - TCS49  -3.4240 1.39 1541  -2.457  0.6177
## 
## Results are averaged over the levels of: semana, bloque 
## P value adjustment: tukey method for comparing a family of 20 estimates
cld_gen <-multcomp::cld(contrast, alpha = 0.05, Letters = LETTERS, reversed=T)
cld_gen
##  gen   emmean    SE   df lower.CL upper.CL .group  
##  TCS49  15.31 0.985 1541    13.38    17.24  A      
##  TCS46  14.84 0.985 1541    12.90    16.77  A      
##  TCS20  13.51 0.985 1541    11.58    15.44  AB     
##  TCS11  12.24 0.985 1541    10.30    14.17  ABC    
##  TCS10  12.19 0.985 1541    10.26    14.12  ABC    
##  TCS48  11.89 0.985 1541     9.95    13.82  ABCD   
##  TCS03  11.42 0.992 1541     9.48    13.37  ABCD   
##  TCS44  10.61 0.985 1541     8.68    12.54  ABCDE  
##  TCS12  10.56 0.985 1541     8.62    12.49  ABCDE  
##  TCS43  10.39 1.077 1541     8.28    12.51  ABCDE  
##  TCS02   9.66 0.985 1541     7.73    11.59   BCDEF 
##  TCS05   8.36 0.985 1541     6.42    10.29    CDEFG
##  TCS45   8.26 0.985 1541     6.33    10.19    CDEFG
##  TCS47   8.24 0.985 1541     6.31    10.18    CDEFG
##  TCS01   7.48 1.039 1541     5.45     9.52    CDEFG
##  TCS04   7.05 0.985 1541     5.11     8.98     DEFG
##  TCS40   6.03 1.004 1541     4.06     8.00      EFG
##  TCS42   4.90 0.992 1541     2.95     6.84       FG
##  TCS41   4.69 0.985 1541     2.76     6.62        G
##  TCS08   3.69 1.231 1541     1.28     6.11        G
## 
## Results are averaged over the levels of: semana, bloque 
## Confidence level used: 0.95 
## P value adjustment: tukey method for comparing a family of 20 estimates 
## significance level used: alpha = 0.05 
## NOTE: Compact letter displays can be misleading
##       because they show NON-findings rather than findings.
##       Consider using 'pairs()', 'pwpp()', or 'pwpm()' instead.
detach(datos5)