M1_M2
Juan Sebastian Herrera Reales
2023-06-05
Muestreo 1 analisis de datos
datos_m <- read_excel("C:/Users/JuanSebH2/Downloads/Universidad/Fisiologia/FisioVeg/DocumentoFinal_muestreos.xlsx",
sheet = "M1")
datos_m=data.frame(datos_m)
datos_m1=datos_m
datos_m1$BLOQUES=as.character(datos_m1$BLOQUES)
collapsibleTree(datos_m1, hierarchy = c("MUESTREO", "BLOQUES", "TRATAMIENTO"), collapsed = F)datos_m1## MUESTREO TRATAMIENTO BLOQUES TEMPERATURA N_HOJAS AREA_DOSEL LONGITUD_P_AEREA
## 1 1 (-DH-Si ) 1 18.40 6 103.67 10.1
## 2 1 (-DH-Si ) 1 18.90 6 74.88 12.6
## 3 1 (-DH-Si ) 2 20.36 6 181.53 14.5
## 4 1 (-DH-Si ) 2 19.00 5 131.15 13.2
## 5 1 (+DH-Si ) 1 19.70 5 69.70 10.4
## 6 1 (+DH-Si ) 1 16.90 6 78.62 10.4
## 7 1 (+DH-Si ) 2 18.90 6 104.22 10.7
## 8 1 (+DH-Si ) 2 20.00 5 96.56 9.8
## 9 1 (+DH+Si1) 1 18.40 6 76.55 9.1
## 10 1 (+DH+Si1) 1 17.00 6 85.95 9.7
## 11 1 (+DH+Si1) 2 17.90 6 81.50 9.8
## 12 1 (+DH+Si1) 2 20.50 6 82.17 10.1
## 13 1 (+DH+Si2) 1 18.70 5 76.14 9.4
## 14 1 (+DH+Si2) 1 19.70 5 77.65 10.7
## 15 1 (+DH+Si2) 2 21.20 6 132.78 10.9
## 16 1 (+DH+Si2) 2 20.90 5 97.48 12.7
## AREA_FOLIAR LONGITUD_RAIZ DIAMETRO_RAIZ PESO_F_RAIZ PESO_F_AEREA PESO_S_RAIZ
## 1 106.90 5.5 23.11 8.04 4.11 0.44
## 2 92.40 3.5 21.52 9.10 3.33 0.45
## 3 150.60 3.5 15.93 4.14 6.84 0.27
## 4 105.20 4.0 21.88 8.69 4.08 0.49
## 5 105.80 3.0 16.84 4.06 4.23 0.26
## 6 74.00 2.5 19.10 4.52 2.92 0.26
## 7 64.60 3.0 16.10 3.03 2.51 0.21
## 8 72.00 3.0 5.81 0.54 2.57 0.03
## 9 79.20 3.5 20.58 5.19 2.98 0.29
## 10 80.00 2.5 13.20 2.51 3.64 0.15
## 11 72.10 2.5 22.49 0.70 2.66 0.06
## 12 81.90 3.4 7.03 7.72 3.02 0.53
## 13 74.02 3.1 18.87 4.98 2.76 0.32
## 14 64.10 2.4 15.10 2.91 2.44 0.14
## 15 117.40 3.5 24.16 13.55 4.23 0.77
## 16 115.90 4.5 27.16 12.70 3.98 0.71
## PESO_S_AEREA SPAD PESO_F_BOCADOS PESO_SAT_H_BOCADOS PESO_S_BOCADOS CRA
## 1 0.33 NA 0.065 0.088 0.010 70.512821
## 2 0.28 NA 0.070 0.098 0.008 68.888889
## 3 0.62 NA 0.085 0.100 0.010 83.333333
## 4 0.32 NA 0.073 0.112 0.009 62.135922
## 5 0.33 NA 0.074 0.096 0.001 76.842105
## 6 0.24 NA 0.065 0.088 0.007 71.604938
## 7 0.25 NA 0.076 0.940 0.004 7.692308
## 8 0.24 NA 0.060 0.980 0.005 5.641026
## 9 0.32 NA 0.087 0.121 0.001 71.666667
## 10 0.30 NA 0.075 0.112 0.004 65.740741
## 11 0.25 NA 0.082 0.122 0.003 66.386555
## 12 0.24 NA 0.060 0.089 0.001 67.045455
## 13 0.24 NA 0.088 0.106 0.008 81.632653
## 14 0.13 NA 0.066 0.084 0.039 60.000000
## 15 0.39 NA 0.067 0.081 0.009 80.555556
## 16 0.35 NA 0.064 0.080 0.020 73.333333
## N_ESTOMAS_ADAX N_ESTOMAS_ABA N_ESTOMAS_ABIERTOS N_ESTOMAS_P_ABIERTOS
## 1 26 39 18 22
## 2 16 29 5 30
## 3 11 20 15 9
## 4 18 34 19 21
## 5 40 45 16 26
## 6 40 21 5 26
## 7 20 26 16 11
## 8 11 27 4 19
## 9 21 41 14 23
## 10 16 14 4 13
## 11 12 16 4 12
## 12 23 28 8 24
## 13 22 24 5 24
## 14 22 26 7 11
## 15 30 48 13 41
## 16 29 55 13 46
## N_ESTOMAS_CERRADOS
## 1 25
## 2 10
## 3 7
## 4 12
## 5 43
## 6 30
## 7 19
## 8 15
## 9 25
## 10 13
## 11 12
## 12 19
## 13 17
## 14 30
## 15 24
## 16 25#Analisis descriptivo Y ESTADISTICO
ggplot(datos_m1)+
aes(TRATAMIENTO,TEMPERATURA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=TEMPERATURA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_T=aov(TEMPERATURA~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_T)## Analysis of Variance Table
##
## Response: TEMPERATURA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 7.6452 7.6452 5.8362 0.04212 *
## TRATAMIENTO 3 6.0657 2.0219 1.5435 0.27677
## BLOQUES:TRATAMIENTO 3 0.4107 0.1369 0.1045 0.95512
## Residuals 8 10.4798 1.3100
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_T)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TEMPERATURA ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 1.3825 0.0628414 2.702159 0.0421243
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.290 -2.881705 2.301705 0.9831182
## (+DH+Si1)-(-DH-Si ) -0.715 -3.306705 1.876705 0.8135346
## (+DH+Si2)-(-DH-Si ) 0.960 -1.631705 3.551705 0.6513763
## (+DH+Si1)-(+DH-Si ) -0.425 -3.016705 2.166705 0.9505884
## (+DH+Si2)-(+DH-Si ) 1.250 -1.341705 3.841705 0.4573683
## (+DH+Si2)-(+DH+Si1) 1.675 -0.916705 4.266705 0.2411400
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 1.030000e+00 -3.499061 5.559061 0.9776138
## 1:(+DH-Si )-1:(-DH-Si ) -3.500000e-01 -4.879061 4.179061 0.9999724
## 2:(+DH-Si )-1:(-DH-Si ) 8.000000e-01 -3.729061 5.329061 0.9945050
## 1:(+DH+Si1)-1:(-DH-Si ) -9.500000e-01 -5.479061 3.579061 0.9854728
## 2:(+DH+Si1)-1:(-DH-Si ) 5.500000e-01 -3.979061 5.079061 0.9994528
## 1:(+DH+Si2)-1:(-DH-Si ) 5.500000e-01 -3.979061 5.079061 0.9994528
## 2:(+DH+Si2)-1:(-DH-Si ) 2.400000e+00 -2.129061 6.929061 0.4840288
## 1:(+DH-Si )-2:(-DH-Si ) -1.380000e+00 -5.909061 3.149061 0.9093774
## 2:(+DH-Si )-2:(-DH-Si ) -2.300000e-01 -4.759061 4.299061 0.9999984
## 1:(+DH+Si1)-2:(-DH-Si ) -1.980000e+00 -6.509061 2.549061 0.6735981
## 2:(+DH+Si1)-2:(-DH-Si ) -4.800000e-01 -5.009061 4.049061 0.9997735
## 1:(+DH+Si2)-2:(-DH-Si ) -4.800000e-01 -5.009061 4.049061 0.9997735
## 2:(+DH+Si2)-2:(-DH-Si ) 1.370000e+00 -3.159061 5.899061 0.9121605
## 2:(+DH-Si )-1:(+DH-Si ) 1.150000e+00 -3.379061 5.679061 0.9608295
## 1:(+DH+Si1)-1:(+DH-Si ) -6.000000e-01 -5.129061 3.929061 0.9990482
## 2:(+DH+Si1)-1:(+DH-Si ) 9.000000e-01 -3.629061 5.429061 0.9892224
## 1:(+DH+Si2)-1:(+DH-Si ) 9.000000e-01 -3.629061 5.429061 0.9892224
## 2:(+DH+Si2)-1:(+DH-Si ) 2.750000e+00 -1.779061 7.279061 0.3485156
## 1:(+DH+Si1)-2:(+DH-Si ) -1.750000e+00 -6.279061 2.779061 0.7761588
## 2:(+DH+Si1)-2:(+DH-Si ) -2.500000e-01 -4.779061 4.279061 0.9999972
## 1:(+DH+Si2)-2:(+DH-Si ) -2.500000e-01 -4.779061 4.279061 0.9999972
## 2:(+DH+Si2)-2:(+DH-Si ) 1.600000e+00 -2.929061 6.129061 0.8365204
## 2:(+DH+Si1)-1:(+DH+Si1) 1.500000e+00 -3.029061 6.029061 0.8722729
## 1:(+DH+Si2)-1:(+DH+Si1) 1.500000e+00 -3.029061 6.029061 0.8722729
## 2:(+DH+Si2)-1:(+DH+Si1) 3.350000e+00 -1.179061 7.879061 0.1857278
## 1:(+DH+Si2)-2:(+DH+Si1) 3.552714e-15 -4.529061 4.529061 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) 1.850000e+00 -2.679061 6.379061 0.7325928
## 2:(+DH+Si2)-1:(+DH+Si2) 1.850000e+00 -2.679061 6.379061 0.7325928ggplot(datos_m1)+
aes(TRATAMIENTO,N_HOJAS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=N_HOJAS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_NH=aov(N_HOJAS~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_NH)## Analysis of Variance Table
##
## Response: N_HOJAS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.00 0.00000 0.0000 1.0000
## TRATAMIENTO 3 1.25 0.41667 1.6667 0.2503
## BLOQUES:TRATAMIENTO 3 0.50 0.16667 0.6667 0.5957
## Residuals 8 2.00 0.25000TukeyHSD(mod1_NH)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_HOJAS ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0 -0.576501 0.576501 1
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.25 -1.3822024 0.8822024 0.8915309
## (+DH+Si1)-(-DH-Si ) 0.25 -0.8822024 1.3822024 0.8915309
## (+DH+Si2)-(-DH-Si ) -0.50 -1.6322024 0.6322024 0.5252407
## (+DH+Si1)-(+DH-Si ) 0.50 -0.6322024 1.6322024 0.5252407
## (+DH+Si2)-(+DH-Si ) -0.25 -1.3822024 0.8822024 0.8915309
## (+DH+Si2)-(+DH+Si1) -0.75 -1.8822024 0.3822024 0.2252476
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.5 -2.478549 1.4785485 0.961735
## 1:(+DH-Si )-1:(-DH-Si ) -0.5 -2.478549 1.4785485 0.961735
## 2:(+DH-Si )-1:(-DH-Si ) -0.5 -2.478549 1.4785485 0.961735
## 1:(+DH+Si1)-1:(-DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 2:(+DH+Si1)-1:(-DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 1:(+DH+Si2)-1:(-DH-Si ) -1.0 -2.978549 0.9785485 0.532305
## 2:(+DH+Si2)-1:(-DH-Si ) -0.5 -2.478549 1.4785485 0.961735
## 1:(+DH-Si )-2:(-DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 2:(+DH-Si )-2:(-DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 1:(+DH+Si1)-2:(-DH-Si ) 0.5 -1.478549 2.4785485 0.961735
## 2:(+DH+Si1)-2:(-DH-Si ) 0.5 -1.478549 2.4785485 0.961735
## 1:(+DH+Si2)-2:(-DH-Si ) -0.5 -2.478549 1.4785485 0.961735
## 2:(+DH+Si2)-2:(-DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 2:(+DH-Si )-1:(+DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 1:(+DH+Si1)-1:(+DH-Si ) 0.5 -1.478549 2.4785485 0.961735
## 2:(+DH+Si1)-1:(+DH-Si ) 0.5 -1.478549 2.4785485 0.961735
## 1:(+DH+Si2)-1:(+DH-Si ) -0.5 -2.478549 1.4785485 0.961735
## 2:(+DH+Si2)-1:(+DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 1:(+DH+Si1)-2:(+DH-Si ) 0.5 -1.478549 2.4785485 0.961735
## 2:(+DH+Si1)-2:(+DH-Si ) 0.5 -1.478549 2.4785485 0.961735
## 1:(+DH+Si2)-2:(+DH-Si ) -0.5 -2.478549 1.4785485 0.961735
## 2:(+DH+Si2)-2:(+DH-Si ) 0.0 -1.978549 1.9785485 1.000000
## 2:(+DH+Si1)-1:(+DH+Si1) 0.0 -1.978549 1.9785485 1.000000
## 1:(+DH+Si2)-1:(+DH+Si1) -1.0 -2.978549 0.9785485 0.532305
## 2:(+DH+Si2)-1:(+DH+Si1) -0.5 -2.478549 1.4785485 0.961735
## 1:(+DH+Si2)-2:(+DH+Si1) -1.0 -2.978549 0.9785485 0.532305
## 2:(+DH+Si2)-2:(+DH+Si1) -0.5 -2.478549 1.4785485 0.961735
## 2:(+DH+Si2)-1:(+DH+Si2) 0.5 -1.478549 2.4785485 0.961735ggplot(datos_m1)+
aes(TRATAMIENTO,LONGITUD_P_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=LONGITUD_P_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PA=aov(LONGITUD_P_AEREA~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PA)## Analysis of Variance Table
##
## Response: LONGITUD_P_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 5.4056 5.4056 6.1210 0.03847 *
## TRATAMIENTO 3 18.8819 6.2940 7.1269 0.01196 *
## BLOQUES:TRATAMIENTO 3 4.2319 1.4106 1.5973 0.26482
## Residuals 8 7.0650 0.8831
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_PA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = LONGITUD_P_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 1.1625 0.07896941 2.246031 0.0384656
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -2.275 -4.4029683 -0.1470317 0.0366713
## (+DH+Si1)-(-DH-Si ) -2.925 -5.0529683 -0.7970317 0.0098074
## (+DH+Si2)-(-DH-Si ) -1.675 -3.8029683 0.4529683 0.1305992
## (+DH+Si1)-(+DH-Si ) -0.650 -2.7779683 1.4779683 0.7653838
## (+DH+Si2)-(+DH-Si ) 0.600 -1.5279683 2.7279683 0.8039300
## (+DH+Si2)-(+DH+Si1) 1.250 -0.8779683 3.3779683 0.3072291
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 2.50 -1.218671 6.21867128 0.2577524
## 1:(+DH-Si )-1:(-DH-Si ) -0.95 -4.668671 2.76867128 0.9596427
## 2:(+DH-Si )-1:(-DH-Si ) -1.10 -4.818671 2.61867128 0.9202886
## 1:(+DH+Si1)-1:(-DH-Si ) -1.95 -5.668671 1.76867128 0.4947462
## 2:(+DH+Si1)-1:(-DH-Si ) -1.40 -5.118671 2.31867128 0.7949526
## 1:(+DH+Si2)-1:(-DH-Si ) -1.30 -5.018671 2.41867128 0.8427684
## 2:(+DH+Si2)-1:(-DH-Si ) 0.45 -3.268671 4.16867128 0.9994650
## 1:(+DH-Si )-2:(-DH-Si ) -3.45 -7.168671 0.26867128 0.0719662
## 2:(+DH-Si )-2:(-DH-Si ) -3.60 -7.318671 0.11867128 0.0587000
## 1:(+DH+Si1)-2:(-DH-Si ) -4.45 -8.168671 -0.73132872 0.0190799
## 2:(+DH+Si1)-2:(-DH-Si ) -3.90 -7.618671 -0.18132872 0.0392013
## 1:(+DH+Si2)-2:(-DH-Si ) -3.80 -7.518671 -0.08132872 0.0448174
## 2:(+DH+Si2)-2:(-DH-Si ) -2.05 -5.768671 1.66867128 0.4437897
## 2:(+DH-Si )-1:(+DH-Si ) -0.15 -3.868671 3.56867128 0.9999997
## 1:(+DH+Si1)-1:(+DH-Si ) -1.00 -4.718671 2.71867128 0.9483619
## 2:(+DH+Si1)-1:(+DH-Si ) -0.45 -4.168671 3.26867128 0.9994650
## 1:(+DH+Si2)-1:(+DH-Si ) -0.35 -4.068671 3.36867128 0.9998965
## 2:(+DH+Si2)-1:(+DH-Si ) 1.40 -2.318671 5.11867128 0.7949526
## 1:(+DH+Si1)-2:(+DH-Si ) -0.85 -4.568671 2.86867128 0.9770120
## 2:(+DH+Si1)-2:(+DH-Si ) -0.30 -4.018671 3.41867128 0.9999631
## 1:(+DH+Si2)-2:(+DH-Si ) -0.20 -3.918671 3.51867128 0.9999977
## 2:(+DH+Si2)-2:(+DH-Si ) 1.55 -2.168671 5.26867128 0.7156532
## 2:(+DH+Si1)-1:(+DH+Si1) 0.55 -3.168671 4.26867128 0.9981086
## 1:(+DH+Si2)-1:(+DH+Si1) 0.65 -3.068671 4.36867128 0.9948323
## 2:(+DH+Si2)-1:(+DH+Si1) 2.40 -1.318671 6.11867128 0.2926012
## 1:(+DH+Si2)-2:(+DH+Si1) 0.10 -3.618671 3.81867128 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) 1.85 -1.868671 5.56867128 0.5483455
## 2:(+DH+Si2)-1:(+DH+Si2) 1.75 -1.968671 5.46867128 0.6037760ggplot(datos_m1)+
aes(TRATAMIENTO,AREA_FOLIAR, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=AREA_FOLIAR, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_AF=aov(AREA_FOLIAR~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_AF)## Analysis of Variance Table
##
## Response: AREA_FOLIAR
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 666.7 666.67 3.0177 0.12057
## TRATAMIENTO 3 3300.2 1100.06 4.9794 0.03089 *
## BLOQUES:TRATAMIENTO 3 2869.5 956.51 4.3296 0.04326 *
## Residuals 8 1767.4 220.92
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_AF)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AREA_FOLIAR ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 12.91 -4.227569 30.04757 0.1205656
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -34.675 -68.33183 -1.01817 0.0436386
## (+DH+Si1)-(-DH-Si ) -35.475 -69.13183 -1.81817 0.0392277
## (+DH+Si2)-(-DH-Si ) -20.920 -54.57683 12.73683 0.2673194
## (+DH+Si1)-(+DH-Si ) -0.800 -34.45683 32.85683 0.9998278
## (+DH+Si2)-(+DH-Si ) 13.755 -19.90183 47.41183 0.5828782
## (+DH+Si2)-(+DH+Si1) 14.555 -19.10183 48.21183 0.5410836
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 28.25 -30.56605 87.06604792 0.5836054
## 1:(+DH-Si )-1:(-DH-Si ) -9.75 -68.56605 49.06604792 0.9962232
## 2:(+DH-Si )-1:(-DH-Si ) -31.35 -90.16605 27.46604792 0.4780598
## 1:(+DH+Si1)-1:(-DH-Si ) -20.05 -78.86605 38.76604792 0.8570710
## 2:(+DH+Si1)-1:(-DH-Si ) -22.65 -81.46605 36.16604792 0.7786432
## 1:(+DH+Si2)-1:(-DH-Si ) -30.59 -89.40605 28.22604792 0.5031274
## 2:(+DH+Si2)-1:(-DH-Si ) 17.00 -41.81605 75.81604792 0.9280521
## 1:(+DH-Si )-2:(-DH-Si ) -38.00 -96.81605 20.81604792 0.2916604
## 2:(+DH-Si )-2:(-DH-Si ) -59.60 -118.41605 -0.78395208 0.0467717
## 1:(+DH+Si1)-2:(-DH-Si ) -48.30 -107.11605 10.51604792 0.1234704
## 2:(+DH+Si1)-2:(-DH-Si ) -50.90 -109.71605 7.91604792 0.0987110
## 1:(+DH+Si2)-2:(-DH-Si ) -58.84 -117.65605 -0.02395208 0.0498980
## 2:(+DH+Si2)-2:(-DH-Si ) -11.25 -70.06605 47.56604792 0.9913071
## 2:(+DH-Si )-1:(+DH-Si ) -21.60 -80.41605 37.21604792 0.8119269
## 1:(+DH+Si1)-1:(+DH-Si ) -10.30 -69.11605 48.51604792 0.9947751
## 2:(+DH+Si1)-1:(+DH-Si ) -12.90 -71.71605 45.91604792 0.9815173
## 1:(+DH+Si2)-1:(+DH-Si ) -20.84 -79.65605 37.97604792 0.8347201
## 2:(+DH+Si2)-1:(+DH-Si ) 26.75 -32.06605 85.56604792 0.6367610
## 1:(+DH+Si1)-2:(+DH-Si ) 11.30 -47.51605 70.11604792 0.9910857
## 2:(+DH+Si1)-2:(+DH-Si ) 8.70 -50.11605 67.51604792 0.9981073
## 1:(+DH+Si2)-2:(+DH-Si ) 0.76 -58.05605 59.57604792 1.0000000
## 2:(+DH+Si2)-2:(+DH-Si ) 48.35 -10.46605 107.16604792 0.1229413
## 2:(+DH+Si1)-1:(+DH+Si1) -2.60 -61.41605 56.21604792 0.9999994
## 1:(+DH+Si2)-1:(+DH+Si1) -10.54 -69.35605 48.27604792 0.9940231
## 2:(+DH+Si2)-1:(+DH+Si1) 37.05 -21.76605 95.86604792 0.3143171
## 1:(+DH+Si2)-2:(+DH+Si1) -7.94 -66.75605 50.87604792 0.9989283
## 2:(+DH+Si2)-2:(+DH+Si1) 39.65 -19.16605 98.46604792 0.2554856
## 2:(+DH+Si2)-1:(+DH+Si2) 47.59 -11.22605 106.40604792 0.1312268ggplot(datos_m1)+
aes(TRATAMIENTO,LONGITUD_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=LONGITUD_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_LR=aov(LONGITUD_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_LR)## Analysis of Variance Table
##
## Response: LONGITUD_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.1225 0.12250 0.2513 0.6297
## TRATAMIENTO 3 3.8675 1.28917 2.6444 0.1207
## BLOQUES:TRATAMIENTO 3 2.0675 0.68917 1.4137 0.3083
## Residuals 8 3.9000 0.48750TukeyHSD(mod1_LR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = LONGITUD_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.175 -0.6300399 0.9800399 0.6296742
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -1.25 -2.831035 0.3310346 0.1285823
## (+DH+Si1)-(-DH-Si ) -1.15 -2.731035 0.4310346 0.1701474
## (+DH+Si2)-(-DH-Si ) -0.75 -2.331035 0.8310346 0.4702462
## (+DH+Si1)-(+DH-Si ) 0.10 -1.481035 1.6810346 0.9968132
## (+DH+Si2)-(+DH-Si ) 0.50 -1.081035 2.0810346 0.7470395
## (+DH+Si2)-(+DH+Si1) 0.40 -1.181035 1.9810346 0.8480631
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -7.500000e-01 -3.512893 2.012893 0.9460255
## 1:(+DH-Si )-1:(-DH-Si ) -1.750000e+00 -4.512893 1.012893 0.3093440
## 2:(+DH-Si )-1:(-DH-Si ) -1.500000e+00 -4.262893 1.262893 0.4593145
## 1:(+DH+Si1)-1:(-DH-Si ) -1.500000e+00 -4.262893 1.262893 0.4593145
## 2:(+DH+Si1)-1:(-DH-Si ) -1.550000e+00 -4.312893 1.212893 0.4261135
## 1:(+DH+Si2)-1:(-DH-Si ) -1.750000e+00 -4.512893 1.012893 0.3093440
## 2:(+DH+Si2)-1:(-DH-Si ) -5.000000e-01 -3.262893 2.262893 0.9936727
## 1:(+DH-Si )-2:(-DH-Si ) -1.000000e+00 -3.762893 1.762893 0.8214320
## 2:(+DH-Si )-2:(-DH-Si ) -7.500000e-01 -3.512893 2.012893 0.9460255
## 1:(+DH+Si1)-2:(-DH-Si ) -7.500000e-01 -3.512893 2.012893 0.9460255
## 2:(+DH+Si1)-2:(-DH-Si ) -8.000000e-01 -3.562893 1.962893 0.9274780
## 1:(+DH+Si2)-2:(-DH-Si ) -1.000000e+00 -3.762893 1.762893 0.8214320
## 2:(+DH+Si2)-2:(-DH-Si ) 2.500000e-01 -2.512893 3.012893 0.9999203
## 2:(+DH-Si )-1:(+DH-Si ) 2.500000e-01 -2.512893 3.012893 0.9999203
## 1:(+DH+Si1)-1:(+DH-Si ) 2.500000e-01 -2.512893 3.012893 0.9999203
## 2:(+DH+Si1)-1:(+DH-Si ) 2.000000e-01 -2.562893 2.962893 0.9999822
## 1:(+DH+Si2)-1:(+DH-Si ) 8.881784e-16 -2.762893 2.762893 1.0000000
## 2:(+DH+Si2)-1:(+DH-Si ) 1.250000e+00 -1.512893 4.012893 0.6417479
## 1:(+DH+Si1)-2:(+DH-Si ) 0.000000e+00 -2.762893 2.762893 1.0000000
## 2:(+DH+Si1)-2:(+DH-Si ) -5.000000e-02 -2.812893 2.712893 1.0000000
## 1:(+DH+Si2)-2:(+DH-Si ) -2.500000e-01 -3.012893 2.512893 0.9999203
## 2:(+DH+Si2)-2:(+DH-Si ) 1.000000e+00 -1.762893 3.762893 0.8214320
## 2:(+DH+Si1)-1:(+DH+Si1) -5.000000e-02 -2.812893 2.712893 1.0000000
## 1:(+DH+Si2)-1:(+DH+Si1) -2.500000e-01 -3.012893 2.512893 0.9999203
## 2:(+DH+Si2)-1:(+DH+Si1) 1.000000e+00 -1.762893 3.762893 0.8214320
## 1:(+DH+Si2)-2:(+DH+Si1) -2.000000e-01 -2.962893 2.562893 0.9999822
## 2:(+DH+Si2)-2:(+DH+Si1) 1.050000e+00 -1.712893 3.812893 0.7882720
## 2:(+DH+Si2)-1:(+DH+Si2) 1.250000e+00 -1.512893 4.012893 0.6417479ggplot(datos_m1)+
aes(TRATAMIENTO,DIAMETRO_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=DIAMETRO_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_DR=aov(DIAMETRO_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_DR)## Analysis of Variance Table
##
## Response: DIAMETRO_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 3.764 3.764 0.1293 0.7284
## TRATAMIENTO 3 140.334 46.778 1.6075 0.2626
## BLOQUES:TRATAMIENTO 3 136.867 45.622 1.5677 0.2713
## Residuals 8 232.806 29.101TukeyHSD(mod1_DR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = DIAMETRO_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.97 -7.189877 5.249877 0.7284357
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -6.1475 -18.362846 6.067846 0.4243124
## (+DH+Si1)-(-DH-Si ) -4.7850 -17.000346 7.430346 0.6131436
## (+DH+Si2)-(-DH-Si ) 0.7125 -11.502846 12.927846 0.9974928
## (+DH+Si1)-(+DH-Si ) 1.3625 -10.852846 13.577846 0.9832717
## (+DH+Si2)-(+DH-Si ) 6.8600 -5.355346 19.075346 0.3403502
## (+DH+Si2)-(+DH+Si1) 5.4975 -6.717846 17.712846 0.5108262
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -3.410 -24.756586 17.936586 0.9969755
## 1:(+DH-Si )-1:(-DH-Si ) -4.345 -25.691586 17.001586 0.9876836
## 2:(+DH-Si )-1:(-DH-Si ) -11.360 -32.706586 9.986586 0.4796861
## 1:(+DH+Si1)-1:(-DH-Si ) -5.425 -26.771586 15.921586 0.9606602
## 2:(+DH+Si1)-1:(-DH-Si ) -7.555 -28.901586 13.791586 0.8354137
## 1:(+DH+Si2)-1:(-DH-Si ) -5.330 -26.676586 16.016586 0.9639438
## 2:(+DH+Si2)-1:(-DH-Si ) 3.345 -18.001586 24.691586 0.9973078
## 1:(+DH-Si )-2:(-DH-Si ) -0.935 -22.281586 20.411586 0.9999994
## 2:(+DH-Si )-2:(-DH-Si ) -7.950 -29.296586 13.396586 0.8024720
## 1:(+DH+Si1)-2:(-DH-Si ) -2.015 -23.361586 19.331586 0.9998944
## 2:(+DH+Si1)-2:(-DH-Si ) -4.145 -25.491586 17.201586 0.9905340
## 1:(+DH+Si2)-2:(-DH-Si ) -1.920 -23.266586 19.426586 0.9999234
## 2:(+DH+Si2)-2:(-DH-Si ) 6.755 -14.591586 28.101586 0.8937599
## 2:(+DH-Si )-1:(+DH-Si ) -7.015 -28.361586 14.331586 0.8761622
## 1:(+DH+Si1)-1:(+DH-Si ) -1.080 -22.426586 20.266586 0.9999985
## 2:(+DH+Si1)-1:(+DH-Si ) -3.210 -24.556586 18.136586 0.9979058
## 1:(+DH+Si2)-1:(+DH-Si ) -0.985 -22.331586 20.361586 0.9999992
## 2:(+DH+Si2)-1:(+DH-Si ) 7.690 -13.656586 29.036586 0.8244256
## 1:(+DH+Si1)-2:(+DH-Si ) 5.935 -15.411586 27.281586 0.9396745
## 2:(+DH+Si1)-2:(+DH-Si ) 3.805 -17.541586 25.151586 0.9942059
## 1:(+DH+Si2)-2:(+DH-Si ) 6.030 -15.316586 27.376586 0.9351248
## 2:(+DH+Si2)-2:(+DH-Si ) 14.705 -6.641586 36.051586 0.2380580
## 2:(+DH+Si1)-1:(+DH+Si1) -2.130 -23.476586 19.216586 0.9998476
## 1:(+DH+Si2)-1:(+DH+Si1) 0.095 -21.251586 21.441586 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si1) 8.770 -12.576586 30.116586 0.7278119
## 1:(+DH+Si2)-2:(+DH+Si1) 2.225 -19.121586 23.571586 0.9997970
## 2:(+DH+Si2)-2:(+DH+Si1) 10.900 -10.446586 32.246586 0.5218885
## 2:(+DH+Si2)-1:(+DH+Si2) 8.675 -12.671586 30.021586 0.7367913ggplot(datos_m1)+
aes(TRATAMIENTO,PESO_F_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=PESO_F_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PFR=aov(PESO_F_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PFR)## Analysis of Variance Table
##
## Response: PESO_F_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 5.954 5.9536 1.0619 0.33294
## TRATAMIENTO 3 84.425 28.1418 5.0193 0.03029 *
## BLOQUES:TRATAMIENTO 3 89.367 29.7892 5.3131 0.02625 *
## Residuals 8 44.854 5.6067
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_PFR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 1.22 -1.510145 3.950145 0.3329351
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -4.4550 -9.8167883 0.9067883 0.1073442
## (+DH+Si1)-(-DH-Si ) -3.4625 -8.8242883 1.8992883 0.2416670
## (+DH+Si2)-(-DH-Si ) 1.0425 -4.3192883 6.4042883 0.9219437
## (+DH+Si1)-(+DH-Si ) 0.9925 -4.3692883 6.3542883 0.9314715
## (+DH+Si2)-(+DH-Si ) 5.4975 0.1357117 10.8592883 0.0446179
## (+DH+Si2)-(+DH+Si1) 4.5050 -0.8567883 9.8667883 0.1029307
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -2.155 -11.52484261 7.214843 0.9762590
## 1:(+DH-Si )-1:(-DH-Si ) -4.280 -13.64984261 5.089843 0.6326253
## 2:(+DH-Si )-1:(-DH-Si ) -6.785 -16.15484261 2.584843 0.2005099
## 1:(+DH+Si1)-1:(-DH-Si ) -4.720 -14.08984261 4.649843 0.5356824
## 2:(+DH+Si1)-1:(-DH-Si ) -4.360 -13.72984261 5.009843 0.6147783
## 1:(+DH+Si2)-1:(-DH-Si ) -4.625 -13.99484261 4.744843 0.5562620
## 2:(+DH+Si2)-1:(-DH-Si ) 4.555 -4.81484261 13.924843 0.5715846
## 1:(+DH-Si )-2:(-DH-Si ) -2.125 -11.49484261 7.244843 0.9779363
## 2:(+DH-Si )-2:(-DH-Si ) -4.630 -13.99984261 4.739843 0.5551723
## 1:(+DH+Si1)-2:(-DH-Si ) -2.565 -11.93484261 6.804843 0.9438582
## 2:(+DH+Si1)-2:(-DH-Si ) -2.205 -11.57484261 7.164843 0.9732678
## 1:(+DH+Si2)-2:(-DH-Si ) -2.470 -11.83984261 6.899843 0.9530256
## 2:(+DH+Si2)-2:(-DH-Si ) 6.710 -2.65984261 16.079843 0.2085355
## 2:(+DH-Si )-1:(+DH-Si ) -2.505 -11.87484261 6.864843 0.9497692
## 1:(+DH+Si1)-1:(+DH-Si ) -0.440 -9.80984261 8.929843 0.9999991
## 2:(+DH+Si1)-1:(+DH-Si ) -0.080 -9.44984261 9.289843 1.0000000
## 1:(+DH+Si2)-1:(+DH-Si ) -0.345 -9.71484261 9.024843 0.9999998
## 2:(+DH+Si2)-1:(+DH-Si ) 8.835 -0.53484261 18.204843 0.0666512
## 1:(+DH+Si1)-2:(+DH-Si ) 2.065 -7.30484261 11.434843 0.9810361
## 2:(+DH+Si1)-2:(+DH-Si ) 2.425 -6.94484261 11.794843 0.9570072
## 1:(+DH+Si2)-2:(+DH-Si ) 2.160 -7.20984261 11.529843 0.9759710
## 2:(+DH+Si2)-2:(+DH-Si ) 11.340 1.97015739 20.709843 0.0178856
## 2:(+DH+Si1)-1:(+DH+Si1) 0.360 -9.00984261 9.729843 0.9999998
## 1:(+DH+Si2)-1:(+DH+Si1) 0.095 -9.27484261 9.464843 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si1) 9.275 -0.09484261 18.644843 0.0526052
## 1:(+DH+Si2)-2:(+DH+Si1) -0.265 -9.63484261 9.104843 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) 8.915 -0.45484261 18.284843 0.0638372
## 2:(+DH+Si2)-1:(+DH+Si2) 9.180 -0.18984261 18.549843 0.0553556ggplot(datos_m1)+
aes(TRATAMIENTO,PESO_F_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=PESO_F_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PFA=aov(PESO_F_AEREA~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PFA)## Analysis of Variance Table
##
## Response: PESO_F_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.7569 0.75690 1.1344 0.31793
## TRATAMIENTO 3 6.3395 2.11318 3.1671 0.08543 .
## BLOQUES:TRATAMIENTO 3 5.8279 1.94262 2.9114 0.10082
## Residuals 8 5.3379 0.66724
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_PFA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.435 -0.5068252 1.376825 0.3179302
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -1.5325 -3.38217 0.3171702 0.1084867
## (+DH+Si1)-(-DH-Si ) -1.5150 -3.36467 0.3346702 0.1131995
## (+DH+Si2)-(-DH-Si ) -1.2375 -3.08717 0.6121702 0.2189980
## (+DH+Si1)-(+DH-Si ) 0.0175 -1.83217 1.8671702 0.9999891
## (+DH+Si2)-(+DH-Si ) 0.2950 -1.55467 2.1446702 0.9542120
## (+DH+Si2)-(+DH+Si1) 0.2775 -1.57217 2.1271702 0.9613236
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 1.740 -1.492339 4.9723392 0.4679807
## 1:(+DH-Si )-1:(-DH-Si ) -0.145 -3.377339 3.0873392 0.9999993
## 2:(+DH-Si )-1:(-DH-Si ) -1.180 -4.412339 2.0523392 0.8158653
## 1:(+DH+Si1)-1:(-DH-Si ) -0.410 -3.642339 2.8223392 0.9992773
## 2:(+DH+Si1)-1:(-DH-Si ) -0.880 -4.112339 2.3523392 0.9452821
## 1:(+DH+Si2)-1:(-DH-Si ) -1.120 -4.352339 2.1123392 0.8479206
## 2:(+DH+Si2)-1:(-DH-Si ) 0.385 -2.847339 3.6173392 0.9995168
## 1:(+DH-Si )-2:(-DH-Si ) -1.885 -5.117339 1.3473392 0.3875443
## 2:(+DH-Si )-2:(-DH-Si ) -2.920 -6.152339 0.3123392 0.0814267
## 1:(+DH+Si1)-2:(-DH-Si ) -2.150 -5.382339 1.0823392 0.2666363
## 2:(+DH+Si1)-2:(-DH-Si ) -2.620 -5.852339 0.6123392 0.1302868
## 1:(+DH+Si2)-2:(-DH-Si ) -2.860 -6.092339 0.3723392 0.0894631
## 2:(+DH+Si2)-2:(-DH-Si ) -1.355 -4.587339 1.8773392 0.7107721
## 2:(+DH-Si )-1:(+DH-Si ) -1.035 -4.267339 2.1973392 0.8884801
## 1:(+DH+Si1)-1:(+DH-Si ) -0.265 -3.497339 2.9673392 0.9999588
## 2:(+DH+Si1)-1:(+DH-Si ) -0.735 -3.967339 2.4973392 0.9776296
## 1:(+DH+Si2)-1:(+DH-Si ) -0.975 -4.207339 2.2573392 0.9132211
## 2:(+DH+Si2)-1:(+DH-Si ) 0.530 -2.702339 3.7623392 0.9964621
## 1:(+DH+Si1)-2:(+DH-Si ) 0.770 -2.462339 4.0023392 0.9715442
## 2:(+DH+Si1)-2:(+DH-Si ) 0.300 -2.932339 3.5323392 0.9999056
## 1:(+DH+Si2)-2:(+DH-Si ) 0.060 -3.172339 3.2923392 1.0000000
## 2:(+DH+Si2)-2:(+DH-Si ) 1.565 -1.667339 4.7973392 0.5756333
## 2:(+DH+Si1)-1:(+DH+Si1) -0.470 -3.702339 2.7623392 0.9982972
## 1:(+DH+Si2)-1:(+DH+Si1) -0.710 -3.942339 2.5223392 0.9813697
## 2:(+DH+Si2)-1:(+DH+Si1) 0.795 -2.437339 4.0273392 0.9665545
## 1:(+DH+Si2)-2:(+DH+Si1) -0.240 -3.472339 2.9923392 0.9999789
## 2:(+DH+Si2)-2:(+DH+Si1) 1.265 -1.967339 4.4973392 0.7665537
## 2:(+DH+Si2)-1:(+DH+Si2) 1.505 -1.727339 4.7373392 0.6141845ggplot(datos_m1)+
aes(TRATAMIENTO,PESO_S_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=PESO_S_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PSR=aov(PESO_S_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PSR)## Analysis of Variance Table
##
## Response: PESO_S_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.03610 0.036100 1.6161 0.23935
## TRATAMIENTO 3 0.22212 0.074042 3.3147 0.07786 .
## BLOQUES:TRATAMIENTO 3 0.25345 0.084483 3.7821 0.05879 .
## Residuals 8 0.17870 0.022338
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_PSR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.095 -0.07732464 0.2673246 0.2393485
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.2225 -0.56093195 0.1159319 0.2300569
## (+DH+Si1)-(-DH-Si ) -0.1550 -0.49343195 0.1834319 0.4973977
## (+DH+Si2)-(-DH-Si ) 0.0725 -0.26593195 0.4109319 0.8995885
## (+DH+Si1)-(+DH-Si ) 0.0675 -0.27093195 0.4059319 0.9165468
## (+DH+Si2)-(+DH-Si ) 0.2950 -0.04343195 0.6334319 0.0893079
## (+DH+Si2)-(+DH+Si1) 0.2275 -0.11093195 0.5659319 0.2160466
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.065 -0.65641724 0.5264172 0.9997129
## 1:(+DH-Si )-1:(-DH-Si ) -0.185 -0.77641724 0.4064172 0.8987396
## 2:(+DH-Si )-1:(-DH-Si ) -0.325 -0.91641724 0.2664172 0.4470033
## 1:(+DH+Si1)-1:(-DH-Si ) -0.225 -0.81641724 0.3664172 0.7875058
## 2:(+DH+Si1)-1:(-DH-Si ) -0.150 -0.74141724 0.4414172 0.9610469
## 1:(+DH+Si2)-1:(-DH-Si ) -0.215 -0.80641724 0.3764172 0.8185975
## 2:(+DH+Si2)-1:(-DH-Si ) 0.295 -0.29641724 0.8864172 0.5456798
## 1:(+DH-Si )-2:(-DH-Si ) -0.120 -0.71141724 0.4714172 0.9878974
## 2:(+DH-Si )-2:(-DH-Si ) -0.260 -0.85141724 0.3314172 0.6685020
## 1:(+DH+Si1)-2:(-DH-Si ) -0.160 -0.75141724 0.4314172 0.9468771
## 2:(+DH+Si1)-2:(-DH-Si ) -0.085 -0.67641724 0.5064172 0.9984156
## 1:(+DH+Si2)-2:(-DH-Si ) -0.150 -0.74141724 0.4414172 0.9610469
## 2:(+DH+Si2)-2:(-DH-Si ) 0.360 -0.23141724 0.9514172 0.3461497
## 2:(+DH-Si )-1:(+DH-Si ) -0.140 -0.73141724 0.4514172 0.9724475
## 1:(+DH+Si1)-1:(+DH-Si ) -0.040 -0.63141724 0.5514172 0.9999888
## 2:(+DH+Si1)-1:(+DH-Si ) 0.035 -0.55641724 0.6264172 0.9999955
## 1:(+DH+Si2)-1:(+DH-Si ) -0.030 -0.62141724 0.5614172 0.9999984
## 2:(+DH+Si2)-1:(+DH-Si ) 0.480 -0.11141724 1.0714172 0.1295978
## 1:(+DH+Si1)-2:(+DH-Si ) 0.100 -0.49141724 0.6914172 0.9957513
## 2:(+DH+Si1)-2:(+DH-Si ) 0.175 -0.41641724 0.7664172 0.9201765
## 1:(+DH+Si2)-2:(+DH-Si ) 0.110 -0.48141724 0.7014172 0.9925905
## 2:(+DH+Si2)-2:(+DH-Si ) 0.620 0.02858276 1.2114172 0.0392855
## 2:(+DH+Si1)-1:(+DH+Si1) 0.075 -0.51641724 0.6664172 0.9992784
## 1:(+DH+Si2)-1:(+DH+Si1) 0.010 -0.58141724 0.6014172 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si1) 0.520 -0.07141724 1.1114172 0.0920236
## 1:(+DH+Si2)-2:(+DH+Si1) -0.065 -0.65641724 0.5264172 0.9997129
## 2:(+DH+Si2)-2:(+DH+Si1) 0.445 -0.14641724 1.0364172 0.1743406
## 2:(+DH+Si2)-1:(+DH+Si2) 0.510 -0.08141724 1.1014172 0.1002618ggplot(datos_m1)+
aes(TRATAMIENTO,PESO_S_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=PESO_S_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PSA=aov(PESO_S_AEREA~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PSA)## Analysis of Variance Table
##
## Response: PESO_S_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.015006 0.0150062 2.0896 0.1863
## TRATAMIENTO 3 0.039519 0.0131729 1.8343 0.2190
## BLOQUES:TRATAMIENTO 3 0.052269 0.0174229 2.4262 0.1407
## Residuals 8 0.057450 0.0071813TukeyHSD(mod1_PSA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.06125 -0.036458 0.158958 0.1863095
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -1.225000e-01 -0.3143908 0.06939077 0.2492709
## (+DH+Si1)-(-DH-Si ) -1.100000e-01 -0.3018908 0.08189077 0.3250897
## (+DH+Si2)-(-DH-Si ) -1.100000e-01 -0.3018908 0.08189077 0.3250897
## (+DH+Si1)-(+DH-Si ) 1.250000e-02 -0.1793908 0.20439077 0.9965231
## (+DH+Si2)-(+DH-Si ) 1.250000e-02 -0.1793908 0.20439077 0.9965231
## (+DH+Si2)-(+DH+Si1) -5.551115e-17 -0.1918908 0.19189077 1.0000000
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.165 -0.1703333 0.50033333 0.5594455
## 1:(+DH-Si )-1:(-DH-Si ) -0.020 -0.3553333 0.31533333 0.9999953
## 2:(+DH-Si )-1:(-DH-Si ) -0.060 -0.3953333 0.27533333 0.9940765
## 1:(+DH+Si1)-1:(-DH-Si ) 0.005 -0.3303333 0.34033333 1.0000000
## 2:(+DH+Si1)-1:(-DH-Si ) -0.060 -0.3953333 0.27533333 0.9940765
## 1:(+DH+Si2)-1:(-DH-Si ) -0.120 -0.4553333 0.21533333 0.8286165
## 2:(+DH+Si2)-1:(-DH-Si ) 0.065 -0.2703333 0.40033333 0.9906270
## 1:(+DH-Si )-2:(-DH-Si ) -0.185 -0.5203333 0.15033333 0.4430223
## 2:(+DH-Si )-2:(-DH-Si ) -0.225 -0.5603333 0.11033333 0.2593644
## 1:(+DH+Si1)-2:(-DH-Si ) -0.160 -0.4953333 0.17533333 0.5901820
## 2:(+DH+Si1)-2:(-DH-Si ) -0.225 -0.5603333 0.11033333 0.2593644
## 1:(+DH+Si2)-2:(-DH-Si ) -0.285 -0.6203333 0.05033333 0.1067837
## 2:(+DH+Si2)-2:(-DH-Si ) -0.100 -0.4353333 0.23533333 0.9174158
## 2:(+DH-Si )-1:(+DH-Si ) -0.040 -0.3753333 0.29533333 0.9995122
## 1:(+DH+Si1)-1:(+DH-Si ) 0.025 -0.3103333 0.36033333 0.9999783
## 2:(+DH+Si1)-1:(+DH-Si ) -0.040 -0.3753333 0.29533333 0.9995122
## 1:(+DH+Si2)-1:(+DH-Si ) -0.100 -0.4353333 0.23533333 0.9174158
## 2:(+DH+Si2)-1:(+DH-Si ) 0.085 -0.2503333 0.42033333 0.9611591
## 1:(+DH+Si1)-2:(+DH-Si ) 0.065 -0.2703333 0.40033333 0.9906270
## 2:(+DH+Si1)-2:(+DH-Si ) 0.000 -0.3353333 0.33533333 1.0000000
## 1:(+DH+Si2)-2:(+DH-Si ) -0.060 -0.3953333 0.27533333 0.9940765
## 2:(+DH+Si2)-2:(+DH-Si ) 0.125 -0.2103333 0.46033333 0.8018479
## 2:(+DH+Si1)-1:(+DH+Si1) -0.065 -0.4003333 0.27033333 0.9906270
## 1:(+DH+Si2)-1:(+DH+Si1) -0.125 -0.4603333 0.21033333 0.8018479
## 2:(+DH+Si2)-1:(+DH+Si1) 0.060 -0.2753333 0.39533333 0.9940765
## 1:(+DH+Si2)-2:(+DH+Si1) -0.060 -0.3953333 0.27533333 0.9940765
## 2:(+DH+Si2)-2:(+DH+Si1) 0.125 -0.2103333 0.46033333 0.8018479
## 2:(+DH+Si2)-1:(+DH+Si2) 0.185 -0.1503333 0.52033333 0.4430223#ggplot(datos_m1)+
# aes(TRATAMIENTO,SPAD, fill=TRATAMIENTO)+
# geom_boxplot()+
# facet_wrap(~BLOQUES)
#ggplot(datos_m1)+
# aes(x=TRATAMIENTO, y=SPAD, fill=TRATAMIENTO)+
# geom_bar(stat = "identity")+
# facet_wrap(~BLOQUES)
#mod1_SPAD=aov(SPAD~BLOQUES*TRATAMIENTO, data=datos_m1)
#anova(mod1_SPAD)
#TukeyHSD(mod1_SPAD)#INTERACCIÓN
ggplot(datos_m1)+
aes(TRATAMIENTO,PESO_S_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=PESO_S_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PSB=aov(PESO_S_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PSB)## Analysis of Variance Table
##
## Response: PESO_S_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.00001806 1.8063e-05 0.2542 0.62774
## TRATAMIENTO 3 0.00067119 2.2373e-04 3.1483 0.08646 .
## BLOQUES:TRATAMIENTO 3 0.00006369 2.1229e-05 0.2987 0.82555
## Residuals 8 0.00056850 7.1063e-05
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_PSB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.002125 -0.01184464 0.007594643 0.6277377
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.00500 -0.02408861 0.01408861 0.8349002
## (+DH+Si1)-(-DH-Si ) -0.00700 -0.02608861 0.01208861 0.6580175
## (+DH+Si2)-(-DH-Si ) 0.00975 -0.00933861 0.02883861 0.4128098
## (+DH+Si1)-(+DH-Si ) -0.00200 -0.02108861 0.01708861 0.9860299
## (+DH+Si2)-(+DH-Si ) 0.01475 -0.00433861 0.03383861 0.1392676
## (+DH+Si2)-(+DH+Si1) 0.01675 -0.00233861 0.03583861 0.0869898
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0005 -0.03285776 0.03385776 1.0000000
## 1:(+DH-Si )-1:(-DH-Si ) -0.0050 -0.03835776 0.02835776 0.9979469
## 2:(+DH-Si )-1:(-DH-Si ) -0.0045 -0.03785776 0.02885776 0.9989331
## 1:(+DH+Si1)-1:(-DH-Si ) -0.0065 -0.03985776 0.02685776 0.9903454
## 2:(+DH+Si1)-1:(-DH-Si ) -0.0070 -0.04035776 0.02635776 0.9854389
## 1:(+DH+Si2)-1:(-DH-Si ) 0.0145 -0.01885776 0.04785776 0.6787911
## 2:(+DH+Si2)-1:(-DH-Si ) 0.0055 -0.02785776 0.03885776 0.9963429
## 1:(+DH-Si )-2:(-DH-Si ) -0.0055 -0.03885776 0.02785776 0.9963429
## 2:(+DH-Si )-2:(-DH-Si ) -0.0050 -0.03835776 0.02835776 0.9979469
## 1:(+DH+Si1)-2:(-DH-Si ) -0.0070 -0.04035776 0.02635776 0.9854389
## 2:(+DH+Si1)-2:(-DH-Si ) -0.0075 -0.04085776 0.02585776 0.9789188
## 1:(+DH+Si2)-2:(-DH-Si ) 0.0140 -0.01935776 0.04735776 0.7097650
## 2:(+DH+Si2)-2:(-DH-Si ) 0.0050 -0.02835776 0.03835776 0.9979469
## 2:(+DH-Si )-1:(+DH-Si ) 0.0005 -0.03285776 0.03385776 1.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) -0.0015 -0.03485776 0.03185776 0.9999993
## 2:(+DH+Si1)-1:(+DH-Si ) -0.0020 -0.03535776 0.03135776 0.9999951
## 1:(+DH+Si2)-1:(+DH-Si ) 0.0195 -0.01385776 0.05285776 0.3851848
## 2:(+DH+Si2)-1:(+DH-Si ) 0.0105 -0.02285776 0.04385776 0.8960694
## 1:(+DH+Si1)-2:(+DH-Si ) -0.0020 -0.03535776 0.03135776 0.9999951
## 2:(+DH+Si1)-2:(+DH-Si ) -0.0025 -0.03585776 0.03085776 0.9999775
## 1:(+DH+Si2)-2:(+DH-Si ) 0.0190 -0.01435776 0.05235776 0.4109060
## 2:(+DH+Si2)-2:(+DH-Si ) 0.0100 -0.02335776 0.04335776 0.9155102
## 2:(+DH+Si1)-1:(+DH+Si1) -0.0005 -0.03385776 0.03285776 1.0000000
## 1:(+DH+Si2)-1:(+DH+Si1) 0.0210 -0.01235776 0.05435776 0.3148831
## 2:(+DH+Si2)-1:(+DH+Si1) 0.0120 -0.02135776 0.04535776 0.8253202
## 1:(+DH+Si2)-2:(+DH+Si1) 0.0215 -0.01185776 0.05485776 0.2937823
## 2:(+DH+Si2)-2:(+DH+Si1) 0.0125 -0.02085776 0.04585776 0.7982171
## 2:(+DH+Si2)-1:(+DH+Si2) -0.0090 -0.04235776 0.02435776 0.9475521ggplot(datos_m1)+
aes(TRATAMIENTO,PESO_F_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=PESO_F_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PFB=aov(PESO_F_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PFB)## Analysis of Variance Table
##
## Response: PESO_F_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.00003306 3.3062e-05 0.3251 0.5842
## TRATAMIENTO 3 0.00011319 3.7729e-05 0.3710 0.7763
## BLOQUES:TRATAMIENTO 3 0.00033369 1.1123e-04 1.0938 0.4060
## Residuals 8 0.00081350 1.0169e-04TukeyHSD(mod1_PFB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.002875 -0.0145019 0.008751898 0.5841904
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.00450 -0.02733431 0.01833431 0.9191270
## (+DH+Si1)-(-DH-Si ) 0.00275 -0.02008431 0.02558431 0.9791647
## (+DH+Si2)-(-DH-Si ) -0.00200 -0.02483431 0.02083431 0.9916978
## (+DH+Si1)-(+DH-Si ) 0.00725 -0.01558431 0.03008431 0.7448827
## (+DH+Si2)-(+DH-Si ) 0.00250 -0.02033431 0.02533431 0.9841426
## (+DH+Si2)-(+DH+Si1) -0.00475 -0.02758431 0.01808431 0.9069000
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0115 -0.02840345 0.05140345 0.9289842
## 1:(+DH-Si )-1:(-DH-Si ) 0.0020 -0.03790345 0.04190345 0.9999986
## 2:(+DH-Si )-1:(-DH-Si ) 0.0005 -0.03940345 0.04040345 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 0.0135 -0.02640345 0.05340345 0.8611893
## 2:(+DH+Si1)-1:(-DH-Si ) 0.0035 -0.03640345 0.04340345 0.9999353
## 1:(+DH+Si2)-1:(-DH-Si ) 0.0095 -0.03040345 0.04940345 0.9716313
## 2:(+DH+Si2)-1:(-DH-Si ) -0.0020 -0.04190345 0.03790345 0.9999986
## 1:(+DH-Si )-2:(-DH-Si ) -0.0095 -0.04940345 0.03040345 0.9716313
## 2:(+DH-Si )-2:(-DH-Si ) -0.0110 -0.05090345 0.02890345 0.9420096
## 1:(+DH+Si1)-2:(-DH-Si ) 0.0020 -0.03790345 0.04190345 0.9999986
## 2:(+DH+Si1)-2:(-DH-Si ) -0.0080 -0.04790345 0.03190345 0.9886768
## 1:(+DH+Si2)-2:(-DH-Si ) -0.0020 -0.04190345 0.03790345 0.9999986
## 2:(+DH+Si2)-2:(-DH-Si ) -0.0135 -0.05340345 0.02640345 0.8611893
## 2:(+DH-Si )-1:(+DH-Si ) -0.0015 -0.04140345 0.03840345 0.9999998
## 1:(+DH+Si1)-1:(+DH-Si ) 0.0115 -0.02840345 0.05140345 0.9289842
## 2:(+DH+Si1)-1:(+DH-Si ) 0.0015 -0.03840345 0.04140345 0.9999998
## 1:(+DH+Si2)-1:(+DH-Si ) 0.0075 -0.03240345 0.04740345 0.9921322
## 2:(+DH+Si2)-1:(+DH-Si ) -0.0040 -0.04390345 0.03590345 0.9998429
## 1:(+DH+Si1)-2:(+DH-Si ) 0.0130 -0.02690345 0.05290345 0.8803859
## 2:(+DH+Si1)-2:(+DH-Si ) 0.0030 -0.03690345 0.04290345 0.9999770
## 1:(+DH+Si2)-2:(+DH-Si ) 0.0090 -0.03090345 0.04890345 0.9785663
## 2:(+DH+Si2)-2:(+DH-Si ) -0.0025 -0.04240345 0.03740345 0.9999934
## 2:(+DH+Si1)-1:(+DH+Si1) -0.0100 -0.04990345 0.02990345 0.9632834
## 1:(+DH+Si2)-1:(+DH+Si1) -0.0040 -0.04390345 0.03590345 0.9998429
## 2:(+DH+Si2)-1:(+DH+Si1) -0.0155 -0.05540345 0.02440345 0.7722185
## 1:(+DH+Si2)-2:(+DH+Si1) 0.0060 -0.03390345 0.04590345 0.9979069
## 2:(+DH+Si2)-2:(+DH+Si1) -0.0055 -0.04540345 0.03440345 0.9987791
## 2:(+DH+Si2)-1:(+DH+Si2) -0.0115 -0.05140345 0.02840345 0.9289842ggplot(datos_m1)+
aes(TRATAMIENTO,PESO_SAT_H_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=PESO_SAT_H_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_PSHB=aov(PESO_SAT_H_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_PSHB)## Analysis of Variance Table
##
## Response: PESO_SAT_H_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.18297 0.182970 821.64 2.373e-09 ***
## TRATAMIENTO 3 0.54700 0.182334 818.79 2.729e-10 ***
## BLOQUES:TRATAMIENTO 3 0.57095 0.190318 854.64 2.301e-10 ***
## Residuals 8 0.00178 0.000223
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_PSHB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_SAT_H_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.213875 0.1966691 0.2310809 0
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.42650 0.39270893 0.46029107 0.0000000
## (+DH+Si1)-(-DH-Si ) 0.01150 -0.02229107 0.04529107 0.7050581
## (+DH+Si2)-(-DH-Si ) -0.01175 -0.04554107 0.02204107 0.6919431
## (+DH+Si1)-(+DH-Si ) -0.41500 -0.44879107 -0.38120893 0.0000000
## (+DH+Si2)-(+DH-Si ) -0.43825 -0.47204107 -0.40445893 0.0000000
## (+DH+Si2)-(+DH+Si1) -0.02325 -0.05704107 0.01054107 0.2018605
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0130 -0.04605064 0.07205064 0.9811449
## 1:(+DH-Si )-1:(-DH-Si ) -0.0010 -0.06005064 0.05805064 1.0000000
## 2:(+DH-Si )-1:(-DH-Si ) 0.8670 0.80794936 0.92605064 0.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 0.0235 -0.03555064 0.08255064 0.7535900
## 2:(+DH+Si1)-1:(-DH-Si ) 0.0125 -0.04655064 0.07155064 0.9847317
## 1:(+DH+Si2)-1:(-DH-Si ) 0.0020 -0.05705064 0.06105064 0.9999999
## 2:(+DH+Si2)-1:(-DH-Si ) -0.0125 -0.07155064 0.04655064 0.9847317
## 1:(+DH-Si )-2:(-DH-Si ) -0.0140 -0.07305064 0.04505064 0.9722292
## 2:(+DH-Si )-2:(-DH-Si ) 0.8540 0.79494936 0.91305064 0.0000000
## 1:(+DH+Si1)-2:(-DH-Si ) 0.0105 -0.04855064 0.06955064 0.9942879
## 2:(+DH+Si1)-2:(-DH-Si ) -0.0005 -0.05955064 0.05855064 1.0000000
## 1:(+DH+Si2)-2:(-DH-Si ) -0.0110 -0.07005064 0.04805064 0.9925247
## 2:(+DH+Si2)-2:(-DH-Si ) -0.0255 -0.08455064 0.03355064 0.6847091
## 2:(+DH-Si )-1:(+DH-Si ) 0.8680 0.80894936 0.92705064 0.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) 0.0245 -0.03455064 0.08355064 0.7195652
## 2:(+DH+Si1)-1:(+DH-Si ) 0.0135 -0.04555064 0.07255064 0.9769904
## 1:(+DH+Si2)-1:(+DH-Si ) 0.0030 -0.05605064 0.06205064 0.9999984
## 2:(+DH+Si2)-1:(+DH-Si ) -0.0115 -0.07055064 0.04755064 0.9903758
## 1:(+DH+Si1)-2:(+DH-Si ) -0.8435 -0.90255064 -0.78444936 0.0000000
## 2:(+DH+Si1)-2:(+DH-Si ) -0.8545 -0.91355064 -0.79544936 0.0000000
## 1:(+DH+Si2)-2:(+DH-Si ) -0.8650 -0.92405064 -0.80594936 0.0000000
## 2:(+DH+Si2)-2:(+DH-Si ) -0.8795 -0.93855064 -0.82044936 0.0000000
## 2:(+DH+Si1)-1:(+DH+Si1) -0.0110 -0.07005064 0.04805064 0.9925247
## 1:(+DH+Si2)-1:(+DH+Si1) -0.0215 -0.08055064 0.03755064 0.8175962
## 2:(+DH+Si2)-1:(+DH+Si1) -0.0360 -0.09505064 0.02305064 0.3446918
## 1:(+DH+Si2)-2:(+DH+Si1) -0.0105 -0.06955064 0.04855064 0.9942879
## 2:(+DH+Si2)-2:(+DH+Si1) -0.0250 -0.08405064 0.03405064 0.7022171
## 2:(+DH+Si2)-1:(+DH+Si2) -0.0145 -0.07355064 0.04455064 0.9668262ggplot(datos_m1)+
aes(TRATAMIENTO,N_ESTOMAS_ADAX, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ADAX, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_NEAx=aov(N_ESTOMAS_ADAX~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_NEAx)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ADAX
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 150.06 150.062 6.3687 0.035607 *
## TRATAMIENTO 3 323.19 107.729 4.5721 0.038030 *
## BLOQUES:TRATAMIENTO 3 549.69 183.229 7.7763 0.009331 **
## Residuals 8 188.50 23.563
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_NEAx)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ADAX ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -6.125 -11.72181 -0.5281873 0.035607
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 10.00 -0.9916972 20.991697 0.0751688
## (+DH+Si1)-(-DH-Si ) 0.25 -10.7416972 11.241697 0.9998491
## (+DH+Si2)-(-DH-Si ) 8.00 -2.9916972 18.991697 0.1698119
## (+DH+Si1)-(+DH-Si ) -9.75 -20.7416972 1.241697 0.0833150
## (+DH+Si2)-(+DH-Si ) -2.00 -12.9916972 8.991697 0.9345413
## (+DH+Si2)-(+DH+Si1) 7.75 -3.2416972 18.741697 0.1875296
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -6.5 -25.7082316 12.708232 0.8610623
## 1:(+DH-Si )-1:(-DH-Si ) 19.0 -0.2082316 38.208232 0.0527955
## 2:(+DH-Si )-1:(-DH-Si ) -5.5 -24.7082316 13.708232 0.9310196
## 1:(+DH+Si1)-1:(-DH-Si ) -2.5 -21.7082316 16.708232 0.9991490
## 2:(+DH+Si1)-1:(-DH-Si ) -3.5 -22.7082316 15.708232 0.9934169
## 1:(+DH+Si2)-1:(-DH-Si ) 1.0 -18.2082316 20.208232 0.9999981
## 2:(+DH+Si2)-1:(-DH-Si ) 8.5 -10.7082316 27.708232 0.6624576
## 1:(+DH-Si )-2:(-DH-Si ) 25.5 6.2917684 44.708232 0.0103690
## 2:(+DH-Si )-2:(-DH-Si ) 1.0 -18.2082316 20.208232 0.9999981
## 1:(+DH+Si1)-2:(-DH-Si ) 4.0 -15.2082316 23.208232 0.9860360
## 2:(+DH+Si1)-2:(-DH-Si ) 3.0 -16.2082316 22.208232 0.9973612
## 1:(+DH+Si2)-2:(-DH-Si ) 7.5 -11.7082316 26.708232 0.7683039
## 2:(+DH+Si2)-2:(-DH-Si ) 15.0 -4.2082316 34.208232 0.1512317
## 2:(+DH-Si )-1:(+DH-Si ) -24.5 -43.7082316 -5.291768 0.0131714
## 1:(+DH+Si1)-1:(+DH-Si ) -21.5 -40.7082316 -2.291768 0.0277021
## 2:(+DH+Si1)-1:(+DH-Si ) -22.5 -41.7082316 -3.291768 0.0215330
## 1:(+DH+Si2)-1:(+DH-Si ) -18.0 -37.2082316 1.208232 0.0686418
## 2:(+DH+Si2)-1:(+DH-Si ) -10.5 -29.7082316 8.708232 0.4523499
## 1:(+DH+Si1)-2:(+DH-Si ) 3.0 -16.2082316 22.208232 0.9973612
## 2:(+DH+Si1)-2:(+DH-Si ) 2.0 -17.2082316 21.208232 0.9997984
## 1:(+DH+Si2)-2:(+DH-Si ) 6.5 -12.7082316 25.708232 0.8610623
## 2:(+DH+Si2)-2:(+DH-Si ) 14.0 -5.2082316 33.208232 0.1959070
## 2:(+DH+Si1)-1:(+DH+Si1) -1.0 -20.2082316 18.208232 0.9999981
## 1:(+DH+Si2)-1:(+DH+Si1) 3.5 -15.7082316 22.708232 0.9934169
## 2:(+DH+Si2)-1:(+DH+Si1) 11.0 -8.2082316 30.208232 0.4055156
## 1:(+DH+Si2)-2:(+DH+Si1) 4.5 -14.7082316 23.708232 0.9738789
## 2:(+DH+Si2)-2:(+DH+Si1) 12.0 -7.2082316 31.208232 0.3218966
## 2:(+DH+Si2)-1:(+DH+Si2) 7.5 -11.7082316 26.708232 0.7683039ggplot(datos_m1)+
aes(TRATAMIENTO,N_ESTOMAS_ABA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ABA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_NEAbx=aov(N_ESTOMAS_ABA~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_NEAbx)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ABA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 14.06 14.062 0.1251 0.7327
## TRATAMIENTO 3 373.19 124.396 1.1064 0.4015
## BLOQUES:TRATAMIENTO 3 809.69 269.896 2.4004 0.1433
## Residuals 8 899.50 112.438TukeyHSD(mod1_NEAbx)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ABA ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 1.875 -10.35104 14.10104 0.732739
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.75 -24.76097 23.26097 0.9996100
## (+DH+Si1)-(-DH-Si ) -5.75 -29.76097 18.26097 0.8671147
## (+DH+Si2)-(-DH-Si ) 7.75 -16.26097 31.76097 0.7357988
## (+DH+Si1)-(+DH-Si ) -5.00 -29.01097 19.01097 0.9066478
## (+DH+Si2)-(+DH-Si ) 8.50 -15.51097 32.51097 0.6807564
## (+DH+Si2)-(+DH+Si1) 13.50 -10.51097 37.51097 0.3394796
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -7.0 -48.95969 34.95969 0.9960775
## 1:(+DH-Si )-1:(-DH-Si ) -1.0 -42.95969 40.95969 1.0000000
## 2:(+DH-Si )-1:(-DH-Si ) -7.5 -49.45969 34.45969 0.9941117
## 1:(+DH+Si1)-1:(-DH-Si ) -6.5 -48.45969 35.45969 0.9974891
## 2:(+DH+Si1)-1:(-DH-Si ) -12.0 -53.95969 29.95969 0.9313952
## 1:(+DH+Si2)-1:(-DH-Si ) -9.0 -50.95969 32.95969 0.9836022
## 2:(+DH+Si2)-1:(-DH-Si ) 17.5 -24.45969 59.45969 0.7151395
## 1:(+DH-Si )-2:(-DH-Si ) 6.0 -35.95969 47.95969 0.9984647
## 2:(+DH-Si )-2:(-DH-Si ) -0.5 -42.45969 41.45969 1.0000000
## 1:(+DH+Si1)-2:(-DH-Si ) 0.5 -41.45969 42.45969 1.0000000
## 2:(+DH+Si1)-2:(-DH-Si ) -5.0 -46.95969 36.95969 0.9995154
## 1:(+DH+Si2)-2:(-DH-Si ) -2.0 -43.95969 39.95969 0.9999990
## 2:(+DH+Si2)-2:(-DH-Si ) 24.5 -17.45969 66.45969 0.3863239
## 2:(+DH-Si )-1:(+DH-Si ) -6.5 -48.45969 35.45969 0.9974891
## 1:(+DH+Si1)-1:(+DH-Si ) -5.5 -47.45969 36.45969 0.9991099
## 2:(+DH+Si1)-1:(+DH-Si ) -11.0 -52.95969 30.95969 0.9542564
## 1:(+DH+Si2)-1:(+DH-Si ) -8.0 -49.95969 33.95969 0.9914647
## 2:(+DH+Si2)-1:(+DH-Si ) 18.5 -23.45969 60.45969 0.6658375
## 1:(+DH+Si1)-2:(+DH-Si ) 1.0 -40.95969 42.95969 1.0000000
## 2:(+DH+Si1)-2:(+DH-Si ) -4.5 -46.45969 37.45969 0.9997553
## 1:(+DH+Si2)-2:(+DH-Si ) -1.5 -43.45969 40.45969 0.9999999
## 2:(+DH+Si2)-2:(+DH-Si ) 25.0 -16.95969 66.95969 0.3666374
## 2:(+DH+Si1)-1:(+DH+Si1) -5.5 -47.45969 36.45969 0.9991099
## 1:(+DH+Si2)-1:(+DH+Si1) -2.5 -44.45969 39.45969 0.9999953
## 2:(+DH+Si2)-1:(+DH+Si1) 24.0 -17.95969 65.95969 0.4067222
## 1:(+DH+Si2)-2:(+DH+Si1) 3.0 -38.95969 44.95969 0.9999837
## 2:(+DH+Si2)-2:(+DH+Si1) 29.5 -12.45969 71.45969 0.2222621
## 2:(+DH+Si2)-1:(+DH+Si2) 26.5 -15.45969 68.45969 0.3119715ggplot(datos_m1)+
aes(TRATAMIENTO,N_ESTOMAS_ABIERTOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ABIERTOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_NEA=aov(N_ESTOMAS_ABIERTOS~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_NEA)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ABIERTOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 20.25 20.250 0.5684 0.4725
## TRATAMIENTO 3 96.25 32.083 0.9006 0.4821
## BLOQUES:TRATAMIENTO 3 68.25 22.750 0.6386 0.6111
## Residuals 8 285.00 35.625TukeyHSD(mod1_NEA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ABIERTOS ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 2.25 -4.631887 9.131887 0.4724919
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -4.00 -17.51548 9.515481 0.7812175
## (+DH+Si1)-(-DH-Si ) -6.75 -20.26548 6.765481 0.4302528
## (+DH+Si2)-(-DH-Si ) -4.75 -18.26548 8.765481 0.6853149
## (+DH+Si1)-(+DH-Si ) -2.75 -16.26548 10.765481 0.9120877
## (+DH+Si2)-(+DH-Si ) -0.75 -14.26548 12.765481 0.9978374
## (+DH+Si2)-(+DH+Si1) 2.00 -11.51548 15.515481 0.9627721
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 5.500000e+00 -18.1186 29.1186 0.9746788
## 1:(+DH-Si )-1:(-DH-Si ) -1.000000e+00 -24.6186 22.6186 0.9999996
## 2:(+DH-Si )-1:(-DH-Si ) -1.500000e+00 -25.1186 22.1186 0.9999927
## 1:(+DH+Si1)-1:(-DH-Si ) -2.500000e+00 -26.1186 21.1186 0.9997754
## 2:(+DH+Si1)-1:(-DH-Si ) -5.500000e+00 -29.1186 18.1186 0.9746788
## 1:(+DH+Si2)-1:(-DH-Si ) -5.500000e+00 -29.1186 18.1186 0.9746788
## 2:(+DH+Si2)-1:(-DH-Si ) 1.500000e+00 -22.1186 25.1186 0.9999927
## 1:(+DH-Si )-2:(-DH-Si ) -6.500000e+00 -30.1186 17.1186 0.9424557
## 2:(+DH-Si )-2:(-DH-Si ) -7.000000e+00 -30.6186 16.6186 0.9196109
## 1:(+DH+Si1)-2:(-DH-Si ) -8.000000e+00 -31.6186 15.6186 0.8605545
## 2:(+DH+Si1)-2:(-DH-Si ) -1.100000e+01 -34.6186 12.6186 0.6139182
## 1:(+DH+Si2)-2:(-DH-Si ) -1.100000e+01 -34.6186 12.6186 0.6139182
## 2:(+DH+Si2)-2:(-DH-Si ) -4.000000e+00 -27.6186 19.6186 0.9957105
## 2:(+DH-Si )-1:(+DH-Si ) -5.000000e-01 -24.1186 23.1186 1.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) -1.500000e+00 -25.1186 22.1186 0.9999927
## 2:(+DH+Si1)-1:(+DH-Si ) -4.500000e+00 -28.1186 19.1186 0.9914980
## 1:(+DH+Si2)-1:(+DH-Si ) -4.500000e+00 -28.1186 19.1186 0.9914980
## 2:(+DH+Si2)-1:(+DH-Si ) 2.500000e+00 -21.1186 26.1186 0.9997754
## 1:(+DH+Si1)-2:(+DH-Si ) -1.000000e+00 -24.6186 22.6186 0.9999996
## 2:(+DH+Si1)-2:(+DH-Si ) -4.000000e+00 -27.6186 19.6186 0.9957105
## 1:(+DH+Si2)-2:(+DH-Si ) -4.000000e+00 -27.6186 19.6186 0.9957105
## 2:(+DH+Si2)-2:(+DH-Si ) 3.000000e+00 -20.6186 26.6186 0.9992709
## 2:(+DH+Si1)-1:(+DH+Si1) -3.000000e+00 -26.6186 20.6186 0.9992709
## 1:(+DH+Si2)-1:(+DH+Si1) -3.000000e+00 -26.6186 20.6186 0.9992709
## 2:(+DH+Si2)-1:(+DH+Si1) 4.000000e+00 -19.6186 27.6186 0.9957105
## 1:(+DH+Si2)-2:(+DH+Si1) -3.552714e-15 -23.6186 23.6186 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) 7.000000e+00 -16.6186 30.6186 0.9196109
## 2:(+DH+Si2)-1:(+DH+Si2) 7.000000e+00 -16.6186 30.6186 0.9196109ggplot(datos_m1)+
aes(TRATAMIENTO,N_ESTOMAS_P_ABIERTOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_P_ABIERTOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_NEPA=aov(N_ESTOMAS_P_ABIERTOS~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_NEPA)## Analysis of Variance Table
##
## Response: N_ESTOMAS_P_ABIERTOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 4.00 4.000 0.0901 0.77165
## TRATAMIENTO 3 368.75 122.917 2.7700 0.11081
## BLOQUES:TRATAMIENTO 3 914.00 304.667 6.8657 0.01327 *
## Residuals 8 355.00 44.375
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_NEPA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_P_ABIERTOS ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 1 -6.680673 8.680673 0.7716519
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.0 -15.084234 15.08423 1.0000000
## (+DH+Si1)-(-DH-Si ) -2.5 -17.584234 12.58423 0.9491346
## (+DH+Si2)-(-DH-Si ) 10.0 -5.084234 25.08423 0.2247525
## (+DH+Si1)-(+DH-Si ) -2.5 -17.584234 12.58423 0.9491346
## (+DH+Si2)-(+DH-Si ) 10.0 -5.084234 25.08423 0.2247525
## (+DH+Si2)-(+DH+Si1) 12.5 -2.584234 27.58423 0.1084117
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -11.0 -37.360029 15.36003 0.7146662
## 1:(+DH-Si )-1:(-DH-Si ) 0.0 -26.360029 26.36003 1.0000000
## 2:(+DH-Si )-1:(-DH-Si ) -11.0 -37.360029 15.36003 0.7146662
## 1:(+DH+Si1)-1:(-DH-Si ) -8.0 -34.360029 18.36003 0.9109075
## 2:(+DH+Si1)-1:(-DH-Si ) -8.0 -34.360029 18.36003 0.9109075
## 1:(+DH+Si2)-1:(-DH-Si ) -8.5 -34.860029 17.86003 0.8852451
## 2:(+DH+Si2)-1:(-DH-Si ) 17.5 -8.860029 43.86003 0.2682442
## 1:(+DH-Si )-2:(-DH-Si ) 11.0 -15.360029 37.36003 0.7146662
## 2:(+DH-Si )-2:(-DH-Si ) 0.0 -26.360029 26.36003 1.0000000
## 1:(+DH+Si1)-2:(-DH-Si ) 3.0 -23.360029 29.36003 0.9996398
## 2:(+DH+Si1)-2:(-DH-Si ) 3.0 -23.360029 29.36003 0.9996398
## 1:(+DH+Si2)-2:(-DH-Si ) 2.5 -23.860029 28.86003 0.9998911
## 2:(+DH+Si2)-2:(-DH-Si ) 28.5 2.139971 54.86003 0.0333942
## 2:(+DH-Si )-1:(+DH-Si ) -11.0 -37.360029 15.36003 0.7146662
## 1:(+DH+Si1)-1:(+DH-Si ) -8.0 -34.360029 18.36003 0.9109075
## 2:(+DH+Si1)-1:(+DH-Si ) -8.0 -34.360029 18.36003 0.9109075
## 1:(+DH+Si2)-1:(+DH-Si ) -8.5 -34.860029 17.86003 0.8852451
## 2:(+DH+Si2)-1:(+DH-Si ) 17.5 -8.860029 43.86003 0.2682442
## 1:(+DH+Si1)-2:(+DH-Si ) 3.0 -23.360029 29.36003 0.9996398
## 2:(+DH+Si1)-2:(+DH-Si ) 3.0 -23.360029 29.36003 0.9996398
## 1:(+DH+Si2)-2:(+DH-Si ) 2.5 -23.860029 28.86003 0.9998911
## 2:(+DH+Si2)-2:(+DH-Si ) 28.5 2.139971 54.86003 0.0333942
## 2:(+DH+Si1)-1:(+DH+Si1) 0.0 -26.360029 26.36003 1.0000000
## 1:(+DH+Si2)-1:(+DH+Si1) -0.5 -26.860029 25.86003 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si1) 25.5 -0.860029 51.86003 0.0589116
## 1:(+DH+Si2)-2:(+DH+Si1) -0.5 -26.860029 25.86003 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) 25.5 -0.860029 51.86003 0.0589116
## 2:(+DH+Si2)-1:(+DH+Si2) 26.0 -0.360029 52.36003 0.0535481ggplot(datos_m1)+
aes(TRATAMIENTO,N_ESTOMAS_CERRADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_CERRADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_NEC=aov(N_ESTOMAS_CERRADOS~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_NEC)## Analysis of Variance Table
##
## Response: N_ESTOMAS_CERRADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 225.00 225.000 4.5113 0.06641 .
## TRATAMIENTO 3 443.25 147.750 2.9624 0.09749 .
## BLOQUES:TRATAMIENTO 3 232.50 77.500 1.5539 0.27441
## Residuals 8 399.00 49.875
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_NEC)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_CERRADOS ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -7.5 -15.64276 0.6427583 0.0664128
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 13.25 -2.741733 29.241733 0.1084735
## (+DH+Si1)-(-DH-Si ) 3.75 -12.241733 19.741733 0.8738465
## (+DH+Si2)-(-DH-Si ) 10.50 -5.491733 26.491733 0.2308907
## (+DH+Si1)-(+DH-Si ) -9.50 -25.491733 6.491733 0.2991214
## (+DH+Si2)-(+DH-Si ) -2.75 -18.741733 13.241733 0.9437719
## (+DH+Si2)-(+DH+Si1) 6.75 -9.241733 22.741733 0.5591889
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -8.0 -35.9459035 19.945903 0.9310936
## 1:(+DH-Si )-1:(-DH-Si ) 19.0 -8.9459035 46.945903 0.2485588
## 2:(+DH-Si )-1:(-DH-Si ) -0.5 -28.4459035 27.445903 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 1.5 -26.4459035 29.445903 0.9999977
## 2:(+DH+Si1)-1:(-DH-Si ) -2.0 -29.9459035 25.945903 0.9999836
## 1:(+DH+Si2)-1:(-DH-Si ) 6.0 -21.9459035 33.945903 0.9835158
## 2:(+DH+Si2)-1:(-DH-Si ) 7.0 -20.9459035 34.945903 0.9633710
## 1:(+DH-Si )-2:(-DH-Si ) 27.0 -0.9459035 54.945903 0.0592754
## 2:(+DH-Si )-2:(-DH-Si ) 7.5 -20.4459035 35.445903 0.9488478
## 1:(+DH+Si1)-2:(-DH-Si ) 9.5 -18.4459035 37.445903 0.8585973
## 2:(+DH+Si1)-2:(-DH-Si ) 6.0 -21.9459035 33.945903 0.9835158
## 1:(+DH+Si2)-2:(-DH-Si ) 14.0 -13.9459035 41.945903 0.5412897
## 2:(+DH+Si2)-2:(-DH-Si ) 15.0 -12.9459035 42.945903 0.4709373
## 2:(+DH-Si )-1:(+DH-Si ) -19.5 -47.4459035 8.445903 0.2280289
## 1:(+DH+Si1)-1:(+DH-Si ) -17.5 -45.4459035 10.445903 0.3197275
## 2:(+DH+Si1)-1:(+DH-Si ) -21.0 -48.9459035 6.945903 0.1751899
## 1:(+DH+Si2)-1:(+DH-Si ) -13.0 -40.9459035 14.945903 0.6150668
## 2:(+DH+Si2)-1:(+DH-Si ) -12.0 -39.9459035 15.945903 0.6897508
## 1:(+DH+Si1)-2:(+DH-Si ) 2.0 -25.9459035 29.945903 0.9999836
## 2:(+DH+Si1)-2:(+DH-Si ) -1.5 -29.4459035 26.445903 0.9999977
## 1:(+DH+Si2)-2:(+DH-Si ) 6.5 -21.4459035 34.445903 0.9748331
## 2:(+DH+Si2)-2:(+DH-Si ) 7.5 -20.4459035 35.445903 0.9488478
## 2:(+DH+Si1)-1:(+DH+Si1) -3.5 -31.4459035 24.445903 0.9993335
## 1:(+DH+Si2)-1:(+DH+Si1) 4.5 -23.4459035 32.445903 0.9968264
## 2:(+DH+Si2)-1:(+DH+Si1) 5.5 -22.4459035 33.445903 0.9897890
## 1:(+DH+Si2)-2:(+DH+Si1) 8.0 -19.9459035 35.945903 0.9310936
## 2:(+DH+Si2)-2:(+DH+Si1) 9.0 -18.9459035 36.945903 0.8858324
## 2:(+DH+Si2)-1:(+DH+Si2) 1.0 -26.9459035 28.945903 0.9999999ggplot(datos_m1)+
aes(TRATAMIENTO,CRA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=CRA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_CRA=aov(CRA~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_CRA)## Analysis of Variance Table
##
## Response: CRA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 911.5 911.52 14.033 0.0056560 **
## TRATAMIENTO 3 2865.7 955.24 14.706 0.0012794 **
## BLOQUES:TRATAMIENTO 3 3703.1 1234.37 19.003 0.0005358 ***
## Residuals 8 519.6 64.96
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_CRA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = CRA ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -15.09567 -24.38827 -5.803057 0.0056561
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -30.772647 -49.022594 -12.52270 0.0028695
## (+DH+Si1)-(-DH-Si ) -3.507887 -21.757834 14.74206 0.9242692
## (+DH+Si2)-(-DH-Si ) 2.662644 -15.587303 20.91259 0.9642098
## (+DH+Si1)-(+DH-Si ) 27.264760 9.014813 45.51471 0.0060339
## (+DH+Si2)-(+DH-Si ) 33.435291 15.185344 51.68524 0.0016844
## (+DH+Si2)-(+DH+Si1) 6.170531 -12.079416 24.42048 0.7089673
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 3.033773 -28.85841 34.92596 0.9998889
## 1:(+DH-Si )-1:(-DH-Si ) 4.522667 -27.36952 36.41485 0.9985419
## 2:(+DH-Si )-1:(-DH-Si ) -63.034188 -94.92637 -31.14201 0.0007591
## 1:(+DH+Si1)-1:(-DH-Si ) -0.997151 -32.88933 30.89503 0.9999999
## 2:(+DH+Si1)-1:(-DH-Si ) -2.984850 -34.87703 28.90733 0.9999003
## 1:(+DH+Si2)-1:(-DH-Si ) 1.115472 -30.77671 33.00765 0.9999999
## 2:(+DH+Si2)-1:(-DH-Si ) 7.243590 -24.64859 39.13577 0.9777644
## 1:(+DH-Si )-2:(-DH-Si ) 1.488894 -30.40329 33.38108 0.9999991
## 2:(+DH-Si )-2:(-DH-Si ) -66.067961 -97.96014 -34.17578 0.0005450
## 1:(+DH+Si1)-2:(-DH-Si ) -4.030924 -35.92311 27.86126 0.9992935
## 2:(+DH+Si1)-2:(-DH-Si ) -6.018623 -37.91081 25.87356 0.9919478
## 1:(+DH+Si2)-2:(-DH-Si ) -1.918301 -33.81048 29.97388 0.9999950
## 2:(+DH+Si2)-2:(-DH-Si ) 4.209817 -27.68237 36.10200 0.9990696
## 2:(+DH-Si )-1:(+DH-Si ) -67.556855 -99.44904 -35.66467 0.0004652
## 1:(+DH+Si1)-1:(+DH-Si ) -5.519818 -37.41200 26.37236 0.9951234
## 2:(+DH+Si1)-1:(+DH-Si ) -7.507517 -39.39970 24.38467 0.9732248
## 1:(+DH+Si2)-1:(+DH-Si ) -3.407195 -35.29938 28.48499 0.9997614
## 2:(+DH+Si2)-1:(+DH-Si ) 2.720923 -29.17126 34.61311 0.9999462
## 1:(+DH+Si1)-2:(+DH-Si ) 62.037037 30.14485 93.92922 0.0008486
## 2:(+DH+Si1)-2:(+DH-Si ) 60.049338 28.15716 91.94152 0.0010640
## 1:(+DH+Si2)-2:(+DH-Si ) 64.149660 32.25748 96.04184 0.0006711
## 2:(+DH+Si2)-2:(+DH-Si ) 70.277778 38.38560 102.16996 0.0003508
## 2:(+DH+Si1)-1:(+DH+Si1) -1.987699 -33.87988 29.90448 0.9999936
## 1:(+DH+Si2)-1:(+DH+Si1) 2.112623 -29.77956 34.00481 0.9999903
## 2:(+DH+Si2)-1:(+DH+Si1) 8.240741 -23.65144 40.13292 0.9573407
## 1:(+DH+Si2)-2:(+DH+Si1) 4.100322 -27.79186 35.99250 0.9992125
## 2:(+DH+Si2)-2:(+DH+Si1) 10.228440 -21.66374 42.12062 0.8877420
## 2:(+DH+Si2)-1:(+DH+Si2) 6.128118 -25.76406 38.02030 0.9910787ggplot(datos_m1)+
aes(TRATAMIENTO,AREA_DOSEL, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m1)+
aes(x=TRATAMIENTO, y=TEMPERATURA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod1_ABD=aov(AREA_DOSEL~BLOQUES*TRATAMIENTO, data=datos_m1)
anova(mod1_ABD)## Analysis of Variance Table
##
## Response: AREA_DOSEL
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 4363.6 4363.6 14.4179 0.005258 **
## TRATAMIENTO 3 4001.9 1334.0 4.4076 0.041486 *
## BLOQUES:TRATAMIENTO 3 2284.4 761.5 2.5160 0.131993
## Residuals 8 2421.2 302.7
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod1_ABD)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AREA_DOSEL ~ BLOQUES * TRATAMIENTO, data = datos_m1)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 33.02875 12.9701 53.0874 0.0052584
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -35.5325 -74.92611 3.861106 0.0778659
## (+DH+Si1)-(-DH-Si ) -41.2650 -80.65861 -1.871394 0.0403890
## (+DH+Si2)-(-DH-Si ) -26.7950 -66.18861 12.598606 0.2088016
## (+DH+Si1)-(+DH-Si ) -5.7325 -45.12611 33.661106 0.9644690
## (+DH+Si2)-(+DH-Si ) 8.7375 -30.65611 48.131106 0.8902920
## (+DH+Si2)-(+DH+Si1) 14.4700 -24.92361 53.863606 0.6569276
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 67.065 -1.776189 135.906189 0.0569208
## 1:(+DH-Si )-1:(-DH-Si ) -15.115 -83.956189 53.726189 0.9814112
## 2:(+DH-Si )-1:(-DH-Si ) 11.115 -57.726189 79.956189 0.9967747
## 1:(+DH+Si1)-1:(-DH-Si ) -8.025 -76.866189 60.816189 0.9995794
## 2:(+DH+Si1)-1:(-DH-Si ) -7.440 -76.281189 61.401189 0.9997427
## 1:(+DH+Si2)-1:(-DH-Si ) -12.380 -81.221189 56.461189 0.9938996
## 2:(+DH+Si2)-1:(-DH-Si ) 25.855 -42.986189 94.696189 0.7966381
## 1:(+DH-Si )-2:(-DH-Si ) -82.180 -151.021189 -13.338811 0.0193456
## 2:(+DH-Si )-2:(-DH-Si ) -55.950 -124.791189 12.891189 0.1288609
## 1:(+DH+Si1)-2:(-DH-Si ) -75.090 -143.931189 -6.248811 0.0318554
## 2:(+DH+Si1)-2:(-DH-Si ) -74.505 -143.346189 -5.663811 0.0332152
## 1:(+DH+Si2)-2:(-DH-Si ) -79.445 -148.286189 -10.603811 0.0234075
## 2:(+DH+Si2)-2:(-DH-Si ) -41.210 -110.051189 27.631189 0.3620938
## 2:(+DH-Si )-1:(+DH-Si ) 26.230 -42.611189 95.071189 0.7864088
## 1:(+DH+Si1)-1:(+DH-Si ) 7.090 -61.751189 75.931189 0.9998123
## 2:(+DH+Si1)-1:(+DH-Si ) 7.675 -61.166189 76.516189 0.9996849
## 1:(+DH+Si2)-1:(+DH-Si ) 2.735 -66.106189 71.576189 0.9999997
## 2:(+DH+Si2)-1:(+DH-Si ) 40.970 -27.871189 109.811189 0.3677284
## 1:(+DH+Si1)-2:(+DH-Si ) -19.140 -87.981189 49.701189 0.9396737
## 2:(+DH+Si1)-2:(+DH-Si ) -18.555 -87.396189 50.286189 0.9477979
## 1:(+DH+Si2)-2:(+DH-Si ) -23.495 -92.336189 45.346189 0.8564275
## 2:(+DH+Si2)-2:(+DH-Si ) 14.740 -54.101189 83.581189 0.9837566
## 2:(+DH+Si1)-1:(+DH+Si1) 0.585 -68.256189 69.426189 1.0000000
## 1:(+DH+Si2)-1:(+DH+Si1) -4.355 -73.196189 64.486189 0.9999929
## 2:(+DH+Si2)-1:(+DH+Si1) 33.880 -34.961189 102.721189 0.5592419
## 1:(+DH+Si2)-2:(+DH+Si1) -4.940 -73.781189 63.901189 0.9999833
## 2:(+DH+Si2)-2:(+DH+Si1) 33.295 -35.546189 102.136189 0.5767063
## 2:(+DH+Si2)-1:(+DH+Si2) 38.235 -30.606189 107.076189 0.4362251#Muetreo 2 analisis de datos
#install.packages("plotly")
#library(car)
#library(ggplot2)
#library(agricolae)
#library(plotly)
#library(dplyr)
#library(tidyverse)
#library(collapsibleTree)
#library(readxl)
#library(multcomp)
datos_m <- read_excel("C:/Users/JuanSebH2/Downloads/Universidad/Fisiologia/FisioVeg/DocumentoFinal_muestreos.xlsx",
sheet = "M2")
datos_m=data.frame(datos_m)
datos_m2=datos_m
datos_m2$BLOQUES=as.character(datos_m2$BLOQUES)
collapsibleTree(datos_m2, hierarchy = c("MUESTREO", "BLOQUES", "TRATAMIENTO"), collapsed = F)datos_m2## MUESTREO TRATAMIENTO BLOQUES TEMPERATURA N_HOJAS AREA_DOSEL LONGITUD_P_AEREA
## 1 2 (-DH-Si ) 1 20.42 6 163.38 11.0
## 2 2 (-DH-Si ) 1 20.88 7 208.79 12.0
## 3 2 (-DH-Si ) 2 20.56 6 174.16 14.0
## 4 2 (-DH-Si ) 2 17.84 5 114.81 13.0
## 5 2 (+DH-Si ) 1 22.04 5 72.22 10.0
## 6 2 (+DH-Si ) 1 20.66 7 125.67 11.0
## 7 2 (+DH-Si ) 2 20.82 6 97.77 12.0
## 8 2 (+DH-Si ) 2 20.60 6 100.71 10.0
## 9 2 (+DH+Si1) 1 21.08 6 132.76 13.3
## 10 2 (+DH+Si1) 1 20.88 6 106.49 12.5
## 11 2 (+DH+Si1) 2 21.38 6 103.79 12.5
## 12 2 (+DH+Si1) 2 20.82 7 140.38 11.0
## 13 2 (+DH+Si2) 1 19.98 6 107.21 11.0
## 14 2 (+DH+Si2) 1 21.22 6 84.88 12.5
## 15 2 (+DH+Si2) 2 21.26 6 108.69 13.5
## 16 2 (+DH+Si2) 2 20.68 7 165.57 9.5
## AREA_FOLIAR LONGITUD_RAIZ DIAMETRO_RAIZ PESO_F_RAIZ PESO_F_AEREA PESO_S_RAIZ
## 1 156.5 5.0 34.92 27.93 5.83 0.821
## 2 90.9 3.5 28.45 14.55 3.53 0.373
## 3 109.2 3.5 28.88 14.77 4.03 0.566
## 4 161.3 3.0 34.27 18.77 5.87 0.656
## 5 113.8 3.5 28.21 12.29 3.78 0.640
## 6 106.1 3.0 23.21 8.60 3.99 0.482
## 7 97.1 3.0 25.59 10.20 3.55 0.352
## 8 100.2 3.0 27.36 8.54 4.24 0.325
## 9 111.2 3.0 25.97 10.27 4.79 0.920
## 10 125.2 2.5 21.73 4.63 5.64 0.409
## 11 116.3 4.0 23.24 10.18 4.50 0.670
## 12 93.2 3.0 23.04 8.83 3.24 0.419
## 13 100.0 3.5 27.65 13.52 3.80 0.638
## 14 104.2 3.0 22.14 7.80 3.56 0.365
## 15 115.3 2.5 21.89 7.34 4.64 0.417
## 16 104.8 3.0 28.05 11.07 3.65 0.627
## PESO_S_AEREA SPAD PESO_F_BOCADOS PESO_SAT_H_BOCADOS PESO_S_BOCADOS CRA
## 1 0.517 39.10 0.145 0.210 0.0125 67.08861
## 2 0.286 37.50 0.129 0.220 0.0099 56.68729
## 3 0.368 33.20 0.134 0.245 0.0120 52.36052
## 4 0.540 32.60 0.137 0.226 0.0108 58.64312
## 5 0.400 43.60 0.140 0.185 0.0117 74.03347
## 6 0.490 38.40 0.135 0.229 0.0142 56.23836
## 7 0.320 37.40 0.149 0.207 0.0112 70.37794
## 8 0.377 39.20 0.158 0.220 0.0107 70.37745
## 9 0.636 43.30 0.155 0.225 0.0194 65.95331
## 10 0.555 53.08 0.168 0.210 0.0174 78.19315
## 11 0.570 49.00 0.135 0.216 0.0179 59.11156
## 12 0.339 40.70 0.136 0.176 0.0153 75.10890
## 13 0.444 52.94 0.123 0.220 0.0145 52.79805
## 14 0.454 48.06 0.135 0.258 0.0161 49.15254
## 15 0.507 46.10 0.112 0.165 0.0146 64.76064
## 16 0.388 46.40 0.125 0.183 0.0133 65.82204
## N_ESTOMAS_ADAX N_ESTOMAS_ABA N_ESTOMAS_ABIERTOS N_ESTOMAS_P_ABIERTOS
## 1 31 46 8 26
## 2 31 29 5 36
## 3 16 37 12 27
## 4 40 52 14 47
## 5 16 40 2 12
## 6 33 27 5 24
## 7 14 30 4 24
## 8 16 32 8 23
## 9 60 61 21 70
## 10 58 65 19 80
## 11 51 78 34 65
## 12 46 62 18 70
## 13 20 42 7 22
## 14 26 41 4 23
## 15 36 31 9 28
## 16 20 28 16 17
## N_ESTOMAS_CERRADOS
## 1 43
## 2 19
## 3 14
## 4 31
## 5 42
## 6 32
## 7 16
## 8 17
## 9 30
## 10 24
## 11 30
## 12 20
## 13 33
## 14 34
## 15 30
## 16 15#Analisis descriptivo Y ESTADISTICO
ggplot(datos_m2)+
aes(TRATAMIENTO,TEMPERATURA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=TEMPERATURA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_T=aov(TEMPERATURA~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_T)## Analysis of Variance Table
##
## Response: TEMPERATURA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.6400 0.64000 0.8685 0.3786
## TRATAMIENTO 3 3.3290 1.10967 1.5059 0.2855
## BLOQUES:TRATAMIENTO 3 2.0234 0.67447 0.9153 0.4758
## Residuals 8 5.8952 0.73690TukeyHSD(mod2_T)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TEMPERATURA ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.4 -1.38977 0.5897702 0.3786366
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 1.105 -0.8388303 3.04883 0.3312832
## (+DH+Si1)-(-DH-Si ) 1.115 -0.8288303 3.05883 0.3246175
## (+DH+Si2)-(-DH-Si ) 0.860 -1.0838303 2.80383 0.5238518
## (+DH+Si1)-(+DH-Si ) 0.010 -1.9338303 1.95383 0.9999982
## (+DH+Si2)-(+DH-Si ) -0.245 -2.1888303 1.69883 0.9762850
## (+DH+Si2)-(+DH+Si1) -0.255 -2.1988303 1.68883 0.9734406
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -1.45 -4.846886 1.946886 0.6950048
## 1:(+DH-Si )-1:(-DH-Si ) 0.70 -2.696886 4.096886 0.9868169
## 2:(+DH-Si )-1:(-DH-Si ) 0.06 -3.336886 3.456886 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 0.33 -3.066886 3.726886 0.9998723
## 2:(+DH+Si1)-1:(-DH-Si ) 0.45 -2.946886 3.846886 0.9990483
## 1:(+DH+Si2)-1:(-DH-Si ) -0.05 -3.446886 3.346886 1.0000000
## 2:(+DH+Si2)-1:(-DH-Si ) 0.32 -3.076886 3.716886 0.9998958
## 1:(+DH-Si )-2:(-DH-Si ) 2.15 -1.246886 5.546886 0.3100030
## 2:(+DH-Si )-2:(-DH-Si ) 1.51 -1.886886 4.906886 0.6582674
## 1:(+DH+Si1)-2:(-DH-Si ) 1.78 -1.616886 5.176886 0.4954702
## 2:(+DH+Si1)-2:(-DH-Si ) 1.90 -1.496886 5.296886 0.4291159
## 1:(+DH+Si2)-2:(-DH-Si ) 1.40 -1.996886 4.796886 0.7251692
## 2:(+DH+Si2)-2:(-DH-Si ) 1.77 -1.626886 5.166886 0.5012256
## 2:(+DH-Si )-1:(+DH-Si ) -0.64 -4.036886 2.756886 0.9920230
## 1:(+DH+Si1)-1:(+DH-Si ) -0.37 -3.766886 3.026886 0.9997292
## 2:(+DH+Si1)-1:(+DH-Si ) -0.25 -3.646886 3.146886 0.9999801
## 1:(+DH+Si2)-1:(+DH-Si ) -0.75 -4.146886 2.646886 0.9808510
## 2:(+DH+Si2)-1:(+DH-Si ) -0.38 -3.776886 3.016886 0.9996779
## 1:(+DH+Si1)-2:(+DH-Si ) 0.27 -3.126886 3.666886 0.9999666
## 2:(+DH+Si1)-2:(+DH-Si ) 0.39 -3.006886 3.786886 0.9996188
## 1:(+DH+Si2)-2:(+DH-Si ) -0.11 -3.506886 3.286886 0.9999999
## 2:(+DH+Si2)-2:(+DH-Si ) 0.26 -3.136886 3.656886 0.9999741
## 2:(+DH+Si1)-1:(+DH+Si1) 0.12 -3.276886 3.516886 0.9999999
## 1:(+DH+Si2)-1:(+DH+Si1) -0.38 -3.776886 3.016886 0.9996779
## 2:(+DH+Si2)-1:(+DH+Si1) -0.01 -3.406886 3.386886 1.0000000
## 1:(+DH+Si2)-2:(+DH+Si1) -0.50 -3.896886 2.896886 0.9981638
## 2:(+DH+Si2)-2:(+DH+Si1) -0.13 -3.526886 3.266886 0.9999998
## 2:(+DH+Si2)-1:(+DH+Si2) 0.37 -3.026886 3.766886 0.9997292ggplot(datos_m2)+
aes(TRATAMIENTO,N_HOJAS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=N_HOJAS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_NH=aov(N_HOJAS~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_NH)## Analysis of Variance Table
##
## Response: N_HOJAS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.00 0.00000 0.0000 1.0000
## TRATAMIENTO 3 0.25 0.08333 0.1667 0.9159
## BLOQUES:TRATAMIENTO 3 1.50 0.50000 1.0000 0.4411
## Residuals 8 4.00 0.50000TukeyHSD(mod2_NH)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_HOJAS ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0 -0.8152956 0.8152956 1
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.00 -1.601176 1.601176 1.0000000
## (+DH+Si1)-(-DH-Si ) 0.25 -1.351176 1.851176 0.9568114
## (+DH+Si2)-(-DH-Si ) 0.25 -1.351176 1.851176 0.9568114
## (+DH+Si1)-(+DH-Si ) 0.25 -1.351176 1.851176 0.9568114
## (+DH+Si2)-(+DH-Si ) 0.25 -1.351176 1.851176 0.9568114
## (+DH+Si2)-(+DH+Si1) 0.00 -1.601176 1.601176 1.0000000
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -1.000000e+00 -3.79809 1.79809 0.8294285
## 1:(+DH-Si )-1:(-DH-Si ) -5.000000e-01 -3.29809 2.29809 0.9941212
## 2:(+DH-Si )-1:(-DH-Si ) -5.000000e-01 -3.29809 2.29809 0.9941212
## 1:(+DH+Si1)-1:(-DH-Si ) -5.000000e-01 -3.29809 2.29809 0.9941212
## 2:(+DH+Si1)-1:(-DH-Si ) 8.881784e-16 -2.79809 2.79809 1.0000000
## 1:(+DH+Si2)-1:(-DH-Si ) -5.000000e-01 -3.29809 2.29809 0.9941212
## 2:(+DH+Si2)-1:(-DH-Si ) 8.881784e-16 -2.79809 2.79809 1.0000000
## 1:(+DH-Si )-2:(-DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 2:(+DH-Si )-2:(-DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 1:(+DH+Si1)-2:(-DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 2:(+DH+Si1)-2:(-DH-Si ) 1.000000e+00 -1.79809 3.79809 0.8294285
## 1:(+DH+Si2)-2:(-DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 2:(+DH+Si2)-2:(-DH-Si ) 1.000000e+00 -1.79809 3.79809 0.8294285
## 2:(+DH-Si )-1:(+DH-Si ) 0.000000e+00 -2.79809 2.79809 1.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) 0.000000e+00 -2.79809 2.79809 1.0000000
## 2:(+DH+Si1)-1:(+DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 1:(+DH+Si2)-1:(+DH-Si ) 0.000000e+00 -2.79809 2.79809 1.0000000
## 2:(+DH+Si2)-1:(+DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 1:(+DH+Si1)-2:(+DH-Si ) 0.000000e+00 -2.79809 2.79809 1.0000000
## 2:(+DH+Si1)-2:(+DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 1:(+DH+Si2)-2:(+DH-Si ) 0.000000e+00 -2.79809 2.79809 1.0000000
## 2:(+DH+Si2)-2:(+DH-Si ) 5.000000e-01 -2.29809 3.29809 0.9941212
## 2:(+DH+Si1)-1:(+DH+Si1) 5.000000e-01 -2.29809 3.29809 0.9941212
## 1:(+DH+Si2)-1:(+DH+Si1) 0.000000e+00 -2.79809 2.79809 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si1) 5.000000e-01 -2.29809 3.29809 0.9941212
## 1:(+DH+Si2)-2:(+DH+Si1) -5.000000e-01 -3.29809 2.29809 0.9941212
## 2:(+DH+Si2)-2:(+DH+Si1) 0.000000e+00 -2.79809 2.79809 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si2) 5.000000e-01 -2.29809 3.29809 0.9941212ggplot(datos_m2)+
aes(TRATAMIENTO,LONGITUD_P_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=LONGITUD_P_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PA=aov(LONGITUD_P_AEREA~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PA)## Analysis of Variance Table
##
## Response: LONGITUD_P_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.3025 0.3025 0.1720 0.6892
## TRATAMIENTO 3 7.5950 2.5317 1.4395 0.3017
## BLOQUES:TRATAMIENTO 3 5.3325 1.7775 1.0107 0.4369
## Residuals 8 14.0700 1.7588TukeyHSD(mod2_PA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = LONGITUD_P_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.275 -1.254087 1.804087 0.689241
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -1.750 -4.753005 1.253005 0.3130240
## (+DH+Si1)-(-DH-Si ) -0.175 -3.178005 2.828005 0.9974996
## (+DH+Si2)-(-DH-Si ) -0.875 -3.878005 2.128005 0.7887386
## (+DH+Si1)-(+DH-Si ) 1.575 -1.428005 4.578005 0.3923490
## (+DH+Si2)-(+DH-Si ) 0.875 -2.128005 3.878005 0.7887386
## (+DH+Si2)-(+DH+Si1) -0.700 -3.703005 2.303005 0.8757068
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 2.00 -3.247818 7.247818 0.7862405
## 1:(+DH-Si )-1:(-DH-Si ) -1.00 -6.247818 4.247818 0.9914910
## 2:(+DH-Si )-1:(-DH-Si ) -0.50 -5.747818 4.747818 0.9998877
## 1:(+DH+Si1)-1:(-DH-Si ) 1.40 -3.847818 6.647818 0.9502748
## 2:(+DH+Si1)-1:(-DH-Si ) 0.25 -4.997818 5.497818 0.9999990
## 1:(+DH+Si2)-1:(-DH-Si ) 0.25 -4.997818 5.497818 0.9999990
## 2:(+DH+Si2)-1:(-DH-Si ) 0.00 -5.247818 5.247818 1.0000000
## 1:(+DH-Si )-2:(-DH-Si ) -3.00 -8.247818 2.247818 0.4072645
## 2:(+DH-Si )-2:(-DH-Si ) -2.50 -7.747818 2.747818 0.5917355
## 1:(+DH+Si1)-2:(-DH-Si ) -0.60 -5.847818 4.647818 0.9996290
## 2:(+DH+Si1)-2:(-DH-Si ) -1.75 -6.997818 3.497818 0.8687772
## 1:(+DH+Si2)-2:(-DH-Si ) -1.75 -6.997818 3.497818 0.8687772
## 2:(+DH+Si2)-2:(-DH-Si ) -2.00 -7.247818 3.247818 0.7862405
## 2:(+DH-Si )-1:(+DH-Si ) 0.50 -4.747818 5.747818 0.9998877
## 1:(+DH+Si1)-1:(+DH-Si ) 2.40 -2.847818 7.647818 0.6314783
## 2:(+DH+Si1)-1:(+DH-Si ) 1.25 -3.997818 6.497818 0.9715584
## 1:(+DH+Si2)-1:(+DH-Si ) 1.25 -3.997818 6.497818 0.9715584
## 2:(+DH+Si2)-1:(+DH-Si ) 1.00 -4.247818 6.247818 0.9914910
## 1:(+DH+Si1)-2:(+DH-Si ) 1.90 -3.347818 7.147818 0.8212268
## 2:(+DH+Si1)-2:(+DH-Si ) 0.75 -4.497818 5.997818 0.9984699
## 1:(+DH+Si2)-2:(+DH-Si ) 0.75 -4.497818 5.997818 0.9984699
## 2:(+DH+Si2)-2:(+DH-Si ) 0.50 -4.747818 5.747818 0.9998877
## 2:(+DH+Si1)-1:(+DH+Si1) -1.15 -6.397818 4.097818 0.9816023
## 1:(+DH+Si2)-1:(+DH+Si1) -1.15 -6.397818 4.097818 0.9816023
## 2:(+DH+Si2)-1:(+DH+Si1) -1.40 -6.647818 3.847818 0.9502748
## 1:(+DH+Si2)-2:(+DH+Si1) 0.00 -5.247818 5.247818 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) -0.25 -5.497818 4.997818 0.9999990
## 2:(+DH+Si2)-1:(+DH+Si2) -0.25 -5.497818 4.997818 0.9999990ggplot(datos_m2)+
aes(TRATAMIENTO,AREA_FOLIAR, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=AREA_FOLIAR, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_AF=aov(AREA_FOLIAR~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_AF)## Analysis of Variance Table
##
## Response: AREA_FOLIAR
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 6.9 6.89 0.0139 0.9091
## TRATAMIENTO 3 1589.1 529.71 1.0669 0.4157
## BLOQUES:TRATAMIENTO 3 498.3 166.10 0.3345 0.8010
## Residuals 8 3972.1 496.51TukeyHSD(mod2_AF)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AREA_FOLIAR ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -1.3125 -27.00429 24.37929 0.909126
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -25.175 -75.63164 25.28164 0.4310065
## (+DH+Si1)-(-DH-Si ) -18.000 -68.45664 32.45664 0.6758682
## (+DH+Si2)-(-DH-Si ) -23.400 -73.85664 27.05664 0.4877485
## (+DH+Si1)-(+DH-Si ) 7.175 -43.28164 57.63164 0.9666839
## (+DH+Si2)-(+DH-Si ) 1.775 -48.68164 52.23164 0.9994437
## (+DH+Si2)-(+DH+Si1) -5.400 -55.85664 45.05664 0.9851477
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 11.55 -76.62409 99.72409 0.9991137
## 1:(+DH-Si )-1:(-DH-Si ) -13.75 -101.92409 74.42409 0.9973859
## 2:(+DH-Si )-1:(-DH-Si ) -25.05 -113.22409 63.12409 0.9334221
## 1:(+DH+Si1)-1:(-DH-Si ) -5.50 -93.67409 82.67409 0.9999936
## 2:(+DH+Si1)-1:(-DH-Si ) -18.95 -107.12409 69.22409 0.9834265
## 1:(+DH+Si2)-1:(-DH-Si ) -21.60 -109.77409 66.57409 0.9672187
## 2:(+DH+Si2)-1:(-DH-Si ) -13.65 -101.82409 74.52409 0.9974993
## 1:(+DH-Si )-2:(-DH-Si ) -25.30 -113.47409 62.87409 0.9303704
## 2:(+DH-Si )-2:(-DH-Si ) -36.60 -124.77409 51.57409 0.7191784
## 1:(+DH+Si1)-2:(-DH-Si ) -17.05 -105.22409 71.12409 0.9907545
## 2:(+DH+Si1)-2:(-DH-Si ) -30.50 -118.67409 57.67409 0.8488988
## 1:(+DH+Si2)-2:(-DH-Si ) -33.15 -121.32409 55.02409 0.7959186
## 2:(+DH+Si2)-2:(-DH-Si ) -25.20 -113.37409 62.97409 0.9316010
## 2:(+DH-Si )-1:(+DH-Si ) -11.30 -99.47409 76.87409 0.9992284
## 1:(+DH+Si1)-1:(+DH-Si ) 8.25 -79.92409 96.42409 0.9999004
## 2:(+DH+Si1)-1:(+DH-Si ) -5.20 -93.37409 82.97409 0.9999956
## 1:(+DH+Si2)-1:(+DH-Si ) -7.85 -96.02409 80.32409 0.9999285
## 2:(+DH+Si2)-1:(+DH-Si ) 0.10 -88.07409 88.27409 1.0000000
## 1:(+DH+Si1)-2:(+DH-Si ) 19.55 -68.62409 107.72409 0.9804184
## 2:(+DH+Si1)-2:(+DH-Si ) 6.10 -82.07409 94.27409 0.9999869
## 1:(+DH+Si2)-2:(+DH-Si ) 3.45 -84.72409 91.62409 0.9999997
## 2:(+DH+Si2)-2:(+DH-Si ) 11.40 -76.77409 99.57409 0.9991841
## 2:(+DH+Si1)-1:(+DH+Si1) -13.45 -101.62409 74.72409 0.9977145
## 1:(+DH+Si2)-1:(+DH+Si1) -16.10 -104.27409 72.07409 0.9933372
## 2:(+DH+Si2)-1:(+DH+Si1) -8.15 -96.32409 80.02409 0.9999082
## 1:(+DH+Si2)-2:(+DH+Si1) -2.65 -90.82409 85.52409 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) 5.30 -82.87409 93.47409 0.9999950
## 2:(+DH+Si2)-1:(+DH+Si2) 7.95 -80.22409 96.12409 0.9999222ggplot(datos_m2)+
aes(TRATAMIENTO,LONGITUD_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=LONGITUD_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_LR=aov(LONGITUD_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_LR)## Analysis of Variance Table
##
## Response: LONGITUD_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.250 0.25000 0.8889 0.3734
## TRATAMIENTO 3 1.375 0.45833 1.6296 0.2579
## BLOQUES:TRATAMIENTO 3 1.625 0.54167 1.9259 0.2039
## Residuals 8 2.250 0.28125TukeyHSD(mod2_LR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = LONGITUD_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.25 -0.8614717 0.3614717 0.3733749
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.625 -1.825882 0.575882 0.3982901
## (+DH+Si1)-(-DH-Si ) -0.625 -1.825882 0.575882 0.3982901
## (+DH+Si2)-(-DH-Si ) -0.750 -1.950882 0.450882 0.2640524
## (+DH+Si1)-(+DH-Si ) 0.000 -1.200882 1.200882 1.0000000
## (+DH+Si2)-(+DH-Si ) -0.125 -1.325882 1.075882 0.9862918
## (+DH+Si2)-(+DH+Si1) -0.125 -1.325882 1.075882 0.9862918
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -1.00 -3.098568 1.0985676 0.5914719
## 1:(+DH-Si )-1:(-DH-Si ) -1.00 -3.098568 1.0985676 0.5914719
## 2:(+DH-Si )-1:(-DH-Si ) -1.25 -3.348568 0.8485676 0.3669059
## 1:(+DH+Si1)-1:(-DH-Si ) -1.50 -3.598568 0.5985676 0.2099214
## 2:(+DH+Si1)-1:(-DH-Si ) -0.75 -2.848568 1.3485676 0.8294285
## 1:(+DH+Si2)-1:(-DH-Si ) -1.00 -3.098568 1.0985676 0.5914719
## 2:(+DH+Si2)-1:(-DH-Si ) -1.50 -3.598568 0.5985676 0.2099214
## 1:(+DH-Si )-2:(-DH-Si ) 0.00 -2.098568 2.0985676 1.0000000
## 2:(+DH-Si )-2:(-DH-Si ) -0.25 -2.348568 1.8485676 0.9995163
## 1:(+DH+Si1)-2:(-DH-Si ) -0.50 -2.598568 1.5985676 0.9715197
## 2:(+DH+Si1)-2:(-DH-Si ) 0.25 -1.848568 2.3485676 0.9995163
## 1:(+DH+Si2)-2:(-DH-Si ) 0.00 -2.098568 2.0985676 1.0000000
## 2:(+DH+Si2)-2:(-DH-Si ) -0.50 -2.598568 1.5985676 0.9715197
## 2:(+DH-Si )-1:(+DH-Si ) -0.25 -2.348568 1.8485676 0.9995163
## 1:(+DH+Si1)-1:(+DH-Si ) -0.50 -2.598568 1.5985676 0.9715197
## 2:(+DH+Si1)-1:(+DH-Si ) 0.25 -1.848568 2.3485676 0.9995163
## 1:(+DH+Si2)-1:(+DH-Si ) 0.00 -2.098568 2.0985676 1.0000000
## 2:(+DH+Si2)-1:(+DH-Si ) -0.50 -2.598568 1.5985676 0.9715197
## 1:(+DH+Si1)-2:(+DH-Si ) -0.25 -2.348568 1.8485676 0.9995163
## 2:(+DH+Si1)-2:(+DH-Si ) 0.50 -1.598568 2.5985676 0.9715197
## 1:(+DH+Si2)-2:(+DH-Si ) 0.25 -1.848568 2.3485676 0.9995163
## 2:(+DH+Si2)-2:(+DH-Si ) -0.25 -2.348568 1.8485676 0.9995163
## 2:(+DH+Si1)-1:(+DH+Si1) 0.75 -1.348568 2.8485676 0.8294285
## 1:(+DH+Si2)-1:(+DH+Si1) 0.50 -1.598568 2.5985676 0.9715197
## 2:(+DH+Si2)-1:(+DH+Si1) 0.00 -2.098568 2.0985676 1.0000000
## 1:(+DH+Si2)-2:(+DH+Si1) -0.25 -2.348568 1.8485676 0.9995163
## 2:(+DH+Si2)-2:(+DH+Si1) -0.75 -2.848568 1.3485676 0.8294285
## 2:(+DH+Si2)-1:(+DH+Si2) -0.50 -2.598568 1.5985676 0.9715197ggplot(datos_m2)+
aes(TRATAMIENTO,DIAMETRO_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=DIAMETRO_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_DR=aov(DIAMETRO_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_DR)## Analysis of Variance Table
##
## Response: DIAMETRO_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.000 0.000 0.0000 0.99773
## TRATAMIENTO 3 151.858 50.619 4.3692 0.04235 *
## BLOQUES:TRATAMIENTO 3 1.107 0.369 0.0318 0.99176
## Residuals 8 92.685 11.586
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_DR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = DIAMETRO_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.005 -3.919541 3.929541 0.9977278
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -5.5375 -13.244988 2.1699882 0.1769148
## (+DH+Si1)-(-DH-Si ) -8.1350 -15.842488 -0.4275118 0.0389746
## (+DH+Si2)-(-DH-Si ) -6.6975 -14.404988 1.0099882 0.0904068
## (+DH+Si1)-(+DH-Si ) -2.5975 -10.304988 5.1099882 0.7109103
## (+DH+Si2)-(+DH-Si ) -1.1600 -8.867488 6.5474882 0.9609822
## (+DH+Si2)-(+DH+Si1) 1.4375 -6.269988 9.1449882 0.9300823
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.110 -13.579 13.359004 1.0000000
## 1:(+DH-Si )-1:(-DH-Si ) -5.975 -19.444 7.494004 0.6601761
## 2:(+DH-Si )-1:(-DH-Si ) -5.210 -18.679 8.259004 0.7753499
## 1:(+DH+Si1)-1:(-DH-Si ) -7.835 -21.304 5.634004 0.3900161
## 2:(+DH+Si1)-1:(-DH-Si ) -8.545 -22.014 4.924004 0.3078839
## 1:(+DH+Si2)-1:(-DH-Si ) -6.790 -20.259 6.679004 0.5349233
## 2:(+DH+Si2)-1:(-DH-Si ) -6.715 -20.184 6.754004 0.5461868
## 1:(+DH-Si )-2:(-DH-Si ) -5.865 -19.334 7.604004 0.6772045
## 2:(+DH-Si )-2:(-DH-Si ) -5.100 -18.569 8.369004 0.7908895
## 1:(+DH+Si1)-2:(-DH-Si ) -7.725 -21.194 5.744004 0.4040093
## 2:(+DH+Si1)-2:(-DH-Si ) -8.435 -21.904 5.034004 0.3196649
## 1:(+DH+Si2)-2:(-DH-Si ) -6.680 -20.149 6.789004 0.5514713
## 2:(+DH+Si2)-2:(-DH-Si ) -6.605 -20.074 6.864004 0.5628507
## 2:(+DH-Si )-1:(+DH-Si ) 0.765 -12.704 14.234004 0.9999966
## 1:(+DH+Si1)-1:(+DH-Si ) -1.860 -15.329 11.609004 0.9987645
## 2:(+DH+Si1)-1:(+DH-Si ) -2.570 -16.039 10.899004 0.9914266
## 1:(+DH+Si2)-1:(+DH-Si ) -0.815 -14.284 12.654004 0.9999948
## 2:(+DH+Si2)-1:(+DH-Si ) -0.740 -14.209 12.729004 0.9999973
## 1:(+DH+Si1)-2:(+DH-Si ) -2.625 -16.094 10.844004 0.9903357
## 2:(+DH+Si1)-2:(+DH-Si ) -3.335 -16.804 10.134004 0.9654155
## 1:(+DH+Si2)-2:(+DH-Si ) -1.580 -15.049 11.889004 0.9995620
## 2:(+DH+Si2)-2:(+DH-Si ) -1.505 -14.974 11.964004 0.9996803
## 2:(+DH+Si1)-1:(+DH+Si1) -0.710 -14.179 12.759004 0.9999980
## 1:(+DH+Si2)-1:(+DH+Si1) 1.045 -12.424 14.514004 0.9999716
## 2:(+DH+Si2)-1:(+DH+Si1) 1.120 -12.349 14.589004 0.9999547
## 1:(+DH+Si2)-2:(+DH+Si1) 1.755 -11.714 15.224004 0.9991429
## 2:(+DH+Si2)-2:(+DH+Si1) 1.830 -11.639 15.299004 0.9988842
## 2:(+DH+Si2)-1:(+DH+Si2) 0.075 -13.394 13.544004 1.0000000ggplot(datos_m2)+
aes(TRATAMIENTO,PESO_F_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=PESO_F_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PFR=aov(PESO_F_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PFR)## Analysis of Variance Table
##
## Response: PESO_F_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 6.113 6.113 0.3354 0.5785
## TRATAMIENTO 3 280.066 93.355 5.1213 0.0288 *
## BLOQUES:TRATAMIENTO 3 21.363 7.121 0.3907 0.7631
## Residuals 8 145.830 18.229
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_PFR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -1.23625 -6.159004 3.686504 0.5784545
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -9.0975 -18.7654 0.5703998 0.0652655
## (+DH+Si1)-(-DH-Si ) -10.5275 -20.1954 -0.8596002 0.0335705
## (+DH+Si2)-(-DH-Si ) -9.0725 -18.7404 0.5953998 0.0660334
## (+DH+Si1)-(+DH-Si ) -1.4300 -11.0979 8.2378998 0.9628185
## (+DH+Si2)-(+DH-Si ) 0.0250 -9.6429 9.6928998 0.9999998
## (+DH+Si2)-(+DH+Si1) 1.4550 -8.2129 11.1228998 0.9609858
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -4.470 -21.36487 12.424867 0.9521972
## 1:(+DH-Si )-1:(-DH-Si ) -10.795 -27.68987 6.099867 0.3014960
## 2:(+DH-Si )-1:(-DH-Si ) -11.870 -28.76487 5.024867 0.2227761
## 1:(+DH+Si1)-1:(-DH-Si ) -13.790 -30.68487 3.104867 0.1266149
## 2:(+DH+Si1)-1:(-DH-Si ) -11.735 -28.62987 5.159867 0.2315738
## 1:(+DH+Si2)-1:(-DH-Si ) -10.580 -27.47487 6.314867 0.3197043
## 2:(+DH+Si2)-1:(-DH-Si ) -12.035 -28.92987 4.859867 0.2124257
## 1:(+DH-Si )-2:(-DH-Si ) -6.325 -23.21987 10.569867 0.7988689
## 2:(+DH-Si )-2:(-DH-Si ) -7.400 -24.29487 9.494867 0.6718770
## 1:(+DH+Si1)-2:(-DH-Si ) -9.320 -26.21487 7.574867 0.4431013
## 2:(+DH+Si1)-2:(-DH-Si ) -7.265 -24.15987 9.629867 0.6884845
## 1:(+DH+Si2)-2:(-DH-Si ) -6.110 -23.00487 10.784867 0.8219479
## 2:(+DH+Si2)-2:(-DH-Si ) -7.565 -24.45987 9.329867 0.6514857
## 2:(+DH-Si )-1:(+DH-Si ) -1.075 -17.96987 15.819867 0.9999926
## 1:(+DH+Si1)-1:(+DH-Si ) -2.995 -19.88987 13.899867 0.9943885
## 2:(+DH+Si1)-1:(+DH-Si ) -0.940 -17.83487 15.954867 0.9999971
## 1:(+DH+Si2)-1:(+DH-Si ) 0.215 -16.67987 17.109867 1.0000000
## 2:(+DH+Si2)-1:(+DH-Si ) -1.240 -18.13487 15.654867 0.9999805
## 1:(+DH+Si1)-2:(+DH-Si ) -1.920 -18.81487 14.974867 0.9996432
## 2:(+DH+Si1)-2:(+DH-Si ) 0.135 -16.75987 17.029867 1.0000000
## 1:(+DH+Si2)-2:(+DH-Si ) 1.290 -15.60487 18.184867 0.9999745
## 2:(+DH+Si2)-2:(+DH-Si ) -0.165 -17.05987 16.729867 1.0000000
## 2:(+DH+Si1)-1:(+DH+Si1) 2.055 -14.83987 18.949867 0.9994471
## 1:(+DH+Si2)-1:(+DH+Si1) 3.210 -13.68487 20.104867 0.9916315
## 2:(+DH+Si2)-1:(+DH+Si1) 1.755 -15.13987 18.649867 0.9998014
## 1:(+DH+Si2)-2:(+DH+Si1) 1.155 -15.73987 18.049867 0.9999880
## 2:(+DH+Si2)-2:(+DH+Si1) -0.300 -17.19487 16.594867 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si2) -1.455 -18.34987 15.439867 0.9999427ggplot(datos_m2)+
aes(TRATAMIENTO,PESO_F_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=PESO_F_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PFA=aov(PESO_F_AEREA~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PFA)## Analysis of Variance Table
##
## Response: PESO_F_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.0900 0.09000 0.1148 0.7435
## TRATAMIENTO 3 2.5675 0.85585 1.0917 0.4067
## BLOQUES:TRATAMIENTO 3 2.0082 0.66942 0.8539 0.5028
## Residuals 8 6.2718 0.78397TukeyHSD(mod2_PFA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.15 -1.170895 0.8708953 0.7434664
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.9250 -2.929958 1.079958 0.4917649
## (+DH+Si1)-(-DH-Si ) -0.2725 -2.277458 1.732458 0.9706525
## (+DH+Si2)-(-DH-Si ) -0.9025 -2.907458 1.102458 0.5106807
## (+DH+Si1)-(+DH-Si ) 0.6525 -1.352458 2.657458 0.7311585
## (+DH+Si2)-(+DH-Si ) 0.0225 -1.982458 2.027458 0.9999818
## (+DH+Si2)-(+DH+Si1) -0.6300 -2.634958 1.374958 0.7505116
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.270 -3.233707 3.773707 0.9999729
## 1:(+DH-Si )-1:(-DH-Si ) -0.795 -4.298707 2.708707 0.9778796
## 2:(+DH-Si )-1:(-DH-Si ) -0.785 -4.288707 2.718707 0.9793051
## 1:(+DH+Si1)-1:(-DH-Si ) 0.535 -2.968707 4.038707 0.9977001
## 2:(+DH+Si1)-1:(-DH-Si ) -0.810 -4.313707 2.693707 0.9756131
## 1:(+DH+Si2)-1:(-DH-Si ) -1.000 -4.503707 2.503707 0.9320171
## 2:(+DH+Si2)-1:(-DH-Si ) -0.535 -4.038707 2.968707 0.9977001
## 1:(+DH-Si )-2:(-DH-Si ) -1.065 -4.568707 2.438707 0.9103095
## 2:(+DH-Si )-2:(-DH-Si ) -1.055 -4.558707 2.448707 0.9138767
## 1:(+DH+Si1)-2:(-DH-Si ) 0.265 -3.238707 3.768707 0.9999761
## 2:(+DH+Si1)-2:(-DH-Si ) -1.080 -4.583707 2.423707 0.9048056
## 1:(+DH+Si2)-2:(-DH-Si ) -1.270 -4.773707 2.233707 0.8204850
## 2:(+DH+Si2)-2:(-DH-Si ) -0.805 -4.308707 2.698707 0.9763859
## 2:(+DH-Si )-1:(+DH-Si ) 0.010 -3.493707 3.513707 1.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) 1.330 -2.173707 4.833707 0.7890994
## 2:(+DH+Si1)-1:(+DH-Si ) -0.015 -3.518707 3.488707 1.0000000
## 1:(+DH+Si2)-1:(+DH-Si ) -0.205 -3.708707 3.298707 0.9999958
## 2:(+DH+Si2)-1:(+DH-Si ) 0.260 -3.243707 3.763707 0.9999790
## 1:(+DH+Si1)-2:(+DH-Si ) 1.320 -2.183707 4.823707 0.7944571
## 2:(+DH+Si1)-2:(+DH-Si ) -0.025 -3.528707 3.478707 1.0000000
## 1:(+DH+Si2)-2:(+DH-Si ) -0.215 -3.718707 3.288707 0.9999942
## 2:(+DH+Si2)-2:(+DH-Si ) 0.250 -3.253707 3.753707 0.9999839
## 2:(+DH+Si1)-1:(+DH+Si1) -1.345 -4.848707 2.158707 0.7809761
## 1:(+DH+Si2)-1:(+DH+Si1) -1.535 -5.038707 1.968707 0.6716592
## 2:(+DH+Si2)-1:(+DH+Si1) -1.070 -4.573707 2.433707 0.9084952
## 1:(+DH+Si2)-2:(+DH+Si1) -0.190 -3.693707 3.313707 0.9999975
## 2:(+DH+Si2)-2:(+DH+Si1) 0.275 -3.228707 3.778707 0.9999693
## 2:(+DH+Si2)-1:(+DH+Si2) 0.465 -3.038707 3.968707 0.9990373ggplot(datos_m2)+
aes(TRATAMIENTO,PESO_S_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=PESO_S_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PSR=aov(PESO_S_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PSR)## Analysis of Variance Table
##
## Response: PESO_S_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.02372 0.023716 0.5603 0.4756
## TRATAMIENTO 3 0.06870 0.022899 0.5410 0.6676
## BLOQUES:TRATAMIENTO 3 0.04081 0.013602 0.3214 0.8100
## Residuals 8 0.33862 0.042328TukeyHSD(mod2_PSR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.077 -0.314216 0.160216 0.4755714
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.15425 -0.6201234 0.3116234 0.7212924
## (+DH+Si1)-(-DH-Si ) 0.00050 -0.4653734 0.4663734 1.0000000
## (+DH+Si2)-(-DH-Si ) -0.09225 -0.5581234 0.3736234 0.9181083
## (+DH+Si1)-(+DH-Si ) 0.15475 -0.3111234 0.6206234 0.7194124
## (+DH+Si2)-(+DH-Si ) 0.06200 -0.4038734 0.5278734 0.9723410
## (+DH+Si2)-(+DH+Si1) -0.09275 -0.5586234 0.3731234 0.9169411
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0140 -0.800124 0.828124 1.0000000
## 1:(+DH-Si )-1:(-DH-Si ) -0.0360 -0.850124 0.778124 0.9999994
## 2:(+DH-Si )-1:(-DH-Si ) -0.2585 -1.072624 0.555624 0.8922630
## 1:(+DH+Si1)-1:(-DH-Si ) 0.0675 -0.746624 0.881624 0.9999556
## 2:(+DH+Si1)-1:(-DH-Si ) -0.0525 -0.866624 0.761624 0.9999919
## 1:(+DH+Si2)-1:(-DH-Si ) -0.0955 -0.909624 0.718624 0.9995621
## 2:(+DH+Si2)-1:(-DH-Si ) -0.0750 -0.889124 0.739124 0.9999102
## 1:(+DH-Si )-2:(-DH-Si ) -0.0500 -0.864124 0.764124 0.9999942
## 2:(+DH-Si )-2:(-DH-Si ) -0.2725 -1.086624 0.541624 0.8668490
## 1:(+DH+Si1)-2:(-DH-Si ) 0.0535 -0.760624 0.867624 0.9999908
## 2:(+DH+Si1)-2:(-DH-Si ) -0.0665 -0.880624 0.747624 0.9999598
## 1:(+DH+Si2)-2:(-DH-Si ) -0.1095 -0.923624 0.704624 0.9989529
## 2:(+DH+Si2)-2:(-DH-Si ) -0.0890 -0.903124 0.725124 0.9997227
## 2:(+DH-Si )-1:(+DH-Si ) -0.2225 -1.036624 0.591624 0.9442879
## 1:(+DH+Si1)-1:(+DH-Si ) 0.1035 -0.710624 0.917624 0.9992668
## 2:(+DH+Si1)-1:(+DH-Si ) -0.0165 -0.830624 0.797624 1.0000000
## 1:(+DH+Si2)-1:(+DH-Si ) -0.0595 -0.873624 0.754624 0.9999810
## 2:(+DH+Si2)-1:(+DH-Si ) -0.0390 -0.853124 0.775124 0.9999989
## 1:(+DH+Si1)-2:(+DH-Si ) 0.3260 -0.488124 1.140124 0.7487001
## 2:(+DH+Si1)-2:(+DH-Si ) 0.2060 -0.608124 1.020124 0.9614943
## 1:(+DH+Si2)-2:(+DH-Si ) 0.1630 -0.651124 0.977124 0.9887611
## 2:(+DH+Si2)-2:(+DH-Si ) 0.1835 -0.630624 0.997624 0.9786406
## 2:(+DH+Si1)-1:(+DH+Si1) -0.1200 -0.934124 0.694124 0.9981480
## 1:(+DH+Si2)-1:(+DH+Si1) -0.1630 -0.977124 0.651124 0.9887611
## 2:(+DH+Si2)-1:(+DH+Si1) -0.1425 -0.956624 0.671624 0.9947904
## 1:(+DH+Si2)-2:(+DH+Si1) -0.0430 -0.857124 0.771124 0.9999979
## 2:(+DH+Si2)-2:(+DH+Si1) -0.0225 -0.836624 0.791624 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si2) 0.0205 -0.793624 0.834624 1.0000000ggplot(datos_m2)+
aes(TRATAMIENTO,PESO_S_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=PESO_S_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PSA=aov(PESO_S_AEREA~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PSA)## Analysis of Variance Table
##
## Response: PESO_S_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.008696 0.0086956 0.8258 0.3900
## TRATAMIENTO 3 0.035830 0.0119432 1.1342 0.3918
## BLOQUES:TRATAMIENTO 3 0.023256 0.0077521 0.7362 0.5593
## Residuals 8 0.084239 0.0105298TukeyHSD(mod2_PSA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.046625 -0.1649402 0.07169016 0.3900434
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.03100 -0.2633616 0.2013616 0.9721475
## (+DH+Si1)-(-DH-Si ) 0.09725 -0.1351116 0.3296116 0.5654604
## (+DH+Si2)-(-DH-Si ) 0.02050 -0.2118616 0.2528616 0.9915203
## (+DH+Si1)-(+DH-Si ) 0.12825 -0.1041116 0.3606116 0.3533724
## (+DH+Si2)-(+DH-Si ) 0.05150 -0.1808616 0.2838616 0.8904957
## (+DH+Si2)-(+DH+Si1) -0.07675 -0.3091116 0.1556116 0.7226824
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0525 -0.353557 0.458557 0.9991840
## 1:(+DH-Si )-1:(-DH-Si ) 0.0435 -0.362557 0.449557 0.9997571
## 2:(+DH-Si )-1:(-DH-Si ) -0.0530 -0.459057 0.353057 0.9991335
## 1:(+DH+Si1)-1:(-DH-Si ) 0.1940 -0.212057 0.600057 0.5888790
## 2:(+DH+Si1)-1:(-DH-Si ) 0.0530 -0.353057 0.459057 0.9991335
## 1:(+DH+Si2)-1:(-DH-Si ) 0.0475 -0.358557 0.453557 0.9995698
## 2:(+DH+Si2)-1:(-DH-Si ) 0.0460 -0.360057 0.452057 0.9996505
## 1:(+DH-Si )-2:(-DH-Si ) -0.0090 -0.415057 0.397057 1.0000000
## 2:(+DH-Si )-2:(-DH-Si ) -0.1055 -0.511557 0.300557 0.9561916
## 1:(+DH+Si1)-2:(-DH-Si ) 0.1415 -0.264557 0.547557 0.8446353
## 2:(+DH+Si1)-2:(-DH-Si ) 0.0005 -0.405557 0.406557 1.0000000
## 1:(+DH+Si2)-2:(-DH-Si ) -0.0050 -0.411057 0.401057 1.0000000
## 2:(+DH+Si2)-2:(-DH-Si ) -0.0065 -0.412557 0.399557 1.0000000
## 2:(+DH-Si )-1:(+DH-Si ) -0.0965 -0.502557 0.309557 0.9718880
## 1:(+DH+Si1)-1:(+DH-Si ) 0.1505 -0.255557 0.556557 0.8057550
## 2:(+DH+Si1)-1:(+DH-Si ) 0.0095 -0.396557 0.415557 1.0000000
## 1:(+DH+Si2)-1:(+DH-Si ) 0.0040 -0.402057 0.410057 1.0000000
## 2:(+DH+Si2)-1:(+DH-Si ) 0.0025 -0.403557 0.408557 1.0000000
## 1:(+DH+Si1)-2:(+DH-Si ) 0.2470 -0.159057 0.653057 0.3468007
## 2:(+DH+Si1)-2:(+DH-Si ) 0.1060 -0.300057 0.512057 0.9551803
## 1:(+DH+Si2)-2:(+DH-Si ) 0.1005 -0.305557 0.506557 0.9654875
## 2:(+DH+Si2)-2:(+DH-Si ) 0.0990 -0.307057 0.505057 0.9679929
## 2:(+DH+Si1)-1:(+DH+Si1) -0.1410 -0.547057 0.265057 0.8466873
## 1:(+DH+Si2)-1:(+DH+Si1) -0.1465 -0.552557 0.259557 0.8234707
## 2:(+DH+Si2)-1:(+DH+Si1) -0.1480 -0.554057 0.258057 0.8169046
## 1:(+DH+Si2)-2:(+DH+Si1) -0.0055 -0.411557 0.400557 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) -0.0070 -0.413057 0.399057 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si2) -0.0015 -0.407557 0.404557 1.0000000ggplot(datos_m2)+
aes(TRATAMIENTO,SPAD, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=SPAD, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_SPAD=aov(SPAD~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_SPAD)## Analysis of Variance Table
##
## Response: SPAD
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 61.54 61.544 4.4428 0.068108 .
## TRATAMIENTO 3 425.61 141.871 10.2414 0.004095 **
## BLOQUES:TRATAMIENTO 3 4.12 1.375 0.0992 0.958209
## Residuals 8 110.82 13.853
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_SPAD)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = SPAD ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -3.9225 -8.213879 0.3688794 0.0681077
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 4.050 -4.3779295 12.47793 0.4602130
## (+DH+Si1)-(-DH-Si ) 10.920 2.4920705 19.34793 0.0136499
## (+DH+Si2)-(-DH-Si ) 12.775 4.3470705 21.20293 0.0055321
## (+DH+Si1)-(+DH-Si ) 6.870 -1.5579295 15.29793 0.1152085
## (+DH+Si2)-(+DH-Si ) 8.725 0.2970705 17.15293 0.0426707
## (+DH+Si2)-(+DH+Si1) 1.855 -6.5729295 10.28293 0.8924104
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -5.400000e+00 -20.1279914 9.327991 0.8130074
## 1:(+DH-Si )-1:(-DH-Si ) 2.700000e+00 -12.0279914 17.427991 0.9931818
## 2:(+DH-Si )-1:(-DH-Si ) -7.105427e-15 -14.7279914 14.727991 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 9.890000e+00 -4.8379914 24.617991 0.2587029
## 2:(+DH+Si1)-1:(-DH-Si ) 6.550000e+00 -8.1779914 21.277991 0.6578357
## 1:(+DH+Si2)-1:(-DH-Si ) 1.220000e+01 -2.5279914 26.927991 0.1190882
## 2:(+DH+Si2)-1:(-DH-Si ) 7.950000e+00 -6.7779914 22.677991 0.4651850
## 1:(+DH-Si )-2:(-DH-Si ) 8.100000e+00 -6.6279914 22.827991 0.4461810
## 2:(+DH-Si )-2:(-DH-Si ) 5.400000e+00 -9.3279914 20.127991 0.8130074
## 1:(+DH+Si1)-2:(-DH-Si ) 1.529000e+01 0.5620086 30.017991 0.0413224
## 2:(+DH+Si1)-2:(-DH-Si ) 1.195000e+01 -2.7779914 26.677991 0.1297480
## 1:(+DH+Si2)-2:(-DH-Si ) 1.760000e+01 2.8720086 32.327991 0.0192315
## 2:(+DH+Si2)-2:(-DH-Si ) 1.335000e+01 -1.3779914 28.077991 0.0801707
## 2:(+DH-Si )-1:(+DH-Si ) -2.700000e+00 -17.4279914 12.027991 0.9931818
## 1:(+DH+Si1)-1:(+DH-Si ) 7.190000e+00 -7.5379914 21.917991 0.5673641
## 2:(+DH+Si1)-1:(+DH-Si ) 3.850000e+00 -10.8779914 18.577991 0.9548819
## 1:(+DH+Si2)-1:(+DH-Si ) 9.500000e+00 -5.2279914 24.227991 0.2930949
## 2:(+DH+Si2)-1:(+DH-Si ) 5.250000e+00 -9.4779914 19.977991 0.8310304
## 1:(+DH+Si1)-2:(+DH-Si ) 9.890000e+00 -4.8379914 24.617991 0.2587029
## 2:(+DH+Si1)-2:(+DH-Si ) 6.550000e+00 -8.1779914 21.277991 0.6578357
## 1:(+DH+Si2)-2:(+DH-Si ) 1.220000e+01 -2.5279914 26.927991 0.1190882
## 2:(+DH+Si2)-2:(+DH-Si ) 7.950000e+00 -6.7779914 22.677991 0.4651850
## 2:(+DH+Si1)-1:(+DH+Si1) -3.340000e+00 -18.0679914 11.387991 0.9779426
## 1:(+DH+Si2)-1:(+DH+Si1) 2.310000e+00 -12.4179914 17.037991 0.9972927
## 2:(+DH+Si2)-1:(+DH+Si1) -1.940000e+00 -16.6679914 12.787991 0.9990820
## 1:(+DH+Si2)-2:(+DH+Si1) 5.650000e+00 -9.0779914 20.377991 0.7814644
## 2:(+DH+Si2)-2:(+DH+Si1) 1.400000e+00 -13.3279914 16.127991 0.9998894
## 2:(+DH+Si2)-1:(+DH+Si2) -4.250000e+00 -18.9779914 10.477991 0.9285745#INTERACCIÓN
ggplot(datos_m2)+
aes(TRATAMIENTO,PESO_S_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=PESO_S_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PSB=aov(PESO_S_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PSB)## Analysis of Variance Table
##
## Response: PESO_S_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 6.1260e-06 6.1260e-06 3.2989 0.106858
## TRATAMIENTO 3 9.6142e-05 3.2047e-05 17.2587 0.000746 ***
## BLOQUES:TRATAMIENTO 3 2.9770e-06 9.9200e-07 0.5344 0.671538
## Residuals 8 1.4855e-05 1.8570e-06
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_PSB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.0012375 -0.002808664 0.0003336636 0.1068581
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.000650 -0.0024356411 0.0037356411 0.9038300
## (+DH+Si1)-(-DH-Si ) 0.006200 0.0031143589 0.0092856411 0.0009131
## (+DH+Si2)-(-DH-Si ) 0.003325 0.0002393589 0.0064106411 0.0353101
## (+DH+Si1)-(+DH-Si ) 0.005550 0.0024643589 0.0086356411 0.0018986
## (+DH+Si2)-(+DH-Si ) 0.002675 -0.0004106411 0.0057606411 0.0912444
## (+DH+Si2)-(+DH+Si1) -0.002875 -0.0059606411 0.0002106411 0.0680594
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.00020 -5.192225e-03 0.0055922254 0.9999998
## 1:(+DH-Si )-1:(-DH-Si ) 0.00175 -3.642225e-03 0.0071422254 0.8822156
## 2:(+DH-Si )-1:(-DH-Si ) -0.00025 -5.642225e-03 0.0051422254 0.9999992
## 1:(+DH+Si1)-1:(-DH-Si ) 0.00720 1.807775e-03 0.0125922254 0.0100123
## 2:(+DH+Si1)-1:(-DH-Si ) 0.00540 7.774638e-06 0.0107922254 0.0496399
## 1:(+DH+Si2)-1:(-DH-Si ) 0.00410 -1.292225e-03 0.0094922254 0.1675999
## 2:(+DH+Si2)-1:(-DH-Si ) 0.00275 -2.642225e-03 0.0081422254 0.5231398
## 1:(+DH-Si )-2:(-DH-Si ) 0.00155 -3.842225e-03 0.0069422254 0.9298034
## 2:(+DH-Si )-2:(-DH-Si ) -0.00045 -5.842225e-03 0.0049422254 0.9999536
## 1:(+DH+Si1)-2:(-DH-Si ) 0.00700 1.607775e-03 0.0123922254 0.0118618
## 2:(+DH+Si1)-2:(-DH-Si ) 0.00520 -1.922254e-04 0.0105922254 0.0598160
## 1:(+DH+Si2)-2:(-DH-Si ) 0.00390 -1.492225e-03 0.0092922254 0.2013667
## 2:(+DH+Si2)-2:(-DH-Si ) 0.00255 -2.842225e-03 0.0079422254 0.5989781
## 2:(+DH-Si )-1:(+DH-Si ) -0.00200 -7.392225e-03 0.0033922254 0.8052662
## 1:(+DH+Si1)-1:(+DH-Si ) 0.00545 5.777464e-05 0.0108422254 0.0473875
## 2:(+DH+Si1)-1:(+DH-Si ) 0.00365 -1.742225e-03 0.0090422254 0.2521343
## 1:(+DH+Si2)-1:(+DH-Si ) 0.00235 -3.042225e-03 0.0077422254 0.6764378
## 2:(+DH+Si2)-1:(+DH-Si ) 0.00100 -4.392225e-03 0.0063922254 0.9927135
## 1:(+DH+Si1)-2:(+DH-Si ) 0.00745 2.057775e-03 0.0128422254 0.0081266
## 2:(+DH+Si1)-2:(+DH-Si ) 0.00565 2.577746e-04 0.0110422254 0.0393880
## 1:(+DH+Si2)-2:(+DH-Si ) 0.00435 -1.042225e-03 0.0097422254 0.1328348
## 2:(+DH+Si2)-2:(+DH-Si ) 0.00300 -2.392225e-03 0.0083922254 0.4345028
## 2:(+DH+Si1)-1:(+DH+Si1) -0.00180 -7.192225e-03 0.0035922254 0.8682483
## 1:(+DH+Si2)-1:(+DH+Si1) -0.00310 -8.492225e-03 0.0022922254 0.4016515
## 2:(+DH+Si2)-1:(+DH+Si1) -0.00445 -9.842225e-03 0.0009422254 0.1209646
## 1:(+DH+Si2)-2:(+DH+Si1) -0.00130 -6.692225e-03 0.0040922254 0.9697567
## 2:(+DH+Si2)-2:(+DH+Si1) -0.00265 -8.042225e-03 0.0027422254 0.5606725
## 2:(+DH+Si2)-1:(+DH+Si2) -0.00135 -6.742225e-03 0.0040422254 0.9634603ggplot(datos_m2)+
aes(TRATAMIENTO,PESO_F_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=PESO_F_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PFB=aov(PESO_F_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PFB)## Analysis of Variance Table
##
## Response: PESO_F_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.0001210 0.00012100 2.2670 0.17058
## TRATAMIENTO 3 0.0014865 0.00049550 9.2834 0.00553 **
## BLOQUES:TRATAMIENTO 3 0.0009235 0.00030783 5.7674 0.02125 *
## Residuals 8 0.0004270 0.00005338
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_PFB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.0055 -0.01392363 0.002923625 0.1705763
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.00925 -0.007293333 0.025793333 0.3436254
## (+DH+Si1)-(-DH-Si ) 0.01225 -0.004293333 0.028793333 0.1606255
## (+DH+Si2)-(-DH-Si ) -0.01250 -0.029043333 0.004043333 0.1502586
## (+DH+Si1)-(+DH-Si ) 0.00300 -0.013543333 0.019543333 0.9351303
## (+DH+Si2)-(+DH-Si ) -0.02175 -0.038293333 -0.005206667 0.0125941
## (+DH+Si2)-(+DH+Si1) -0.02475 -0.041293333 -0.008206667 0.0059833
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.0015 -0.030409837 0.027409837 0.9999982
## 1:(+DH-Si )-1:(-DH-Si ) 0.0005 -0.028409837 0.029409837 1.0000000
## 2:(+DH-Si )-1:(-DH-Si ) 0.0165 -0.012409837 0.045409837 0.4088806
## 1:(+DH+Si1)-1:(-DH-Si ) 0.0245 -0.004409837 0.053409837 0.1081108
## 2:(+DH+Si1)-1:(-DH-Si ) -0.0015 -0.030409837 0.027409837 0.9999982
## 1:(+DH+Si2)-1:(-DH-Si ) -0.0080 -0.036909837 0.020909837 0.9409730
## 2:(+DH+Si2)-1:(-DH-Si ) -0.0185 -0.047409837 0.010409837 0.3001439
## 1:(+DH-Si )-2:(-DH-Si ) 0.0020 -0.026909837 0.030909837 0.9999869
## 2:(+DH-Si )-2:(-DH-Si ) 0.0180 -0.010909837 0.046909837 0.3250087
## 1:(+DH+Si1)-2:(-DH-Si ) 0.0260 -0.002909837 0.054909837 0.0831043
## 2:(+DH+Si1)-2:(-DH-Si ) 0.0000 -0.028909837 0.028909837 1.0000000
## 1:(+DH+Si2)-2:(-DH-Si ) -0.0065 -0.035409837 0.022409837 0.9789180
## 2:(+DH+Si2)-2:(-DH-Si ) -0.0170 -0.045909837 0.011909837 0.3794003
## 2:(+DH-Si )-1:(+DH-Si ) 0.0160 -0.012909837 0.044909837 0.4398069
## 1:(+DH+Si1)-1:(+DH-Si ) 0.0240 -0.004909837 0.052909837 0.1180017
## 2:(+DH+Si1)-1:(+DH-Si ) -0.0020 -0.030909837 0.026909837 0.9999869
## 1:(+DH+Si2)-1:(+DH-Si ) -0.0085 -0.037409837 0.020409837 0.9223840
## 2:(+DH+Si2)-1:(+DH-Si ) -0.0190 -0.047909837 0.009909837 0.2768264
## 1:(+DH+Si1)-2:(+DH-Si ) 0.0080 -0.020909837 0.036909837 0.9409730
## 2:(+DH+Si1)-2:(+DH-Si ) -0.0180 -0.046909837 0.010909837 0.3250087
## 1:(+DH+Si2)-2:(+DH-Si ) -0.0245 -0.053409837 0.004409837 0.1081108
## 2:(+DH+Si2)-2:(+DH-Si ) -0.0350 -0.063909837 -0.006090163 0.0178521
## 2:(+DH+Si1)-1:(+DH+Si1) -0.0260 -0.054909837 0.002909837 0.0831043
## 1:(+DH+Si2)-1:(+DH+Si1) -0.0325 -0.061409837 -0.003590163 0.0270521
## 2:(+DH+Si2)-1:(+DH+Si1) -0.0430 -0.071909837 -0.014090163 0.0051170
## 1:(+DH+Si2)-2:(+DH+Si1) -0.0065 -0.035409837 0.022409837 0.9789180
## 2:(+DH+Si2)-2:(+DH+Si1) -0.0170 -0.045909837 0.011909837 0.3794003
## 2:(+DH+Si2)-1:(+DH+Si2) -0.0105 -0.039409837 0.018409837 0.8191953ggplot(datos_m2)+
aes(TRATAMIENTO,PESO_SAT_H_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=PESO_SAT_H_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_PSHB=aov(PESO_SAT_H_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_PSHB)## Analysis of Variance Table
##
## Response: PESO_SAT_H_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.0008851 0.00088506 2.2992 0.16791
## TRATAMIENTO 3 0.0009452 0.00031506 0.8185 0.51911
## BLOQUES:TRATAMIENTO 3 0.0042647 0.00142156 3.6930 0.06194 .
## Residuals 8 0.0030795 0.00038494
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_PSHB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_SAT_H_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.014875 -0.0374967 0.007746698 0.1679089
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.01500 -0.05942722 0.02942722 0.7098162
## (+DH+Si1)-(-DH-Si ) -0.01850 -0.06292722 0.02592722 0.5691950
## (+DH+Si2)-(-DH-Si ) -0.01875 -0.06317722 0.02567722 0.5592848
## (+DH+Si1)-(+DH-Si ) -0.00350 -0.04792722 0.04092722 0.9939126
## (+DH+Si2)-(+DH-Si ) -0.00375 -0.04817722 0.04067722 0.9925485
## (+DH+Si2)-(+DH+Si1) -0.00025 -0.04467722 0.04417722 0.9999977
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0205 -0.05713755 0.09813755 0.9526529
## 1:(+DH-Si )-1:(-DH-Si ) -0.0080 -0.08563755 0.06963755 0.9998117
## 2:(+DH-Si )-1:(-DH-Si ) -0.0015 -0.07913755 0.07613755 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 0.0025 -0.07513755 0.08013755 0.9999999
## 2:(+DH+Si1)-1:(-DH-Si ) -0.0190 -0.09663755 0.05863755 0.9673817
## 1:(+DH+Si2)-1:(-DH-Si ) 0.0240 -0.05363755 0.10163755 0.9036467
## 2:(+DH+Si2)-1:(-DH-Si ) -0.0410 -0.11863755 0.03663755 0.4875356
## 1:(+DH-Si )-2:(-DH-Si ) -0.0285 -0.10613755 0.04913755 0.8122095
## 2:(+DH-Si )-2:(-DH-Si ) -0.0220 -0.09963755 0.05563755 0.9341935
## 1:(+DH+Si1)-2:(-DH-Si ) -0.0180 -0.09563755 0.05963755 0.9752482
## 2:(+DH+Si1)-2:(-DH-Si ) -0.0395 -0.11713755 0.03813755 0.5255771
## 1:(+DH+Si2)-2:(-DH-Si ) 0.0035 -0.07413755 0.08113755 0.9999993
## 2:(+DH+Si2)-2:(-DH-Si ) -0.0615 -0.13913755 0.01613755 0.1429453
## 2:(+DH-Si )-1:(+DH-Si ) 0.0065 -0.07113755 0.08413755 0.9999526
## 1:(+DH+Si1)-1:(+DH-Si ) 0.0105 -0.06713755 0.08813755 0.9989159
## 2:(+DH+Si1)-1:(+DH-Si ) -0.0110 -0.08863755 0.06663755 0.9985501
## 1:(+DH+Si2)-1:(+DH-Si ) 0.0320 -0.04563755 0.10963755 0.7251093
## 2:(+DH+Si2)-1:(+DH-Si ) -0.0330 -0.11063755 0.04463755 0.6987409
## 1:(+DH+Si1)-2:(+DH-Si ) 0.0040 -0.07363755 0.08163755 0.9999983
## 2:(+DH+Si1)-2:(+DH-Si ) -0.0175 -0.09513755 0.06013755 0.9786352
## 1:(+DH+Si2)-2:(+DH-Si ) 0.0255 -0.05213755 0.10313755 0.8764243
## 2:(+DH+Si2)-2:(+DH-Si ) -0.0395 -0.11713755 0.03813755 0.5255771
## 2:(+DH+Si1)-1:(+DH+Si1) -0.0215 -0.09913755 0.05613755 0.9407702
## 1:(+DH+Si2)-1:(+DH+Si1) 0.0215 -0.05613755 0.09913755 0.9407702
## 2:(+DH+Si2)-1:(+DH+Si1) -0.0435 -0.12113755 0.03413755 0.4273887
## 1:(+DH+Si2)-2:(+DH+Si1) 0.0430 -0.03463755 0.12063755 0.4390561
## 2:(+DH+Si2)-2:(+DH+Si1) -0.0220 -0.09963755 0.05563755 0.9341935
## 2:(+DH+Si2)-1:(+DH+Si2) -0.0650 -0.14263755 0.01263755 0.1138633ggplot(datos_m2)+
aes(TRATAMIENTO,N_ESTOMAS_ADAX, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ADAX, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_NEAx=aov(N_ESTOMAS_ADAX~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_NEAx)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ADAX
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 81.0 81.00 1.0891 0.32718
## TRATAMIENTO 3 2686.2 895.42 12.0392 0.00246 **
## BLOQUES:TRATAMIENTO 3 153.5 51.17 0.6880 0.58431
## Residuals 8 595.0 74.37
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_NEAx)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ADAX ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -4.5 -14.4436 5.443598 0.3271826
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -9.75 -29.27844 9.77844 0.4304948
## (+DH+Si1)-(-DH-Si ) 24.25 4.72156 43.77844 0.0171846
## (+DH+Si2)-(-DH-Si ) -4.00 -23.52844 15.52844 0.9105538
## (+DH+Si1)-(+DH-Si ) 34.00 14.47156 53.52844 0.0023412
## (+DH+Si2)-(+DH-Si ) 5.75 -13.77844 25.27844 0.7837123
## (+DH+Si2)-(+DH+Si1) -28.25 -47.77844 -8.72156 0.0072997
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -3.0 -37.1263774 31.126377 0.9999343
## 1:(+DH-Si )-1:(-DH-Si ) -6.5 -40.6263774 27.626377 0.9915131
## 2:(+DH-Si )-1:(-DH-Si ) -16.0 -50.1263774 18.126377 0.6074283
## 1:(+DH+Si1)-1:(-DH-Si ) 28.0 -6.1263774 62.126377 0.1239232
## 2:(+DH+Si1)-1:(-DH-Si ) 17.5 -16.6263774 51.626377 0.5175470
## 1:(+DH+Si2)-1:(-DH-Si ) -8.0 -42.1263774 26.126377 0.9737936
## 2:(+DH+Si2)-1:(-DH-Si ) -3.0 -37.1263774 31.126377 0.9999343
## 1:(+DH-Si )-2:(-DH-Si ) -3.5 -37.6263774 30.626377 0.9998174
## 2:(+DH-Si )-2:(-DH-Si ) -13.0 -47.1263774 21.126377 0.7865695
## 1:(+DH+Si1)-2:(-DH-Si ) 31.0 -3.1263774 65.126377 0.0793819
## 2:(+DH+Si1)-2:(-DH-Si ) 20.5 -13.6263774 54.626377 0.3587551
## 1:(+DH+Si2)-2:(-DH-Si ) -5.0 -39.1263774 29.126377 0.9982154
## 2:(+DH+Si2)-2:(-DH-Si ) 0.0 -34.1263774 34.126377 1.0000000
## 2:(+DH-Si )-1:(+DH-Si ) -9.5 -43.6263774 24.626377 0.9393274
## 1:(+DH+Si1)-1:(+DH-Si ) 34.5 0.3736226 68.626377 0.0473321
## 2:(+DH+Si1)-1:(+DH-Si ) 24.0 -10.1263774 58.126377 0.2220325
## 1:(+DH+Si2)-1:(+DH-Si ) -1.5 -35.6263774 32.626377 0.9999994
## 2:(+DH+Si2)-1:(+DH-Si ) 3.5 -30.6263774 37.626377 0.9998174
## 1:(+DH+Si1)-2:(+DH-Si ) 44.0 9.8736226 78.126377 0.0123550
## 2:(+DH+Si1)-2:(+DH-Si ) 33.5 -0.6263774 67.626377 0.0548274
## 1:(+DH+Si2)-2:(+DH-Si ) 8.0 -26.1263774 42.126377 0.9737936
## 2:(+DH+Si2)-2:(+DH-Si ) 13.0 -21.1263774 47.126377 0.7865695
## 2:(+DH+Si1)-1:(+DH+Si1) -10.5 -44.6263774 23.626377 0.9055425
## 1:(+DH+Si2)-1:(+DH+Si1) -36.0 -70.1263774 -1.873623 0.0380251
## 2:(+DH+Si2)-1:(+DH+Si1) -31.0 -65.1263774 3.126377 0.0793819
## 1:(+DH+Si2)-2:(+DH+Si1) -25.5 -59.6263774 8.626377 0.1789065
## 2:(+DH+Si2)-2:(+DH+Si1) -20.5 -54.6263774 13.626377 0.3587551
## 2:(+DH+Si2)-1:(+DH+Si2) 5.0 -29.1263774 39.126377 0.9982154ggplot(datos_m2)+
aes(TRATAMIENTO,N_ESTOMAS_ABA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ABA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_NEAbx=aov(N_ESTOMAS_ABA~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_NEAbx)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ABA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.06 0.06 0.001 0.9751597
## TRATAMIENTO 3 2901.69 967.23 15.971 0.0009705 ***
## BLOQUES:TRATAMIENTO 3 248.19 82.73 1.366 0.3209649
## Residuals 8 484.50 60.56
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_NEAbx)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ABA ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.125 -9.097883 8.847883 0.9751597
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -8.75 -26.372031 8.872031 0.4347636
## (+DH+Si1)-(-DH-Si ) 25.50 7.877969 43.122031 0.0072866
## (+DH+Si2)-(-DH-Si ) -5.50 -23.122031 12.122031 0.7541157
## (+DH+Si1)-(+DH-Si ) 34.25 16.627969 51.872031 0.0011409
## (+DH+Si2)-(+DH-Si ) 3.25 -14.372031 20.872031 0.9321384
## (+DH+Si2)-(+DH+Si1) -31.00 -48.622031 -13.377969 0.0021909
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 7.0 -23.794886 37.794886 0.9776703
## 1:(+DH-Si )-1:(-DH-Si ) -4.0 -34.794886 26.794886 0.9991597
## 2:(+DH-Si )-1:(-DH-Si ) -6.5 -37.294886 24.294886 0.9849688
## 1:(+DH+Si1)-1:(-DH-Si ) 25.5 -5.294886 56.294886 0.1192660
## 2:(+DH+Si1)-1:(-DH-Si ) 32.5 1.705114 63.294886 0.0379371
## 1:(+DH+Si2)-1:(-DH-Si ) 4.0 -26.794886 34.794886 0.9991597
## 2:(+DH+Si2)-1:(-DH-Si ) -8.0 -38.794886 22.794886 0.9562180
## 1:(+DH-Si )-2:(-DH-Si ) -11.0 -41.794886 19.794886 0.8297493
## 2:(+DH-Si )-2:(-DH-Si ) -13.5 -44.294886 17.294886 0.6710818
## 1:(+DH+Si1)-2:(-DH-Si ) 18.5 -12.294886 49.294886 0.3586900
## 2:(+DH+Si1)-2:(-DH-Si ) 25.5 -5.294886 56.294886 0.1192660
## 1:(+DH+Si2)-2:(-DH-Si ) -3.0 -33.794886 27.794886 0.9998699
## 2:(+DH+Si2)-2:(-DH-Si ) -15.0 -45.794886 15.794886 0.5696094
## 2:(+DH-Si )-1:(+DH-Si ) -2.5 -33.294886 28.294886 0.9999615
## 1:(+DH+Si1)-1:(+DH-Si ) 29.5 -1.294886 60.294886 0.0617784
## 2:(+DH+Si1)-1:(+DH-Si ) 36.5 5.705114 67.294886 0.0201468
## 1:(+DH+Si2)-1:(+DH-Si ) 8.0 -22.794886 38.794886 0.9562180
## 2:(+DH+Si2)-1:(+DH-Si ) -4.0 -34.794886 26.794886 0.9991597
## 1:(+DH+Si1)-2:(+DH-Si ) 32.0 1.205114 62.794886 0.0411226
## 2:(+DH+Si1)-2:(+DH-Si ) 39.0 8.205114 69.794886 0.0137370
## 1:(+DH+Si2)-2:(+DH-Si ) 10.5 -20.294886 41.294886 0.8569545
## 2:(+DH+Si2)-2:(+DH-Si ) -1.5 -32.294886 29.294886 0.9999988
## 2:(+DH+Si1)-1:(+DH+Si1) 7.0 -23.794886 37.794886 0.9776703
## 1:(+DH+Si2)-1:(+DH+Si1) -21.5 -52.294886 9.294886 0.2275971
## 2:(+DH+Si2)-1:(+DH+Si1) -33.5 -64.294886 -2.705114 0.0323172
## 1:(+DH+Si2)-2:(+DH+Si1) -28.5 -59.294886 2.294886 0.0727988
## 2:(+DH+Si2)-2:(+DH+Si1) -40.5 -71.294886 -9.705114 0.0109726
## 2:(+DH+Si2)-1:(+DH+Si2) -12.0 -42.794886 18.794886 0.7698227ggplot(datos_m2)+
aes(TRATAMIENTO,N_ESTOMAS_ABIERTOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ABIERTOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_NEA=aov(N_ESTOMAS_ABIERTOS~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_NEA)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ABIERTOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 121.00 121.000 5.4382 0.048012 *
## TRATAMIENTO 3 748.25 249.417 11.2097 0.003087 **
## BLOQUES:TRATAMIENTO 3 12.50 4.167 0.1873 0.902124
## Residuals 8 178.00 22.250
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_NEA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ABIERTOS ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 5.5 0.06130011 10.9387 0.0480119
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -5.00 -15.681176 5.681176 0.4805354
## (+DH+Si1)-(-DH-Si ) 13.25 2.568824 23.931176 0.0172794
## (+DH+Si2)-(-DH-Si ) -0.75 -11.431176 9.931176 0.9956609
## (+DH+Si1)-(+DH-Si ) 18.25 7.568824 28.931176 0.0026392
## (+DH+Si2)-(+DH-Si ) 4.25 -6.431176 14.931176 0.6020957
## (+DH+Si2)-(+DH+Si1) -14.00 -24.681176 -3.318824 0.0128091
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 6.5 -12.1655894 25.165589 0.8450343
## 1:(+DH-Si )-1:(-DH-Si ) -3.0 -21.6655894 15.665589 0.9968621
## 2:(+DH-Si )-1:(-DH-Si ) -0.5 -19.1655894 18.165589 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 13.5 -5.1655894 32.165589 0.2013741
## 2:(+DH+Si1)-1:(-DH-Si ) 19.5 0.8344106 38.165589 0.0399995
## 1:(+DH+Si2)-1:(-DH-Si ) -1.0 -19.6655894 17.665589 0.9999977
## 2:(+DH+Si2)-1:(-DH-Si ) 6.0 -12.6655894 24.665589 0.8867003
## 1:(+DH-Si )-2:(-DH-Si ) -9.5 -28.1655894 9.165589 0.5252100
## 2:(+DH-Si )-2:(-DH-Si ) -7.0 -25.6655894 11.665589 0.7976667
## 1:(+DH+Si1)-2:(-DH-Si ) 7.0 -11.6655894 25.665589 0.7976667
## 2:(+DH+Si1)-2:(-DH-Si ) 13.0 -5.6655894 31.665589 0.2294779
## 1:(+DH+Si2)-2:(-DH-Si ) -7.5 -26.1655894 11.165589 0.7459568
## 2:(+DH+Si2)-2:(-DH-Si ) -0.5 -19.1655894 18.165589 1.0000000
## 2:(+DH-Si )-1:(+DH-Si ) 2.5 -16.1655894 21.165589 0.9989802
## 1:(+DH+Si1)-1:(+DH-Si ) 16.5 -2.1655894 35.165589 0.0898398
## 2:(+DH+Si1)-1:(+DH-Si ) 22.5 3.8344106 41.165589 0.0183011
## 1:(+DH+Si2)-1:(+DH-Si ) 2.0 -16.6655894 20.665589 0.9997567
## 2:(+DH+Si2)-1:(+DH-Si ) 9.0 -9.6655894 27.665589 0.5797599
## 1:(+DH+Si1)-2:(+DH-Si ) 14.0 -4.6655894 32.665589 0.1764197
## 2:(+DH+Si1)-2:(+DH-Si ) 20.0 1.3344106 38.665589 0.0350291
## 1:(+DH+Si2)-2:(+DH-Si ) -0.5 -19.1655894 18.165589 1.0000000
## 2:(+DH+Si2)-2:(+DH-Si ) 6.5 -12.1655894 25.165589 0.8450343
## 2:(+DH+Si1)-1:(+DH+Si1) 6.0 -12.6655894 24.665589 0.8867003
## 1:(+DH+Si2)-1:(+DH+Si1) -14.5 -33.1655894 4.165589 0.1543582
## 2:(+DH+Si2)-1:(+DH+Si1) -7.5 -26.1655894 11.165589 0.7459568
## 1:(+DH+Si2)-2:(+DH+Si1) -20.5 -39.1655894 -1.834411 0.0307030
## 2:(+DH+Si2)-2:(+DH+Si1) -13.5 -32.1655894 5.165589 0.2013741
## 2:(+DH+Si2)-1:(+DH+Si2) 7.0 -11.6655894 25.665589 0.7976667ggplot(datos_m2)+
aes(TRATAMIENTO,N_ESTOMAS_P_ABIERTOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_P_ABIERTOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_NEPA=aov(N_ESTOMAS_P_ABIERTOS~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_NEPA)## Analysis of Variance Table
##
## Response: N_ESTOMAS_P_ABIERTOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 4.0 4.00 0.0717 0.7956
## TRATAMIENTO 3 6625.2 2208.42 39.6129 3.795e-05 ***
## BLOQUES:TRATAMIENTO 3 118.5 39.50 0.7085 0.5735
## Residuals 8 446.0 55.75
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_NEPA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_P_ABIERTOS ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 1 -7.608996 9.608996 0.795584
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -13.25 -30.15739 3.657387 0.1326277
## (+DH+Si1)-(-DH-Si ) 37.25 20.34261 54.157387 0.0004868
## (+DH+Si2)-(-DH-Si ) -11.50 -28.40739 5.407387 0.2088110
## (+DH+Si1)-(+DH-Si ) 50.50 33.59261 67.407387 0.0000550
## (+DH+Si2)-(+DH-Si ) 1.75 -15.15739 18.657387 0.9865132
## (+DH+Si2)-(+DH+Si1) -48.75 -65.65739 -31.842613 0.0000713
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 6.000000e+00 -23.5460301 35.54603 0.9878409
## 1:(+DH-Si )-1:(-DH-Si ) -1.300000e+01 -42.5460301 16.54603 0.6677315
## 2:(+DH-Si )-1:(-DH-Si ) -7.500000e+00 -37.0460301 22.04603 0.9608859
## 1:(+DH+Si1)-1:(-DH-Si ) 4.400000e+01 14.4539699 73.54603 0.0050773
## 2:(+DH+Si1)-1:(-DH-Si ) 3.650000e+01 6.9539699 66.04603 0.0158868
## 1:(+DH+Si2)-1:(-DH-Si ) -8.500000e+00 -38.0460301 21.04603 0.9295446
## 2:(+DH+Si2)-1:(-DH-Si ) -8.500000e+00 -38.0460301 21.04603 0.9295446
## 1:(+DH-Si )-2:(-DH-Si ) -1.900000e+01 -48.5460301 10.54603 0.2957909
## 2:(+DH-Si )-2:(-DH-Si ) -1.350000e+01 -43.0460301 16.04603 0.6323543
## 1:(+DH+Si1)-2:(-DH-Si ) 3.800000e+01 8.4539699 67.54603 0.0125387
## 2:(+DH+Si1)-2:(-DH-Si ) 3.050000e+01 0.9539699 60.04603 0.0425483
## 1:(+DH+Si2)-2:(-DH-Si ) -1.450000e+01 -44.0460301 15.04603 0.5620831
## 2:(+DH+Si2)-2:(-DH-Si ) -1.450000e+01 -44.0460301 15.04603 0.5620831
## 2:(+DH-Si )-1:(+DH-Si ) 5.500000e+00 -24.0460301 35.04603 0.9925547
## 1:(+DH+Si1)-1:(+DH-Si ) 5.700000e+01 27.4539699 86.54603 0.0008990
## 2:(+DH+Si1)-1:(+DH-Si ) 4.950000e+01 19.9539699 79.04603 0.0023534
## 1:(+DH+Si2)-1:(+DH-Si ) 4.500000e+00 -25.0460301 34.04603 0.9977358
## 2:(+DH+Si2)-1:(+DH-Si ) 4.500000e+00 -25.0460301 34.04603 0.9977358
## 1:(+DH+Si1)-2:(+DH-Si ) 5.150000e+01 21.9539699 81.04603 0.0018038
## 2:(+DH+Si1)-2:(+DH-Si ) 4.400000e+01 14.4539699 73.54603 0.0050773
## 1:(+DH+Si2)-2:(+DH-Si ) -1.000000e+00 -30.5460301 28.54603 0.9999999
## 2:(+DH+Si2)-2:(+DH-Si ) -1.000000e+00 -30.5460301 28.54603 0.9999999
## 2:(+DH+Si1)-1:(+DH+Si1) -7.500000e+00 -37.0460301 22.04603 0.9608859
## 1:(+DH+Si2)-1:(+DH+Si1) -5.250000e+01 -82.0460301 -22.95397 0.0015833
## 2:(+DH+Si2)-1:(+DH+Si1) -5.250000e+01 -82.0460301 -22.95397 0.0015833
## 1:(+DH+Si2)-2:(+DH+Si1) -4.500000e+01 -74.5460301 -15.45397 0.0043964
## 2:(+DH+Si2)-2:(+DH+Si1) -4.500000e+01 -74.5460301 -15.45397 0.0043964
## 2:(+DH+Si2)-1:(+DH+Si2) -7.105427e-15 -29.5460301 29.54603 1.0000000ggplot(datos_m2)+
aes(TRATAMIENTO,N_ESTOMAS_CERRADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_CERRADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_NEC=aov(N_ESTOMAS_CERRADOS~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_NEC)## Analysis of Variance Table
##
## Response: N_ESTOMAS_CERRADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 441.00 441.00 5.3133 0.05007 .
## TRATAMIENTO 3 8.25 2.75 0.0331 0.99127
## BLOQUES:TRATAMIENTO 3 176.50 58.83 0.7088 0.57334
## Residuals 8 664.00 83.00
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_NEC)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_CERRADOS ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -10.5 -21.00435 0.004348785 0.0500746
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.00 -20.62971 20.62971 1.0000000
## (+DH+Si1)-(-DH-Si ) -0.75 -21.37971 19.87971 0.9993862
## (+DH+Si2)-(-DH-Si ) 1.25 -19.37971 21.87971 0.9971934
## (+DH+Si1)-(+DH-Si ) -0.75 -21.37971 19.87971 0.9993862
## (+DH+Si2)-(+DH-Si ) 1.25 -19.37971 21.87971 0.9971934
## (+DH+Si2)-(+DH+Si1) 2.00 -18.62971 22.62971 0.9888420
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -8.500000e+00 -44.55087 27.55087 0.9730050
## 1:(+DH-Si )-1:(-DH-Si ) 6.000000e+00 -30.05087 42.05087 0.9961326
## 2:(+DH-Si )-1:(-DH-Si ) -1.450000e+01 -50.55087 21.55087 0.7451583
## 1:(+DH+Si1)-1:(-DH-Si ) -4.000000e+00 -40.05087 32.05087 0.9996946
## 2:(+DH+Si1)-1:(-DH-Si ) -6.000000e+00 -42.05087 30.05087 0.9961326
## 1:(+DH+Si2)-1:(-DH-Si ) 2.500000e+00 -33.55087 38.55087 0.9999867
## 2:(+DH+Si2)-1:(-DH-Si ) -8.500000e+00 -44.55087 27.55087 0.9730050
## 1:(+DH-Si )-2:(-DH-Si ) 1.450000e+01 -21.55087 50.55087 0.7451583
## 2:(+DH-Si )-2:(-DH-Si ) -6.000000e+00 -42.05087 30.05087 0.9961326
## 1:(+DH+Si1)-2:(-DH-Si ) 4.500000e+00 -31.55087 40.55087 0.9993476
## 2:(+DH+Si1)-2:(-DH-Si ) 2.500000e+00 -33.55087 38.55087 0.9999867
## 1:(+DH+Si2)-2:(-DH-Si ) 1.100000e+01 -25.05087 47.05087 0.9088354
## 2:(+DH+Si2)-2:(-DH-Si ) 7.105427e-15 -36.05087 36.05087 1.0000000
## 2:(+DH-Si )-1:(+DH-Si ) -2.050000e+01 -56.55087 15.55087 0.4125585
## 1:(+DH+Si1)-1:(+DH-Si ) -1.000000e+01 -46.05087 26.05087 0.9403166
## 2:(+DH+Si1)-1:(+DH-Si ) -1.200000e+01 -48.05087 24.05087 0.8697164
## 1:(+DH+Si2)-1:(+DH-Si ) -3.500000e+00 -39.55087 32.55087 0.9998728
## 2:(+DH+Si2)-1:(+DH-Si ) -1.450000e+01 -50.55087 21.55087 0.7451583
## 1:(+DH+Si1)-2:(+DH-Si ) 1.050000e+01 -25.55087 46.55087 0.9255590
## 2:(+DH+Si1)-2:(+DH-Si ) 8.500000e+00 -27.55087 44.55087 0.9730050
## 1:(+DH+Si2)-2:(+DH-Si ) 1.700000e+01 -19.05087 53.05087 0.6017817
## 2:(+DH+Si2)-2:(+DH-Si ) 6.000000e+00 -30.05087 42.05087 0.9961326
## 2:(+DH+Si1)-1:(+DH+Si1) -2.000000e+00 -38.05087 34.05087 0.9999971
## 1:(+DH+Si2)-1:(+DH+Si1) 6.500000e+00 -29.55087 42.55087 0.9938072
## 2:(+DH+Si2)-1:(+DH+Si1) -4.500000e+00 -40.55087 31.55087 0.9993476
## 1:(+DH+Si2)-2:(+DH+Si1) 8.500000e+00 -27.55087 44.55087 0.9730050
## 2:(+DH+Si2)-2:(+DH+Si1) -2.500000e+00 -38.55087 33.55087 0.9999867
## 2:(+DH+Si2)-1:(+DH+Si2) -1.100000e+01 -47.05087 25.05087 0.9088354ggplot(datos_m2)+
aes(TRATAMIENTO,CRA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=CRA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_CRA=aov(CRA~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_CRA)## Analysis of Variance Table
##
## Response: CRA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 16.85 16.846 0.3047 0.5960
## TRATAMIENTO 3 428.45 142.816 2.5835 0.1259
## BLOQUES:TRATAMIENTO 3 280.99 93.665 1.6944 0.2448
## Residuals 8 442.23 55.279TukeyHSD(mod2_CRA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = CRA ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 2.052173 -6.520406 10.62475 0.5960081
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 9.061919 -7.773946 25.897784 0.3724518
## (+DH+Si1)-(-DH-Si ) 10.896844 -5.939021 27.732709 0.2401813
## (+DH+Si2)-(-DH-Si ) -0.561566 -17.397431 16.274299 0.9995255
## (+DH+Si1)-(+DH-Si ) 1.834924 -15.000940 18.670789 0.9843482
## (+DH+Si2)-(+DH-Si ) -9.623485 -26.459350 7.212380 0.3272015
## (+DH+Si2)-(+DH+Si1) -11.458410 -28.294275 5.377455 0.2084358
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -6.386131 -35.80717 23.034912 0.9825174
## 1:(+DH-Si )-1:(-DH-Si ) 3.247965 -26.17308 32.669008 0.9997045
## 2:(+DH-Si )-1:(-DH-Si ) 8.489743 -20.93130 37.910786 0.9285810
## 1:(+DH+Si1)-1:(-DH-Si ) 10.185277 -19.23577 39.606320 0.8484303
## 2:(+DH+Si1)-1:(-DH-Si ) 5.222280 -24.19876 34.643322 0.9943461
## 1:(+DH+Si2)-1:(-DH-Si ) -10.912652 -40.33369 18.508391 0.8052501
## 2:(+DH+Si2)-1:(-DH-Si ) 3.403389 -26.01765 32.824432 0.9995998
## 1:(+DH-Si )-2:(-DH-Si ) 9.634096 -19.78695 39.055139 0.8779097
## 2:(+DH-Si )-2:(-DH-Si ) 14.875874 -14.54517 44.296917 0.5319044
## 1:(+DH+Si1)-2:(-DH-Si ) 16.571408 -12.84963 45.992451 0.4221009
## 2:(+DH+Si1)-2:(-DH-Si ) 11.608410 -17.81263 41.029453 0.7602907
## 1:(+DH+Si2)-2:(-DH-Si ) -4.526521 -33.94756 24.894522 0.9975917
## 2:(+DH+Si2)-2:(-DH-Si ) 9.789520 -19.63152 39.210563 0.8699077
## 2:(+DH-Si )-1:(+DH-Si ) 5.241778 -24.17926 34.662821 0.9942219
## 1:(+DH+Si1)-1:(+DH-Si ) 6.937312 -22.48373 36.358355 0.9729955
## 2:(+DH+Si1)-1:(+DH-Si ) 1.974315 -27.44673 31.395357 0.9999894
## 1:(+DH+Si2)-1:(+DH-Si ) -14.160617 -43.58166 15.260426 0.5815408
## 2:(+DH+Si2)-1:(+DH-Si ) 0.155424 -29.26562 29.576467 1.0000000
## 1:(+DH+Si1)-2:(+DH-Si ) 1.695534 -27.72551 31.116577 0.9999963
## 2:(+DH+Si1)-2:(+DH-Si ) -3.267463 -32.68851 26.153579 0.9996927
## 1:(+DH+Si2)-2:(+DH-Si ) -19.402395 -48.82344 10.018648 0.2738955
## 2:(+DH+Si2)-2:(+DH-Si ) -5.086354 -34.50740 24.334689 0.9951558
## 2:(+DH+Si1)-1:(+DH+Si1) -4.962998 -34.38404 24.458045 0.9958101
## 1:(+DH+Si2)-1:(+DH+Si1) -21.097929 -50.51897 8.323114 0.2075416
## 2:(+DH+Si2)-1:(+DH+Si1) -6.781888 -36.20293 22.639154 0.9759793
## 1:(+DH+Si2)-2:(+DH+Si1) -16.134931 -45.55597 13.286112 0.4490599
## 2:(+DH+Si2)-2:(+DH+Si1) -1.818891 -31.23993 27.602152 0.9999939
## 2:(+DH+Si2)-1:(+DH+Si2) 14.316041 -15.10500 43.737083 0.5706423ggplot(datos_m2)+
aes(TRATAMIENTO,AREA_DOSEL, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m2)+
aes(x=TRATAMIENTO, y=TEMPERATURA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod2_ABD=aov(AREA_DOSEL~BLOQUES*TRATAMIENTO, data=datos_m2)
anova(mod2_ABD)## Analysis of Variance Table
##
## Response: AREA_DOSEL
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 1.3 1.3 0.0014 0.97094
## TRATAMIENTO 3 9524.8 3174.9 3.5741 0.06647 .
## BLOQUES:TRATAMIENTO 3 3423.4 1141.1 1.2846 0.34400
## Residuals 8 7106.5 888.3
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod2_ABD)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AREA_DOSEL ~ BLOQUES * TRATAMIENTO, data = datos_m2)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.56 -33.80466 34.92466 0.9709448
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -66.1925 -133.68198 1.296976 0.0545265
## (+DH+Si1)-(-DH-Si ) -44.4300 -111.91948 23.059476 0.2292001
## (+DH+Si2)-(-DH-Si ) -48.6975 -116.18698 18.791976 0.1745344
## (+DH+Si1)-(+DH-Si ) 21.7625 -45.72698 89.251976 0.7363382
## (+DH+Si2)-(+DH-Si ) 17.4950 -49.99448 84.984476 0.8389384
## (+DH+Si2)-(+DH+Si1) -4.2675 -71.75698 63.221976 0.9968158
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -41.600 -159.53934 76.33934 0.8374598
## 1:(+DH-Si )-1:(-DH-Si ) -87.140 -205.07934 30.79934 0.1864761
## 2:(+DH-Si )-1:(-DH-Si ) -86.845 -204.78434 31.09434 0.1887952
## 1:(+DH+Si1)-1:(-DH-Si ) -66.460 -184.39934 51.47934 0.4216377
## 2:(+DH+Si1)-1:(-DH-Si ) -64.000 -181.93934 53.93934 0.4597972
## 1:(+DH+Si2)-1:(-DH-Si ) -90.040 -207.97934 27.89934 0.1650385
## 2:(+DH+Si2)-1:(-DH-Si ) -48.955 -166.89434 68.98434 0.7191818
## 1:(+DH-Si )-2:(-DH-Si ) -45.540 -163.47934 72.39934 0.7766639
## 2:(+DH-Si )-2:(-DH-Si ) -45.245 -163.18434 72.69434 0.7814526
## 1:(+DH+Si1)-2:(-DH-Si ) -24.860 -142.79934 93.07934 0.9850800
## 2:(+DH+Si1)-2:(-DH-Si ) -22.400 -140.33934 95.53934 0.9916491
## 1:(+DH+Si2)-2:(-DH-Si ) -48.440 -166.37934 69.49934 0.7280527
## 2:(+DH+Si2)-2:(-DH-Si ) -7.355 -125.29434 110.58434 0.9999936
## 2:(+DH-Si )-1:(+DH-Si ) 0.295 -117.64434 118.23434 1.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) 20.680 -97.25934 138.61934 0.9947362
## 2:(+DH+Si1)-1:(+DH-Si ) 23.140 -94.79934 141.07934 0.9899643
## 1:(+DH+Si2)-1:(+DH-Si ) -2.900 -120.83934 115.03934 1.0000000
## 2:(+DH+Si2)-1:(+DH-Si ) 38.185 -79.75434 156.12434 0.8833446
## 1:(+DH+Si1)-2:(+DH-Si ) 20.385 -97.55434 138.32434 0.9951621
## 2:(+DH+Si1)-2:(+DH-Si ) 22.845 -95.09434 140.78434 0.9906639
## 1:(+DH+Si2)-2:(+DH-Si ) -3.195 -121.13434 114.74434 1.0000000
## 2:(+DH+Si2)-2:(+DH-Si ) 37.890 -80.04934 155.82934 0.8869588
## 2:(+DH+Si1)-1:(+DH+Si1) 2.460 -115.47934 120.39934 1.0000000
## 1:(+DH+Si2)-1:(+DH+Si1) -23.580 -141.51934 94.35934 0.9888490
## 2:(+DH+Si2)-1:(+DH+Si1) 17.505 -100.43434 135.44434 0.9980672
## 1:(+DH+Si2)-2:(+DH+Si1) -26.040 -143.97934 91.89934 0.9808505
## 2:(+DH+Si2)-2:(+DH+Si1) 15.045 -102.89434 132.98434 0.9992507
## 2:(+DH+Si2)-1:(+DH+Si2) 41.085 -76.85434 159.02434 0.8448294