Asignacion : Construir un diagrama que permita visualizar el comportamiento de la floracion, la dosis y varieda. Usando la librerÃa ggpubr
library(readxl)
dff <- read_excel("floracion.xlsx", sheet = "FLORACION")
yf=cbind(floreada=dff$floreada,no.floreada=dff$total-dff$floreada)
head(yf)
## floreada no.floreada
## [1,] 0 12
## [2,] 0 17
## [3,] 4 6
## [4,] 9 2
## [5,] 10 0
## [6,] 0 17
pf<-dff$floreada/dff$total; pf #proporcion de floracion
## [1] 0.00000000 0.00000000 0.40000000 0.81818182 1.00000000 0.00000000
## [7] 0.20000000 0.50000000 0.90000000 0.50000000 0.14285714 0.06666667
## [13] 0.17647059 0.25000000 1.00000000 0.11111111 0.15789474 0.53571429
## [19] 0.73076923 0.77777778 0.00000000 0.00000000 0.15789474 0.75000000
## [25] 1.00000000 0.00000000 0.08333333 0.00000000 0.06666667 0.00000000
boxplot(pf)
hist(pf, breaks = 20)
pfc<-split(pf,dff$variedad);pfc
## $A
## [1] 0.0000000 0.0000000 0.4000000 0.8181818 1.0000000 0.0000000
##
## $B
## [1] 0.00000000 0.20000000 0.50000000 0.90000000 0.50000000 0.08333333
##
## $C
## [1] 0.14285714 0.06666667 0.17647059 0.25000000 1.00000000 0.00000000
##
## $D
## [1] 0.11111111 0.15789474 0.53571429 0.73076923 0.77777778 0.06666667
##
## $E
## [1] 0.0000000 0.0000000 0.1578947 0.7500000 1.0000000 0.0000000
class(pfc)
## [1] "list"
ball <- data.frame(round(pf, digits = 2), dff$dosis, dff$variedad)
colnames(ball) <- c("prop.flor", "dosis", "variedad")
library(ggpubr)
## Loading required package: ggplot2
my_cols <- c("#0D0887FF", "#6A00A8FF", "#B12A90FF",
"#E16462FF", "#FCA636FF", "#F0F921FF")
ggballoonplot(data = ball, y = "dosis", x = "variedad", size = "prop.flor",
fill = "prop.flor") + scale_fill_gradientn(colors = my_cols) +
guides(size = FALSE) + xlab("Dosis") + ylab("Variedad")
Dado que con la anterior librerÃa no se logró adicionar los nombres de los ejes X y Y, se utilizó la librerÃa ggplot2 que también permite la contrucción de gráficos de balón.
library(ggplot2)
hec_sp <- ggplot(ball, aes(x = dosis, y = variedad)) +
geom_point(aes(size = prop.flor), shape = 21, colour = "black", fill = "cornsilk") +
scale_size_area(max_size = 20, guide = FALSE) +
geom_text(aes(
y = as.numeric(as.factor(variedad)) - sqrt(prop.flor)/34, label = prop.flor),
vjust = 1.3,
colour = "black",
size = 2
)
hec_sp
hec_sp +
geom_point(aes(y = as.numeric(as.factor(variedad)) - sqrt(prop.flor)/34), colour = "red", size = 1)