library(agricolae)
library(readxl)
P5 <- read_excel("C:/Users/Johan Forigua/Documents/Universidad Nacional/Semestre 10/Poscosecha/Practica 5/P5.xlsx")
P5$Trt <- as.factor(P5$Trt)
MedpH<- tapply(P5$pH, P5$Trt, mean);MedpH
## AIM AVAP C IIM IVAP
## 4.620000 4.760000 4.756667 4.550000 18.366667
SdpH <- tapply(P5$pH, P5$Trt, sd);SdpH
## AIM AVAP C IIM IVAP
## 0.01000000 0.05291503 0.04041452 0.03000000 1.15902258
MedTex <- tapply(P5$Textura, P5$Trt, mean);MedTex
## AIM AVAP C IIM IVAP
## 9.966667 10.633333 14.200000 14.166667 18.366667
SdTex <- tapply(P5$Textura, P5$Trt, sd);SdTex
## AIM AVAP C IIM IVAP
## 0.8736895 1.6165808 1.7349352 2.3671361 1.1590226
MedRm <- tapply(P5$RM, P5$Trt, mean);MedRm
## AIM AVAP C IIM IVAP
## 11.173333 15.596667 10.956667 6.873333 9.940000
Sdrm <- tapply(P5$RM, P5$Trt, sd);Sdrm
## AIM AVAP C IIM IVAP
## 0.2830783 1.4779829 0.4080850 0.5173329 0.4503332
MedAc <- tapply(P5$Aci, P5$Trt, mean);MedAc
## AIM AVAP C IIM IVAP
## 0.3966667 0.2966667 0.3833333 0.5833333 0.5033333
SdAc <- tapply(P5$Aci, P5$Trt, sd);SdAc
## AIM AVAP C IIM IVAP
## 0.005773503 0.028867513 0.005773503 0.035118846 0.015275252
MedBr <- tapply(P5$brix, P5$Trt, mean);MedBr
## AIM AVAP C IIM IVAP
## 4.400000 4.600000 4.200000 4.000000 4.966667
SdBr <- tapply(P5$brix, P5$Trt, sd);SdBr
## AIM AVAP C IIM IVAP
## 0.1000000 0.1000000 0.1000000 0.1000000 0.1154701
boxplot(P5$pH ~ P5$Trt, col=c("red","pink","darkgreen", "yellow", "orange"),
ylim=c(4.4,22), xlab = "Tratamientos", ylab = "")
MdpH <- aov(pH~Trt,data = P5)
summary(MdpH)
## Df Sum Sq Mean Sq F value Pr(>F)
## Trt 4 450.2 112.56 417.3 4.47e-11 ***
## Residuals 10 2.7 0.27
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Sp <- HSD.test(MdpH,"Trt");Sp
## $statistics
## MSerror Df Mean CV MSD
## 0.2697533 10 7.410667 7.008517 1.39565
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey Trt 5 4.654293 0.05
##
## $means
## pH std r Min Max Q25 Q50 Q75
## AIM 4.620000 0.01000000 3 4.61 4.63 4.615 4.62 4.625
## AVAP 4.760000 0.05291503 3 4.70 4.80 4.740 4.78 4.790
## C 4.756667 0.04041452 3 4.72 4.80 4.735 4.75 4.775
## IIM 4.550000 0.03000000 3 4.52 4.58 4.535 4.55 4.565
## IVAP 18.366667 1.15902258 3 17.60 19.70 17.700 17.80 18.750
##
## $comparison
## NULL
##
## $groups
## pH groups
## IVAP 18.366667 a
## AVAP 4.760000 b
## C 4.756667 b
## AIM 4.620000 b
## IIM 4.550000 b
##
## attr(,"class")
## [1] "group"
text(c(1:5),c(MedpH+SdpH+1.5), labels = c("b","b","b","b","a") )
boxplot(P5$Textura~P5$Trt, col=c("red","pink","darkgreen", "yellow", "orange"),
ylim=c(9.01,20.5), xlab = "Tratamientos",ylab = "N")
Md <- aov(Textura~Trt,data = P5)
summary(Md)
## Df Sum Sq Mean Sq F value Pr(>F)
## Trt 4 135.95 33.99 12.74 0.000614 ***
## Residuals 10 26.67 2.67
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Sp <- HSD.test(Md,"Trt");Sp
## $statistics
## MSerror Df Mean CV MSD
## 2.666667 10 13.46667 12.12619 4.38811
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey Trt 5 4.654293 0.05
##
## $means
## Textura std r Min Max Q25 Q50 Q75
## AIM 9.966667 0.8736895 3 9.0 10.7 9.6 10.2 10.45
## AVAP 10.633333 1.6165808 3 9.7 12.5 9.7 9.7 11.10
## C 14.200000 1.7349352 3 13.1 16.2 13.2 13.3 14.75
## IIM 14.166667 2.3671361 3 12.8 16.9 12.8 12.8 14.85
## IVAP 18.366667 1.1590226 3 17.6 19.7 17.7 17.8 18.75
##
## $comparison
## NULL
##
## $groups
## Textura groups
## IVAP 18.366667 a
## C 14.200000 ab
## IIM 14.166667 ab
## AVAP 10.633333 b
## AIM 9.966667 b
##
## attr(,"class")
## [1] "group"
text(c(1:5),c(MedTex+SdTex+1), labels = c("b","b","ab","ab","a") )
boxplot(P5$RM~P5$Trt, col=c("red","pink","darkgreen", "yellow", "orange"),
ylim=c(6.36,18.07), xlab = "Tratamientos",ylab = "")
Md <- aov(RM~Trt,data = P5)
summary(Md)
## Df Sum Sq Mean Sq F value Pr(>F)
## Trt 4 117.8 29.45 50.76 1.31e-06 ***
## Residuals 10 5.8 0.58
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Sp <- HSD.test(Md,"Trt");Sp
## $statistics
## MSerror Df Mean CV MSD
## 0.5803067 10 10.908 6.983669 2.047019
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey Trt 5 4.654293 0.05
##
## $means
## RM std r Min Max Q25 Q50 Q75
## AIM 11.173333 0.2830783 3 11.00 11.50 11.010 11.02 11.260
## AVAP 15.596667 1.4779829 3 13.94 16.78 15.005 16.07 16.425
## C 10.956667 0.4080850 3 10.51 11.31 10.780 11.05 11.180
## IIM 6.873333 0.5173329 3 6.45 7.45 6.585 6.72 7.085
## IVAP 9.940000 0.4503332 3 9.42 10.20 9.810 10.20 10.200
##
## $comparison
## NULL
##
## $groups
## RM groups
## AVAP 15.596667 a
## AIM 11.173333 b
## C 10.956667 b
## IVAP 9.940000 b
## IIM 6.873333 c
##
## attr(,"class")
## [1] "group"
text(c(1:5),c(MedRm+Sdrm+.7), labels = c("b","a","b","c","b") )
boxplot(P5$Aci~P5$Trt, col=c("red","pink","darkgreen", "yellow", "orange"),
ylim=c(0.267,0.681), xlab = "Tratamientos",ylab = "")
Md <- aov(Aci~Trt,data = P5)
summary(Md)
## Df Sum Sq Mean Sq F value Pr(>F)
## Trt 4 0.14976 0.03744 79.1 1.58e-07 ***
## Residuals 10 0.00473 0.00047
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Sp <- HSD.test(Md,"Trt");Sp
## $statistics
## MSerror Df Mean CV MSD
## 0.0004733333 10 0.4326667 5.028403 0.0584624
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey Trt 5 4.654293 0.05
##
## $means
## Aci std r Min Max Q25 Q50 Q75
## AIM 0.3966667 0.005773503 3 0.39 0.40 0.395 0.40 0.400
## AVAP 0.2966667 0.028867513 3 0.28 0.33 0.280 0.28 0.305
## C 0.3833333 0.005773503 3 0.38 0.39 0.380 0.38 0.385
## IIM 0.5833333 0.035118846 3 0.55 0.62 0.565 0.58 0.600
## IVAP 0.5033333 0.015275252 3 0.49 0.52 0.495 0.50 0.510
##
## $comparison
## NULL
##
## $groups
## Aci groups
## IIM 0.5833333 a
## IVAP 0.5033333 b
## AIM 0.3966667 c
## C 0.3833333 c
## AVAP 0.2966667 d
##
## attr(,"class")
## [1] "group"
text(c(1:5),c(MedAc+SdAc+0.04), labels = c("c","d","c","a","b"))
boxplot(P5$brix~P5$Trt, col=c("red","pink","darkgreen", "yellow", "orange"),
ylim=c(3.9,5.3), xlab = "Tratamientos",ylab = "")
Md <- aov(brix~Trt,data = P5)
summary(Md)
## Df Sum Sq Mean Sq F value Pr(>F)
## Trt 4 1.6667 0.4167 39.06 4.49e-06 ***
## Residuals 10 0.1067 0.0107
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Sp <- HSD.test(Md,"Trt");Sp
## $statistics
## MSerror Df Mean CV MSD
## 0.01066667 10 4.433333 2.329614 0.2775284
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey Trt 5 4.654293 0.05
##
## $means
## brix std r Min Max Q25 Q50 Q75
## AIM 4.400000 0.1000000 3 4.3 4.5 4.35 4.4 4.45
## AVAP 4.600000 0.1000000 3 4.5 4.7 4.55 4.6 4.65
## C 4.200000 0.1000000 3 4.1 4.3 4.15 4.2 4.25
## IIM 4.000000 0.1000000 3 3.9 4.1 3.95 4.0 4.05
## IVAP 4.966667 0.1154701 3 4.9 5.1 4.90 4.9 5.00
##
## $comparison
## NULL
##
## $groups
## brix groups
## IVAP 4.966667 a
## AVAP 4.600000 b
## AIM 4.400000 bc
## C 4.200000 cd
## IIM 4.000000 d
##
## attr(,"class")
## [1] "group"
text(c(1:5),c(MedBr-SdBr+.35), labels = c("bc","b","cd","d","a"))
