FisioVeg
Juan Sebastian Herrera Reales
2023-06-04
Muestreo 3 analisis de datos
datos_m <- read_excel("C:/Users/JuanSebH2/Downloads/Universidad/Fisiologia/FisioVeg/DocumentoFinal_muestreos.xlsx",
sheet = "M3")
datos_m=data.frame(datos_m)
datos_m3=datos_m
datos_m3$BLOQUES=as.character(datos_m3$BLOQUES)
collapsibleTree(datos_m3, hierarchy = c("MUESTREO", "BLOQUES", "TRATAMIENTO"), collapsed = F)#Analisis descriptivo Y ESTADISTICO
ggplot(datos_m3)+
aes(TRATAMIENTO,TEMPERATURA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=TEMPERATURA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_T=aov(TEMPERATURA~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_T)## Analysis of Variance Table
##
## Response: TEMPERATURA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.6045 0.60451 6.0751 0.0390278 *
## TRATAMIENTO 3 7.9606 2.65354 26.6671 0.0001618 ***
## BLOQUES:TRATAMIENTO 3 0.6368 0.21227 2.1333 0.1741951
## Residuals 8 0.7961 0.09951
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_T)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TEMPERATURA ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.38875 -0.75246 -0.02503998 0.0390278
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 1.8850 1.1707023 2.5992977 0.0001358
## (+DH+Si1)-(-DH-Si ) 1.4925 0.7782023 2.2067977 0.0007005
## (+DH+Si2)-(-DH-Si ) 1.0000 0.2857023 1.7142977 0.0088303
## (+DH+Si1)-(+DH-Si ) -0.3925 -1.1067977 0.3217977 0.3567362
## (+DH+Si2)-(+DH-Si ) -0.8850 -1.5992977 -0.1707023 0.0173929
## (+DH+Si2)-(+DH+Si1) -0.4925 -1.2067977 0.2217977 0.2006146
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.680 -1.92825088 0.56825088 0.4558484
## 1:(+DH-Si )-1:(-DH-Si ) 1.520 0.27174912 2.76825088 0.0172665
## 2:(+DH-Si )-1:(-DH-Si ) 1.570 0.32174912 2.81825088 0.0143043
## 1:(+DH+Si1)-1:(-DH-Si ) 1.175 -0.07325088 2.42325088 0.0671929
## 2:(+DH+Si1)-1:(-DH-Si ) 1.130 -0.11825088 2.37825088 0.0806478
## 1:(+DH+Si2)-1:(-DH-Si ) 1.100 -0.14825088 2.34825088 0.0910999
## 2:(+DH+Si2)-1:(-DH-Si ) 0.220 -1.02825088 1.46825088 0.9945756
## 1:(+DH-Si )-2:(-DH-Si ) 2.200 0.95174912 3.44825088 0.0016734
## 2:(+DH-Si )-2:(-DH-Si ) 2.250 1.00174912 3.49825088 0.0014360
## 1:(+DH+Si1)-2:(-DH-Si ) 1.855 0.60674912 3.10325088 0.0051457
## 2:(+DH+Si1)-2:(-DH-Si ) 1.810 0.56174912 3.05825088 0.0060123
## 1:(+DH+Si2)-2:(-DH-Si ) 1.780 0.53174912 3.02825088 0.0066777
## 2:(+DH+Si2)-2:(-DH-Si ) 0.900 -0.34825088 2.14825088 0.2036072
## 2:(+DH-Si )-1:(+DH-Si ) 0.050 -1.19825088 1.29825088 0.9999997
## 1:(+DH+Si1)-1:(+DH-Si ) -0.345 -1.59325088 0.90325088 0.9413034
## 2:(+DH+Si1)-1:(+DH-Si ) -0.390 -1.63825088 0.85825088 0.8992400
## 1:(+DH+Si2)-1:(+DH-Si ) -0.420 -1.66825088 0.82825088 0.8641044
## 2:(+DH+Si2)-1:(+DH-Si ) -1.300 -2.54825088 -0.05174912 0.0406500
## 1:(+DH+Si1)-2:(+DH-Si ) -0.395 -1.64325088 0.85325088 0.8937616
## 2:(+DH+Si1)-2:(+DH-Si ) -0.440 -1.68825088 0.80825088 0.8378525
## 1:(+DH+Si2)-2:(+DH-Si ) -0.470 -1.71825088 0.77825088 0.7948624
## 2:(+DH+Si2)-2:(+DH-Si ) -1.350 -2.59825088 -0.10174912 0.0333398
## 2:(+DH+Si1)-1:(+DH+Si1) -0.045 -1.29325088 1.20325088 0.9999999
## 1:(+DH+Si2)-1:(+DH+Si1) -0.075 -1.32325088 1.17325088 0.9999950
## 2:(+DH+Si2)-1:(+DH+Si1) -0.955 -2.20325088 0.29325088 0.1637041
## 1:(+DH+Si2)-2:(+DH+Si1) -0.030 -1.27825088 1.21825088 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) -0.910 -2.15825088 0.33825088 0.1957470
## 2:(+DH+Si2)-1:(+DH+Si2) -0.880 -2.12825088 0.36825088 0.2201853ggplot(datos_m3)+
aes(TRATAMIENTO,N_HOJAS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=N_HOJAS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_NH=aov(N_HOJAS~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_NH)## Analysis of Variance Table
##
## Response: N_HOJAS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.25 0.2500 0.4000 0.54474
## TRATAMIENTO 3 6.50 2.1667 3.4667 0.07092 .
## BLOQUES:TRATAMIENTO 3 2.25 0.7500 1.2000 0.37003
## Residuals 8 5.00 0.6250
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_NH)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_HOJAS ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.25 -0.6615282 1.161528 0.5447373
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -1.75 -3.540169 0.04016919 0.0553252
## (+DH+Si1)-(-DH-Si ) -1.00 -2.790169 0.79016919 0.3443382
## (+DH+Si2)-(-DH-Si ) -1.25 -3.040169 0.54016919 0.1931571
## (+DH+Si1)-(+DH-Si ) 0.75 -1.040169 2.54016919 0.5647096
## (+DH+Si2)-(+DH-Si ) 0.50 -1.290169 2.29016919 0.8081467
## (+DH+Si2)-(+DH+Si1) -0.25 -2.040169 1.54016919 0.9683292
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 1.0 -2.12836 4.1283599 0.8892546
## 1:(+DH-Si )-1:(-DH-Si ) -1.0 -4.12836 2.1283599 0.8892546
## 2:(+DH-Si )-1:(-DH-Si ) -1.5 -4.62836 1.6283599 0.5853167
## 1:(+DH+Si1)-1:(-DH-Si ) -1.0 -4.12836 2.1283599 0.8892546
## 2:(+DH+Si1)-1:(-DH-Si ) 0.0 -3.12836 3.1283599 1.0000000
## 1:(+DH+Si2)-1:(-DH-Si ) -0.5 -3.62836 2.6283599 0.9969659
## 2:(+DH+Si2)-1:(-DH-Si ) -1.0 -4.12836 2.1283599 0.8892546
## 1:(+DH-Si )-2:(-DH-Si ) -2.0 -5.12836 1.1283599 0.3009915
## 2:(+DH-Si )-2:(-DH-Si ) -2.5 -5.62836 0.6283599 0.1380016
## 1:(+DH+Si1)-2:(-DH-Si ) -2.0 -5.12836 1.1283599 0.3009915
## 2:(+DH+Si1)-2:(-DH-Si ) -1.0 -4.12836 2.1283599 0.8892546
## 1:(+DH+Si2)-2:(-DH-Si ) -1.5 -4.62836 1.6283599 0.5853167
## 2:(+DH+Si2)-2:(-DH-Si ) -2.0 -5.12836 1.1283599 0.3009915
## 2:(+DH-Si )-1:(+DH-Si ) -0.5 -3.62836 2.6283599 0.9969659
## 1:(+DH+Si1)-1:(+DH-Si ) 0.0 -3.12836 3.1283599 1.0000000
## 2:(+DH+Si1)-1:(+DH-Si ) 1.0 -2.12836 4.1283599 0.8892546
## 1:(+DH+Si2)-1:(+DH-Si ) 0.5 -2.62836 3.6283599 0.9969659
## 2:(+DH+Si2)-1:(+DH-Si ) 0.0 -3.12836 3.1283599 1.0000000
## 1:(+DH+Si1)-2:(+DH-Si ) 0.5 -2.62836 3.6283599 0.9969659
## 2:(+DH+Si1)-2:(+DH-Si ) 1.5 -1.62836 4.6283599 0.5853167
## 1:(+DH+Si2)-2:(+DH-Si ) 1.0 -2.12836 4.1283599 0.8892546
## 2:(+DH+Si2)-2:(+DH-Si ) 0.5 -2.62836 3.6283599 0.9969659
## 2:(+DH+Si1)-1:(+DH+Si1) 1.0 -2.12836 4.1283599 0.8892546
## 1:(+DH+Si2)-1:(+DH+Si1) 0.5 -2.62836 3.6283599 0.9969659
## 2:(+DH+Si2)-1:(+DH+Si1) 0.0 -3.12836 3.1283599 1.0000000
## 1:(+DH+Si2)-2:(+DH+Si1) -0.5 -3.62836 2.6283599 0.9969659
## 2:(+DH+Si2)-2:(+DH+Si1) -1.0 -4.12836 2.1283599 0.8892546
## 2:(+DH+Si2)-1:(+DH+Si2) -0.5 -3.62836 2.6283599 0.9969659ggplot(datos_m3)+
aes(TRATAMIENTO,LONGITUD_P_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=LONGITUD_P_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PA=aov(LONGITUD_P_AEREA~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PA)## Analysis of Variance Table
##
## Response: LONGITUD_P_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 3.423 3.4225 1.6975 0.228870
## TRATAMIENTO 3 72.585 24.1950 12.0000 0.002486 **
## BLOQUES:TRATAMIENTO 3 1.953 0.6508 0.3228 0.809017
## Residuals 8 16.130 2.0163
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = LONGITUD_P_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.925 -2.562202 0.712202 0.2288698
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -5.175 -8.390335 -1.959665 0.0038382
## (+DH+Si1)-(-DH-Si ) -4.875 -8.090335 -1.659665 0.0055238
## (+DH+Si2)-(-DH-Si ) -4.650 -7.865335 -1.434665 0.0073120
## (+DH+Si1)-(+DH-Si ) 0.300 -2.915335 3.515335 0.9900180
## (+DH+Si2)-(+DH-Si ) 0.525 -2.690335 3.740335 0.9511666
## (+DH+Si2)-(+DH+Si1) 0.225 -2.990335 3.440335 0.9957045
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.50 -6.118869 5.11886879 0.9999287
## 1:(+DH-Si )-1:(-DH-Si ) -5.35 -10.968869 0.26886879 0.0636145
## 2:(+DH-Si )-1:(-DH-Si ) -5.50 -11.118869 0.11886879 0.0556052
## 1:(+DH+Si1)-1:(-DH-Si ) -4.60 -10.218869 1.01886879 0.1250610
## 2:(+DH+Si1)-1:(-DH-Si ) -5.65 -11.268869 -0.03113121 0.0486309
## 1:(+DH+Si2)-1:(-DH-Si ) -3.90 -9.518869 1.71886879 0.2321341
## 2:(+DH+Si2)-1:(-DH-Si ) -5.90 -11.518869 -0.28113121 0.0389546
## 1:(+DH-Si )-2:(-DH-Si ) -4.85 -10.468869 0.76886879 0.0998420
## 2:(+DH-Si )-2:(-DH-Si ) -5.00 -10.618869 0.61886879 0.0872007
## 1:(+DH+Si1)-2:(-DH-Si ) -4.10 -9.718869 1.51886879 0.1951079
## 2:(+DH+Si1)-2:(-DH-Si ) -5.15 -10.768869 0.46886879 0.0761642
## 1:(+DH+Si2)-2:(-DH-Si ) -3.40 -9.018869 2.21886879 0.3517871
## 2:(+DH+Si2)-2:(-DH-Si ) -5.40 -11.018869 0.21886879 0.0608209
## 2:(+DH-Si )-1:(+DH-Si ) -0.15 -5.768869 5.46886879 1.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) 0.75 -4.868869 6.36886879 0.9990019
## 2:(+DH+Si1)-1:(+DH-Si ) -0.30 -5.918869 5.31886879 0.9999978
## 1:(+DH+Si2)-1:(+DH-Si ) 1.45 -4.168869 7.06886879 0.9576083
## 2:(+DH+Si2)-1:(+DH-Si ) -0.55 -6.168869 5.06886879 0.9998657
## 1:(+DH+Si1)-2:(+DH-Si ) 0.90 -4.718869 6.51886879 0.9969261
## 2:(+DH+Si1)-2:(+DH-Si ) -0.15 -5.768869 5.46886879 1.0000000
## 1:(+DH+Si2)-2:(+DH-Si ) 1.60 -4.018869 7.21886879 0.9327215
## 2:(+DH+Si2)-2:(+DH-Si ) -0.40 -6.018869 5.21886879 0.9999841
## 2:(+DH+Si1)-1:(+DH+Si1) -1.05 -6.668869 4.56886879 0.9923881
## 1:(+DH+Si2)-1:(+DH+Si1) 0.70 -4.918869 6.31886879 0.9993556
## 2:(+DH+Si2)-1:(+DH+Si1) -1.30 -6.918869 4.31886879 0.9755146
## 1:(+DH+Si2)-2:(+DH+Si1) 1.75 -3.868869 7.36886879 0.9005654
## 2:(+DH+Si2)-2:(+DH+Si1) -0.25 -5.868869 5.36886879 0.9999994
## 2:(+DH+Si2)-1:(+DH+Si2) -2.00 -7.618869 3.61886879 0.8319315ggplot(datos_m3)+
aes(TRATAMIENTO,AREA_FOLIAR, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=AREA_FOLIAR, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_AF=aov(AREA_FOLIAR~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_AF)## Analysis of Variance Table
##
## Response: AREA_FOLIAR
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 416.4 416.4 0.8499 0.3835315
## TRATAMIENTO 3 29179.5 9726.5 19.8515 0.0004605 ***
## BLOQUES:TRATAMIENTO 3 3301.2 1100.4 2.2459 0.1602476
## Residuals 8 3919.7 490.0
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_AF)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AREA_FOLIAR ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -10.20337 -35.72517 15.31842 0.3835315
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -110.77525 -160.89803 -60.65247 0.0004763
## (+DH+Si1)-(-DH-Si ) -87.19100 -137.31378 -37.06822 0.0023541
## (+DH+Si2)-(-DH-Si ) -91.41800 -141.54078 -41.29522 0.0017344
## (+DH+Si1)-(+DH-Si ) 23.58425 -26.53853 73.70703 0.4765504
## (+DH+Si2)-(+DH-Si ) 19.35725 -30.76553 69.48003 0.6230442
## (+DH+Si2)-(+DH+Si1) -4.22700 -54.34978 45.89578 0.9925678
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 24.1500 -63.44066 111.7406603 0.9419618
## 1:(+DH-Si )-1:(-DH-Si ) -86.9215 -174.51216 0.6691603 0.0519536
## 2:(+DH-Si )-1:(-DH-Si ) -110.4790 -198.06966 -22.8883397 0.0140694
## 1:(+DH+Si1)-1:(-DH-Si ) -79.4370 -167.02766 8.1536603 0.0799781
## 2:(+DH+Si1)-1:(-DH-Si ) -70.7950 -158.38566 16.7956603 0.1318134
## 1:(+DH+Si2)-1:(-DH-Si ) -54.3190 -141.90966 33.2716603 0.3287009
## 2:(+DH+Si2)-1:(-DH-Si ) -104.3670 -191.95766 -16.7763397 0.0195525
## 1:(+DH-Si )-2:(-DH-Si ) -111.0715 -198.66216 -23.4808397 0.0136332
## 2:(+DH-Si )-2:(-DH-Si ) -134.6290 -222.21966 -47.0383397 0.0041448
## 1:(+DH+Si1)-2:(-DH-Si ) -103.5870 -191.17766 -15.9963397 0.0204023
## 2:(+DH+Si1)-2:(-DH-Si ) -94.9450 -182.53566 -7.3543397 0.0329401
## 1:(+DH+Si2)-2:(-DH-Si ) -78.4690 -166.05966 9.1216603 0.0845869
## 2:(+DH+Si2)-2:(-DH-Si ) -128.5170 -216.10766 -40.9263397 0.0055801
## 2:(+DH-Si )-1:(+DH-Si ) -23.5575 -111.14816 64.0331603 0.9483287
## 1:(+DH+Si1)-1:(+DH-Si ) 7.4845 -80.10616 95.0751603 0.9999457
## 2:(+DH+Si1)-1:(+DH-Si ) 16.1265 -71.46416 103.7171603 0.9930112
## 1:(+DH+Si2)-1:(+DH-Si ) 32.6025 -54.98816 120.1931603 0.8028598
## 2:(+DH+Si2)-1:(+DH-Si ) -17.4455 -105.03616 70.1451603 0.9890841
## 1:(+DH+Si1)-2:(+DH-Si ) 31.0420 -56.54866 118.6326603 0.8345940
## 2:(+DH+Si1)-2:(+DH-Si ) 39.6840 -47.90666 127.2746603 0.6404138
## 1:(+DH+Si2)-2:(+DH-Si ) 56.1600 -31.43066 143.7506603 0.2984157
## 2:(+DH+Si2)-2:(+DH-Si ) 6.1120 -81.47866 93.7026603 0.9999862
## 2:(+DH+Si1)-1:(+DH+Si1) 8.6420 -78.94866 96.2326603 0.9998585
## 1:(+DH+Si2)-1:(+DH+Si1) 25.1180 -62.47266 112.7086603 0.9305521
## 2:(+DH+Si2)-1:(+DH+Si1) -24.9300 -112.52066 62.6606603 0.9328662
## 1:(+DH+Si2)-2:(+DH+Si1) 16.4760 -71.11466 104.0666603 0.9920966
## 2:(+DH+Si2)-2:(+DH+Si1) -33.5720 -121.16266 54.0186603 0.7821135
## 2:(+DH+Si2)-1:(+DH+Si2) -50.0480 -137.63866 37.5426603 0.4077544ggplot(datos_m3)+
aes(TRATAMIENTO,LONGITUD_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=LONGITUD_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_LR=aov(LONGITUD_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_LR)## Analysis of Variance Table
##
## Response: LONGITUD_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 1.3225 1.3225 1.6480 0.2352
## TRATAMIENTO 3 3.8475 1.2825 1.5981 0.2646
## BLOQUES:TRATAMIENTO 3 4.1275 1.3758 1.7144 0.2409
## Residuals 8 6.4200 0.8025TukeyHSD(mod3_LR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = LONGITUD_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.575 -1.607887 0.4578865 0.2351662
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -1.25 -3.278507 0.7785074 0.2732823
## (+DH+Si1)-(-DH-Si ) -0.95 -2.978507 1.0785074 0.4801861
## (+DH+Si2)-(-DH-Si ) -0.35 -2.378507 1.6785074 0.9432600
## (+DH+Si1)-(+DH-Si ) 0.30 -1.728507 2.3285074 0.9628326
## (+DH+Si2)-(+DH-Si ) 0.90 -1.128507 2.9285074 0.5217059
## (+DH+Si2)-(+DH+Si1) 0.60 -1.428507 2.6285074 0.7815028
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.00 -3.544861 3.544861 1.0000000
## 1:(+DH-Si )-1:(-DH-Si ) -1.65 -5.194861 1.894861 0.6144862
## 2:(+DH-Si )-1:(-DH-Si ) -0.85 -4.394861 2.694861 0.9705781
## 1:(+DH+Si1)-1:(-DH-Si ) -0.20 -3.744861 3.344861 0.9999968
## 2:(+DH+Si1)-1:(-DH-Si ) -1.70 -5.244861 1.844861 0.5851448
## 1:(+DH+Si2)-1:(-DH-Si ) 0.45 -3.094861 3.994861 0.9992736
## 2:(+DH+Si2)-1:(-DH-Si ) -1.15 -4.694861 2.394861 0.8824029
## 1:(+DH-Si )-2:(-DH-Si ) -1.65 -5.194861 1.894861 0.6144862
## 2:(+DH-Si )-2:(-DH-Si ) -0.85 -4.394861 2.694861 0.9705781
## 1:(+DH+Si1)-2:(-DH-Si ) -0.20 -3.744861 3.344861 0.9999968
## 2:(+DH+Si1)-2:(-DH-Si ) -1.70 -5.244861 1.844861 0.5851448
## 1:(+DH+Si2)-2:(-DH-Si ) 0.45 -3.094861 3.994861 0.9992736
## 2:(+DH+Si2)-2:(-DH-Si ) -1.15 -4.694861 2.394861 0.8824029
## 2:(+DH-Si )-1:(+DH-Si ) 0.80 -2.744861 4.344861 0.9784990
## 1:(+DH+Si1)-1:(+DH-Si ) 1.45 -2.094861 4.994861 0.7314431
## 2:(+DH+Si1)-1:(+DH-Si ) -0.05 -3.594861 3.494861 1.0000000
## 1:(+DH+Si2)-1:(+DH-Si ) 2.10 -1.444861 5.644861 0.3721859
## 2:(+DH+Si2)-1:(+DH-Si ) 0.50 -3.044861 4.044861 0.9985900
## 1:(+DH+Si1)-2:(+DH-Si ) 0.65 -2.894861 4.194861 0.9931732
## 2:(+DH+Si1)-2:(+DH-Si ) -0.85 -4.394861 2.694861 0.9705781
## 1:(+DH+Si2)-2:(+DH-Si ) 1.30 -2.244861 4.844861 0.8128642
## 2:(+DH+Si2)-2:(+DH-Si ) -0.30 -3.844861 3.244861 0.9999490
## 2:(+DH+Si1)-1:(+DH+Si1) -1.50 -5.044861 2.044861 0.7026653
## 1:(+DH+Si2)-1:(+DH+Si1) 0.65 -2.894861 4.194861 0.9931732
## 2:(+DH+Si2)-1:(+DH+Si1) -0.95 -4.494861 2.594861 0.9491913
## 1:(+DH+Si2)-2:(+DH+Si1) 2.15 -1.394861 5.694861 0.3495766
## 2:(+DH+Si2)-2:(+DH+Si1) 0.55 -2.994861 4.094861 0.9974650
## 2:(+DH+Si2)-1:(+DH+Si2) -1.60 -5.144861 1.944861 0.6439798ggplot(datos_m3)+
aes(TRATAMIENTO,DIAMETRO_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=DIAMETRO_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_DR=aov(DIAMETRO_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_DR)## Analysis of Variance Table
##
## Response: DIAMETRO_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.397 0.397 0.3084 0.5938
## TRATAMIENTO 3 263.750 87.917 68.3239 4.826e-06 ***
## BLOQUES:TRATAMIENTO 3 0.449 0.150 0.1163 0.9480
## Residuals 8 10.294 1.287
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_DR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = DIAMETRO_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.315 -1.622914 0.9929143 0.5938195
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -11.345 -13.91364 -8.77636 0.0000029
## (+DH+Si1)-(-DH-Si ) -6.775 -9.34364 -4.20636 0.0001363
## (+DH+Si2)-(-DH-Si ) -7.055 -9.62364 -4.48636 0.0001016
## (+DH+Si1)-(+DH-Si ) 4.570 2.00136 7.13864 0.0020372
## (+DH+Si2)-(+DH-Si ) 4.290 1.72136 6.85864 0.0030475
## (+DH+Si2)-(+DH+Si1) -0.280 -2.84864 2.28864 0.9843407
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.150 -4.638755 4.338755 0.9999999
## 1:(+DH-Si )-1:(-DH-Si ) -11.075 -15.563755 -6.586245 0.0001543
## 2:(+DH-Si )-1:(-DH-Si ) -11.765 -16.253755 -7.276245 0.0000987
## 1:(+DH+Si1)-1:(-DH-Si ) -6.925 -11.413755 -2.436245 0.0040465
## 2:(+DH+Si1)-1:(-DH-Si ) -6.775 -11.263755 -2.286245 0.0046596
## 1:(+DH+Si2)-1:(-DH-Si ) -6.845 -11.333755 -2.356245 0.0043617
## 2:(+DH+Si2)-1:(-DH-Si ) -7.415 -11.903755 -2.926245 0.0025848
## 1:(+DH-Si )-2:(-DH-Si ) -10.925 -15.413755 -6.436245 0.0001705
## 2:(+DH-Si )-2:(-DH-Si ) -11.615 -16.103755 -7.126245 0.0001086
## 1:(+DH+Si1)-2:(-DH-Si ) -6.775 -11.263755 -2.286245 0.0046596
## 2:(+DH+Si1)-2:(-DH-Si ) -6.625 -11.113755 -2.136245 0.0053754
## 1:(+DH+Si2)-2:(-DH-Si ) -6.695 -11.183755 -2.206245 0.0050275
## 2:(+DH+Si2)-2:(-DH-Si ) -7.265 -11.753755 -2.776245 0.0029589
## 2:(+DH-Si )-1:(+DH-Si ) -0.690 -5.178755 3.798755 0.9976046
## 1:(+DH+Si1)-1:(+DH-Si ) 4.150 -0.338755 8.638755 0.0731482
## 2:(+DH+Si1)-1:(+DH-Si ) 4.300 -0.188755 8.788755 0.0617790
## 1:(+DH+Si2)-1:(+DH-Si ) 4.230 -0.258755 8.718755 0.0668414
## 2:(+DH+Si2)-1:(+DH-Si ) 3.660 -0.828755 8.148755 0.1271615
## 1:(+DH+Si1)-2:(+DH-Si ) 4.840 0.351245 9.328755 0.0338802
## 2:(+DH+Si1)-2:(+DH-Si ) 4.990 0.501245 9.478755 0.0287553
## 1:(+DH+Si2)-2:(+DH-Si ) 4.920 0.431245 9.408755 0.0310370
## 2:(+DH+Si2)-2:(+DH-Si ) 4.350 -0.138755 8.838755 0.0584048
## 2:(+DH+Si1)-1:(+DH+Si1) 0.150 -4.338755 4.638755 0.9999999
## 1:(+DH+Si2)-1:(+DH+Si1) 0.080 -4.408755 4.568755 1.0000000
## 2:(+DH+Si2)-1:(+DH+Si1) -0.490 -4.978755 3.998755 0.9997253
## 1:(+DH+Si2)-2:(+DH+Si1) -0.070 -4.558755 4.418755 1.0000000
## 2:(+DH+Si2)-2:(+DH+Si1) -0.640 -5.128755 3.848755 0.9984922
## 2:(+DH+Si2)-1:(+DH+Si2) -0.570 -5.058755 3.918755 0.9992722ggplot(datos_m3)+
aes(TRATAMIENTO,PESO_F_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=PESO_F_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PFR=aov(PESO_F_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PFR)## Analysis of Variance Table
##
## Response: PESO_F_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.78 0.78 0.0684 0.8003
## TRATAMIENTO 3 1096.16 365.39 31.8999 8.443e-05 ***
## BLOQUES:TRATAMIENTO 3 13.62 4.54 0.3962 0.7594
## Residuals 8 91.63 11.45
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PFR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 0.4425 -3.459712 4.344712 0.8003175
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -21.1150 -28.778635 -13.451365 0.0000993
## (+DH+Si1)-(-DH-Si ) -16.6075 -24.271135 -8.943865 0.0005456
## (+DH+Si2)-(-DH-Si ) -18.5425 -26.206135 -10.878865 0.0002526
## (+DH+Si1)-(+DH-Si ) 4.5075 -3.156135 12.171135 0.3063019
## (+DH+Si2)-(+DH-Si ) 2.5725 -5.091135 10.236135 0.7132583
## (+DH+Si2)-(+DH+Si1) -1.9350 -9.598635 5.728635 0.8487976
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 3.350 -10.042369 16.7423695 0.9636177
## 1:(+DH-Si )-1:(-DH-Si ) -19.360 -32.752369 -5.9676305 0.0061293
## 2:(+DH-Si )-1:(-DH-Si ) -19.520 -32.912369 -6.1276305 0.0058192
## 1:(+DH+Si1)-1:(-DH-Si ) -14.060 -27.452369 -0.6676305 0.0389898
## 2:(+DH+Si1)-1:(-DH-Si ) -15.805 -29.197369 -2.4126305 0.0206456
## 1:(+DH+Si2)-1:(-DH-Si ) -17.030 -30.422369 -3.6376305 0.0134105
## 2:(+DH+Si2)-1:(-DH-Si ) -16.705 -30.097369 -3.3126305 0.0150175
## 1:(+DH-Si )-2:(-DH-Si ) -22.710 -36.102369 -9.3176305 0.0021706
## 2:(+DH-Si )-2:(-DH-Si ) -22.870 -36.262369 -9.4776305 0.0020708
## 1:(+DH+Si1)-2:(-DH-Si ) -17.410 -30.802369 -4.0176305 0.0117625
## 2:(+DH+Si1)-2:(-DH-Si ) -19.155 -32.547369 -5.7626305 0.0065531
## 1:(+DH+Si2)-2:(-DH-Si ) -20.380 -33.772369 -6.9876305 0.0044203
## 2:(+DH+Si2)-2:(-DH-Si ) -20.055 -33.447369 -6.6626305 0.0049004
## 2:(+DH-Si )-1:(+DH-Si ) -0.160 -13.552369 13.2323695 1.0000000
## 1:(+DH+Si1)-1:(+DH-Si ) 5.300 -8.092369 18.6923695 0.7579600
## 2:(+DH+Si1)-1:(+DH-Si ) 3.555 -9.837369 16.9473695 0.9514408
## 1:(+DH+Si2)-1:(+DH-Si ) 2.330 -11.062369 15.7223695 0.9949720
## 2:(+DH+Si2)-1:(+DH-Si ) 2.655 -10.737369 16.0473695 0.9893639
## 1:(+DH+Si1)-2:(+DH-Si ) 5.460 -7.932369 18.8523695 0.7341619
## 2:(+DH+Si1)-2:(+DH-Si ) 3.715 -9.677369 17.1073695 0.9403057
## 1:(+DH+Si2)-2:(+DH-Si ) 2.490 -10.902369 15.8823695 0.9926058
## 2:(+DH+Si2)-2:(+DH-Si ) 2.815 -10.577369 16.2073695 0.9853067
## 2:(+DH+Si1)-1:(+DH+Si1) -1.745 -15.137369 11.6473695 0.9991429
## 1:(+DH+Si2)-1:(+DH+Si1) -2.970 -16.362369 10.4223695 0.9803961
## 2:(+DH+Si2)-1:(+DH+Si1) -2.645 -16.037369 10.7473695 0.9895862
## 1:(+DH+Si2)-2:(+DH+Si1) -1.225 -14.617369 12.1673695 0.9999144
## 2:(+DH+Si2)-2:(+DH+Si1) -0.900 -14.292369 12.4923695 0.9999893
## 2:(+DH+Si2)-1:(+DH+Si2) 0.325 -13.067369 13.7173695 1.0000000ggplot(datos_m3)+
aes(TRATAMIENTO,PESO_F_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=PESO_F_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PFA=aov(PESO_F_AEREA~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PFA)## Analysis of Variance Table
##
## Response: PESO_F_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 1.607 1.6066 0.8876 0.373705
## TRATAMIENTO 3 50.374 16.7913 9.2769 0.005541 **
## BLOQUES:TRATAMIENTO 3 3.842 1.2806 0.7075 0.574037
## Residuals 8 14.480 1.8100
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PFA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.63375 -2.184958 0.9174583 0.3737049
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -4.6750 -7.72145 -1.6285496 0.0051365
## (+DH+Si1)-(-DH-Si ) -3.6750 -6.72145 -0.6285496 0.0200303
## (+DH+Si2)-(-DH-Si ) -3.5775 -6.62395 -0.5310496 0.0230278
## (+DH+Si1)-(+DH-Si ) 1.0000 -2.04645 4.0464504 0.7262690
## (+DH+Si2)-(+DH-Si ) 1.0975 -1.94895 4.1439504 0.6695844
## (+DH+Si2)-(+DH+Si1) 0.0975 -2.94895 3.1439504 0.9995806
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.615 -4.708739 5.9387388 0.9996033
## 1:(+DH-Si )-1:(-DH-Si ) -3.780 -9.103739 1.5437388 0.2148441
## 2:(+DH-Si )-1:(-DH-Si ) -4.955 -10.278739 0.3687388 0.0708867
## 1:(+DH+Si1)-1:(-DH-Si ) -3.340 -8.663739 1.9837388 0.3180200
## 2:(+DH+Si1)-1:(-DH-Si ) -3.395 -8.718739 1.9287388 0.3032357
## 1:(+DH+Si2)-1:(-DH-Si ) -2.310 -7.633739 3.0137388 0.6804040
## 2:(+DH+Si2)-1:(-DH-Si ) -4.230 -9.553739 1.0937388 0.1412226
## 1:(+DH-Si )-2:(-DH-Si ) -4.395 -9.718739 0.9287388 0.1207902
## 2:(+DH-Si )-2:(-DH-Si ) -5.570 -10.893739 -0.2462612 0.0396919
## 1:(+DH+Si1)-2:(-DH-Si ) -3.955 -9.278739 1.3687388 0.1827814
## 2:(+DH+Si1)-2:(-DH-Si ) -4.010 -9.333739 1.3137388 0.1736486
## 1:(+DH+Si2)-2:(-DH-Si ) -2.925 -8.248739 2.3987388 0.4471907
## 2:(+DH+Si2)-2:(-DH-Si ) -4.845 -10.168739 0.4787388 0.0787063
## 2:(+DH-Si )-1:(+DH-Si ) -1.175 -6.498739 4.1487388 0.9808883
## 1:(+DH+Si1)-1:(+DH-Si ) 0.440 -4.883739 5.7637388 0.9999565
## 2:(+DH+Si1)-1:(+DH-Si ) 0.385 -4.938739 5.7087388 0.9999824
## 1:(+DH+Si2)-1:(+DH-Si ) 1.470 -3.853739 6.7937388 0.9415635
## 2:(+DH+Si2)-1:(+DH-Si ) -0.450 -5.773739 4.8737388 0.9999495
## 1:(+DH+Si1)-2:(+DH-Si ) 1.615 -3.708739 6.9387388 0.9110731
## 2:(+DH+Si1)-2:(+DH-Si ) 1.560 -3.763739 6.8837388 0.9235313
## 1:(+DH+Si2)-2:(+DH-Si ) 2.645 -2.678739 7.9687388 0.5496843
## 2:(+DH+Si2)-2:(+DH-Si ) 0.725 -4.598739 6.0487388 0.9988678
## 2:(+DH+Si1)-1:(+DH+Si1) -0.055 -5.378739 5.2687388 1.0000000
## 1:(+DH+Si2)-1:(+DH+Si1) 1.030 -4.293739 6.3537388 0.9907259
## 2:(+DH+Si2)-1:(+DH+Si1) -0.890 -6.213739 4.4337388 0.9960285
## 1:(+DH+Si2)-2:(+DH+Si1) 1.085 -4.238739 6.4087388 0.9875968
## 2:(+DH+Si2)-2:(+DH+Si1) -0.835 -6.158739 4.4887388 0.9972926
## 2:(+DH+Si2)-1:(+DH+Si2) -1.920 -7.243739 3.4037388 0.8237141ggplot(datos_m3)+
aes(TRATAMIENTO,PESO_S_RAIZ, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=PESO_S_RAIZ, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PSR=aov(PESO_S_RAIZ~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PSR)## Analysis of Variance Table
##
## Response: PESO_S_RAIZ
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.042436 0.042436 1.8273 0.21343
## TRATAMIENTO 3 0.307491 0.102497 4.4134 0.04136 *
## BLOQUES:TRATAMIENTO 3 0.076790 0.025597 1.1022 0.40301
## Residuals 8 0.185791 0.023224
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PSR)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_RAIZ ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.103 -0.2787104 0.07271038 0.2134265
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.3820 -0.72708129 -0.03691871 0.0309802
## (+DH+Si1)-(-DH-Si ) -0.1690 -0.51408129 0.17608129 0.4454822
## (+DH+Si2)-(-DH-Si ) -0.1145 -0.45958129 0.23058129 0.7200536
## (+DH+Si1)-(+DH-Si ) 0.2130 -0.13208129 0.55808129 0.2721271
## (+DH+Si2)-(+DH-Si ) 0.2675 -0.07758129 0.61258129 0.1377472
## (+DH+Si2)-(+DH+Si1) 0.0545 -0.29058129 0.39958129 0.9554278
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -0.2110 -0.8140371 0.3920371 0.8422500
## 1:(+DH-Si )-1:(-DH-Si ) -0.5025 -1.1055371 0.1005371 0.1161571
## 2:(+DH-Si )-1:(-DH-Si ) -0.4725 -1.0755371 0.1305371 0.1492590
## 1:(+DH+Si1)-1:(-DH-Si ) -0.1400 -0.7430371 0.4630371 0.9750753
## 2:(+DH+Si1)-1:(-DH-Si ) -0.4090 -1.0120371 0.1940371 0.2505319
## 1:(+DH+Si2)-1:(-DH-Si ) -0.2390 -0.8420371 0.3640371 0.7568193
## 2:(+DH+Si2)-1:(-DH-Si ) -0.2010 -0.8040371 0.4020371 0.8690230
## 1:(+DH-Si )-2:(-DH-Si ) -0.2915 -0.8945371 0.3115371 0.5772489
## 2:(+DH-Si )-2:(-DH-Si ) -0.2615 -0.8645371 0.3415371 0.6809595
## 1:(+DH+Si1)-2:(-DH-Si ) 0.0710 -0.5320371 0.6740371 0.9995515
## 2:(+DH+Si1)-2:(-DH-Si ) -0.1980 -0.8010371 0.4050371 0.8765912
## 1:(+DH+Si2)-2:(-DH-Si ) -0.0280 -0.6310371 0.5750371 0.9999992
## 2:(+DH+Si2)-2:(-DH-Si ) 0.0100 -0.5930371 0.6130371 1.0000000
## 2:(+DH-Si )-1:(+DH-Si ) 0.0300 -0.5730371 0.6330371 0.9999986
## 1:(+DH+Si1)-1:(+DH-Si ) 0.3625 -0.2405371 0.9655371 0.3580921
## 2:(+DH+Si1)-1:(+DH-Si ) 0.0935 -0.5095371 0.6965371 0.9974755
## 1:(+DH+Si2)-1:(+DH-Si ) 0.2635 -0.3395371 0.8665371 0.6740600
## 2:(+DH+Si2)-1:(+DH-Si ) 0.3015 -0.3015371 0.9045371 0.5433122
## 1:(+DH+Si1)-2:(+DH-Si ) 0.3325 -0.2705371 0.9355371 0.4435994
## 2:(+DH+Si1)-2:(+DH-Si ) 0.0635 -0.5395371 0.6665371 0.9997829
## 1:(+DH+Si2)-2:(+DH-Si ) 0.2335 -0.3695371 0.8365371 0.7745938
## 2:(+DH+Si2)-2:(+DH-Si ) 0.2715 -0.3315371 0.8745371 0.6463575
## 2:(+DH+Si1)-1:(+DH+Si1) -0.2690 -0.8720371 0.3340371 0.6550269
## 1:(+DH+Si2)-1:(+DH+Si1) -0.0990 -0.7020371 0.5040371 0.9964361
## 2:(+DH+Si2)-1:(+DH+Si1) -0.0610 -0.6640371 0.5420371 0.9998333
## 1:(+DH+Si2)-2:(+DH+Si1) 0.1700 -0.4330371 0.7730371 0.9357235
## 2:(+DH+Si2)-2:(+DH+Si1) 0.2080 -0.3950371 0.8110371 0.8505223
## 2:(+DH+Si2)-1:(+DH+Si2) 0.0380 -0.5650371 0.6410371 0.9999931ggplot(datos_m3)+
aes(TRATAMIENTO,PESO_S_AEREA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=PESO_S_AEREA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PSA=aov(PESO_S_AEREA~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PSA)## Analysis of Variance Table
##
## Response: PESO_S_AEREA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.015563 0.015563 0.5219 0.49062
## TRATAMIENTO 3 0.266386 0.088795 2.9777 0.09652 .
## BLOQUES:TRATAMIENTO 3 0.039019 0.013006 0.4362 0.73319
## Residuals 8 0.238563 0.029820
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PSA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_AEREA ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.062375 -0.2614823 0.1367323 0.4906182
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.35025 -0.7412809 0.0407809 0.0801089
## (+DH+Si1)-(-DH-Si ) -0.26200 -0.6530309 0.1290309 0.2180795
## (+DH+Si2)-(-DH-Si ) -0.22150 -0.6125309 0.1695309 0.3339199
## (+DH+Si1)-(+DH-Si ) 0.08825 -0.3027809 0.4792809 0.8853723
## (+DH+Si2)-(+DH-Si ) 0.12875 -0.2622809 0.5197809 0.7245111
## (+DH+Si2)-(+DH+Si1) 0.04050 -0.3505309 0.4315309 0.9864878
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0360 -0.6473351 0.7193351 0.9999980
## 1:(+DH-Si )-1:(-DH-Si ) -0.3200 -1.0033351 0.3633351 0.6085835
## 2:(+DH-Si )-1:(-DH-Si ) -0.3445 -1.0278351 0.3388351 0.5348736
## 1:(+DH+Si1)-1:(-DH-Si ) -0.2270 -0.9103351 0.4563351 0.8707440
## 2:(+DH+Si1)-1:(-DH-Si ) -0.2610 -0.9443351 0.4223351 0.7846486
## 1:(+DH+Si2)-1:(-DH-Si ) -0.0900 -0.7733351 0.5933351 0.9990826
## 2:(+DH+Si2)-1:(-DH-Si ) -0.3170 -1.0003351 0.3663351 0.6177466
## 1:(+DH-Si )-2:(-DH-Si ) -0.3560 -1.0393351 0.3273351 0.5014043
## 2:(+DH-Si )-2:(-DH-Si ) -0.3805 -1.0638351 0.3028351 0.4336479
## 1:(+DH+Si1)-2:(-DH-Si ) -0.2630 -0.9463351 0.4203351 0.7790681
## 2:(+DH+Si1)-2:(-DH-Si ) -0.2970 -0.9803351 0.3863351 0.6788917
## 1:(+DH+Si2)-2:(-DH-Si ) -0.1260 -0.8093351 0.5573351 0.9929503
## 2:(+DH+Si2)-2:(-DH-Si ) -0.3530 -1.0363351 0.3303351 0.5100457
## 2:(+DH-Si )-1:(+DH-Si ) -0.0245 -0.7078351 0.6588351 0.9999999
## 1:(+DH+Si1)-1:(+DH-Si ) 0.0930 -0.5903351 0.7763351 0.9988723
## 2:(+DH+Si1)-1:(+DH-Si ) 0.0590 -0.6243351 0.7423351 0.9999417
## 1:(+DH+Si2)-1:(+DH-Si ) 0.2300 -0.4533351 0.9133351 0.8639257
## 2:(+DH+Si2)-1:(+DH-Si ) 0.0030 -0.6803351 0.6863351 1.0000000
## 1:(+DH+Si1)-2:(+DH-Si ) 0.1175 -0.5658351 0.8008351 0.9953074
## 2:(+DH+Si1)-2:(+DH-Si ) 0.0835 -0.5998351 0.7668351 0.9994306
## 1:(+DH+Si2)-2:(+DH-Si ) 0.2545 -0.4288351 0.9378351 0.8024477
## 2:(+DH+Si2)-2:(+DH-Si ) 0.0275 -0.6558351 0.7108351 0.9999997
## 2:(+DH+Si1)-1:(+DH+Si1) -0.0340 -0.7173351 0.6493351 0.9999986
## 1:(+DH+Si2)-1:(+DH+Si1) 0.1370 -0.5463351 0.8203351 0.9886758
## 2:(+DH+Si2)-1:(+DH+Si1) -0.0900 -0.7733351 0.5933351 0.9990826
## 1:(+DH+Si2)-2:(+DH+Si1) 0.1710 -0.5123351 0.8543351 0.9635450
## 2:(+DH+Si2)-2:(+DH+Si1) -0.0560 -0.7393351 0.6273351 0.9999589
## 2:(+DH+Si2)-1:(+DH+Si2) -0.2270 -0.9103351 0.4563351 0.8707440ggplot(datos_m3)+
aes(TRATAMIENTO,SPAD, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=SPAD, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_SPAD=aov(SPAD~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_SPAD)## Analysis of Variance Table
##
## Response: SPAD
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 11.902 11.9025 0.6756 0.4349
## TRATAMIENTO 3 4.788 1.5958 0.0906 0.9632
## BLOQUES:TRATAMIENTO 3 58.228 19.4092 1.1017 0.4032
## Residuals 8 140.940 17.6175TukeyHSD(mod3_SPAD)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = SPAD ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -1.725 -6.564519 3.114519 0.4349148
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.10 -9.604433 9.404433 0.9999850
## (+DH+Si1)-(-DH-Si ) 0.30 -9.204433 9.804433 0.9995976
## (+DH+Si2)-(-DH-Si ) -1.15 -10.654433 8.354433 0.9788856
## (+DH+Si1)-(+DH-Si ) 0.40 -9.104433 9.904433 0.9990501
## (+DH+Si2)-(+DH-Si ) -1.05 -10.554433 8.454433 0.9837266
## (+DH+Si2)-(+DH+Si1) -1.45 -10.954433 8.054433 0.9594871
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -2.40 -19.0092 14.209205 0.9983619
## 1:(+DH-Si )-1:(-DH-Si ) -3.10 -19.7092 13.509205 0.9924409
## 2:(+DH-Si )-1:(-DH-Si ) 0.50 -16.1092 17.109205 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 2.65 -13.9592 19.259205 0.9969976
## 2:(+DH+Si1)-1:(-DH-Si ) -4.45 -21.0592 12.159205 0.9492539
## 1:(+DH+Si2)-1:(-DH-Si ) -1.85 -18.4592 14.759205 0.9996868
## 2:(+DH+Si2)-1:(-DH-Si ) -2.85 -19.4592 13.759205 0.9953649
## 1:(+DH-Si )-2:(-DH-Si ) -0.70 -17.3092 15.909205 0.9999996
## 2:(+DH-Si )-2:(-DH-Si ) 2.90 -13.7092 19.509205 0.9948655
## 1:(+DH+Si1)-2:(-DH-Si ) 5.05 -11.5592 21.659205 0.9102026
## 2:(+DH+Si1)-2:(-DH-Si ) -2.05 -18.6592 14.559205 0.9993929
## 1:(+DH+Si2)-2:(-DH-Si ) 0.55 -16.0592 17.159205 0.9999999
## 2:(+DH+Si2)-2:(-DH-Si ) -0.45 -17.0592 16.159205 1.0000000
## 2:(+DH-Si )-1:(+DH-Si ) 3.60 -13.0092 20.209205 0.9826522
## 1:(+DH+Si1)-1:(+DH-Si ) 5.75 -10.8592 22.359205 0.8484247
## 2:(+DH+Si1)-1:(+DH-Si ) -1.35 -17.9592 15.259205 0.9999611
## 1:(+DH+Si2)-1:(+DH-Si ) 1.25 -15.3592 17.859205 0.9999769
## 2:(+DH+Si2)-1:(+DH-Si ) 0.25 -16.3592 16.859205 1.0000000
## 1:(+DH+Si1)-2:(+DH-Si ) 2.15 -14.4592 18.759205 0.9991778
## 2:(+DH+Si1)-2:(+DH-Si ) -4.95 -21.5592 11.659205 0.9176366
## 1:(+DH+Si2)-2:(+DH-Si ) -2.35 -18.9592 14.259205 0.9985625
## 2:(+DH+Si2)-2:(+DH-Si ) -3.35 -19.9592 13.259205 0.9882912
## 2:(+DH+Si1)-1:(+DH+Si1) -7.10 -23.7092 9.509205 0.6937392
## 1:(+DH+Si2)-1:(+DH+Si1) -4.50 -21.1092 12.109205 0.9465094
## 2:(+DH+Si2)-1:(+DH+Si1) -5.50 -22.1092 11.109205 0.8723533
## 1:(+DH+Si2)-2:(+DH+Si1) 2.60 -14.0092 19.209205 0.9973244
## 2:(+DH+Si2)-2:(+DH+Si1) 1.60 -15.0092 18.209205 0.9998792
## 2:(+DH+Si2)-1:(+DH+Si2) -1.00 -17.6092 15.609205 0.9999949#INTERACCIÓN
ggplot(datos_m3)+
aes(TRATAMIENTO,PESO_S_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=PESO_S_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PSB=aov(PESO_S_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PSB)## Analysis of Variance Table
##
## Response: PESO_S_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 6.3000e-08 6.3000e-08 0.1111 0.74745
## TRATAMIENTO 3 9.8188e-05 3.2729e-05 58.1852 8.915e-06 ***
## BLOQUES:TRATAMIENTO 3 2.2687e-05 7.5620e-06 13.4444 0.00172 **
## Residuals 8 4.5000e-06 5.6200e-07
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PSB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_S_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.000125 -0.0009897516 0.0007397516 0.747451
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.00100 -0.0006983036 0.002698304 0.3054963
## (+DH+Si1)-(-DH-Si ) -0.00475 -0.0064483036 -0.003051696 0.0000890
## (+DH+Si2)-(-DH-Si ) 0.00150 -0.0001983036 0.003198304 0.0847586
## (+DH+Si1)-(+DH-Si ) -0.00575 -0.0074483036 -0.004051696 0.0000216
## (+DH+Si2)-(+DH-Si ) 0.00050 -0.0011983036 0.002198304 0.7837632
## (+DH+Si2)-(+DH+Si1) 0.00625 0.0045516964 0.007948304 0.0000116
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 5.000000e-04 -0.0024678228 3.467823e-03 0.9958417
## 1:(+DH-Si )-1:(-DH-Si ) -6.938894e-18 -0.0029678228 2.967823e-03 1.0000000
## 2:(+DH-Si )-1:(-DH-Si ) 2.500000e-03 -0.0004678228 5.467823e-03 0.1109362
## 1:(+DH+Si1)-1:(-DH-Si ) -2.500000e-03 -0.0054678228 4.678228e-04 0.1109362
## 2:(+DH+Si1)-1:(-DH-Si ) -6.500000e-03 -0.0094678228 -3.532177e-03 0.0003665
## 1:(+DH+Si2)-1:(-DH-Si ) 1.500000e-03 -0.0014678228 4.467823e-03 0.5323050
## 2:(+DH+Si2)-1:(-DH-Si ) 2.000000e-03 -0.0009678228 4.967823e-03 0.2557796
## 1:(+DH-Si )-2:(-DH-Si ) -5.000000e-04 -0.0034678228 2.467823e-03 0.9958417
## 2:(+DH-Si )-2:(-DH-Si ) 2.000000e-03 -0.0009678228 4.967823e-03 0.2557796
## 1:(+DH+Si1)-2:(-DH-Si ) -3.000000e-03 -0.0059678228 -3.217722e-05 0.0473572
## 2:(+DH+Si1)-2:(-DH-Si ) -7.000000e-03 -0.0099678228 -4.032177e-03 0.0002145
## 1:(+DH+Si2)-2:(-DH-Si ) 1.000000e-03 -0.0019678228 3.967823e-03 0.8633546
## 2:(+DH+Si2)-2:(-DH-Si ) 1.500000e-03 -0.0014678228 4.467823e-03 0.5323050
## 2:(+DH-Si )-1:(+DH-Si ) 2.500000e-03 -0.0004678228 5.467823e-03 0.1109362
## 1:(+DH+Si1)-1:(+DH-Si ) -2.500000e-03 -0.0054678228 4.678228e-04 0.1109362
## 2:(+DH+Si1)-1:(+DH-Si ) -6.500000e-03 -0.0094678228 -3.532177e-03 0.0003665
## 1:(+DH+Si2)-1:(+DH-Si ) 1.500000e-03 -0.0014678228 4.467823e-03 0.5323050
## 2:(+DH+Si2)-1:(+DH-Si ) 2.000000e-03 -0.0009678228 4.967823e-03 0.2557796
## 1:(+DH+Si1)-2:(+DH-Si ) -5.000000e-03 -0.0079678228 -2.032177e-03 0.0022672
## 2:(+DH+Si1)-2:(+DH-Si ) -9.000000e-03 -0.0119678228 -6.032177e-03 0.0000330
## 1:(+DH+Si2)-2:(+DH-Si ) -1.000000e-03 -0.0039678228 1.967823e-03 0.8633546
## 2:(+DH+Si2)-2:(+DH-Si ) -5.000000e-04 -0.0034678228 2.467823e-03 0.9958417
## 2:(+DH+Si1)-1:(+DH+Si1) -4.000000e-03 -0.0069678228 -1.032177e-03 0.0094598
## 1:(+DH+Si2)-1:(+DH+Si1) 4.000000e-03 0.0010321772 6.967823e-03 0.0094598
## 2:(+DH+Si2)-1:(+DH+Si1) 4.500000e-03 0.0015321772 7.467823e-03 0.0045244
## 1:(+DH+Si2)-2:(+DH+Si1) 8.000000e-03 0.0050321772 1.096782e-02 0.0000801
## 2:(+DH+Si2)-2:(+DH+Si1) 8.500000e-03 0.0055321772 1.146782e-02 0.0000509
## 2:(+DH+Si2)-1:(+DH+Si2) 5.000000e-04 -0.0024678228 3.467823e-03 0.9958417ggplot(datos_m3)+
aes(TRATAMIENTO,PESO_F_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=PESO_F_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PFB=aov(PESO_F_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PFB)## Analysis of Variance Table
##
## Response: PESO_F_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.00012488 1.2488e-04 2.8523 0.12972
## TRATAMIENTO 3 0.00078765 2.6255e-04 5.9966 0.01917 *
## BLOQUES:TRATAMIENTO 3 0.00040418 1.3473e-04 3.0772 0.09049 .
## Residuals 8 0.00035026 4.3783e-05
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PFB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_F_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.0055875 -0.01321678 0.002041778 0.12972
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.005700 -0.009283299 0.0206833 0.6334040
## (+DH+Si1)-(-DH-Si ) 0.005100 -0.009883299 0.0200833 0.7049655
## (+DH+Si2)-(-DH-Si ) 0.018975 0.003991701 0.0339583 0.0154620
## (+DH+Si1)-(+DH-Si ) -0.000600 -0.015583299 0.0143833 0.9991810
## (+DH+Si2)-(+DH-Si ) 0.013275 -0.001708299 0.0282583 0.0837103
## (+DH+Si2)-(+DH+Si1) 0.013875 -0.001108299 0.0288583 0.0698387
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.00280 -0.0233836434 0.02898364 0.9997599
## 1:(+DH-Si )-1:(-DH-Si ) 0.01480 -0.0113836434 0.04098364 0.4185654
## 2:(+DH-Si )-1:(-DH-Si ) -0.00060 -0.0267836434 0.02558364 1.0000000
## 1:(+DH+Si1)-1:(-DH-Si ) 0.00350 -0.0226836434 0.02968364 0.9989928
## 2:(+DH+Si1)-1:(-DH-Si ) 0.00950 -0.0166836434 0.03568364 0.8198651
## 1:(+DH+Si2)-1:(-DH-Si ) 0.02825 0.0020663566 0.05443364 0.0337688
## 2:(+DH+Si2)-1:(-DH-Si ) 0.01250 -0.0136836434 0.03868364 0.5896424
## 1:(+DH-Si )-2:(-DH-Si ) 0.01200 -0.0141836434 0.03818364 0.6294526
## 2:(+DH-Si )-2:(-DH-Si ) -0.00340 -0.0295836434 0.02278364 0.9991613
## 1:(+DH+Si1)-2:(-DH-Si ) 0.00070 -0.0254836434 0.02688364 1.0000000
## 2:(+DH+Si1)-2:(-DH-Si ) 0.00670 -0.0194836434 0.03288364 0.9593195
## 1:(+DH+Si2)-2:(-DH-Si ) 0.02545 -0.0007336434 0.05163364 0.0575600
## 2:(+DH+Si2)-2:(-DH-Si ) 0.00970 -0.0164836434 0.03588364 0.8060796
## 2:(+DH-Si )-1:(+DH-Si ) -0.01540 -0.0415836434 0.01078364 0.3792038
## 1:(+DH+Si1)-1:(+DH-Si ) -0.01130 -0.0374836434 0.01488364 0.6852603
## 2:(+DH+Si1)-1:(+DH-Si ) -0.00530 -0.0314836434 0.02088364 0.9880571
## 1:(+DH+Si2)-1:(+DH-Si ) 0.01345 -0.0127336434 0.03963364 0.5157991
## 2:(+DH+Si2)-1:(+DH-Si ) -0.00230 -0.0284836434 0.02388364 0.9999346
## 1:(+DH+Si1)-2:(+DH-Si ) 0.00410 -0.0220836434 0.03028364 0.9973196
## 2:(+DH+Si1)-2:(+DH-Si ) 0.01010 -0.0160836434 0.03628364 0.7774198
## 1:(+DH+Si2)-2:(+DH-Si ) 0.02885 0.0026663566 0.05503364 0.0301726
## 2:(+DH+Si2)-2:(+DH-Si ) 0.01310 -0.0130836434 0.03928364 0.5426170
## 2:(+DH+Si1)-1:(+DH+Si1) 0.00600 -0.0201836434 0.03218364 0.9767088
## 1:(+DH+Si2)-1:(+DH+Si1) 0.02475 -0.0014336434 0.05093364 0.0658716
## 2:(+DH+Si2)-1:(+DH+Si1) 0.00900 -0.0171836434 0.03518364 0.8524789
## 1:(+DH+Si2)-2:(+DH+Si1) 0.01875 -0.0074336434 0.04493364 0.2085676
## 2:(+DH+Si2)-2:(+DH+Si1) 0.00300 -0.0231836434 0.02918364 0.9996238
## 2:(+DH+Si2)-1:(+DH+Si2) -0.01575 -0.0419336434 0.01043364 0.3574600ggplot(datos_m3)+
aes(TRATAMIENTO,PESO_SAT_H_BOCADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=PESO_SAT_H_BOCADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_PSHB=aov(PESO_SAT_H_BOCADOS~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_PSHB)## Analysis of Variance Table
##
## Response: PESO_SAT_H_BOCADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 0.0001000 0.00010000 0.7968 0.39810
## TRATAMIENTO 3 0.0018322 0.00061075 4.8665 0.03269 *
## BLOQUES:TRATAMIENTO 3 0.0005335 0.00017783 1.4170 0.30746
## Residuals 8 0.0010040 0.00012550
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_PSHB)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = PESO_SAT_H_BOCADOS ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -0.005 -0.01791671 0.007916711 0.3981029
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.01125 -0.014117399 0.0366174 0.5220375
## (+DH+Si1)-(-DH-Si ) 0.00800 -0.017367399 0.0333674 0.7485534
## (+DH+Si2)-(-DH-Si ) 0.02925 0.003882601 0.0546174 0.0252777
## (+DH+Si1)-(+DH-Si ) -0.00325 -0.028617399 0.0221174 0.9751599
## (+DH+Si2)-(+DH-Si ) 0.01800 -0.007367399 0.0433674 0.1839015
## (+DH+Si2)-(+DH+Si1) 0.02125 -0.004117399 0.0466174 0.1041035
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) 0.0050 -0.039330085 0.04933008 0.9996602
## 1:(+DH-Si )-1:(-DH-Si ) 0.0220 -0.022330085 0.06633008 0.5508126
## 2:(+DH-Si )-1:(-DH-Si ) 0.0055 -0.038830085 0.04983008 0.9993723
## 1:(+DH+Si1)-1:(-DH-Si ) 0.0065 -0.037830085 0.05083008 0.9982069
## 2:(+DH+Si1)-1:(-DH-Si ) 0.0145 -0.029830085 0.05883008 0.8784522
## 1:(+DH+Si2)-1:(-DH-Si ) 0.0400 -0.004330085 0.08433008 0.0818607
## 2:(+DH+Si2)-1:(-DH-Si ) 0.0235 -0.020830085 0.06783008 0.4836378
## 1:(+DH-Si )-2:(-DH-Si ) 0.0170 -0.027330085 0.06133008 0.7817248
## 2:(+DH-Si )-2:(-DH-Si ) 0.0005 -0.043830085 0.04483008 1.0000000
## 1:(+DH+Si1)-2:(-DH-Si ) 0.0015 -0.042830085 0.04583008 0.9999999
## 2:(+DH+Si1)-2:(-DH-Si ) 0.0095 -0.034830085 0.05383008 0.9836805
## 1:(+DH+Si2)-2:(-DH-Si ) 0.0350 -0.009330085 0.07933008 0.1448328
## 2:(+DH+Si2)-2:(-DH-Si ) 0.0185 -0.025830085 0.06283008 0.7146145
## 2:(+DH-Si )-1:(+DH-Si ) -0.0165 -0.060830085 0.02783008 0.8028720
## 1:(+DH+Si1)-1:(+DH-Si ) -0.0155 -0.059830085 0.02883008 0.8426635
## 2:(+DH+Si1)-1:(+DH-Si ) -0.0075 -0.051830085 0.03683008 0.9957364
## 1:(+DH+Si2)-1:(+DH-Si ) 0.0180 -0.026330085 0.06233008 0.7374813
## 2:(+DH+Si2)-1:(+DH-Si ) 0.0015 -0.042830085 0.04583008 0.9999999
## 1:(+DH+Si1)-2:(+DH-Si ) 0.0010 -0.043330085 0.04533008 1.0000000
## 2:(+DH+Si1)-2:(+DH-Si ) 0.0090 -0.035330085 0.05333008 0.9878577
## 1:(+DH+Si2)-2:(+DH-Si ) 0.0345 -0.009830085 0.07883008 0.1532629
## 2:(+DH+Si2)-2:(+DH-Si ) 0.0180 -0.026330085 0.06233008 0.7374813
## 2:(+DH+Si1)-1:(+DH+Si1) 0.0080 -0.036330085 0.05233008 0.9937745
## 1:(+DH+Si2)-1:(+DH+Si1) 0.0335 -0.010830085 0.07783008 0.1715335
## 2:(+DH+Si2)-1:(+DH+Si1) 0.0170 -0.027330085 0.06133008 0.7817248
## 1:(+DH+Si2)-2:(+DH+Si1) 0.0255 -0.018830085 0.06983008 0.4010848
## 2:(+DH+Si2)-2:(+DH+Si1) 0.0090 -0.035330085 0.05333008 0.9878577
## 2:(+DH+Si2)-1:(+DH+Si2) -0.0165 -0.060830085 0.02783008 0.8028720ggplot(datos_m3)+
aes(TRATAMIENTO,N_ESTOMAS_ADAX, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ADAX, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_NEAx=aov(N_ESTOMAS_ADAX~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_NEAx)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ADAX
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 27.56 27.562 0.2176 0.6533
## TRATAMIENTO 3 74.19 24.729 0.1952 0.8968
## BLOQUES:TRATAMIENTO 3 150.69 50.229 0.3965 0.7592
## Residuals 8 1013.50 126.688TukeyHSD(mod3_NEAx)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ADAX ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 2.625 -10.35268 15.60268 0.6533411
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 5.75 -19.73713 31.23713 0.8854765
## (+DH+Si1)-(-DH-Si ) 4.00 -21.48713 29.48713 0.9561939
## (+DH+Si2)-(-DH-Si ) 2.00 -23.48713 27.48713 0.9939825
## (+DH+Si1)-(+DH-Si ) -1.75 -27.23713 23.73713 0.9959383
## (+DH+Si2)-(+DH-Si ) -3.75 -29.23713 21.73713 0.9633605
## (+DH+Si2)-(+DH+Si1) -2.00 -27.48713 23.48713 0.9939825
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -6.0 -50.53932 38.53932 0.9989424
## 1:(+DH-Si )-1:(-DH-Si ) 0.5 -44.03932 45.03932 1.0000000
## 2:(+DH-Si )-1:(-DH-Si ) 5.0 -39.53932 49.53932 0.9996705
## 1:(+DH+Si1)-1:(-DH-Si ) -4.5 -49.03932 40.03932 0.9998346
## 2:(+DH+Si1)-1:(-DH-Si ) 6.5 -38.03932 51.03932 0.9982583
## 1:(+DH+Si2)-1:(-DH-Si ) -1.5 -46.03932 43.03932 0.9999999
## 2:(+DH+Si2)-1:(-DH-Si ) -0.5 -45.03932 44.03932 1.0000000
## 1:(+DH-Si )-2:(-DH-Si ) 6.5 -38.03932 51.03932 0.9982583
## 2:(+DH-Si )-2:(-DH-Si ) 11.0 -33.53932 55.53932 0.9658544
## 1:(+DH+Si1)-2:(-DH-Si ) 1.5 -43.03932 46.03932 0.9999999
## 2:(+DH+Si1)-2:(-DH-Si ) 12.5 -32.03932 57.03932 0.9370201
## 1:(+DH+Si2)-2:(-DH-Si ) 4.5 -40.03932 49.03932 0.9998346
## 2:(+DH+Si2)-2:(-DH-Si ) 5.5 -39.03932 50.03932 0.9993909
## 2:(+DH-Si )-1:(+DH-Si ) 4.5 -40.03932 49.03932 0.9998346
## 1:(+DH+Si1)-1:(+DH-Si ) -5.0 -49.53932 39.53932 0.9996705
## 2:(+DH+Si1)-1:(+DH-Si ) 6.0 -38.53932 50.53932 0.9989424
## 1:(+DH+Si2)-1:(+DH-Si ) -2.0 -46.53932 42.53932 0.9999993
## 2:(+DH+Si2)-1:(+DH-Si ) -1.0 -45.53932 43.53932 1.0000000
## 1:(+DH+Si1)-2:(+DH-Si ) -9.5 -54.03932 35.03932 0.9840907
## 2:(+DH+Si1)-2:(+DH-Si ) 1.5 -43.03932 46.03932 0.9999999
## 1:(+DH+Si2)-2:(+DH-Si ) -6.5 -51.03932 38.03932 0.9982583
## 2:(+DH+Si2)-2:(+DH-Si ) -5.5 -50.03932 39.03932 0.9993909
## 2:(+DH+Si1)-1:(+DH+Si1) 11.0 -33.53932 55.53932 0.9658544
## 1:(+DH+Si2)-1:(+DH+Si1) 3.0 -41.53932 47.53932 0.9999891
## 2:(+DH+Si2)-1:(+DH+Si1) 4.0 -40.53932 48.53932 0.9999242
## 1:(+DH+Si2)-2:(+DH+Si1) -8.0 -52.53932 36.53932 0.9939424
## 2:(+DH+Si2)-2:(+DH+Si1) -7.0 -51.53932 37.53932 0.9972590
## 2:(+DH+Si2)-1:(+DH+Si2) 1.0 -43.53932 45.53932 1.0000000ggplot(datos_m3)+
aes(TRATAMIENTO,N_ESTOMAS_ABA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ABA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_NEAbx=aov(N_ESTOMAS_ABA~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_NEAbx)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ABA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 162.56 162.562 1.4670 0.2604
## TRATAMIENTO 3 130.69 43.562 0.3931 0.7615
## BLOQUES:TRATAMIENTO 3 448.19 149.396 1.3482 0.3259
## Residuals 8 886.50 110.812TukeyHSD(mod3_NEAbx)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ABA ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 6.375 -5.762366 18.51237 0.2603805
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.50 -23.33683 24.33683 0.9998816
## (+DH+Si1)-(-DH-Si ) 6.75 -17.08683 30.58683 0.8019887
## (+DH+Si2)-(-DH-Si ) 0.00 -23.83683 23.83683 1.0000000
## (+DH+Si1)-(+DH-Si ) 6.25 -17.58683 30.08683 0.8345063
## (+DH+Si2)-(+DH-Si ) -0.50 -24.33683 23.33683 0.9998816
## (+DH+Si2)-(+DH+Si1) -6.75 -30.58683 17.08683 0.8019887
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -9.500000e+00 -51.15538 32.15538 0.9772812
## 1:(+DH-Si )-1:(-DH-Si ) -7.000000e+00 -48.65538 34.65538 0.9959040
## 2:(+DH-Si )-1:(-DH-Si ) -1.500000e+00 -43.15538 40.15538 0.9999999
## 1:(+DH+Si1)-1:(-DH-Si ) -8.000000e+00 -49.65538 33.65538 0.9911048
## 2:(+DH+Si1)-1:(-DH-Si ) 1.200000e+01 -29.65538 53.65538 0.9291139
## 1:(+DH+Si2)-1:(-DH-Si ) -9.500000e+00 -51.15538 32.15538 0.9772812
## 2:(+DH+Si2)-1:(-DH-Si ) 0.000000e+00 -41.65538 41.65538 1.0000000
## 1:(+DH-Si )-2:(-DH-Si ) 2.500000e+00 -39.15538 44.15538 0.9999950
## 2:(+DH-Si )-2:(-DH-Si ) 8.000000e+00 -33.65538 49.65538 0.9911048
## 1:(+DH+Si1)-2:(-DH-Si ) 1.500000e+00 -40.15538 43.15538 0.9999999
## 2:(+DH+Si1)-2:(-DH-Si ) 2.150000e+01 -20.15538 63.15538 0.5109238
## 1:(+DH+Si2)-2:(-DH-Si ) 3.552714e-15 -41.65538 41.65538 1.0000000
## 2:(+DH+Si2)-2:(-DH-Si ) 9.500000e+00 -32.15538 51.15538 0.9772812
## 2:(+DH-Si )-1:(+DH-Si ) 5.500000e+00 -36.15538 47.15538 0.9990681
## 1:(+DH+Si1)-1:(+DH-Si ) -1.000000e+00 -42.65538 40.65538 1.0000000
## 2:(+DH+Si1)-1:(+DH-Si ) 1.900000e+01 -22.65538 60.65538 0.6340083
## 1:(+DH+Si2)-1:(+DH-Si ) -2.500000e+00 -44.15538 39.15538 0.9999950
## 2:(+DH+Si2)-1:(+DH-Si ) 7.000000e+00 -34.65538 48.65538 0.9959040
## 1:(+DH+Si1)-2:(+DH-Si ) -6.500000e+00 -48.15538 35.15538 0.9973756
## 2:(+DH+Si1)-2:(+DH-Si ) 1.350000e+01 -28.15538 55.15538 0.8828785
## 1:(+DH+Si2)-2:(+DH-Si ) -8.000000e+00 -49.65538 33.65538 0.9911048
## 2:(+DH+Si2)-2:(+DH-Si ) 1.500000e+00 -40.15538 43.15538 0.9999999
## 2:(+DH+Si1)-1:(+DH+Si1) 2.000000e+01 -21.65538 61.65538 0.5839800
## 1:(+DH+Si2)-1:(+DH+Si1) -1.500000e+00 -43.15538 40.15538 0.9999999
## 2:(+DH+Si2)-1:(+DH+Si1) 8.000000e+00 -33.65538 49.65538 0.9911048
## 1:(+DH+Si2)-2:(+DH+Si1) -2.150000e+01 -63.15538 20.15538 0.5109238
## 2:(+DH+Si2)-2:(+DH+Si1) -1.200000e+01 -53.65538 29.65538 0.9291139
## 2:(+DH+Si2)-1:(+DH+Si2) 9.500000e+00 -32.15538 51.15538 0.9772812ggplot(datos_m3)+
aes(TRATAMIENTO,N_ESTOMAS_ABIERTOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_ABIERTOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_NEA=aov(N_ESTOMAS_ABIERTOS~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_NEA)## Analysis of Variance Table
##
## Response: N_ESTOMAS_ABIERTOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 10.562 10.562 0.6190 0.45408
## TRATAMIENTO 3 242.188 80.729 4.7314 0.03502 *
## BLOQUES:TRATAMIENTO 3 53.187 17.729 1.0391 0.42607
## Residuals 8 136.500 17.063
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_NEA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_ABIERTOS ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 1.625 -3.13768 6.38768 0.4540759
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -10.50 -19.853527 -1.146473 0.0289138
## (+DH+Si1)-(-DH-Si ) -7.25 -16.603527 2.103527 0.1377907
## (+DH+Si2)-(-DH-Si ) -4.00 -13.353527 5.353527 0.5494603
## (+DH+Si1)-(+DH-Si ) 3.25 -6.103527 12.603527 0.6924097
## (+DH+Si2)-(+DH-Si ) 6.50 -2.853527 15.853527 0.1959619
## (+DH+Si2)-(+DH+Si1) 3.25 -6.103527 12.603527 0.6924097
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -4.5 -20.845493 11.845493 0.9423601
## 1:(+DH-Si )-1:(-DH-Si ) -14.0 -30.345493 2.345493 0.1032691
## 2:(+DH-Si )-1:(-DH-Si ) -11.5 -27.845493 4.845493 0.2217203
## 1:(+DH+Si1)-1:(-DH-Si ) -12.0 -28.345493 4.345493 0.1908623
## 2:(+DH+Si1)-1:(-DH-Si ) -7.0 -23.345493 9.345493 0.6921283
## 1:(+DH+Si2)-1:(-DH-Si ) -8.0 -24.345493 8.345493 0.5648066
## 2:(+DH+Si2)-1:(-DH-Si ) -4.5 -20.845493 11.845493 0.9423601
## 1:(+DH-Si )-2:(-DH-Si ) -9.5 -25.845493 6.845493 0.3908735
## 2:(+DH-Si )-2:(-DH-Si ) -7.0 -23.345493 9.345493 0.6921283
## 1:(+DH+Si1)-2:(-DH-Si ) -7.5 -23.845493 8.845493 0.6283217
## 2:(+DH+Si1)-2:(-DH-Si ) -2.5 -18.845493 13.845493 0.9976769
## 1:(+DH+Si2)-2:(-DH-Si ) -3.5 -19.845493 12.845493 0.9837524
## 2:(+DH+Si2)-2:(-DH-Si ) 0.0 -16.345493 16.345493 1.0000000
## 2:(+DH-Si )-1:(+DH-Si ) 2.5 -13.845493 18.845493 0.9976769
## 1:(+DH+Si1)-1:(+DH-Si ) 2.0 -14.345493 18.345493 0.9994257
## 2:(+DH+Si1)-1:(+DH-Si ) 7.0 -9.345493 23.345493 0.6921283
## 1:(+DH+Si2)-1:(+DH-Si ) 6.0 -10.345493 22.345493 0.8122398
## 2:(+DH+Si2)-1:(+DH-Si ) 9.5 -6.845493 25.845493 0.3908735
## 1:(+DH+Si1)-2:(+DH-Si ) -0.5 -16.845493 15.845493 1.0000000
## 2:(+DH+Si1)-2:(+DH-Si ) 4.5 -11.845493 20.845493 0.9423601
## 1:(+DH+Si2)-2:(+DH-Si ) 3.5 -12.845493 19.845493 0.9837524
## 2:(+DH+Si2)-2:(+DH-Si ) 7.0 -9.345493 23.345493 0.6921283
## 2:(+DH+Si1)-1:(+DH+Si1) 5.0 -11.345493 21.345493 0.9078497
## 1:(+DH+Si2)-1:(+DH+Si1) 4.0 -12.345493 20.345493 0.9673892
## 2:(+DH+Si2)-1:(+DH+Si1) 7.5 -8.845493 23.845493 0.6283217
## 1:(+DH+Si2)-2:(+DH+Si1) -1.0 -17.345493 15.345493 0.9999944
## 2:(+DH+Si2)-2:(+DH+Si1) 2.5 -13.845493 18.845493 0.9976769
## 2:(+DH+Si2)-1:(+DH+Si2) 3.5 -12.845493 19.845493 0.9837524ggplot(datos_m3)+
aes(TRATAMIENTO,N_ESTOMAS_P_ABIERTOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_P_ABIERTOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_NEPA=aov(N_ESTOMAS_P_ABIERTOS~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_NEPA)## Analysis of Variance Table
##
## Response: N_ESTOMAS_P_ABIERTOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 30.25 30.250 0.3617 0.5642
## TRATAMIENTO 3 360.25 120.083 1.4360 0.3026
## BLOQUES:TRATAMIENTO 3 122.25 40.750 0.4873 0.7006
## Residuals 8 669.00 83.625TukeyHSD(mod3_NEPA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_P_ABIERTOS ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 2.75 -7.793824 13.29382 0.5641949
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 3.25 -17.457236 23.95724 0.9561878
## (+DH+Si1)-(-DH-Si ) 9.25 -11.457236 29.95724 0.5165250
## (+DH+Si2)-(-DH-Si ) 12.00 -8.707236 32.70724 0.3170999
## (+DH+Si1)-(+DH-Si ) 6.00 -14.707236 26.70724 0.7913798
## (+DH+Si2)-(+DH-Si ) 8.75 -11.957236 29.45724 0.5583732
## (+DH+Si2)-(+DH+Si1) 2.75 -17.957236 23.45724 0.9725045
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -1.5 -37.68635 34.68635 0.9999996
## 1:(+DH-Si )-1:(-DH-Si ) 1.5 -34.68635 37.68635 0.9999996
## 2:(+DH-Si )-1:(-DH-Si ) 3.5 -32.68635 39.68635 0.9998759
## 1:(+DH+Si1)-1:(-DH-Si ) 2.5 -33.68635 38.68635 0.9999871
## 2:(+DH+Si1)-1:(-DH-Si ) 14.5 -21.68635 50.68635 0.7481570
## 1:(+DH+Si2)-1:(-DH-Si ) 12.0 -24.18635 48.18635 0.8716282
## 2:(+DH+Si2)-1:(-DH-Si ) 10.5 -25.68635 46.68635 0.9267908
## 1:(+DH-Si )-2:(-DH-Si ) 3.0 -33.18635 39.18635 0.9999556
## 2:(+DH-Si )-2:(-DH-Si ) 5.0 -31.18635 41.18635 0.9987601
## 1:(+DH+Si1)-2:(-DH-Si ) 4.0 -32.18635 40.18635 0.9997020
## 2:(+DH+Si1)-2:(-DH-Si ) 16.0 -20.18635 52.18635 0.6632153
## 1:(+DH+Si2)-2:(-DH-Si ) 13.5 -22.68635 49.68635 0.8012842
## 2:(+DH+Si2)-2:(-DH-Si ) 12.0 -24.18635 48.18635 0.8716282
## 2:(+DH-Si )-1:(+DH-Si ) 2.0 -34.18635 38.18635 0.9999972
## 1:(+DH+Si1)-1:(+DH-Si ) 1.0 -35.18635 37.18635 1.0000000
## 2:(+DH+Si1)-1:(+DH-Si ) 13.0 -23.18635 49.18635 0.8261712
## 1:(+DH+Si2)-1:(+DH-Si ) 10.5 -25.68635 46.68635 0.9267908
## 2:(+DH+Si2)-1:(+DH-Si ) 9.0 -27.18635 45.18635 0.9646385
## 1:(+DH+Si1)-2:(+DH-Si ) -1.0 -37.18635 35.18635 1.0000000
## 2:(+DH+Si1)-2:(+DH-Si ) 11.0 -25.18635 47.18635 0.9102864
## 1:(+DH+Si2)-2:(+DH-Si ) 8.5 -27.68635 44.68635 0.9735205
## 2:(+DH+Si2)-2:(+DH-Si ) 7.0 -29.18635 43.18635 0.9907340
## 2:(+DH+Si1)-1:(+DH+Si1) 12.0 -24.18635 48.18635 0.8716282
## 1:(+DH+Si2)-1:(+DH+Si1) 9.5 -26.68635 45.68635 0.9539425
## 2:(+DH+Si2)-1:(+DH+Si1) 8.0 -28.18635 44.18635 0.9807181
## 1:(+DH+Si2)-2:(+DH+Si1) -2.5 -38.68635 33.68635 0.9999871
## 2:(+DH+Si2)-2:(+DH+Si1) -4.0 -40.18635 32.18635 0.9997020
## 2:(+DH+Si2)-1:(+DH+Si2) -1.5 -37.68635 34.68635 0.9999996ggplot(datos_m3)+
aes(TRATAMIENTO,N_ESTOMAS_CERRADOS, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=N_ESTOMAS_CERRADOS, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_NEC=aov(N_ESTOMAS_CERRADOS~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_NEC)## Analysis of Variance Table
##
## Response: N_ESTOMAS_CERRADOS
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 95.06 95.063 0.9070 0.36880
## TRATAMIENTO 3 932.19 310.729 2.9646 0.09735 .
## BLOQUES:TRATAMIENTO 3 275.69 91.896 0.8768 0.49253
## Residuals 8 838.50 104.813
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1TukeyHSD(mod3_NEC)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = N_ESTOMAS_CERRADOS ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 4.875 -6.929202 16.6792 0.3688014
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 14.00 -9.182518 37.182518 0.2874111
## (+DH+Si1)-(-DH-Si ) 9.25 -13.932518 32.432518 0.6001108
## (+DH+Si2)-(-DH-Si ) -5.50 -28.682518 17.682518 0.8701445
## (+DH+Si1)-(+DH-Si ) -4.75 -27.932518 18.432518 0.9104786
## (+DH+Si2)-(+DH-Si ) -19.50 -42.682518 3.682518 0.1024928
## (+DH+Si2)-(+DH+Si1) -14.75 -37.932518 8.432518 0.2514982
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -8.5 -49.01196 32.01196 0.9854509
## 1:(+DH-Si )-1:(-DH-Si ) 7.0 -33.51196 47.51196 0.9951711
## 2:(+DH-Si )-1:(-DH-Si ) 12.5 -28.01196 53.01196 0.9044053
## 1:(+DH+Si1)-1:(-DH-Si ) -2.0 -42.51196 38.51196 0.9999987
## 2:(+DH+Si1)-1:(-DH-Si ) 12.0 -28.51196 52.51196 0.9198100
## 1:(+DH+Si2)-1:(-DH-Si ) -14.0 -54.51196 26.51196 0.8494425
## 2:(+DH+Si2)-1:(-DH-Si ) -5.5 -46.01196 35.01196 0.9988896
## 1:(+DH-Si )-2:(-DH-Si ) 15.5 -25.01196 56.01196 0.7834079
## 2:(+DH-Si )-2:(-DH-Si ) 21.0 -19.51196 61.51196 0.5065313
## 1:(+DH+Si1)-2:(-DH-Si ) 6.5 -34.01196 47.01196 0.9968945
## 2:(+DH+Si1)-2:(-DH-Si ) 20.5 -20.01196 61.01196 0.5310939
## 1:(+DH+Si2)-2:(-DH-Si ) -5.5 -46.01196 35.01196 0.9988896
## 2:(+DH+Si2)-2:(-DH-Si ) 3.0 -37.51196 43.51196 0.9999793
## 2:(+DH-Si )-1:(+DH-Si ) 5.5 -35.01196 46.01196 0.9988896
## 1:(+DH+Si1)-1:(+DH-Si ) -9.0 -49.51196 31.51196 0.9802132
## 2:(+DH+Si1)-1:(+DH-Si ) 5.0 -35.51196 45.51196 0.9993930
## 1:(+DH+Si2)-1:(+DH-Si ) -21.0 -61.51196 19.51196 0.5065313
## 2:(+DH+Si2)-1:(+DH-Si ) -12.5 -53.01196 28.01196 0.9044053
## 1:(+DH+Si1)-2:(+DH-Si ) -14.5 -55.01196 26.01196 0.8285012
## 2:(+DH+Si1)-2:(+DH-Si ) -0.5 -41.01196 40.01196 1.0000000
## 1:(+DH+Si2)-2:(+DH-Si ) -26.5 -67.01196 14.01196 0.2808771
## 2:(+DH+Si2)-2:(+DH-Si ) -18.0 -58.51196 22.51196 0.6587093
## 2:(+DH+Si1)-1:(+DH+Si1) 14.0 -26.51196 54.51196 0.8494425
## 1:(+DH+Si2)-1:(+DH+Si1) -12.0 -52.51196 28.51196 0.9198100
## 2:(+DH+Si2)-1:(+DH+Si1) -3.5 -44.01196 37.01196 0.9999415
## 1:(+DH+Si2)-2:(+DH+Si1) -26.0 -66.51196 14.51196 0.2975551
## 2:(+DH+Si2)-2:(+DH+Si1) -17.5 -58.01196 23.01196 0.6844217
## 2:(+DH+Si2)-1:(+DH+Si2) 8.5 -32.01196 49.01196 0.9854509ggplot(datos_m3)+
aes(TRATAMIENTO,CRA, fill=TRATAMIENTO)+
geom_boxplot()+
facet_wrap(~BLOQUES)
ggplot(datos_m3)+
aes(x=TRATAMIENTO, y=CRA, fill=TRATAMIENTO)+
geom_bar(stat = "identity")+
facet_wrap(~BLOQUES)
mod3_CRA=aov(CRA~BLOQUES*TRATAMIENTO, data=datos_m3)
anova(mod3_CRA)## Analysis of Variance Table
##
## Response: CRA
## Df Sum Sq Mean Sq F value Pr(>F)
## BLOQUES 1 7.07 7.069 0.1761 0.6858
## TRATAMIENTO 3 39.32 13.106 0.3265 0.8065
## BLOQUES:TRATAMIENTO 3 4.43 1.478 0.0368 0.9898
## Residuals 8 321.16 40.145TukeyHSD(mod3_CRA)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = CRA ~ BLOQUES * TRATAMIENTO, data = datos_m3)
##
## $BLOQUES
## diff lwr upr p adj
## 2-1 -1.32941 -8.634874 5.976054 0.6857948
##
## $TRATAMIENTO
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -2.8755892 -17.22294 11.47176 0.9154741
## (+DH+Si1)-(-DH-Si ) -0.4980160 -14.84537 13.84934 0.9994655
## (+DH+Si2)-(-DH-Si ) -3.7338596 -18.08121 10.61349 0.8374210
## (+DH+Si1)-(+DH-Si ) 2.3775733 -11.96978 16.72493 0.9491521
## (+DH+Si2)-(+DH-Si ) -0.8582703 -15.20562 13.48908 0.9972979
## (+DH+Si2)-(+DH+Si1) -3.2358436 -17.58320 11.11151 0.8855636
##
## $`BLOQUES:TRATAMIENTO`
## diff lwr upr p adj
## 2:(-DH-Si )-1:(-DH-Si ) -1.39197673 -26.46429 23.68034 0.9999971
## 1:(+DH-Si )-1:(-DH-Si ) -2.80871855 -27.88103 22.26360 0.9996749
## 2:(+DH-Si )-1:(-DH-Si ) -4.33443663 -29.40675 20.73788 0.9951565
## 1:(+DH+Si1)-1:(-DH-Si ) -1.33209696 -26.40441 23.74022 0.9999979
## 2:(+DH+Si1)-1:(-DH-Si ) -1.05591168 -26.12823 24.01640 0.9999996
## 1:(+DH+Si2)-1:(-DH-Si ) -3.09178322 -28.16410 21.98053 0.9993963
## 2:(+DH+Si2)-1:(-DH-Si ) -5.76791264 -30.84023 19.30440 0.9762279
## 1:(+DH-Si )-2:(-DH-Si ) -1.41674181 -26.48906 23.65557 0.9999967
## 2:(+DH-Si )-2:(-DH-Si ) -2.94245990 -28.01477 22.12985 0.9995607
## 1:(+DH+Si1)-2:(-DH-Si ) 0.05987978 -25.01243 25.13219 1.0000000
## 2:(+DH+Si1)-2:(-DH-Si ) 0.33606505 -24.73625 25.40838 1.0000000
## 1:(+DH+Si2)-2:(-DH-Si ) -1.69980648 -26.77212 23.37251 0.9999886
## 2:(+DH+Si2)-2:(-DH-Si ) -4.37593591 -29.44825 20.69638 0.9948775
## 2:(+DH-Si )-1:(+DH-Si ) -1.52571808 -26.59803 23.54660 0.9999946
## 1:(+DH+Si1)-1:(+DH-Si ) 1.47662159 -23.59569 26.54894 0.9999957
## 2:(+DH+Si1)-1:(+DH-Si ) 1.75280687 -23.31951 26.82512 0.9999860
## 1:(+DH+Si2)-1:(+DH-Si ) -0.28306467 -25.35538 24.78925 1.0000000
## 2:(+DH+Si2)-1:(+DH-Si ) -2.95919409 -28.03151 22.11312 0.9995444
## 1:(+DH+Si1)-2:(+DH-Si ) 3.00233967 -22.06998 28.07465 0.9994999
## 2:(+DH+Si1)-2:(+DH-Si ) 3.27852495 -21.79379 28.35084 0.9991234
## 1:(+DH+Si2)-2:(+DH-Si ) 1.24265341 -23.82966 26.31497 0.9999987
## 2:(+DH+Si2)-2:(+DH-Si ) -1.43347601 -26.50579 23.63884 0.9999965
## 2:(+DH+Si1)-1:(+DH+Si1) 0.27618528 -24.79613 25.34850 1.0000000
## 1:(+DH+Si2)-1:(+DH+Si1) -1.75968626 -26.83200 23.31263 0.9999856
## 2:(+DH+Si2)-1:(+DH+Si1) -4.43581569 -29.50813 20.63650 0.9944532
## 1:(+DH+Si2)-2:(+DH+Si1) -2.03587154 -27.10819 23.03644 0.9999614
## 2:(+DH+Si2)-2:(+DH+Si1) -4.71200096 -29.78432 20.36031 0.9921364
## 2:(+DH+Si2)-1:(+DH+Si2) -2.67612943 -27.74844 22.39619 0.9997628#Comprobación de supuestos
#Comprobación de supuestos
shapiro.test(mod3_AF$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_AF$residuals
## W = 0.89984, p-value = 0.07991shapiro.test(mod3_DR$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_DR$residuals
## W = 0.94496, p-value = 0.4143shapiro.test(mod3_LR$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_LR$residuals
## W = 0.98644, p-value = 0.9949shapiro.test(mod3_NEA$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_NEA$residuals
## W = 0.96095, p-value = 0.679shapiro.test(mod3_NEAbx$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_NEAbx$residuals
## W = 0.99238, p-value = 0.9999shapiro.test(mod3_NEAx$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_NEAx$residuals
## W = 0.97108, p-value = 0.8558shapiro.test(mod3_NEC$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_NEC$residuals
## W = 0.9898, p-value = 0.9991shapiro.test(mod3_NEPA$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_NEPA$residuals
## W = 0.97527, p-value = 0.915shapiro.test(mod3_NH$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_NH$residuals
## W = 0.87436, p-value = 0.03173shapiro.test(mod3_PA$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PA$residuals
## W = 0.95101, p-value = 0.5058shapiro.test(mod3_PFA$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PFA$residuals
## W = 0.88968, p-value = 0.05505shapiro.test(mod3_PFB$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PFB$residuals
## W = 0.97929, p-value = 0.9577shapiro.test(mod3_PFR$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PFR$residuals
## W = 0.97212, p-value = 0.8717shapiro.test(mod3_PSA$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PSA$residuals
## W = 0.97987, p-value = 0.9625shapiro.test(mod3_PSB$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PSB$residuals
## W = 0.91708, p-value = 0.1514shapiro.test(mod3_PSHB$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PSHB$residuals
## W = 0.97693, p-value = 0.9346shapiro.test(mod3_PSR$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_PSR$residuals
## W = 0.92175, p-value = 0.1799shapiro.test(mod3_SPAD$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_SPAD$residuals
## W = 0.9741, p-value = 0.8999shapiro.test(mod3_T$residuals)##
## Shapiro-Wilk normality test
##
## data: mod3_T$residuals
## W = 0.96305, p-value = 0.7172#Para todas las pruebas de normalidad de acepta la H_o de que existe normalidadbartlett.test(datos_m3$TEMPERATURA~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$TEMPERATURA by TRATAMIENTO
## Bartlett's K-squared = 11.074, df = 3, p-value = 0.01133bartlett.test(datos_m3$N_HOJAS~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$N_HOJAS by TRATAMIENTO
## Bartlett's K-squared = 1.0648, df = 3, p-value = 0.7856bartlett.test(datos_m3$LONGITUD_P_AEREA~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$LONGITUD_P_AEREA by TRATAMIENTO
## Bartlett's K-squared = 3.2935, df = 3, p-value = 0.3485bartlett.test(datos_m3$AREA_FOLIAR~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$AREA_FOLIAR by TRATAMIENTO
## Bartlett's K-squared = 0.62666, df = 3, p-value = 0.8903bartlett.test(datos_m3$LONGITUD_RAIZ~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$LONGITUD_RAIZ by TRATAMIENTO
## Bartlett's K-squared = 7.9351, df = 3, p-value = 0.04737bartlett.test(datos_m3$DIAMETRO_RAIZ~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$DIAMETRO_RAIZ by TRATAMIENTO
## Bartlett's K-squared = 2.8279, df = 3, p-value = 0.4189bartlett.test(datos_m3$PESO_F_RAIZ~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$PESO_F_RAIZ by TRATAMIENTO
## Bartlett's K-squared = 5.1275, df = 3, p-value = 0.1627bartlett.test(datos_m3$PESO_F_AEREA~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$PESO_F_AEREA by TRATAMIENTO
## Bartlett's K-squared = 0.1399, df = 3, p-value = 0.9867bartlett.test(datos_m3$PESO_S_RAIZ~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$PESO_S_RAIZ by TRATAMIENTO
## Bartlett's K-squared = 1.0391, df = 3, p-value = 0.7918bartlett.test(datos_m3$PESO_S_AEREA~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$PESO_S_AEREA by TRATAMIENTO
## Bartlett's K-squared = 4.2125, df = 3, p-value = 0.2394bartlett.test(datos_m3$SPAD~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$SPAD by TRATAMIENTO
## Bartlett's K-squared = 3.4638, df = 3, p-value = 0.3255bartlett.test(datos_m3$PESO_F_BOCADOS~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$PESO_F_BOCADOS by TRATAMIENTO
## Bartlett's K-squared = 2.602, df = 3, p-value = 0.4571bartlett.test(datos_m3$PESO_SAT_H_BOCADOS~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$PESO_SAT_H_BOCADOS by TRATAMIENTO
## Bartlett's K-squared = 4.5748, df = 3, p-value = 0.2057bartlett.test(datos_m3$PESO_S_BOCADOS~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$PESO_S_BOCADOS by TRATAMIENTO
## Bartlett's K-squared = 5.8171, df = 3, p-value = 0.1209bartlett.test(datos_m3$CRA~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$CRA by TRATAMIENTO
## Bartlett's K-squared = 5.8279, df = 3, p-value = 0.1203bartlett.test(datos_m3$N_ESTOMAS_ADAX~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$N_ESTOMAS_ADAX by TRATAMIENTO
## Bartlett's K-squared = 7.4806, df = 3, p-value = 0.05806bartlett.test(datos_m3$N_ESTOMAS_ABA~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$N_ESTOMAS_ABA by TRATAMIENTO
## Bartlett's K-squared = 2.7366, df = 3, p-value = 0.434bartlett.test(datos_m3$N_ESTOMAS_ABIERTOS~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$N_ESTOMAS_ABIERTOS by TRATAMIENTO
## Bartlett's K-squared = 1.6845, df = 3, p-value = 0.6404bartlett.test(datos_m3$N_ESTOMAS_P_ABIERTOS~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$N_ESTOMAS_P_ABIERTOS by TRATAMIENTO
## Bartlett's K-squared = 2.3654, df = 3, p-value = 0.5001bartlett.test(datos_m3$N_ESTOMAS_CERRADOS~TRATAMIENTO, data = datos_m3)##
## Bartlett test of homogeneity of variances
##
## data: datos_m3$N_ESTOMAS_CERRADOS by TRATAMIENTO
## Bartlett's K-squared = 3.7196, df = 3, p-value = 0.2934#Se comprueba el suspuesto de homocedasticidad, los datos son homocedasticos. Homogeneidad de varianzasIntento_grafico_combinado <- read_excel("C:/Users/JuanSebH2/Downloads/Universidad/Fisiologia/FisioVeg/Intento_grafico_combinado.xlsx")
datos_c=data.frame(Intento_grafico_combinado);datos_c## MUESTREO TRATAMIENTO BLOQUES TEMPERATURA
## 1 1 (-DH-Si ) 1 18.40
## 2 1 (-DH-Si ) 1 18.90
## 3 1 (-DH-Si ) 2 20.36
## 4 1 (-DH-Si ) 2 19.00
## 5 1 (+DH-Si ) 1 19.70
## 6 1 (+DH-Si ) 1 16.90
## 7 1 (+DH-Si ) 2 18.90
## 8 1 (+DH-Si ) 2 20.00
## 9 1 (+DH+Si1) 1 18.40
## 10 1 (+DH+Si1) 1 17.00
## 11 1 (+DH+Si1) 2 17.90
## 12 1 (+DH+Si1) 2 20.50
## 13 1 (+DH+Si2) 1 18.70
## 14 1 (+DH+Si2) 1 19.70
## 15 1 (+DH+Si2) 2 21.20
## 16 1 (+DH+Si2) 2 20.90
## 17 2 (-DH-Si ) 1 20.42
## 18 2 (-DH-Si ) 1 20.88
## 19 2 (-DH-Si ) 2 20.56
## 20 2 (-DH-Si ) 2 17.84
## 21 2 (+DH-Si ) 1 22.04
## 22 2 (+DH-Si ) 1 20.66
## 23 2 (+DH-Si ) 2 20.82
## 24 2 (+DH-Si ) 2 20.60
## 25 2 (+DH+Si1) 1 21.08
## 26 2 (+DH+Si1) 1 20.88
## 27 2 (+DH+Si1) 2 21.38
## 28 2 (+DH+Si1) 2 20.82
## 29 2 (+DH+Si2) 1 19.98
## 30 2 (+DH+Si2) 1 21.22
## 31 2 (+DH+Si2) 2 21.26
## 32 2 (+DH+Si2) 2 20.68
## 33 3 (-DH-Si ) 1 17.24
## 34 3 (-DH-Si ) 1 16.48
## 35 3 (-DH-Si ) 2 16.28
## 36 3 (-DH-Si ) 2 16.08
## 37 3 (+DH-Si ) 1 18.46
## 38 3 (+DH-Si ) 1 18.30
## 39 3 (+DH-Si ) 2 18.48
## 40 3 (+DH-Si ) 2 18.38
## 41 3 (+DH+Si1) 1 17.94
## 42 3 (+DH+Si1) 1 18.13
## 43 3 (+DH+Si1) 2 18.12
## 44 3 (+DH+Si1) 2 17.86
## 45 3 (+DH+Si2) 1 18.38
## 46 3 (+DH+Si2) 1 17.54
## 47 3 (+DH+Si2) 2 17.26
## 48 3 (+DH+Si2) 2 16.90#ggplot(datos_c)+
#aes(x=datos_c$MUESTREO, y=datos_c$TEMPERATURA)+
#geom_bar()