datos_indices <-read_excel("C:/Users/JuanSebH2/Downloads/Universidad/Fisiologia/FisioVeg/DocumentoFinal_muestreos.xlsx",
sheet = "TRC", col_types = c("text",
"text", "numeric", "numeric", "numeric",
"numeric"))
datos_indices=data.frame(datos_indices); datos_indices
## TRATAMIENTOS DIAS.DESPUES.DE.LA.SIEMBRA TRC AFE IAF TAN
## 1 (-DH-Si ) 0 0.000000 0.0000 0.000000 0.000000
## 2 (+DH-Si ) 0 0.000000 0.0000 0.000000 0.000000
## 3 (+DH+Si1) 0 0.000000 0.0000 0.000000 0.000000
## 4 (+DH+Si2) 0 0.000000 0.0000 0.000000 0.000000
## 5 (-DH-Si ) 32 0.208810 153.1991 0.568875 0.001037
## 6 (+DH-Si ) 32 0.190034 148.4174 0.395500 0.000782
## 7 (+DH+Si1) 32 0.194576 142.9771 0.391500 0.000931
## 8 (+DH+Si2) 32 0.202723 179.2078 0.464275 0.001154
## 9 (-DH-Si ) 39 0.031721 305.1969 1.216250 0.000248
## 10 (+DH-Si ) 39 0.091648 282.1988 0.917000 0.000650
## 11 (+DH+Si1) 39 0.109602 252.9008 0.948875 0.000893
## 12 (+DH+Si2) 39 0.053480 298.2403 0.994650 0.000356
## 13 (-DH-Si ) 46 0.190000 303.7055 1.813095 0.002378
## 14 (+DH-Si ) 46 0.160000 268.2735 1.133344 0.001863
## 15 (+DH+Si1) 46 0.150000 243.5577 1.287140 0.002252
## 16 (+DH+Si2) 46 0.180000 242.1115 1.239005 0.002655
datos_indices$TRC=round(datos_indices$TRC,2)
datos_indices$AFE=round(datos_indices$AFE,2)
datos_indices$IAF=round(datos_indices$IAF,2)
datos_indices
## TRATAMIENTOS DIAS.DESPUES.DE.LA.SIEMBRA TRC AFE IAF TAN
## 1 (-DH-Si ) 0 0.00 0.00 0.00 0.000000
## 2 (+DH-Si ) 0 0.00 0.00 0.00 0.000000
## 3 (+DH+Si1) 0 0.00 0.00 0.00 0.000000
## 4 (+DH+Si2) 0 0.00 0.00 0.00 0.000000
## 5 (-DH-Si ) 32 0.21 153.20 0.57 0.001037
## 6 (+DH-Si ) 32 0.19 148.42 0.40 0.000782
## 7 (+DH+Si1) 32 0.19 142.98 0.39 0.000931
## 8 (+DH+Si2) 32 0.20 179.21 0.46 0.001154
## 9 (-DH-Si ) 39 0.03 305.20 1.22 0.000248
## 10 (+DH-Si ) 39 0.09 282.20 0.92 0.000650
## 11 (+DH+Si1) 39 0.11 252.90 0.95 0.000893
## 12 (+DH+Si2) 39 0.05 298.24 0.99 0.000356
## 13 (-DH-Si ) 46 0.19 303.71 1.81 0.002378
## 14 (+DH-Si ) 46 0.16 268.27 1.13 0.001863
## 15 (+DH+Si1) 46 0.15 243.56 1.29 0.002252
## 16 (+DH+Si2) 46 0.18 242.11 1.24 0.002655
tratamiento_colors <- c("#9E3D22", "#418E4D","#F5C35C", "#2F2F2F")
tratamiento_shapes <- c(16, 17, 18, 19)
ggplot(datos_indices) +
aes(DIAS.DESPUES.DE.LA.SIEMBRA, TRC, group = TRATAMIENTOS, fill = TRATAMIENTOS, shape = TRATAMIENTOS, color = TRATAMIENTOS) +
geom_line(size = 1) +
geom_point(size = 3) +
scale_fill_manual(values = tratamiento_colors) +
scale_shape_manual(values = tratamiento_shapes) +
scale_color_manual(values = tratamiento_colors) +
labs(title = "TASA RELATIVA DE CRECIMIENTO (f/día)",
x = "Días después de la siembra", y = "TRC") +
theme_update()+
theme(plot.title = element_text(hjust = 0.5, face = "bold", family = "Times New Roman"))
## Warning: Using `size` aesthetic for lines was deprecated in ggplot2 3.4.0.
## ℹ Please use `linewidth` instead.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
## Warning in grid.Call(C_stringMetric, as.graphicsAnnot(x$label)): font family
## not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
ggplot(datos_indices) +
aes(DIAS.DESPUES.DE.LA.SIEMBRA, datos_indices$AFE, group = TRATAMIENTOS, fill = TRATAMIENTOS, shape = TRATAMIENTOS, color = TRATAMIENTOS) +
geom_line(size = 1) +
geom_point(size = 3) +
scale_fill_manual(values = tratamiento_colors) +
scale_shape_manual(values = tratamiento_shapes) +
scale_color_manual(values = tratamiento_colors) +
labs(title = "ÁREA FOLIAR ESPECIFICA (cm^2/g)",
x = "Días después de la siembra", y = "AFE") +
theme_update()+
theme(plot.title = element_text(hjust = 0.5, face = "bold", family = "Times New Roman"))
## Warning: Use of `datos_indices$AFE` is discouraged.
## ℹ Use `AFE` instead.
## Warning: Use of `datos_indices$AFE` is discouraged.
## ℹ Use `AFE` instead.
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
ggplot(datos_indices) +
aes(DIAS.DESPUES.DE.LA.SIEMBRA, datos_indices$IAF, group = TRATAMIENTOS, fill = TRATAMIENTOS, shape = TRATAMIENTOS, color = TRATAMIENTOS) +
geom_line(size = 1) +
geom_point(size = 3) +
scale_fill_manual(values = tratamiento_colors) +
scale_shape_manual(values = tratamiento_shapes) +
scale_color_manual(values = tratamiento_colors) +
labs(title = "INDICE DE ÁREA FOLIAR (IAF)",
x = "Días después de la siembra", y = "IAF") +
theme_update()+
theme(plot.title = element_text(hjust = 0.5, face = "bold", family = "Times New Roman"))
## Warning: Use of `datos_indices$IAF` is discouraged.
## ℹ Use `IAF` instead.
## Warning: Use of `datos_indices$IAF` is discouraged.
## ℹ Use `IAF` instead.
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
ggplot(datos_indices) +
aes(DIAS.DESPUES.DE.LA.SIEMBRA, datos_indices$TAN, group = TRATAMIENTOS, fill = TRATAMIENTOS, shape = TRATAMIENTOS, color = TRATAMIENTOS) +
geom_line(size = 1) +
geom_point(size = 3) +
scale_fill_manual(values = tratamiento_colors) +
scale_shape_manual(values = tratamiento_shapes) +
scale_color_manual(values = tratamiento_colors) +
labs(title = "TASA DE ASIMILACIÓN NETA (g/cm^2*día)",
x = "Días después de la siembra", y = "TAN") +
theme_update()+
theme(plot.title = element_text(hjust = 0.5, face = "bold", family = "Times New Roman"))
## Warning: Use of `datos_indices$TAN` is discouraged.
## ℹ Use `TAN` instead.
## Warning: Use of `datos_indices$TAN` is discouraged.
## ℹ Use `TAN` instead.
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
## Warning in grid.Call(C_textBounds, as.graphicsAnnot(x$label), x$x, x$y, : font
## family not found in Windows font database
datos_indices_2 <- read_excel("C:/Users/JuanSebH2/Downloads/Universidad/Fisiologia/FisioVeg/DocumentoFinal_muestreos.xlsx",
sheet = "indices", col_types = c("text",
"text", "numeric", "numeric", "numeric",
"numeric"))
datos_indices_2=data.frame(datos_indices_2);datos_indices_2
## TRATAMIENTOS DIAS.DESPUES.DE.LA.SIEMBRA TRC AFE AIF
## 1 (-DH-Si ) 0 0.000000 0.0000 0.000000
## 2 (-DH-Si ) 0 0.000000 0.0000 0.000000
## 3 (-DH-Si ) 0 0.000000 0.0000 0.000000
## 4 (-DH-Si ) 0 0.000000 0.0000 0.000000
## 5 (+DH-Si ) 0 0.000000 0.0000 0.000000
## 6 (+DH-Si ) 0 0.000000 0.0000 0.000000
## 7 (+DH-Si ) 0 0.000000 0.0000 0.000000
## 8 (+DH-Si ) 0 0.000000 0.0000 0.000000
## 9 (+DH+Si1) 0 0.000000 0.0000 0.000000
## 10 (+DH+Si1) 0 0.000000 0.0000 0.000000
## 11 (+DH+Si1) 0 0.000000 0.0000 0.000000
## 12 (+DH+Si1) 0 0.000000 0.0000 0.000000
## 13 (+DH+Si2) 0 0.000000 0.0000 0.000000
## 14 (+DH+Si2) 0 0.000000 0.0000 0.000000
## 15 (+DH+Si2) 0 0.000000 0.0000 0.000000
## 16 (+DH+Si2) 0 0.000000 0.0000 0.000000
## 17 (-DH-Si ) 32 0.207700 161.9697 0.534500
## 18 (-DH-Si ) 32 0.206033 165.0000 0.462000
## 19 (-DH-Si ) 32 0.212226 121.4516 0.753000
## 20 (-DH-Si ) 32 0.209282 164.3750 0.526000
## 21 (+DH-Si ) 32 0.199379 160.3030 0.529000
## 22 (+DH-Si ) 32 0.194207 154.1667 0.370000
## 23 (+DH-Si ) 32 0.191601 129.2000 0.323000
## 24 (+DH-Si ) 32 0.174951 150.0000 0.360000
## 25 (+DH+Si1) 32 0.200421 123.7500 0.396000
## 26 (+DH+Si1) 32 0.190914 133.3333 0.400000
## 27 (+DH+Si1) 32 0.179268 144.2000 0.360500
## 28 (+DH+Si1) 32 0.207700 170.6250 0.409500
## 29 (+DH+Si2) 32 0.197748 154.2083 0.370100
## 30 (+DH+Si2) 32 0.174951 246.5385 0.320500
## 31 (+DH+Si2) 32 0.220505 150.5128 0.587000
## 32 (+DH+Si2) 32 0.217688 165.5714 0.579500
## 33 (-DH-Si ) 39 0.078934 313.3237 1.317000
## 34 (-DH-Si ) 39 -0.014617 323.9161 0.916500
## 35 (-DH-Si ) 39 0.006894 269.8212 1.299000
## 36 (-DH-Si ) 39 0.055672 313.7269 1.332500
## 37 (+DH-Si ) 39 0.080979 302.5530 1.098000
## 38 (+DH-Si ) 39 0.094964 262.4320 0.900500
## 39 (+DH-Si ) 39 0.054147 280.9187 0.808500
## 40 (+DH-Si ) 39 0.136502 282.8913 0.861000
## 41 (+DH+Si1) 39 0.133774 211.1714 0.952000
## 42 (+DH+Si1) 39 0.108835 246.1261 1.026000
## 43 (+DH+Si1) 39 0.198042 246.2175 0.942000
## 44 (+DH+Si1) 39 -0.002244 308.0881 0.875500
## 45 (+DH+Si2) 39 0.094090 266.8209 0.870100
## 46 (+DH+Si2) 39 0.158523 361.2962 0.841500
## 47 (+DH+Si2) 39 -0.032495 264.2209 1.163500
## 48 (+DH+Si2) 39 -0.006197 300.6230 1.103500
## 49 (-DH-Si ) 46 0.194392 262.0316 1.960110
## 50 (-DH-Si ) 46 0.210391 357.2040 1.487580
## 51 (-DH-Si ) 46 0.187408 307.5118 1.742750
## 52 (-DH-Si ) 46 0.165371 288.0745 2.061940
## 53 (+DH-Si ) 46 0.147529 297.7431 1.131885
## 54 (+DH-Si ) 46 0.150378 241.8142 1.309090
## 55 (+DH-Si ) 46 0.166403 301.9943 1.131685
## 56 (+DH-Si ) 46 0.157216 231.5425 0.960715
## 57 (+DH+Si1) 46 0.135025 211.9690 1.157565
## 58 (+DH+Si1) 46 0.170102 238.5266 1.440755
## 59 (+DH+Si1) 46 0.112217 222.3248 1.243190
## 60 (+DH+Si1) 46 0.187281 301.4105 1.307050
## 61 (+DH+Si2) 46 0.196165 212.7903 1.345500
## 62 (+DH+Si2) 46 0.151069 261.2818 1.343000
## 63 (+DH+Si2) 46 0.191113 249.0238 1.215190
## 64 (+DH+Si2) 46 0.172530 245.3501 1.052330
## TAN
## 1 0.000000
## 2 0.000000
## 3 0.000000
## 4 0.000000
## 5 0.000000
## 6 0.000000
## 7 0.000000
## 8 0.000000
## 9 0.000000
## 10 0.000000
## 11 0.000000
## 12 0.000000
## 13 0.000000
## 14 0.000000
## 15 0.000000
## 16 0.000000
## 17 0.001052
## 18 0.001117
## 19 0.000926
## 20 0.001120
## 21 0.000812
## 22 0.000909
## 23 0.000928
## 24 0.000501
## 25 0.001052
## 26 0.000770
## 27 0.000575
## 28 0.001294
## 29 0.001018
## 30 0.000548
## 31 0.001471
## 32 0.001358
## 33 0.000624
## 34 -0.000111
## 35 0.000049
## 36 0.000420
## 37 0.000586
## 38 0.000757
## 39 0.000380
## 40 0.000723
## 41 0.001433
## 42 0.000728
## 43 0.001437
## 44 -0.000020
## 45 0.000863
## 46 0.000950
## 47 -0.000290
## 48 -0.000058
## 49 0.002866
## 50 0.002245
## 51 0.002183
## 52 0.001838
## 53 0.002374
## 54 0.002003
## 55 0.001884
## 56 0.002095
## 57 0.003023
## 58 0.002201
## 59 0.001703
## 60 0.002310
## 61 0.003464
## 62 0.001665
## 63 0.003056
## 64 0.003232
names(datos_indices_2)
## [1] "TRATAMIENTOS" "DIAS.DESPUES.DE.LA.SIEMBRA"
## [3] "TRC" "AFE"
## [5] "AIF" "TAN"
m1=datos_indices_2[c(17:32),];m1
## TRATAMIENTOS DIAS.DESPUES.DE.LA.SIEMBRA TRC AFE AIF TAN
## 17 (-DH-Si ) 32 0.207700 161.9697 0.5345 0.001052
## 18 (-DH-Si ) 32 0.206033 165.0000 0.4620 0.001117
## 19 (-DH-Si ) 32 0.212226 121.4516 0.7530 0.000926
## 20 (-DH-Si ) 32 0.209282 164.3750 0.5260 0.001120
## 21 (+DH-Si ) 32 0.199379 160.3030 0.5290 0.000812
## 22 (+DH-Si ) 32 0.194207 154.1667 0.3700 0.000909
## 23 (+DH-Si ) 32 0.191601 129.2000 0.3230 0.000928
## 24 (+DH-Si ) 32 0.174951 150.0000 0.3600 0.000501
## 25 (+DH+Si1) 32 0.200421 123.7500 0.3960 0.001052
## 26 (+DH+Si1) 32 0.190914 133.3333 0.4000 0.000770
## 27 (+DH+Si1) 32 0.179268 144.2000 0.3605 0.000575
## 28 (+DH+Si1) 32 0.207700 170.6250 0.4095 0.001294
## 29 (+DH+Si2) 32 0.197748 154.2083 0.3701 0.001018
## 30 (+DH+Si2) 32 0.174951 246.5385 0.3205 0.000548
## 31 (+DH+Si2) 32 0.220505 150.5128 0.5870 0.001471
## 32 (+DH+Si2) 32 0.217688 165.5714 0.5795 0.001358
m2=datos_indices_2[c(33:48),];m2
## TRATAMIENTOS DIAS.DESPUES.DE.LA.SIEMBRA TRC AFE AIF TAN
## 33 (-DH-Si ) 39 0.078934 313.3237 1.3170 0.000624
## 34 (-DH-Si ) 39 -0.014617 323.9161 0.9165 -0.000111
## 35 (-DH-Si ) 39 0.006894 269.8212 1.2990 0.000049
## 36 (-DH-Si ) 39 0.055672 313.7269 1.3325 0.000420
## 37 (+DH-Si ) 39 0.080979 302.5530 1.0980 0.000586
## 38 (+DH-Si ) 39 0.094964 262.4320 0.9005 0.000757
## 39 (+DH-Si ) 39 0.054147 280.9187 0.8085 0.000380
## 40 (+DH-Si ) 39 0.136502 282.8913 0.8610 0.000723
## 41 (+DH+Si1) 39 0.133774 211.1714 0.9520 0.001433
## 42 (+DH+Si1) 39 0.108835 246.1261 1.0260 0.000728
## 43 (+DH+Si1) 39 0.198042 246.2175 0.9420 0.001437
## 44 (+DH+Si1) 39 -0.002244 308.0881 0.8755 -0.000020
## 45 (+DH+Si2) 39 0.094090 266.8209 0.8701 0.000863
## 46 (+DH+Si2) 39 0.158523 361.2962 0.8415 0.000950
## 47 (+DH+Si2) 39 -0.032495 264.2209 1.1635 -0.000290
## 48 (+DH+Si2) 39 -0.006197 300.6230 1.1035 -0.000058
m3=datos_indices_2[c(49:64),];m3
## TRATAMIENTOS DIAS.DESPUES.DE.LA.SIEMBRA TRC AFE AIF TAN
## 49 (-DH-Si ) 46 0.194392 262.0316 1.960110 0.002866
## 50 (-DH-Si ) 46 0.210391 357.2040 1.487580 0.002245
## 51 (-DH-Si ) 46 0.187408 307.5118 1.742750 0.002183
## 52 (-DH-Si ) 46 0.165371 288.0745 2.061940 0.001838
## 53 (+DH-Si ) 46 0.147529 297.7431 1.131885 0.002374
## 54 (+DH-Si ) 46 0.150378 241.8142 1.309090 0.002003
## 55 (+DH-Si ) 46 0.166403 301.9943 1.131685 0.001884
## 56 (+DH-Si ) 46 0.157216 231.5425 0.960715 0.002095
## 57 (+DH+Si1) 46 0.135025 211.9690 1.157565 0.003023
## 58 (+DH+Si1) 46 0.170102 238.5266 1.440755 0.002201
## 59 (+DH+Si1) 46 0.112217 222.3248 1.243190 0.001703
## 60 (+DH+Si1) 46 0.187281 301.4105 1.307050 0.002310
## 61 (+DH+Si2) 46 0.196165 212.7903 1.345500 0.003464
## 62 (+DH+Si2) 46 0.151069 261.2818 1.343000 0.001665
## 63 (+DH+Si2) 46 0.191113 249.0238 1.215190 0.003056
## 64 (+DH+Si2) 46 0.172530 245.3501 1.052330 0.003232
TRC_m1=aov(TRC~TRATAMIENTOS,data=m1)
summary(TRC_m1)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 0.0008402 0.0002801 1.566 0.249
## Residuals 12 0.0021459 0.0001788
TukeyHSD(TRC_m1)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TRC ~ TRATAMIENTOS, data = m1)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.01877575 -0.04684891 0.00929741 0.2462547
## (+DH+Si1)-(-DH-Si ) -0.01423450 -0.04230766 0.01383866 0.4642548
## (+DH+Si2)-(-DH-Si ) -0.00608725 -0.03416041 0.02198591 0.9157032
## (+DH+Si1)-(+DH-Si ) 0.00454125 -0.02353191 0.03261441 0.9620292
## (+DH+Si2)-(+DH-Si ) 0.01268850 -0.01538466 0.04076166 0.5558259
## (+DH+Si2)-(+DH+Si1) 0.00814725 -0.01992591 0.03622041 0.8241635
TukeyC::TukeyC(TRC_m1)
## Results
## Means G1
## (-DH-Si ) 0.21 a
## (+DH+Si2) 0.20 a
## (+DH+Si1) 0.19 a
## (+DH-Si ) 0.19 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (-DH-Si ) (+DH+Si2) (+DH+Si1) (+DH-Si )
## (-DH-Si ) 0.000 0.006 0.014 0.019
## (+DH+Si2) 0.916 0.000 0.008 0.013
## (+DH+Si1) 0.464 0.824 0.000 0.005
## (+DH-Si ) 0.246 0.556 0.962 0.000
##
## MSD
## (-DH-Si ) (+DH+Si2) (+DH+Si1) (+DH-Si )
## (-DH-Si ) 0.000 0.028 0.028 0.028
## (+DH+Si2) 0.028 0.000 0.028 0.028
## (+DH+Si1) 0.028 0.028 0.000 0.028
## (+DH-Si ) 0.028 0.028 0.028 0.000
TRC_m2=aov(TRC~TRATAMIENTOS,data=m2)
summary(TRC_m2)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 0.01506 0.005020 1.123 0.379
## Residuals 12 0.05366 0.004472
TukeyHSD(TRC_m2)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TRC ~ TRATAMIENTOS, data = m2)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.05992725 -0.08045817 0.20031267 0.5990387
## (+DH+Si1)-(-DH-Si ) 0.07788100 -0.06250442 0.21826642 0.3908802
## (+DH+Si2)-(-DH-Si ) 0.02175950 -0.11862592 0.16214492 0.9663117
## (+DH+Si1)-(+DH-Si ) 0.01795375 -0.12243167 0.15833917 0.9804864
## (+DH+Si2)-(+DH-Si ) -0.03816775 -0.17855317 0.10221767 0.8498601
## (+DH+Si2)-(+DH+Si1) -0.05612150 -0.19650692 0.08426392 0.6459936
TukeyC::TukeyC(TRC_m2)
## Results
## Means G1
## (+DH+Si1) 0.11 a
## (+DH-Si ) 0.09 a
## (+DH+Si2) 0.05 a
## (-DH-Si ) 0.03 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (+DH+Si1) (+DH-Si ) (+DH+Si2) (-DH-Si )
## (+DH+Si1) 0.000 0.018 0.056 0.078
## (+DH-Si ) 0.980 0.000 0.038 0.060
## (+DH+Si2) 0.646 0.850 0.000 0.022
## (-DH-Si ) 0.391 0.599 0.966 0.000
##
## MSD
## (+DH+Si1) (+DH-Si ) (+DH+Si2) (-DH-Si )
## (+DH+Si1) 0.00 0.14 0.14 0.14
## (+DH-Si ) 0.14 0.00 0.14 0.14
## (+DH+Si2) 0.14 0.14 0.00 0.14
## (-DH-Si ) 0.14 0.14 0.14 0.00
TRC_m3=aov(TRC~TRATAMIENTOS,data=m3)
summary(TRC_m3)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 0.003977 0.0013257 2.671 0.0948 .
## Residuals 12 0.005956 0.0004963
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
TukeyHSD(TRC_m3)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TRC ~ TRATAMIENTOS, data = m3)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.03400900 -0.08077768 0.012759683 0.1900561
## (+DH+Si1)-(-DH-Si ) -0.03823425 -0.08500293 0.008534433 0.1244047
## (+DH+Si2)-(-DH-Si ) -0.01167125 -0.05843993 0.035097433 0.8787429
## (+DH+Si1)-(+DH-Si ) -0.00422525 -0.05099393 0.042543433 0.9928886
## (+DH+Si2)-(+DH-Si ) 0.02233775 -0.02443093 0.069106433 0.5124691
## (+DH+Si2)-(+DH+Si1) 0.02656300 -0.02020568 0.073331683 0.3718170
TukeyC::TukeyC(TRC_m3)
## Results
## Means G1
## (-DH-Si ) 0.19 a
## (+DH+Si2) 0.18 a
## (+DH-Si ) 0.16 a
## (+DH+Si1) 0.15 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (-DH-Si ) (+DH+Si2) (+DH-Si ) (+DH+Si1)
## (-DH-Si ) 0.000 0.012 0.034 0.038
## (+DH+Si2) 0.879 0.000 0.022 0.027
## (+DH-Si ) 0.190 0.512 0.000 0.004
## (+DH+Si1) 0.124 0.372 0.993 0.000
##
## MSD
## (-DH-Si ) (+DH+Si2) (+DH-Si ) (+DH+Si1)
## (-DH-Si ) 0.000 0.047 0.047 0.047
## (+DH+Si2) 0.047 0.000 0.047 0.047
## (+DH-Si ) 0.047 0.047 0.000 0.047
## (+DH+Si1) 0.047 0.047 0.047 0.000
AFE_m1=aov(AFE~TRATAMIENTOS,data=m1)
summary(AFE_m1)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 3094 1031.4 1.332 0.31
## Residuals 12 9291 774.3
TukeyHSD(AFE_m1)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AFE ~ TRATAMIENTOS, data = m1)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -4.781652 -63.19767 53.63436 0.9946790
## (+DH+Si1)-(-DH-Si ) -10.221995 -68.63801 48.19402 0.9527762
## (+DH+Si2)-(-DH-Si ) 26.008683 -32.40733 84.42470 0.5673862
## (+DH+Si1)-(+DH-Si ) -5.440343 -63.85636 52.97567 0.9922261
## (+DH+Si2)-(+DH-Si ) 30.790335 -27.62568 89.20635 0.4326410
## (+DH+Si2)-(+DH+Si1) 36.230678 -22.18534 94.64669 0.3022206
TukeyC::TukeyC(AFE_m1)
## Results
## Means G1
## (+DH+Si2) 179.21 a
## (-DH-Si ) 153.20 a
## (+DH-Si ) 148.42 a
## (+DH+Si1) 142.98 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (+DH+Si2) (-DH-Si ) (+DH-Si ) (+DH+Si1)
## (+DH+Si2) 0.000 26.009 30.790 36.231
## (-DH-Si ) 0.567 0.000 4.782 10.222
## (+DH-Si ) 0.433 0.995 0.000 5.440
## (+DH+Si1) 0.302 0.953 0.992 0.000
##
## MSD
## (+DH+Si2) (-DH-Si ) (+DH-Si ) (+DH+Si1)
## (+DH+Si2) 0.000 58.416 58.416 58.416
## (-DH-Si ) 58.416 0.000 58.416 58.416
## (+DH-Si ) 58.416 58.416 0.000 58.416
## (+DH+Si1) 58.416 58.416 58.416 0.000
AFE_m2=aov(AFE~TRATAMIENTOS,data=m2)
summary(AFE_m2)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 6484 2161 1.914 0.181
## Residuals 12 13552 1129
TukeyHSD(AFE_m2)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AFE ~ TRATAMIENTOS, data = m2)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -22.998193 -93.54619 47.54981 0.7697199
## (+DH+Si1)-(-DH-Si ) -52.296148 -122.84415 18.25185 0.1781814
## (+DH+Si2)-(-DH-Si ) -6.956693 -77.50469 63.59131 0.9908113
## (+DH+Si1)-(+DH-Si ) -29.297955 -99.84595 41.25004 0.6191052
## (+DH+Si2)-(+DH-Si ) 16.041500 -54.50650 86.58950 0.9045222
## (+DH+Si2)-(+DH+Si1) 45.339455 -25.20854 115.88745 0.2753006
TukeyC::TukeyC(AFE_m2)
## Results
## Means G1
## (-DH-Si ) 305.20 a
## (+DH+Si2) 298.24 a
## (+DH-Si ) 282.20 a
## (+DH+Si1) 252.90 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (-DH-Si ) (+DH+Si2) (+DH-Si ) (+DH+Si1)
## (-DH-Si ) 0.000 6.957 22.998 52.296
## (+DH+Si2) 0.991 0.000 16.041 45.339
## (+DH-Si ) 0.770 0.905 0.000 29.298
## (+DH+Si1) 0.178 0.275 0.619 0.000
##
## MSD
## (-DH-Si ) (+DH+Si2) (+DH-Si ) (+DH+Si1)
## (-DH-Si ) 0.000 70.548 70.548 70.548
## (+DH+Si2) 70.548 0.000 70.548 70.548
## (+DH-Si ) 70.548 70.548 0.000 70.548
## (+DH+Si1) 70.548 70.548 70.548 0.000
AFE_m3=aov(AFE~TRATAMIENTOS,data=m3)
summary(AFE_m3)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 9964 3321 2.654 0.0961 .
## Residuals 12 15019 1252
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
TukeyHSD(AFE_m3)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AFE ~ TRATAMIENTOS, data = m3)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -35.431950 -109.70114 38.83724 0.5133830
## (+DH+Si1)-(-DH-Si ) -60.147757 -134.41695 14.12143 0.1290704
## (+DH+Si2)-(-DH-Si ) -61.593995 -135.86319 12.67520 0.1174995
## (+DH+Si1)-(+DH-Si ) -24.715808 -98.98500 49.55338 0.7588473
## (+DH+Si2)-(+DH-Si ) -26.162045 -100.43124 48.10715 0.7269378
## (+DH+Si2)-(+DH+Si1) -1.446237 -75.71543 72.82295 0.9999265
TukeyC::TukeyC(AFE_m3)
## Results
## Means G1
## (-DH-Si ) 303.71 a
## (+DH-Si ) 268.27 a
## (+DH+Si1) 243.56 a
## (+DH+Si2) 242.11 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (-DH-Si ) (+DH-Si ) (+DH+Si1) (+DH+Si2)
## (-DH-Si ) 0.000 35.432 60.148 61.594
## (+DH-Si ) 0.513 0.000 24.716 26.162
## (+DH+Si1) 0.129 0.759 0.000 1.446
## (+DH+Si2) 0.117 0.727 1.000 0.000
##
## MSD
## (-DH-Si ) (+DH-Si ) (+DH+Si1) (+DH+Si2)
## (-DH-Si ) 0.000 74.269 74.269 74.269
## (+DH-Si ) 74.269 0.000 74.269 74.269
## (+DH+Si1) 74.269 74.269 0.000 74.269
## (+DH+Si2) 74.269 74.269 74.269 0.000
AIF_m1=aov(AIF~TRATAMIENTOS,data=m1)
summary(AIF_m1)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 0.0825 0.02750 2.489 0.11
## Residuals 12 0.1326 0.01105
TukeyHSD(AIF_m1)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AIF ~ TRATAMIENTOS, data = m1)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.173375 -0.3940454 0.04729537 0.1448428
## (+DH+Si1)-(-DH-Si ) -0.177375 -0.3980454 0.04329537 0.1328732
## (+DH+Si2)-(-DH-Si ) -0.104600 -0.3252704 0.11607037 0.5185108
## (+DH+Si1)-(+DH-Si ) -0.004000 -0.2246704 0.21667037 0.9999407
## (+DH+Si2)-(+DH-Si ) 0.068775 -0.1518954 0.28944537 0.7921348
## (+DH+Si2)-(+DH+Si1) 0.072775 -0.1478954 0.29344537 0.7636610
TukeyC::TukeyC(AIF_m1)
## Results
## Means G1
## (-DH-Si ) 0.57 a
## (+DH+Si2) 0.46 a
## (+DH-Si ) 0.40 a
## (+DH+Si1) 0.39 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (-DH-Si ) (+DH+Si2) (+DH-Si ) (+DH+Si1)
## (-DH-Si ) 0.000 0.105 0.173 0.177
## (+DH+Si2) 0.519 0.000 0.069 0.073
## (+DH-Si ) 0.145 0.792 0.000 0.004
## (+DH+Si1) 0.133 0.764 1.000 0.000
##
## MSD
## (-DH-Si ) (+DH+Si2) (+DH-Si ) (+DH+Si1)
## (-DH-Si ) 0.000 0.221 0.221 0.221
## (+DH+Si2) 0.221 0.000 0.221 0.221
## (+DH-Si ) 0.221 0.221 0.000 0.221
## (+DH+Si1) 0.221 0.221 0.221 0.000
AIF_m2=aov(AIF~TRATAMIENTOS,data=m2)
summary(AIF_m2)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 0.2193 0.07310 3.386 0.054 .
## Residuals 12 0.2590 0.02158
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
TukeyHSD(AIF_m2)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AIF ~ TRATAMIENTOS, data = m2)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.299250 -0.6076796 0.009179634 0.0582129
## (+DH+Si1)-(-DH-Si ) -0.267375 -0.5758046 0.041054634 0.0977874
## (+DH+Si2)-(-DH-Si ) -0.221600 -0.5300296 0.086829634 0.1976974
## (+DH+Si1)-(+DH-Si ) 0.031875 -0.2765546 0.340304634 0.9894636
## (+DH+Si2)-(+DH-Si ) 0.077650 -0.2307796 0.386079634 0.8760142
## (+DH+Si2)-(+DH+Si1) 0.045775 -0.2626546 0.354204634 0.9701891
TukeyC::TukeyC(AIF_m2)
## Results
## Means G1
## (-DH-Si ) 1.22 a
## (+DH+Si2) 0.99 a
## (+DH+Si1) 0.95 a
## (+DH-Si ) 0.92 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (-DH-Si ) (+DH+Si2) (+DH+Si1) (+DH-Si )
## (-DH-Si ) 0.000 0.222 0.267 0.299
## (+DH+Si2) 0.198 0.000 0.046 0.078
## (+DH+Si1) 0.098 0.970 0.000 0.032
## (+DH-Si ) 0.058 0.876 0.989 0.000
##
## MSD
## (-DH-Si ) (+DH+Si2) (+DH+Si1) (+DH-Si )
## (-DH-Si ) 0.000 0.308 0.308 0.308
## (+DH+Si2) 0.308 0.000 0.308 0.308
## (+DH+Si1) 0.308 0.308 0.000 0.308
## (+DH-Si ) 0.308 0.308 0.308 0.000
AIF_m3=aov(AIF~TRATAMIENTOS,data=m3)
summary(AIF_m3)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 1.1054 0.3685 12.44 0.000541 ***
## Residuals 12 0.3554 0.0296
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
TukeyHSD(AIF_m3)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = AIF ~ TRATAMIENTOS, data = m3)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.6797512 -1.0410476 -0.3184549 0.0005937
## (+DH+Si1)-(-DH-Si ) -0.5259550 -0.8872513 -0.1646587 0.0047408
## (+DH+Si2)-(-DH-Si ) -0.5740900 -0.9353863 -0.2127937 0.0024267
## (+DH+Si1)-(+DH-Si ) 0.1537962 -0.2075001 0.5150926 0.6011112
## (+DH+Si2)-(+DH-Si ) 0.1056612 -0.2556351 0.4669576 0.8209195
## (+DH+Si2)-(+DH+Si1) -0.0481350 -0.4094313 0.3131613 0.9780622
TukeyC::TukeyC(AIF_m3)
## Results
## Means G1 G2
## (-DH-Si ) 1.81 a
## (+DH+Si1) 1.29 b
## (+DH+Si2) 1.24 b
## (+DH-Si ) 1.13 b
##
## Sig.level
## 0.05
##
## Diff_Prob
## (-DH-Si ) (+DH+Si1) (+DH+Si2) (+DH-Si )
## (-DH-Si ) 0.000 0.526 0.574 0.680
## (+DH+Si1) 0.005 0.000 0.048 0.154
## (+DH+Si2) 0.002 0.978 0.000 0.106
## (+DH-Si ) 0.001 0.601 0.821 0.000
##
## MSD
## (-DH-Si ) (+DH+Si1) (+DH+Si2) (+DH-Si )
## (-DH-Si ) 0.000 0.361 0.361 0.361
## (+DH+Si1) 0.361 0.000 0.361 0.361
## (+DH+Si2) 0.361 0.361 0.000 0.361
## (+DH-Si ) 0.361 0.361 0.361 0.000
TAN_m1=aov(TAN~TRATAMIENTOS,data=m1)
summary(TAN_m1)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 2.362e-07 7.874e-08 0.988 0.431
## Residuals 12 9.563e-07 7.969e-08
TukeyHSD(TAN_m1)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TAN ~ TRATAMIENTOS, data = m1)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.00026625 -0.0008588908 0.0003263908 0.5604788
## (+DH+Si1)-(-DH-Si ) -0.00013100 -0.0007236408 0.0004616408 0.9113281
## (+DH+Si2)-(-DH-Si ) 0.00004500 -0.0005476408 0.0006376408 0.9957364
## (+DH+Si1)-(+DH-Si ) 0.00013525 -0.0004573908 0.0007278908 0.9036104
## (+DH+Si2)-(+DH-Si ) 0.00031125 -0.0002813908 0.0009038908 0.4355846
## (+DH+Si2)-(+DH+Si1) 0.00017600 -0.0004166408 0.0007686408 0.8142826
TukeyC::TukeyC(TAN_m1)
## Results
## Means G1
## (+DH+Si2) 0.00 a
## (-DH-Si ) 0.00 a
## (+DH+Si1) 0.00 a
## (+DH-Si ) 0.00 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (+DH+Si2) (-DH-Si ) (+DH+Si1) (+DH-Si )
## (+DH+Si2) 0.000 0.000 0.000 0
## (-DH-Si ) 0.996 0.000 0.000 0
## (+DH+Si1) 0.814 0.911 0.000 0
## (+DH-Si ) 0.436 0.560 0.904 0
##
## MSD
## (+DH+Si2) (-DH-Si ) (+DH+Si1) (+DH-Si )
## (+DH+Si2) 0.000 0.001 0.001 0.001
## (-DH-Si ) 0.001 0.000 0.001 0.001
## (+DH+Si1) 0.001 0.001 0.000 0.001
## (+DH-Si ) 0.001 0.001 0.001 0.000
TAN_m2=aov(TAN~TRATAMIENTOS,data=m2)
summary(TAN_m2)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 9.890e-07 3.297e-07 1.287 0.323
## Residuals 12 3.074e-06 2.561e-07
TukeyHSD(TAN_m2)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TAN ~ TRATAMIENTOS, data = m2)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) 0.00036600 -0.0006964905 0.0014284905 0.7398178
## (+DH+Si1)-(-DH-Si ) 0.00064900 -0.0004134905 0.0017114905 0.3140148
## (+DH+Si2)-(-DH-Si ) 0.00012075 -0.0009417405 0.0011832405 0.9861146
## (+DH+Si1)-(+DH-Si ) 0.00028300 -0.0007794905 0.0013454905 0.8572622
## (+DH+Si2)-(+DH-Si ) -0.00024525 -0.0013077405 0.0008172405 0.9007202
## (+DH+Si2)-(+DH+Si1) -0.00052825 -0.0015907405 0.0005342405 0.4802068
TukeyC::TukeyC(TAN_m2)
## Results
## Means G1
## (+DH+Si1) 0.00 a
## (+DH-Si ) 0.00 a
## (+DH+Si2) 0.00 a
## (-DH-Si ) 0.00 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (+DH+Si1) (+DH-Si ) (+DH+Si2) (-DH-Si )
## (+DH+Si1) 0.000 0.000 0.001 0.001
## (+DH-Si ) 0.857 0.000 0.000 0.000
## (+DH+Si2) 0.480 0.901 0.000 0.000
## (-DH-Si ) 0.314 0.740 0.986 0.000
##
## MSD
## (+DH+Si1) (+DH-Si ) (+DH+Si2) (-DH-Si )
## (+DH+Si1) 0.000 0.001 0.001 0.001
## (+DH-Si ) 0.001 0.000 0.001 0.001
## (+DH+Si2) 0.001 0.001 0.000 0.001
## (-DH-Si ) 0.001 0.001 0.001 0.000
TAN_m3=aov(TAN~TRATAMIENTOS,data=m3)
summary(TAN_m3)
## Df Sum Sq Mean Sq F value Pr(>F)
## TRATAMIENTOS 3 1.296e-06 4.319e-07 1.465 0.273
## Residuals 12 3.538e-06 2.949e-07
TukeyHSD(TAN_m3)
## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = TAN ~ TRATAMIENTOS, data = m3)
##
## $TRATAMIENTOS
## diff lwr upr p adj
## (+DH-Si )-(-DH-Si ) -0.00019400 -0.0013339455 0.0009459455 0.9562736
## (+DH+Si1)-(-DH-Si ) 0.00002625 -0.0011136955 0.0011661955 0.9998785
## (+DH+Si2)-(-DH-Si ) 0.00057125 -0.0005686955 0.0017111955 0.4738110
## (+DH+Si1)-(+DH-Si ) 0.00022025 -0.0009196955 0.0013601955 0.9380823
## (+DH+Si2)-(+DH-Si ) 0.00076525 -0.0003746955 0.0019051955 0.2436164
## (+DH+Si2)-(+DH+Si1) 0.00054500 -0.0005949455 0.0016849455 0.5116829
TukeyC::TukeyC(TAN_m3)
## Results
## Means G1
## (+DH+Si2) 0.00 a
## (+DH+Si1) 0.00 a
## (-DH-Si ) 0.00 a
## (+DH-Si ) 0.00 a
##
## Sig.level
## 0.05
##
## Diff_Prob
## (+DH+Si2) (+DH+Si1) (-DH-Si ) (+DH-Si )
## (+DH+Si2) 0.000 0.001 0.001 0.001
## (+DH+Si1) 0.512 0.000 0.000 0.000
## (-DH-Si ) 0.474 1.000 0.000 0.000
## (+DH-Si ) 0.244 0.938 0.956 0.000
##
## MSD
## (+DH+Si2) (+DH+Si1) (-DH-Si ) (+DH-Si )
## (+DH+Si2) 0.000 0.001 0.001 0.001
## (+DH+Si1) 0.001 0.000 0.001 0.001
## (-DH-Si ) 0.001 0.001 0.000 0.001
## (+DH-Si ) 0.001 0.001 0.001 0.000