inferencia_pesos_dos_grupos

mu1= 145.2; sd1= 2.8
mu2= 138.6; sd2= 3.2

hist(n1,prob=T, 
     xlim=c( range(group2)[1],range(group1)[2]), ylim = c( 0,0.17),
     col='white', border = 'white', main = "Distribucion del peso en las dos poblaciones"); 
lines(density(group1), col="cyan", lwd=3); 

lines(density(group2), col="brown3", lwd=3); 
abline(v=c(mean(n1),mean(n2)), col=c("cyan", "brown3"))

El promedio de peso de los hombres: 144.9 +/- 0.3. El promedio de peso de las mujeres: 138.4 +/- 0.3

####################
set.seed(2241)
n1 <- rnorm(100, mean = mu1, sd= sd1)
n2 <- rnorm(100, mean = mu2, sd= sd2)
#
pesos.n <- data.frame(cbind(peso = c(n1,n2),
                            sex = c(rep("b.men", 100), rep("a.women", 100) ) ) )
pesos.n$sex <- as.factor(pesos.n$sex) ; pesos.n$peso <- as.numeric(pesos.n$peso)
#str(pesos.n)

men.q.975 <- mean(n1) + (1.96*sd(n1)); men.q.025 <- mean(n1) - (1.96*sd(n1)) 
women.q.975 <- mean(n2) + (1.96*sd(n2)); women.q.025 <- mean(n2) - (1.96*sd(n2))

##############
library(ggplot2)
plotly::ggplotly(ggplot(data=pesos.N, aes(x= sex, color= sex)) + 
                   geom_point(aes(y= peso), position = position_jitter(width=0.5),
                              shape = 21, size = 2
                              ,
                              fill = "NA", stroke=0.1
                              )
          #       + geom_hline(yintercept = mu1, lty=2, col= "red3" ) 
          #       + geom_hline(yintercept = mu2, lwd=0.8, lty=2, col= "cyan2" ) 
                # +
                 + geom_point(data= pesos.n, col= "black",
                              aes(y= peso), position = position_jitter(width=0.2),
                              shape = 19, size = 1,
                              fill = "NA", stroke=0.1
                              )
                 +  stat_summary(data= pesos.n,
                               aes(y = peso, x= sex), geom = "pointrange", 
                               col= "yellow2",
                                fun.data = "mean_cl_normal")

            #     + geom_hline(yintercept = men.q.025, col= "blue", lty=2 ) 
            #     + geom_hline(yintercept = women.q.975,
            #                  col= "blue", lty=2 ) 

  )