library(readxl)
## Warning: package 'readxl' was built under R version 4.3.3
BP <- read_excel("C:/Users/Asus/Documents/Maestría/2024-1/BP.xlsx")
BP
## # A tibble: 1,380 × 5
## Muestreo Finca Variedad Tama Enraizamiento
## <dbl> <chr> <chr> <dbl> <dbl>
## 1 1 Arenales Golden Amazon 10.8 NA
## 2 1 Arenales Golden Amazon 12.6 NA
## 3 1 Arenales Golden Amazon 8.3 NA
## 4 1 Arenales Golden Amazon 8.9 NA
## 5 1 Arenales Golden Amazon 9.2 NA
## 6 1 Arenales Golden Amazon 9.6 NA
## 7 1 Arenales Golden Amazon 8.6 NA
## 8 1 Arenales Golden Amazon 9.7 NA
## 9 1 Arenales Golden Amazon 8.8 NA
## 10 1 Arenales Golden Amazon 7.8 NA
## # ℹ 1,370 more rows
anova<-aov(Tama~Finca+Variedad, data = BP)
summary(anova)
## Df Sum Sq Mean Sq F value Pr(>F)
## Finca 1 110 110 56.16 1.24e-13 ***
## Variedad 4 11094 2774 1415.40 < 2e-16 ***
## Residuals 1284 2516 2
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 90 observations deleted due to missingness
library(agricolae)
## Warning: package 'agricolae' was built under R version 4.3.3
tukey_dca <-HSD.test(anova, "Finca", group = TRUE)
tukey_dca
## $statistics
## MSerror Df Mean CV
## 1.959499 1284 6.341273 22.07477
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey Finca 2 2.774423 0.05
##
## $means
## Tama std r se Min Max Q25 Q50 Q75
## Arenales 6.028070 2.989106 600 0.05714746 0.3 13.4 3.5 4.5 9.2
## Monserrate 6.613623 3.462004 690 0.05329028 1.5 18.2 3.5 5.5 8.7
##
## $comparison
## NULL
##
## $groups
## Tama groups
## Monserrate 6.613623 a
## Arenales 6.028070 b
##
## attr(,"class")
## [1] "group"
tukey_dca <-HSD.test(anova, "Variedad", group = TRUE)
tukey_dca
## $statistics
## MSerror Df Mean CV
## 1.959499 1284 6.341273 22.07477
##
## $parameters
## test name.t ntr StudentizedRange alpha
## Tukey Variedad 5 3.863146 0.05
##
## $means
## Tama std r se Min Max Q25 Q50 Q75
## Briba Green 4.110741 1.441497 270 0.08519040 0.3 9.5 3.5 3.50 5.00000
## Corine Purple 4.077778 1.241876 270 0.08519040 1.5 9.5 3.5 3.50 4.50000
## Golden Amazon 11.143967 1.830724 240 0.09035806 7.4 16.2 9.7 10.77 12.35000
## Michigan Purple 4.299630 1.500768 270 0.08519040 1.5 9.5 3.5 4.00 5.50000
## Monte de Oro 8.891208 1.278755 240 0.09035806 6.2 18.2 8.1 8.70 9.58125
##
## $comparison
## NULL
##
## $groups
## Tama groups
## Golden Amazon 11.143967 a
## Monte de Oro 8.891208 b
## Michigan Purple 4.299630 c
## Briba Green 4.110741 c
## Corine Purple 4.077778 c
##
## attr(,"class")
## [1] "group"
library(dplyr)
## Warning: package 'dplyr' was built under R version 4.3.3
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
medias <- BP %>%
group_by(Finca, Variedad) %>%
summarise(Media_Tama = mean(Tama, na.rm = TRUE))
## `summarise()` has grouped output by 'Finca'. You can override using the
## `.groups` argument.
print(medias)
## # A tibble: 10 × 3
## # Groups: Finca [2]
## Finca Variedad Media_Tama
## <chr> <chr> <dbl>
## 1 Arenales Briba Green 3.67
## 2 Arenales Corine Purple 3.81
## 3 Arenales Golden Amazon 9.94
## 4 Arenales Michigan Purple 3.82
## 5 Arenales Monte de Oro 8.90
## 6 Monserrate Briba Green 4.47
## 7 Monserrate Corine Purple 4.29
## 8 Monserrate Golden Amazon 12.4
## 9 Monserrate Michigan Purple 4.68
## 10 Monserrate Monte de Oro 8.88
CLP<- c("c","c","a","c","b","c","c","a","c","b")
data<-data.frame(LA=as.numeric(BP$Tama), muestreo=as.factor(BP$Finca),Tratamiento=BP$Variedad)
dt <-data %>%
group_by(muestreo, Tratamiento) %>%
summarise(w = mean(LA,na.rm=TRUE), sd = sd(LA,na.rm=TRUE))
## `summarise()` has grouped output by 'muestreo'. You can override using the
## `.groups` argument.
dt$CLP<-CLP
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.3.3
ggplot(dt, aes(muestreo,w, fill = Tratamiento)) +
geom_bar(position="dodge", stat="identity") +
theme_classic() +
xlab("Finca") + ylab("Tamaño") + coord_cartesian(ylim=c(0,18)) +
scale_fill_manual(values=c("cadetblue","cadetblue1","cadetblue3","cyan4","aquamarine1"),name = "Variedad", labels = c("Briba Green","Corine Purple","Golden Amazon", "Michigan Purple","Monte de Oro")) +
geom_text(aes(label = CLP), position = position_dodge(width=1), vjust=-4)+
geom_errorbar(aes(ymin = w-sd, ymax = w+sd),width = .5, position = position_dodge(0.9))+
geom_text(data = NULL, x = 1, y = 25, label = "**") +
geom_text(data = NULL, x = 2, y = 25, label = "**") +
geom_text(data = NULL, x = 3, y = 25, label = "**") +
geom_text(data = NULL, x = 4, y = 25, label = "**")
