library(readxl)
#install.packages("metan")
library(metan)
## Registered S3 method overwritten by 'GGally':
##   method from   
##   +.gg   ggplot2
## |=========================================================|
## | Multi-Environment Trial Analysis (metan) v1.18.0        |
## | Author: Tiago Olivoto                                   |
## | Type 'citation('metan')' to know how to cite metan      |
## | Type 'vignette('metan_start')' for a short tutorial     |
## | Visit 'https://bit.ly/pkgmetan' for a complete tutorial |
## |=========================================================|
data = read_excel("C:/Users/Cristian/Desktop/copias enayos/Datos Parcial Fitomejoramiento/Datos.xlsx")

View(data)


DF = as.data.frame(data)

head(DF)
##   Genotipos Rendimiento Localidad Repeticion
## 1      22-1         3.5  Monteria          1
## 2   Testigo         3.7  Monteria          1
## 3     22-14         3.6  Monteria          1
## 4     22-31         7.2  Monteria          1
## 5     22-53         6.8  Monteria          1
## 6     22-70         6.0  Monteria          1
str(DF)
## 'data.frame':    135 obs. of  4 variables:
##  $ Genotipos  : chr  "22-1" "Testigo" "22-14" "22-31" ...
##  $ Rendimiento: num  3.5 3.7 3.6 7.2 6.8 6 6.1 6.8 6 3.5 ...
##  $ Localidad  : chr  "Monteria" "Monteria" "Monteria" "Monteria" ...
##  $ Repeticion : num  1 1 1 1 1 1 1 1 1 2 ...
#Transformacion a FACTORES
DF$Loc = as.factor(DF$Localidad)
DF$Gen = as.factor(DF$Genotipos)
DF$Rep = as.factor(DF$Repeticion)

View(DF)
#ANOVA
anova = aov(Rendimiento~Loc + Gen + Loc*Gen, data = DF)
summary(anova)
##             Df Sum Sq Mean Sq F value   Pr(>F)    
## Loc          4  56.80  14.199   47.86  < 2e-16 ***
## Gen          8 196.57  24.571   82.82  < 2e-16 ***
## Loc:Gen     32  30.29   0.946    3.19 8.56e-06 ***
## Residuals   90  26.70   0.297                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Within-environment analysis of variance

#Performs a within-environment analysis of variance in randomized complete block or alpha-lattice designs and returns values such as Mean Squares, p-values, coefficient of variation, heritability, and accuracy of selection.

#For analysis in randomized complete block design:

#MEAN: The grand mean.

#⁠DFG, DFB, and DFE⁠: The degree of freedom for genotype blocks, and error, respectively.

#⁠MSG, MSB, and MSE⁠: The mean squares for genotype blocks, and error, respectively.

#⁠FCG and FCB⁠: The F-calculated for genotype and blocks, respectively.

#⁠PFG and PFB⁠: The P-values for genotype and blocks, respectively.

#CV: coefficient of variation.

#h2: broad-sense heritability.

#AS: accuracy of selection (square root of h2)

indiv_RC = anova_ind(.data = DF, env = Loc, gen = Gen, rep = Rep, resp = Rendimiento)
## Evaluating trait Rendimiento |===================================| 100% 00:00:00 
indiv_RC
## Variable Rendimiento 
## ---------------------------------------------------------------------------
## Within-environment ANOVA results
## ---------------------------------------------------------------------------
## # A tibble: 5 × 15
##   ENV    MEAN   DFG   MSG    FCG      PFG   DFB    MSB   FCB    PFB   DFE    MSE
##   <chr> <dbl> <int> <dbl>  <dbl>    <dbl> <int>  <dbl> <dbl>  <dbl> <int>  <dbl>
## 1 Agua…  5.18     8  5.67 381.   8.92e-17     2 0.0211 1.42  0.270     16 0.0149
## 2 Espi…  5.07     8  6.48   4.75 3.92e- 3     2 1.85   1.36  0.286     16 1.36  
## 3 Mont…  5.54     8  6.62 454.   2.25e-17     2 0.0133 0.914 0.421     16 0.0146
## 4 Neiva  3.7      8  1.73  99    3.70e-12     2 0.0400 2.29  0.134     16 0.0175
## 5 Yopal  5.3      8  7.86 682.   8.76e-19     2 0.0578 5.01  0.0204    16 0.0115
## # ℹ 3 more variables: CV <dbl>, h2 <dbl>, AS <dbl>
## ---------------------------------------------------------------------------
#Joint analysis of variance
#Performs a joint analysis of variance to check for the presence of genotype-vs-environment interactions using both randomized complete block and alpha-lattice designs.
joint_RC = anova_joint(.data = DF, env = Loc, gen = Gen, rep = Rep, resp = Rendimiento)
## variable Rendimiento 
## ---------------------------------------------------------------------------
## Joint ANOVA table
## ---------------------------------------------------------------------------
##     Source     Df Sum Sq Mean Sq F value   Pr(>F)
##        ENV   4.00  56.80  14.199   49.95 5.25e-21
##   REP(ENV)  10.00   3.96   0.396    1.39 1.99e-01
##        GEN   8.00 196.57  24.571   86.44 3.82e-36
##    GEN:ENV  32.00  30.29   0.946    3.33 7.13e-06
##  Residuals  80.00  22.74   0.284      NA       NA
##      CV(%)  10.75     NA      NA      NA       NA
##  MSR+/MSR- 118.22     NA      NA      NA       NA
##     OVmean   4.96     NA      NA      NA       NA
## ---------------------------------------------------------------------------
## 
## All variables with significant (p < 0.05) genotype-vs-environment interaction
## Done!
#Normas de reaccion - Analisis Preliminar
with(DF, interaction.plot(Loc, Gen, Rendimiento))

with(DF, interaction.plot(Gen, Loc, Rendimiento))

library(phia)
## Loading required package: car
## Loading required package: carData
#Test Contrasts of Factor Interactions

#Efectos simples de Localidad

#Calculates and tests different types of contrasts for factor interactions, in linear, generalized and mixed linear models: simple main effects, interaction contrasts, residual effects, and others.


testInteractions(anova, pairwise = "Loc", fixed = "Gen", adjust = "none")
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in rbind(deparse.level, ...): number of columns of result, 6, is not a
## multiple of vector length 5 of arg 2
## F Test: 
## P-value adjustment method: none
##                              Value     SE   Df Sum of Sq       F    Pr(>F)    
##  Aguachica-Espinal :    22-1   0.1  0.445  1.0     0.015  0.0506 0.8225983    
## Aguachica-Monteria :    22-1   0.2  0.445  1.0     0.060  0.2022 0.6539943    
##    Aguachica-Neiva :    22-1   0.7  0.445  1.0     0.735  2.4775 0.1189926    
##    Aguachica-Yopal :    22-1   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   Espinal-Monteria :    22-1   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##      Espinal-Neiva :    22-1   0.6  0.445  1.0     0.540  1.8202 0.1806710    
##      Espinal-Yopal :    22-1   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##     Monteria-Neiva :    22-1   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##     Monteria-Yopal :    22-1  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##        Neiva-Yopal :    22-1  -0.6  0.445  1.0     0.540  1.8202 0.1806710    
##  Aguachica-Espinal :   22-14   0.0  0.445  1.0     0.000  0.0000 1.0000000    
## Aguachica-Monteria :   22-14  -0.5  0.445  1.0     0.375  1.2640 0.2638766    
##    Aguachica-Neiva :   22-14   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##    Aguachica-Yopal :   22-14  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   Espinal-Monteria :   22-14  -0.5  0.445  1.0     0.375  1.2640 0.2638766    
##      Espinal-Neiva :   22-14   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##      Espinal-Yopal :   22-14  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##     Monteria-Neiva :   22-14   1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##     Monteria-Yopal :   22-14   0.4  0.445  1.0     0.240  0.8090 0.3708196    
##        Neiva-Yopal :   22-14  -0.6  0.445  1.0     0.540  1.8202 0.1806710    
##  Aguachica-Espinal :   22-23   0.2  0.445  1.0     0.060  0.2022 0.6539943    
## Aguachica-Monteria :   22-23   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##    Aguachica-Neiva :   22-23   2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##    Aguachica-Yopal :   22-23   0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   Espinal-Monteria :   22-23  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
##      Espinal-Neiva :   22-23   2.0  0.445  1.0     6.000 20.2247 2.047e-05 ***
##      Espinal-Yopal :   22-23   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##     Monteria-Neiva :   22-23   2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##     Monteria-Yopal :   22-23   0.2  0.445  1.0     0.060  0.2022 0.6539943    
##        Neiva-Yopal :   22-23  -2.0  0.445  1.0     6.000 20.2247 2.047e-05 ***
##  Aguachica-Espinal :   22-31  -0.4  0.445  1.0     0.240  0.8090 0.3708196    
## Aguachica-Monteria :   22-31  -0.4  0.445  1.0     0.240  0.8090 0.3708196    
##    Aguachica-Neiva :   22-31   2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##    Aguachica-Yopal :   22-31  -0.4  0.445  1.0     0.240  0.8090 0.3708196    
##   Espinal-Monteria :   22-31   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##      Espinal-Neiva :   22-31   2.6  0.445  1.0    10.140 34.1798 7.956e-08 ***
##      Espinal-Yopal :   22-31   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##     Monteria-Neiva :   22-31   2.6  0.445  1.0    10.140 34.1798 7.956e-08 ***
##     Monteria-Yopal :   22-31   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##        Neiva-Yopal :   22-31  -2.6  0.445  1.0    10.140 34.1798 7.956e-08 ***
##  Aguachica-Espinal :   22-53   2.0  0.445  1.0     6.000 20.2247 2.047e-05 ***
## Aguachica-Monteria :   22-53  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##    Aguachica-Neiva :   22-53   2.5  0.445  1.0     9.375 31.6011 2.101e-07 ***
##    Aguachica-Yopal :   22-53  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   Espinal-Monteria :   22-53  -2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##      Espinal-Neiva :   22-53   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##      Espinal-Yopal :   22-53  -2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##     Monteria-Neiva :   22-53   2.6  0.445  1.0    10.140 34.1798 7.956e-08 ***
##     Monteria-Yopal :   22-53  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##        Neiva-Yopal :   22-53  -2.7  0.445  1.0    10.935 36.8596 2.968e-08 ***
##  Aguachica-Espinal :   22-55  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
## Aguachica-Monteria :   22-55  -0.7  0.445  1.0     0.735  2.4775 0.1189926    
##    Aguachica-Neiva :   22-55   1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##    Aguachica-Yopal :   22-55  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   Espinal-Monteria :   22-55  -0.6  0.445  1.0     0.540  1.8202 0.1806710    
##      Espinal-Neiva :   22-55   1.8  0.445  1.0     4.860 16.3820 0.0001093 ***
##      Espinal-Yopal :   22-55  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##     Monteria-Neiva :   22-55   2.4  0.445  1.0     8.640 29.1236 5.457e-07 ***
##     Monteria-Yopal :   22-55   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##        Neiva-Yopal :   22-55  -1.9  0.445  1.0     5.415 18.2528 4.792e-05 ***
##  Aguachica-Espinal :   22-70  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
## Aguachica-Monteria :   22-70   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##    Aguachica-Neiva :   22-70   1.4  0.445  1.0     2.940  9.9101 0.0022300 ** 
##    Aguachica-Yopal :   22-70  -1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##   Espinal-Monteria :   22-70   1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##      Espinal-Neiva :   22-70   2.3  0.445  1.0     7.935 26.7472 1.392e-06 ***
##      Espinal-Yopal :   22-70  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##     Monteria-Neiva :   22-70   1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
##     Monteria-Yopal :   22-70  -1.1  0.445  1.0     1.815  6.1180 0.0152583 *  
##        Neiva-Yopal :   22-70  -2.4  0.445  1.0     8.640 29.1236 5.457e-07 ***
##  Aguachica-Espinal :   22-99  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
## Aguachica-Monteria :   22-99  -2.0  0.445  1.0     6.000 20.2247 2.047e-05 ***
##    Aguachica-Neiva :   22-99   1.1  0.445  1.0     1.815  6.1180 0.0152583 *  
##    Aguachica-Yopal :   22-99   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   Espinal-Monteria :   22-99  -1.8  0.445  1.0     4.860 16.3820 0.0001093 ***
##      Espinal-Neiva :   22-99   1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
##      Espinal-Yopal :   22-99   0.3  0.445  1.0     0.135  0.4551 0.5016728    
##     Monteria-Neiva :   22-99   3.1  0.445  1.0    14.415 48.5899 5.080e-10 ***
##     Monteria-Yopal :   22-99   2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##        Neiva-Yopal :   22-99  -1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##  Aguachica-Espinal : Testigo   0.3  0.445  1.0     0.135  0.4551 0.5016728    
## Aguachica-Monteria : Testigo   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##    Aguachica-Neiva : Testigo   1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##    Aguachica-Yopal : Testigo   0.4  0.445  1.0     0.240  0.8090 0.3708196    
##   Espinal-Monteria : Testigo  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
##      Espinal-Neiva : Testigo   0.7  0.445  1.0     0.735  2.4775 0.1189926    
##      Espinal-Yopal : Testigo   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##     Monteria-Neiva : Testigo   0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##     Monteria-Yopal : Testigo   0.3  0.445  1.0     0.135  0.4551 0.5016728    
##        Neiva-Yopal : Testigo  -0.6  0.445  1.0     0.540  1.8202 0.1806710    
## Residuals                          90.000 26.7                                
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Efectos simples de Genotipo

testInteractions(anova, pairwise = "Gen", fixed = "Loc", adjust = "none")
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in printHypothesis(L, rhs, names(b)): one or more coefficients in the hypothesis include
##      arithmetic operators in their names;
##   the printed representation of the hypothesis will be omitted
## Warning in rbind(deparse.level, ...): number of columns of result, 6, is not a
## multiple of vector length 5 of arg 2
## F Test: 
## P-value adjustment method: none
##                           Value     SE   Df Sum of Sq       F    Pr(>F)    
##    22-1-22-14 : Aguachica   0.7  0.445  1.0     0.735  2.4775 0.1189926    
##    22-1-22-23 : Aguachica  -2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##    22-1-22-31 : Aguachica  -3.1  0.445  1.0    14.415 48.5899 5.080e-10 ***
##    22-1-22-53 : Aguachica  -3.0  0.445  1.0    13.500 45.5056 1.427e-09 ***
##    22-1-22-55 : Aguachica  -1.4  0.445  1.0     2.940  9.9101 0.0022300 ** 
##    22-1-22-70 : Aguachica  -2.3  0.445  1.0     7.935 26.7472 1.392e-06 ***
##    22-1-22-99 : Aguachica  -1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##  22-1-Testigo : Aguachica  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   22-14-22-23 : Aguachica  -2.9  0.445  1.0    12.615 42.5225 3.971e-09 ***
##   22-14-22-31 : Aguachica  -3.8  0.445  1.0    21.660 73.0112 3.044e-13 ***
##   22-14-22-53 : Aguachica  -3.7  0.445  1.0    20.535 69.2191 8.906e-13 ***
##   22-14-22-55 : Aguachica  -2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##   22-14-22-70 : Aguachica  -3.0  0.445  1.0    13.500 45.5056 1.427e-09 ***
##   22-14-22-99 : Aguachica  -1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
## 22-14-Testigo : Aguachica  -0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
##   22-23-22-31 : Aguachica  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-23-22-53 : Aguachica  -0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
##   22-23-22-55 : Aguachica   0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
##   22-23-22-70 : Aguachica  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   22-23-22-99 : Aguachica   1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
## 22-23-Testigo : Aguachica   2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##   22-31-22-53 : Aguachica   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   22-31-22-55 : Aguachica   1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##   22-31-22-70 : Aguachica   0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
##   22-31-22-99 : Aguachica   2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
## 22-31-Testigo : Aguachica   3.0  0.445  1.0    13.500 45.5056 1.427e-09 ***
##   22-53-22-55 : Aguachica   1.6  0.445  1.0     3.840 12.9438 0.0005242 ***
##   22-53-22-70 : Aguachica   0.7  0.445  1.0     0.735  2.4775 0.1189926    
##   22-53-22-99 : Aguachica   2.0  0.445  1.0     6.000 20.2247 2.047e-05 ***
## 22-53-Testigo : Aguachica   2.9  0.445  1.0    12.615 42.5225 3.971e-09 ***
##   22-55-22-70 : Aguachica  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-55-22-99 : Aguachica   0.4  0.445  1.0     0.240  0.8090 0.3708196    
## 22-55-Testigo : Aguachica   1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
##   22-70-22-99 : Aguachica   1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
## 22-70-Testigo : Aguachica   2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
## 22-99-Testigo : Aguachica   0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##    22-1-22-14 :   Espinal   0.6  0.445  1.0     0.540  1.8202 0.1806710    
##    22-1-22-23 :   Espinal  -2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##    22-1-22-31 :   Espinal  -3.6  0.445  1.0    19.440 65.5281 2.597e-12 ***
##    22-1-22-53 :   Espinal  -1.1  0.445  1.0     1.815  6.1180 0.0152583 *  
##    22-1-22-55 :   Espinal  -1.6  0.445  1.0     3.840 12.9438 0.0005242 ***
##    22-1-22-70 :   Espinal  -3.3  0.445  1.0    16.335 55.0618 6.275e-11 ***
##    22-1-22-99 :   Espinal  -1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
##  22-1-Testigo :   Espinal   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   22-14-22-23 :   Espinal  -2.7  0.445  1.0    10.935 36.8596 2.968e-08 ***
##   22-14-22-31 :   Espinal  -4.2  0.445  1.0    26.460 89.1910 4.094e-15 ***
##   22-14-22-53 :   Espinal  -1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##   22-14-22-55 :   Espinal  -2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##   22-14-22-70 :   Espinal  -3.9  0.445  1.0    22.815 76.9045 1.038e-13 ***
##   22-14-22-99 :   Espinal  -1.9  0.445  1.0     5.415 18.2528 4.792e-05 ***
## 22-14-Testigo :   Espinal  -0.5  0.445  1.0     0.375  1.2640 0.2638766    
##   22-23-22-31 :   Espinal  -1.5  0.445  1.0     3.375 11.3764 0.0010983 ** 
##   22-23-22-53 :   Espinal   1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##   22-23-22-55 :   Espinal   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##   22-23-22-70 :   Espinal  -1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
##   22-23-22-99 :   Espinal   0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
## 22-23-Testigo :   Espinal   2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##   22-31-22-53 :   Espinal   2.5  0.445  1.0     9.375 31.6011 2.101e-07 ***
##   22-31-22-55 :   Espinal   2.0  0.445  1.0     6.000 20.2247 2.047e-05 ***
##   22-31-22-70 :   Espinal   0.3  0.445  1.0     0.135  0.4551 0.5016728    
##   22-31-22-99 :   Espinal   2.3  0.445  1.0     7.935 26.7472 1.392e-06 ***
## 22-31-Testigo :   Espinal   3.7  0.445  1.0    20.535 69.2191 8.906e-13 ***
##   22-53-22-55 :   Espinal  -0.5  0.445  1.0     0.375  1.2640 0.2638766    
##   22-53-22-70 :   Espinal  -2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##   22-53-22-99 :   Espinal  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
## 22-53-Testigo :   Espinal   1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
##   22-55-22-70 :   Espinal  -1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##   22-55-22-99 :   Espinal   0.3  0.445  1.0     0.135  0.4551 0.5016728    
## 22-55-Testigo :   Espinal   1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##   22-70-22-99 :   Espinal   2.0  0.445  1.0     6.000 20.2247 2.047e-05 ***
## 22-70-Testigo :   Espinal   3.4  0.445  1.0    17.340 58.4494 2.182e-11 ***
## 22-99-Testigo :   Espinal   1.4  0.445  1.0     2.940  9.9101 0.0022300 ** 
##    22-1-22-14 :  Monteria   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##    22-1-22-23 :  Monteria  -2.4  0.445  1.0     8.640 29.1236 5.457e-07 ***
##    22-1-22-31 :  Monteria  -3.7  0.445  1.0    20.535 69.2191 8.906e-13 ***
##    22-1-22-53 :  Monteria  -3.3  0.445  1.0    16.335 55.0618 6.275e-11 ***
##    22-1-22-55 :  Monteria  -2.3  0.445  1.0     7.935 26.7472 1.392e-06 ***
##    22-1-22-70 :  Monteria  -2.4  0.445  1.0     8.640 29.1236 5.457e-07 ***
##    22-1-22-99 :  Monteria  -3.2  0.445  1.0    15.360 51.7753 1.792e-10 ***
##  22-1-Testigo :  Monteria  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   22-14-22-23 :  Monteria  -2.4  0.445  1.0     8.640 29.1236 5.457e-07 ***
##   22-14-22-31 :  Monteria  -3.7  0.445  1.0    20.535 69.2191 8.906e-13 ***
##   22-14-22-53 :  Monteria  -3.3  0.445  1.0    16.335 55.0618 6.275e-11 ***
##   22-14-22-55 :  Monteria  -2.3  0.445  1.0     7.935 26.7472 1.392e-06 ***
##   22-14-22-70 :  Monteria  -2.4  0.445  1.0     8.640 29.1236 5.457e-07 ***
##   22-14-22-99 :  Monteria  -3.2  0.445  1.0    15.360 51.7753 1.792e-10 ***
## 22-14-Testigo :  Monteria  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   22-23-22-31 :  Monteria  -1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
##   22-23-22-53 :  Monteria  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-23-22-55 :  Monteria   0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   22-23-22-70 :  Monteria   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##   22-23-22-99 :  Monteria  -0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
## 22-23-Testigo :  Monteria   2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##   22-31-22-53 :  Monteria   0.4  0.445  1.0     0.240  0.8090 0.3708196    
##   22-31-22-55 :  Monteria   1.4  0.445  1.0     2.940  9.9101 0.0022300 ** 
##   22-31-22-70 :  Monteria   1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
##   22-31-22-99 :  Monteria   0.5  0.445  1.0     0.375  1.2640 0.2638766    
## 22-31-Testigo :  Monteria   3.5  0.445  1.0    18.375 61.9382 7.547e-12 ***
##   22-53-22-55 :  Monteria   1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##   22-53-22-70 :  Monteria   0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-53-22-99 :  Monteria   0.1  0.445  1.0     0.015  0.0506 0.8225983    
## 22-53-Testigo :  Monteria   3.1  0.445  1.0    14.415 48.5899 5.080e-10 ***
##   22-55-22-70 :  Monteria  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   22-55-22-99 :  Monteria  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
## 22-55-Testigo :  Monteria   2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##   22-70-22-99 :  Monteria  -0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
## 22-70-Testigo :  Monteria   2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
## 22-99-Testigo :  Monteria   3.0  0.445  1.0    13.500 45.5056 1.427e-09 ***
##    22-1-22-14 :     Neiva   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##    22-1-22-23 :     Neiva  -0.7  0.445  1.0     0.735  2.4775 0.1189926    
##    22-1-22-31 :     Neiva  -1.6  0.445  1.0     3.840 12.9438 0.0005242 ***
##    22-1-22-53 :     Neiva  -1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
##    22-1-22-55 :     Neiva  -0.4  0.445  1.0     0.240  0.8090 0.3708196    
##    22-1-22-70 :     Neiva  -1.6  0.445  1.0     3.840 12.9438 0.0005242 ***
##    22-1-22-99 :     Neiva  -0.6  0.445  1.0     0.540  1.8202 0.1806710    
##  22-1-Testigo :     Neiva   0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   22-14-22-23 :     Neiva  -1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
##   22-14-22-31 :     Neiva  -2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##   22-14-22-53 :     Neiva  -1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##   22-14-22-55 :     Neiva  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-14-22-70 :     Neiva  -2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##   22-14-22-99 :     Neiva  -1.1  0.445  1.0     1.815  6.1180 0.0152583 *  
## 22-14-Testigo :     Neiva  -0.3  0.445  1.0     0.135  0.4551 0.5016728    
##   22-23-22-31 :     Neiva  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-23-22-53 :     Neiva  -0.5  0.445  1.0     0.375  1.2640 0.2638766    
##   22-23-22-55 :     Neiva   0.3  0.445  1.0     0.135  0.4551 0.5016728    
##   22-23-22-70 :     Neiva  -0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-23-22-99 :     Neiva   0.1  0.445  1.0     0.015  0.0506 0.8225983    
## 22-23-Testigo :     Neiva   0.9  0.445  1.0     1.215  4.0955 0.0459630 *  
##   22-31-22-53 :     Neiva   0.4  0.445  1.0     0.240  0.8090 0.3708196    
##   22-31-22-55 :     Neiva   1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
##   22-31-22-70 :     Neiva   0.0  0.445  1.0     0.000  0.0000 1.0000000    
##   22-31-22-99 :     Neiva   1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
## 22-31-Testigo :     Neiva   1.8  0.445  1.0     4.860 16.3820 0.0001093 ***
##   22-53-22-55 :     Neiva   0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
##   22-53-22-70 :     Neiva  -0.4  0.445  1.0     0.240  0.8090 0.3708196    
##   22-53-22-99 :     Neiva   0.6  0.445  1.0     0.540  1.8202 0.1806710    
## 22-53-Testigo :     Neiva   1.4  0.445  1.0     2.940  9.9101 0.0022300 ** 
##   22-55-22-70 :     Neiva  -1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
##   22-55-22-99 :     Neiva  -0.2  0.445  1.0     0.060  0.2022 0.6539943    
## 22-55-Testigo :     Neiva   0.6  0.445  1.0     0.540  1.8202 0.1806710    
##   22-70-22-99 :     Neiva   1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
## 22-70-Testigo :     Neiva   1.8  0.445  1.0     4.860 16.3820 0.0001093 ***
## 22-99-Testigo :     Neiva   0.8  0.445  1.0     0.960  3.2360 0.0753903 .  
##    22-1-22-14 :     Yopal   0.5  0.445  1.0     0.375  1.2640 0.2638766    
##    22-1-22-23 :     Yopal  -2.1  0.445  1.0     6.615 22.2978 8.542e-06 ***
##    22-1-22-31 :     Yopal  -3.6  0.445  1.0    19.440 65.5281 2.597e-12 ***
##    22-1-22-53 :     Yopal  -3.3  0.445  1.0    16.335 55.0618 6.275e-11 ***
##    22-1-22-55 :     Yopal  -1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##    22-1-22-70 :     Yopal  -3.4  0.445  1.0    17.340 58.4494 2.182e-11 ***
##    22-1-22-99 :     Yopal  -1.0  0.445  1.0     1.500  5.0562 0.0269794 *  
##  22-1-Testigo :     Yopal   0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   22-14-22-23 :     Yopal  -2.6  0.445  1.0    10.140 34.1798 7.956e-08 ***
##   22-14-22-31 :     Yopal  -4.1  0.445  1.0    25.215 84.9944 1.203e-14 ***
##   22-14-22-53 :     Yopal  -3.8  0.445  1.0    21.660 73.0112 3.044e-13 ***
##   22-14-22-55 :     Yopal  -2.2  0.445  1.0     7.260 24.4719 3.485e-06 ***
##   22-14-22-70 :     Yopal  -3.9  0.445  1.0    22.815 76.9045 1.038e-13 ***
##   22-14-22-99 :     Yopal  -1.5  0.445  1.0     3.375 11.3764 0.0010983 ** 
## 22-14-Testigo :     Yopal  -0.3  0.445  1.0     0.135  0.4551 0.5016728    
##   22-23-22-31 :     Yopal  -1.5  0.445  1.0     3.375 11.3764 0.0010983 ** 
##   22-23-22-53 :     Yopal  -1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
##   22-23-22-55 :     Yopal   0.4  0.445  1.0     0.240  0.8090 0.3708196    
##   22-23-22-70 :     Yopal  -1.3  0.445  1.0     2.535  8.5449 0.0043818 ** 
##   22-23-22-99 :     Yopal   1.1  0.445  1.0     1.815  6.1180 0.0152583 *  
## 22-23-Testigo :     Yopal   2.3  0.445  1.0     7.935 26.7472 1.392e-06 ***
##   22-31-22-53 :     Yopal   0.3  0.445  1.0     0.135  0.4551 0.5016728    
##   22-31-22-55 :     Yopal   1.9  0.445  1.0     5.415 18.2528 4.792e-05 ***
##   22-31-22-70 :     Yopal   0.2  0.445  1.0     0.060  0.2022 0.6539943    
##   22-31-22-99 :     Yopal   2.6  0.445  1.0    10.140 34.1798 7.956e-08 ***
## 22-31-Testigo :     Yopal   3.8  0.445  1.0    21.660 73.0112 3.044e-13 ***
##   22-53-22-55 :     Yopal   1.6  0.445  1.0     3.840 12.9438 0.0005242 ***
##   22-53-22-70 :     Yopal  -0.1  0.445  1.0     0.015  0.0506 0.8225983    
##   22-53-22-99 :     Yopal   2.3  0.445  1.0     7.935 26.7472 1.392e-06 ***
## 22-53-Testigo :     Yopal   3.5  0.445  1.0    18.375 61.9382 7.547e-12 ***
##   22-55-22-70 :     Yopal  -1.7  0.445  1.0     4.335 14.6124 0.0002428 ***
##   22-55-22-99 :     Yopal   0.7  0.445  1.0     0.735  2.4775 0.1189926    
## 22-55-Testigo :     Yopal   1.9  0.445  1.0     5.415 18.2528 4.792e-05 ***
##   22-70-22-99 :     Yopal   2.4  0.445  1.0     8.640 29.1236 5.457e-07 ***
## 22-70-Testigo :     Yopal   3.6  0.445  1.0    19.440 65.5281 2.597e-12 ***
## 22-99-Testigo :     Yopal   1.2  0.445  1.0     2.160  7.2809 0.0083220 ** 
## Residuals                       90.000 26.7                                
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#Regresion Metodo Eberhart & Russell 

#Regression-based stability analysis using the Eberhart and Russell (1966) model.

Regresion = ge_reg(.data = DF, env = Loc, gen = Gen, rep = Rep, resp = Rendimiento)
## Evaluating trait Rendimiento |===================================| 100% 00:00:00 
Regresion
## Variable Rendimiento 
## ---------------------------------------------------------------------------
## Joint-regression Analysis of variance
## ---------------------------------------------------------------------------
## # A tibble: 16 × 6
##    SV                       Df `Sum Sq` `Mean Sq` `F value`  `Pr(>F)`
##    <chr>                 <dbl>    <dbl>     <dbl>     <dbl>     <dbl>
##  1 "Total"                  44  284.       6.45      NA     NA       
##  2 "GEN"                     8  197.      24.6       33.3    4.92e-12
##  3 "ENV + (GEN x ENV)"      36   87.1      2.42      NA     NA       
##  4 "ENV (linear)"            1   56.8     56.8       NA     NA       
##  5 " GEN x ENV (linear)"     8   10.3      1.29       1.75   1.32e- 1
##  6 "Pooled deviation"       27   19.9      0.739     NA     NA       
##  7 "22-1"                    3    0.189    0.0629     0.221  8.81e- 1
##  8 "22-14"                   3    0.219    0.0731     0.257  8.56e- 1
##  9 "22-23"                   3    0.485    0.162      0.568  6.38e- 1
## 10 "22-31"                   3    0.918    0.306      1.08   3.64e- 1
## 11 "22-53"                   3    7.02     2.34       8.23   7.74e- 5
## 12 "22-55"                   3    0.145    0.0482     0.170  9.17e- 1
## 13 "22-70"                   3    4.70     1.57       5.51   1.72e- 3
## 14 "22-99"                   3    5.95     1.98       6.97   3.16e- 4
## 15 "Testigo"                 3    0.324    0.108      0.380  7.67e- 1
## 16 "Pooled error"           80   22.7      0.284     NA     NA       
## ---------------------------------------------------------------------------
## Regression parameters
## ---------------------------------------------------------------------------
## # A tibble: 9 × 10
##   GEN      b0    b1 `t(b1=1)`[,1] pval_t[,1]     s2di `F(s2di=0)`  pval_f   RMSE
##   <chr> <dbl> <dbl>         <dbl>      <dbl>    <dbl>       <dbl>   <dbl>  <dbl>
## 1 22-1   3.58 0.341       -3.10      0.00264 -0.0738        0.221 8.81e-1 0.112 
## 2 22-14  3.12 0.455       -2.57      0.0120  -0.0704        0.257 8.56e-1 0.121 
## 3 22-23  5.48 1.27         1.28      0.203   -0.0409        0.568 6.38e-1 0.180 
## 4 22-31  6.7  1.51         2.42      0.0180   0.00725       1.08  3.64e-1 0.247 
## 5 22-53  5.96 1.45         2.13      0.0364   0.685         8.23  7.74e-5 0.684 
## 6 22-55  5.06 1.25         1.18      0.243   -0.0787        0.170 9.17e-1 0.0982
## 7 22-70  6.18 1.02         0.0871    0.931    0.427         5.51  1.72e-3 0.560 
## 8 22-99  5    1.21         0.988     0.326    0.566         6.97  3.16e-4 0.630 
## 9 Test…  3.54 0.489       -2.41      0.0184  -0.0587        0.380 7.67e-1 0.147 
## # ℹ 1 more variable: R2 <dbl>
## ---------------------------------------------------------------------------
## Variance of b0: 0.01895 
## Variance of b1: 0.04504278
plot(Regresion)

#Modelo AMMI

Ammi = performs_ammi(DF, env = Loc, gen = Gen, rep = Rep, resp = Rendimiento, verbose = T) #2-axis PC1 y PC2 = >72%
## variable Rendimiento 
## ---------------------------------------------------------------------------
## AMMI analysis table
## ---------------------------------------------------------------------------
##     Source  Df Sum Sq Mean Sq F value   Pr(>F) Proportion Accumulated
##        ENV   4  56.80  14.199   35.86 6.67e-06         NA          NA
##   REP(ENV)  10   3.96   0.396    1.39 1.99e-01         NA          NA
##        GEN   8 196.57  24.571   86.44 3.82e-36         NA          NA
##    GEN:ENV  32  30.29   0.946    3.33 7.13e-06         NA          NA
##        PC1  11  11.82   1.074    3.78 2.00e-04       39.0        39.0
##        PC2   9  10.00   1.111    3.91 4.00e-04       33.0        72.0
##        PC3   7   7.42   1.060    3.73 1.50e-03       24.5        96.5
##        PC4   5   1.05   0.210    0.74 5.96e-01        3.5       100.0
##  Residuals  80  22.74   0.284      NA       NA         NA          NA
##      Total 166 340.64   2.052      NA       NA         NA          NA
## ---------------------------------------------------------------------------
## 
## All variables with significant (p < 0.05) genotype-vs-environment interaction
## Done!
#Estimacion de la variable respuesta basada en ejes simplificativos del IPCA 

#predict is a generic function for predictions from the results of various model fitting functions. The function invokes particular methods which depend on the class of the first argument.
Predicted_RC = predict(Ammi, naxis = 2) 
View(Predicted_RC)
#BiPlot Type 1
a = plot_scores(Ammi)
b = plot_scores(Ammi,
                col.gen = "black",
                col.env = "gray70",
                col.segm.env = "gray70",
                axis.expand = 1.5,
                plot_theme = theme_metan(grid = "both"))

arrange_ggplot(a, b, tag_levels = "a")

#Biplot Type 2

c = plot_scores(Ammi, type = 2,
                polygon = T,
              col.segm.env = "transparent",
                axis.expand = 1.5,
                plot_theme = theme_metan_minimal())
#Exactitud de prediccion de los modelos 

Regresion_C = cor(Predicted_RC$Y, Predicted_RC$Ypred)*100
Ammi_C = cor(Predicted_RC$Y, Predicted_RC$YpredAMMI)*100