Sexo_A <- factor(c("Mujer", "Hombre", "Mujer", "Mujer", "Mujer", "Hombre", "Mujer",  
"Hombre", "Hombre", "Mujer", "Mujer" ,"Hombre", "Hombre", "Mujer", "Mujer", "Hombre", 
 "Mujer", "Mujer", "Mujer", "Mujer" ), levels = c("Mujer", "Hombre"))

Edad_A <- c(25, 30, 28, 20, 23, 22, 22, 22, 21, 21, 22, 20, 22, 29, 29, 21, 30, 21, 22, 23)

Estatura_A <- c(1.82, 1.83, 1.78, 1.79, 1.80, 1.90, 1.79,1.83, NA, 1.65, 1.73, 1.79, 1.80, 7, 
 1.69, 1.75, 1.66, NA, 1.79, 1.80)

Gr_Sang_A <- c("A", "B", "A", "AB", "0", "A", "B", "A", "B", "AB", "A", "B","0", "0", "A",
 "B", "AB", "B", "B", "B")

 GrA<- data.frame (Sexo_A,Edad_A, Estatura_A, Gr_Sang_A)
GrA

Sexo_B <- factor(c("Mujer", "Hombre", " Hombre ", "Mujer", " Hombre ", "Hombre", "Mujer",   "Hombre", "Hombre", "Mujer", "Mujer" ,"Hombre", "Hombre", "Mujer", " Hombre", "Hombre",     "Mujer", " Hombre ", "Mujer", " Hombre " ), levels = c("Mujer", "Hombre"))
Edad_B<- c(25, 30, 38, 26, 23, 22, 32, 26, 25, 28, 22, 30, 22, 39, 29, 25, 30, 23, 32, 23)
Estatura_B <- c(1.72, 1.73, 1.78, 1.79, 1.83, 1.90, 1.79,1.83, NA, 1.75, 1.79, 1.79, 1.85, 1.77, 
                1.79, 1.75, 1.66, NA, 1.79, 1.83)
Gr_Sang_B <- c("B", "B", "A", "AB", "0", "AB", "B", "B", "B", "AB", "A", "B","0", "0", "B",
               "B", "A", "0", "B", "B")
GrB<- data.frame (Sexo_B, Edad_B,  Estatura_B,  Gr_Sang_B)
GrB



pie(table(GrA$Gr_Sang_A), col = c( "blue", "yellow", "pink", "green"), main = "Diagrama de sectores
 para Grupo Sanguineo (GrA)")

pie(table(GrB$Gr_Sang_B), col = c( "blue", "yellow", "pink", "green"), main = "Diagrama de sectores para Grupo Sanguineo (GrB)")

hist(table(GrA$Estatura_A), col = "blue", 
     main = "Histograma para Estatura (GrA)", xlab="Estatura", ylab="Frecuencia")

hist(table(GrB$Estatura_B), col = "red", 
     main = "Histograma para Estatura (GrB)", xlab="Estatura", ylab="Frecuencia")


boxplot(GrA$Edad_A, xlab="Edad", main = "Cajas y bigotes de la Edad (GrA)")

boxplot(GrB$Edad_B, xlab="Edad", main = "Cajas y bigotes de la Edad (GrB)")

quantile(GrA$Estatura_A, na.rm = TRUE, probs = c(0.40))

quantile(GrB$Estatura_B, na.rm = TRUE, probs = c(0.70))

media_Est_A <- mean(GrA$Estatura_A, na.rm = TRUE)
var_Est_A <- 19/20 * var(GrA$Estatura_A, na.rm = TRUE)
dt_Est_A <- sqrt(var_Est_A)
 CV_Est_A <- dt_Est_A/media_Est_A
CV_Est_A

media_Est_B <- mean(GrB$Estatura_B, na.rm = TRUE)
var_Est_B <- 19/20 * var(GrB$Estatura_B, na.rm = TRUE)
dt_Est_B <- sqrt(var_Est_B)
CV_Est_B <- dt_Est_B/media_Est_B
CV_Est_B

mean(GrA$Estatura_A, na.rm = TRUE)
median(GrA$Estatura_A, na.rm = TRUE)

mean(GrB$Estatura_B, na.rm = TRUE)
median(GrB$Estatura_B, na.rm = TRUE)

install.packages("moments") # Es necesario si aśn no se ha instalado
library("moments")
skewness(GrA$Estatura_A, na.rm = TRUE)
kurtosis (GrA$Estatura_A, na.rm = TRUE)