library(readxl)
df1 <- read_excel("C:/Users/ginna/Downloads/practicarepro.xlsx",
sheet = "albahaca")
#View(practicarepro)
Anova para Albahaca
summary(df1)
## tratamiento nraices pesoraiz longraiz
## Length:20 Min. :10.00 Min. :0.0050 Min. : 5.00
## Class :character 1st Qu.:15.00 1st Qu.:0.1700 1st Qu.: 7.25
## Mode :character Median :16.50 Median :0.2800 Median : 8.50
## Mean :25.60 Mean :0.3405 Mean : 8.55
## 3rd Qu.:28.75 3rd Qu.:0.3225 3rd Qu.: 9.75
## Max. :79.00 Max. :1.3200 Max. :14.50
## NA's :10 NA's :10 NA's :10
library(agricolae)
## Warning: package 'agricolae' was built under R version 4.0.5
tabla1=aov(nraices~tratamiento, data=df1)
summary(tabla1)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 1367 455.6 1.176 0.394
## Residuals 6 2326 387.6
## 10 observations deleted due to missingness
out1<-duncan.test(tabla1,"tratamiento")
out1
## $statistics
## MSerror Df Mean CV
## 387.5778 6 25.6 76.90232
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## NULL
##
## $means
## nraices std r Min Max Q25 Q50 Q75
## TRT1 10.00000 NA 1 10 10 10.0 10 10
## TRT2 37.20000 24.1080899 5 17 79 28.0 29 33
## TRT3 14.66667 0.5773503 3 14 15 14.5 15 15
## TRT4 16.00000 NA 1 16 16 16.0 16 16
##
## $comparison
## NULL
##
## $groups
## nraices groups
## TRT2 37.20000 a
## TRT4 16.00000 a
## TRT3 14.66667 a
## TRT1 10.00000 a
##
## attr(,"class")
## [1] "group"
ttest1<- HSD.test(tabla1,"tratamiento")
ttest1
## $statistics
## MSerror Df Mean CV
## 387.5778 6 25.6 76.90232
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 4.895599 0.05
##
## $means
## nraices std r Min Max Q25 Q50 Q75
## TRT1 10.00000 NA 1 10 10 10.0 10 10
## TRT2 37.20000 24.1080899 5 17 79 28.0 29 33
## TRT3 14.66667 0.5773503 3 14 15 14.5 15 15
## TRT4 16.00000 NA 1 16 16 16.0 16 16
##
## $comparison
## NULL
##
## $groups
## nraices groups
## TRT2 37.20000 a
## TRT4 16.00000 a
## TRT3 14.66667 a
## TRT1 10.00000 a
##
## attr(,"class")
## [1] "group"
tabla2=aov(pesoraiz~tratamiento, data=df1)
summary(tabla2)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 0.3385 0.1128 0.808 0.534
## Residuals 6 0.8377 0.1396
## 10 observations deleted due to missingness
out2<-duncan.test(tabla2,"tratamiento")
out2
## $statistics
## MSerror Df Mean CV
## 0.1396161 6 0.3405 109.7364
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## NULL
##
## $means
## pesoraiz std r Min Max Q25 Q50 Q75
## TRT1 0.1300000 NA 1 0.130 0.13 0.1300 0.13 0.13
## TRT2 0.5180000 0.4504664 5 0.260 1.32 0.3000 0.33 0.38
## TRT3 0.1283333 0.1140541 3 0.005 0.23 0.0775 0.15 0.19
## TRT4 0.3000000 NA 1 0.300 0.30 0.3000 0.30 0.30
##
## $comparison
## NULL
##
## $groups
## pesoraiz groups
## TRT2 0.5180000 a
## TRT4 0.3000000 a
## TRT1 0.1300000 a
## TRT3 0.1283333 a
##
## attr(,"class")
## [1] "group"
ttest2<- HSD.test(tabla2,"tratamiento")
ttest2
## $statistics
## MSerror Df Mean CV
## 0.1396161 6 0.3405 109.7364
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 4.895599 0.05
##
## $means
## pesoraiz std r Min Max Q25 Q50 Q75
## TRT1 0.1300000 NA 1 0.130 0.13 0.1300 0.13 0.13
## TRT2 0.5180000 0.4504664 5 0.260 1.32 0.3000 0.33 0.38
## TRT3 0.1283333 0.1140541 3 0.005 0.23 0.0775 0.15 0.19
## TRT4 0.3000000 NA 1 0.300 0.30 0.3000 0.30 0.30
##
## $comparison
## NULL
##
## $groups
## pesoraiz groups
## TRT2 0.5180000 a
## TRT4 0.3000000 a
## TRT1 0.1300000 a
## TRT3 0.1283333 a
##
## attr(,"class")
## [1] "group"
tabla3=aov(longraiz~tratamiento, data=df1)
summary(tabla3)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 7.36 2.453 0.242 0.864
## Residuals 6 60.87 10.144
## 10 observations deleted due to missingness
out3<-duncan.test(tabla3,"tratamiento")
out3
## $statistics
## MSerror Df Mean CV
## 10.14444 6 8.55 37.25186
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## NULL
##
## $means
## longraiz std r Min Max Q25 Q50 Q75
## TRT1 10.000000 NA 1 10 10.0 10.0 10 10.0
## TRT2 8.900000 3.612478 5 5 14.5 7.0 8 10.0
## TRT3 7.333333 2.081666 3 5 9.0 6.5 8 8.5
## TRT4 9.000000 NA 1 9 9.0 9.0 9 9.0
##
## $comparison
## NULL
##
## $groups
## longraiz groups
## TRT1 10.000000 a
## TRT4 9.000000 a
## TRT2 8.900000 a
## TRT3 7.333333 a
##
## attr(,"class")
## [1] "group"
ttest3<- HSD.test(tabla3,"tratamiento")
ttest3
## $statistics
## MSerror Df Mean CV
## 10.14444 6 8.55 37.25186
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 4.895599 0.05
##
## $means
## longraiz std r Min Max Q25 Q50 Q75
## TRT1 10.000000 NA 1 10 10.0 10.0 10 10.0
## TRT2 8.900000 3.612478 5 5 14.5 7.0 8 10.0
## TRT3 7.333333 2.081666 3 5 9.0 6.5 8 8.5
## TRT4 9.000000 NA 1 9 9.0 9.0 9 9.0
##
## $comparison
## NULL
##
## $groups
## longraiz groups
## TRT1 10.000000 a
## TRT4 9.000000 a
## TRT2 8.900000 a
## TRT3 7.333333 a
##
## attr(,"class")
## [1] "group"
Analisis para Ruda
library(readxl)
df2 <- read_excel("C:/Users/ginna/Downloads/practicarepro.xlsx",
sheet = "ruda")
#View(practicarepro)
summary(df2)
## tratamiento nraices pesoraiz longraiz
## Length:20 Min. : 6.00 Min. :0.0160 Min. : 2.400
## Class :character 1st Qu.:15.75 1st Qu.:0.2225 1st Qu.: 5.375
## Mode :character Median :23.00 Median :0.4100 Median : 6.900
## Mean :24.19 Mean :0.8310 Mean : 7.075
## 3rd Qu.:30.00 3rd Qu.:0.5700 3rd Qu.: 9.100
## Max. :55.00 Max. :7.5000 Max. :12.100
## NA's :4 NA's :4 NA's :4
library(agricolae)
tabla4=aov(nraices~tratamiento, data=df2)
summary(tabla4)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 102.5 34.16 0.175 0.911
## Residuals 12 2337.9 194.83
## 4 observations deleted due to missingness
out4<-duncan.test(tabla4,"tratamiento")
out4
## $statistics
## MSerror Df Mean CV
## 194.8292 12 24.1875 57.708
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## NULL
##
## $means
## nraices std r Min Max Q25 Q50 Q75
## TRT1 24.75 7.088723 4 15 30 21.75 27.0 30
## TRT2 25.00 18.248288 3 13 46 14.50 16.0 31
## TRT3 20.00 12.192894 4 6 35 13.50 19.5 26
## TRT4 26.60 16.395121 5 13 55 19.00 23.0 23
##
## $comparison
## NULL
##
## $groups
## nraices groups
## TRT4 26.60 a
## TRT2 25.00 a
## TRT1 24.75 a
## TRT3 20.00 a
##
## attr(,"class")
## [1] "group"
ttest4<- HSD.test(tabla4,"tratamiento")
ttest4
## $statistics
## MSerror Df Mean CV
## 194.8292 12 24.1875 57.708
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 4.19866 0.05
##
## $means
## nraices std r Min Max Q25 Q50 Q75
## TRT1 24.75 7.088723 4 15 30 21.75 27.0 30
## TRT2 25.00 18.248288 3 13 46 14.50 16.0 31
## TRT3 20.00 12.192894 4 6 35 13.50 19.5 26
## TRT4 26.60 16.395121 5 13 55 19.00 23.0 23
##
## $comparison
## NULL
##
## $groups
## nraices groups
## TRT4 26.60 a
## TRT2 25.00 a
## TRT1 24.75 a
## TRT3 20.00 a
##
## attr(,"class")
## [1] "group"
tabla5=aov(pesoraiz~tratamiento, data=df2)
summary(tabla5)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 15.99 5.329 1.996 0.168
## Residuals 12 32.04 2.670
## 4 observations deleted due to missingness
out5<-duncan.test(tabla5,"tratamiento")
out5
## $statistics
## MSerror Df Mean CV
## 2.669692 12 0.831 196.6209
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## NULL
##
## $means
## pesoraiz std r Min Max Q25 Q50 Q75
## TRT1 0.3075 0.1376893 4 0.190 0.50 0.220 0.27 0.3575
## TRT2 2.9100 3.9755880 3 0.550 7.50 0.615 0.68 4.0900
## TRT3 0.3240 0.2844597 4 0.016 0.68 0.154 0.30 0.4700
## TRT4 0.4080 0.1775387 5 0.170 0.63 0.310 0.42 0.5100
##
## $comparison
## NULL
##
## $groups
## pesoraiz groups
## TRT2 2.9100 a
## TRT4 0.4080 a
## TRT3 0.3240 a
## TRT1 0.3075 a
##
## attr(,"class")
## [1] "group"
ttest5<- HSD.test(tabla5,"tratamiento")
ttest5
## $statistics
## MSerror Df Mean CV
## 2.669692 12 0.831 196.6209
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 4.19866 0.05
##
## $means
## pesoraiz std r Min Max Q25 Q50 Q75
## TRT1 0.3075 0.1376893 4 0.190 0.50 0.220 0.27 0.3575
## TRT2 2.9100 3.9755880 3 0.550 7.50 0.615 0.68 4.0900
## TRT3 0.3240 0.2844597 4 0.016 0.68 0.154 0.30 0.4700
## TRT4 0.4080 0.1775387 5 0.170 0.63 0.310 0.42 0.5100
##
## $comparison
## NULL
##
## $groups
## pesoraiz groups
## TRT2 2.9100 a
## TRT4 0.4080 a
## TRT3 0.3240 a
## TRT1 0.3075 a
##
## attr(,"class")
## [1] "group"
tabla6=aov(longraiz~tratamiento, data=df2)
summary(tabla6)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 17.81 5.938 0.811 0.512
## Residuals 12 87.86 7.321
## 4 observations deleted due to missingness
out6<-duncan.test(tabla6,"tratamiento")
out6
## $statistics
## MSerror Df Mean CV
## 7.321347 12 7.075 38.24451
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## NULL
##
## $means
## longraiz std r Min Max Q25 Q50 Q75
## TRT1 7.875000 2.657536 4 5.0 11.0 6.125 7.75 9.500
## TRT2 5.866667 1.401190 3 4.5 7.3 5.150 5.80 6.550
## TRT3 8.350000 3.454948 4 4.0 12.1 6.625 8.65 10.375
## TRT4 6.140000 2.594802 5 2.4 9.4 5.500 5.90 7.500
##
## $comparison
## NULL
##
## $groups
## longraiz groups
## TRT3 8.350000 a
## TRT1 7.875000 a
## TRT4 6.140000 a
## TRT2 5.866667 a
##
## attr(,"class")
## [1] "group"
ttest6<- HSD.test(tabla6,"tratamiento")
ttest6
## $statistics
## MSerror Df Mean CV
## 7.321347 12 7.075 38.24451
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 4.19866 0.05
##
## $means
## longraiz std r Min Max Q25 Q50 Q75
## TRT1 7.875000 2.657536 4 5.0 11.0 6.125 7.75 9.500
## TRT2 5.866667 1.401190 3 4.5 7.3 5.150 5.80 6.550
## TRT3 8.350000 3.454948 4 4.0 12.1 6.625 8.65 10.375
## TRT4 6.140000 2.594802 5 2.4 9.4 5.500 5.90 7.500
##
## $comparison
## NULL
##
## $groups
## longraiz groups
## TRT3 8.350000 a
## TRT1 7.875000 a
## TRT4 6.140000 a
## TRT2 5.866667 a
##
## attr(,"class")
## [1] "group"
ANALISIS PARA ESTACAS DE MANZANO
library(readxl)
df3 <- read_excel("C:/Users/ginna/Downloads/practicarepro.xlsx",
sheet = "frutal")
#View(practicarepro)
tabla7=aov(longitudTallo~tratamiento, data=df3)
summary(tabla7)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 53.59 17.865 27.22 0.00402 **
## Residuals 4 2.63 0.656
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
out7<-duncan.test(tabla7,"tratamiento")
out7
## $statistics
## MSerror Df Mean CV
## 0.65625 4 32.9375 2.459484
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## Table CriticalRange
## 2 3.926503 2.249187
## 3 4.012542 2.298472
## 4 4.033093 2.310244
##
## $means
## longitudTallo std r Min Max Q25 Q50 Q75
## t1 31.00 1.4142136 2 30 32.0 30.500 31.00 31.500
## t2 30.00 0.0000000 2 30 30.0 30.000 30.00 30.000
## t3 34.25 0.3535534 2 34 34.5 34.125 34.25 34.375
## t4 36.50 0.7071068 2 36 37.0 36.250 36.50 36.750
##
## $comparison
## NULL
##
## $groups
## longitudTallo groups
## t4 36.50 a
## t3 34.25 b
## t1 31.00 c
## t2 30.00 c
##
## attr(,"class")
## [1] "group"
ttest7<- HSD.test(tabla7,"tratamiento")
ttest7
## $statistics
## MSerror Df Mean CV MSD
## 0.65625 4 32.9375 2.459484 3.297769
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 5.757058 0.05
##
## $means
## longitudTallo std r Min Max Q25 Q50 Q75
## t1 31.00 1.4142136 2 30 32.0 30.500 31.00 31.500
## t2 30.00 0.0000000 2 30 30.0 30.000 30.00 30.000
## t3 34.25 0.3535534 2 34 34.5 34.125 34.25 34.375
## t4 36.50 0.7071068 2 36 37.0 36.250 36.50 36.750
##
## $comparison
## NULL
##
## $groups
## longitudTallo groups
## t4 36.50 a
## t3 34.25 ab
## t1 31.00 bc
## t2 30.00 c
##
## attr(,"class")
## [1] "group"
tabla8=aov(diametroTallo~tratamiento, data=df3)
summary(tabla8)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 0.465 0.1550 8.857 0.0307 *
## Residuals 4 0.070 0.0175
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
out8<-duncan.test(tabla8,"tratamiento")
out8
## $statistics
## MSerror Df Mean CV
## 0.0175 4 1.125 11.75889
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## Table CriticalRange
## 2 3.926503 0.3672908
## 3 4.012542 0.3753389
## 4 4.033093 0.3772613
##
## $means
## diametroTallo std r Min Max Q25 Q50 Q75
## t1 1.25 0.21213203 2 1.1 1.4 1.175 1.25 1.325
## t2 1.40 0.00000000 2 1.4 1.4 1.400 1.40 1.400
## t3 0.75 0.07071068 2 0.7 0.8 0.725 0.75 0.775
## t4 1.10 0.14142136 2 1.0 1.2 1.050 1.10 1.150
##
## $comparison
## NULL
##
## $groups
## diametroTallo groups
## t2 1.40 a
## t1 1.25 a
## t4 1.10 ab
## t3 0.75 b
##
## attr(,"class")
## [1] "group"
ttest8<- HSD.test(tabla8,"tratamiento")
ttest8
## $statistics
## MSerror Df Mean CV MSD
## 0.0175 4 1.125 11.75889 0.5385235
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 5.757058 0.05
##
## $means
## diametroTallo std r Min Max Q25 Q50 Q75
## t1 1.25 0.21213203 2 1.1 1.4 1.175 1.25 1.325
## t2 1.40 0.00000000 2 1.4 1.4 1.400 1.40 1.400
## t3 0.75 0.07071068 2 0.7 0.8 0.725 0.75 0.775
## t4 1.10 0.14142136 2 1.0 1.2 1.050 1.10 1.150
##
## $comparison
## NULL
##
## $groups
## diametroTallo groups
## t2 1.40 a
## t1 1.25 ab
## t4 1.10 ab
## t3 0.75 b
##
## attr(,"class")
## [1] "group"
tabla9=aov(numeroBrotes~tratamiento, data=df3)
summary(tabla9)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 3 24.38 8.125 5 0.077 .
## Residuals 4 6.50 1.625
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
out9<-duncan.test(tabla9,"tratamiento")
out9
## $statistics
## MSerror Df Mean CV
## 1.625 4 3.125 40.79216
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 4 0.05
##
## $duncan
## Table CriticalRange
## 2 3.926503 3.539302
## 3 4.012542 3.616857
## 4 4.033093 3.635381
##
## $means
## numeroBrotes std r Min Max Q25 Q50 Q75
## t1 3.0 1.4142136 2 2 4 2.50 3.0 3.50
## t2 6.0 1.4142136 2 5 7 5.50 6.0 6.50
## t3 1.5 0.7071068 2 1 2 1.25 1.5 1.75
## t4 2.0 1.4142136 2 1 3 1.50 2.0 2.50
##
## $comparison
## NULL
##
## $groups
## numeroBrotes groups
## t2 6.0 a
## t1 3.0 ab
## t4 2.0 b
## t3 1.5 b
##
## attr(,"class")
## [1] "group"
ttest9<- HSD.test(tabla9,"tratamiento")
ttest9
## $statistics
## MSerror Df Mean CV MSD
## 1.625 4 3.125 40.79216 5.189342
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 4 5.757058 0.05
##
## $means
## numeroBrotes std r Min Max Q25 Q50 Q75
## t1 3.0 1.4142136 2 2 4 2.50 3.0 3.50
## t2 6.0 1.4142136 2 5 7 5.50 6.0 6.50
## t3 1.5 0.7071068 2 1 2 1.25 1.5 1.75
## t4 2.0 1.4142136 2 1 3 1.50 2.0 2.50
##
## $comparison
## NULL
##
## $groups
## numeroBrotes groups
## t2 6.0 a
## t1 3.0 a
## t4 2.0 a
## t3 1.5 a
##
## attr(,"class")
## [1] "group"
tabla10=aov(pesoFrescoBrotes~tratamiento, data=df3)
summary(tabla10)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 2 0.8011 0.4006 1.167 0.461
## Residuals 2 0.6864 0.3432
## 3 observations deleted due to missingness
out10<-duncan.test(tabla10,"tratamiento")
out10
## $statistics
## MSerror Df Mean CV
## 0.343225 2 0.89 65.8263
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 3 0.05
##
## $duncan
## NULL
##
## $means
## pesoFrescoBrotes std r Min Max Q25 Q50 Q75
## t1 1.205 0.7424621 2 0.68 1.73 0.9425 1.205 1.4675
## t2 0.960 0.3676955 2 0.70 1.22 0.8300 0.960 1.0900
## t3 0.120 NA 1 0.12 0.12 0.1200 0.120 0.1200
##
## $comparison
## NULL
##
## $groups
## pesoFrescoBrotes groups
## t1 1.205 a
## t2 0.960 a
## t3 0.120 a
##
## attr(,"class")
## [1] "group"
ttest10<- HSD.test(tabla10,"tratamiento")
ttest10
## $statistics
## MSerror Df Mean CV
## 0.343225 2 0.89 65.8263
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 3 8.330783 0.05
##
## $means
## pesoFrescoBrotes std r Min Max Q25 Q50 Q75
## t1 1.205 0.7424621 2 0.68 1.73 0.9425 1.205 1.4675
## t2 0.960 0.3676955 2 0.70 1.22 0.8300 0.960 1.0900
## t3 0.120 NA 1 0.12 0.12 0.1200 0.120 0.1200
##
## $comparison
## NULL
##
## $groups
## pesoFrescoBrotes groups
## t1 1.205 a
## t2 0.960 a
## t3 0.120 a
##
## attr(,"class")
## [1] "group"
tabla11=aov(longitudBrotes~tratamiento, data=df3)
summary(tabla11)
## Df Sum Sq Mean Sq F value Pr(>F)
## tratamiento 2 0.5139 0.2569 0.273 0.778
## Residuals 3 2.8236 0.9412
## 2 observations deleted due to missingness
out11<-duncan.test(tabla11,"tratamiento")
out11
## $statistics
## MSerror Df Mean CV
## 0.9412167 3 1.125 86.23673
##
## $parameters
## test name.t ntr alpha
## Duncan tratamiento 3 0.05
##
## $duncan
## Table CriticalRange
## 2 4.500659 3.087492
## 3 4.515652 3.097778
##
## $means
## longitudBrotes std r Min Max Q25 Q50 Q75
## t1 0.735 0.0212132 2 0.72 0.75 0.7275 0.735 0.7425
## t2 1.440 1.3576450 2 0.48 2.40 0.9600 1.440 1.9200
## t3 1.200 0.9899495 2 0.50 1.90 0.8500 1.200 1.5500
##
## $comparison
## NULL
##
## $groups
## longitudBrotes groups
## t2 1.440 a
## t3 1.200 a
## t1 0.735 a
##
## attr(,"class")
## [1] "group"
ttest11<- HSD.test(tabla11,"tratamiento")
ttest11
## $statistics
## MSerror Df Mean CV MSD
## 0.9412167 3 1.125 86.23673 4.054082
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey tratamiento 3 5.909663 0.05
##
## $means
## longitudBrotes std r Min Max Q25 Q50 Q75
## t1 0.735 0.0212132 2 0.72 0.75 0.7275 0.735 0.7425
## t2 1.440 1.3576450 2 0.48 2.40 0.9600 1.440 1.9200
## t3 1.200 0.9899495 2 0.50 1.90 0.8500 1.200 1.5500
##
## $comparison
## NULL
##
## $groups
## longitudBrotes groups
## t2 1.440 a
## t3 1.200 a
## t1 0.735 a
##
## attr(,"class")
## [1] "group"