# leer datos
D = read.csv("DatosEquipo7_filtrado.csv")
D$p34_2_bin = ifelse(D$p34_2 %in% c(2,3), "Sí", "No")
D$p34_3_bin = ifelse(D$p34_3 %in% c(2,3), "Sí", "No")
grupo_2 = subset(D, p34_2_bin == "Sí") # vio leer a sus padres
grupo_3 = subset(D, p34_3_bin == "Sí") # le leían de niño
# Nivel de confianza
conf.level = 0.96
# Intervalos de confianza
# p4 en p34_2 = Sí
ic_p4_g2 = t.test(grupo_2$p4, conf.level = conf.level)
# p4 en grupo p34_3 = Sí
ic_p4_g3 = t.test(grupo_3$p4, conf.level = conf.level)
# p26 en p34_2 = Sí
ic_p26_g2 = t.test(grupo_2$p26, conf.level = conf.level)
# p26 en p34_3 = Sí
ic_p26_g3 = t.test(grupo_3$p26, conf.level = conf.level)
# Intervalos de confianza diferencia de medias
# Diferencia de medias p4
ic_diff_p4 = t.test(grupo_2$p4, grupo_3$p4, conf.level = conf.level)
# Diferencia de medias p26
ic_diff_p26 = t.test(grupo_2$p26, grupo_3$p26, conf.level = conf.level)
#resultados
cat("Media de p4 en grupo p34_2 = Sí: \n")
## Media de p4 en grupo p34_2 = Sí:
cat("Media =", ic_p4_g2$estimate, "\n Intervalos de confianza96% = [", ic_p4_g2$conf.int[1], ",", ic_p4_g2$conf.int[2], "]\n")
## Media = 0.8477157
## Intervalos de confianza96% = [ 0.6739368 , 1.021495 ]
cat("Media de p4 en grupo p34_3 = Sí: \n")
## Media de p4 en grupo p34_3 = Sí:
cat("Media =", ic_p4_g3$estimate, "\n Intervalos de confianza96% = [", ic_p4_g3$conf.int[1], ",", ic_p4_g3$conf.int[2], "]\n")
## Media = 0.9708365
## Intervalos de confianza96% = [ 0.8194057 , 1.122267 ]
cat("Media de p26 en grupo p34_2 = Sí:\n")
## Media de p26 en grupo p34_2 = Sí:
cat("Media =", ic_p26_g2$estimate, "\n Intervalos de confianza96% = [", ic_p26_g2$conf.int[1], ",", ic_p26_g2$conf.int[2], "]\n")
## Media = 20.2132
## Intervalos de confianza96% = [ 18.27073 , 22.15566 ]
cat("Media de p26 en grupo p34_3 = Sí: \n")
## Media de p26 en grupo p34_3 = Sí:
cat("Media =", ic_p26_g3$estimate, "\n Intervalos de confianza96% = [", ic_p26_g3$conf.int[1], ",", ic_p26_g3$conf.int[2], "]\n")
## Media = 22.46355
## Intervalos de confianza96% = [ 20.74663 , 24.18047 ]
cat("Diferencia de medias de p4 entre grupos p34_2 = Sí y p34_3 = Sí:\n")
## Diferencia de medias de p4 entre grupos p34_2 = Sí y p34_3 = Sí:
cat("Diferencia =", ic_diff_p4$estimate[1] - ic_diff_p4$estimate[2], "\n Intervalos de confianza96% = [", ic_diff_p4$conf.int[1], ",", ic_diff_p4$conf.int[2], "]\n")
## Diferencia = -0.1231208
## Intervalos de confianza96% = [ -0.3534894 , 0.1072479 ]
cat("Diferencia de medias de p26 entre grupos p34_2 = Sí y p34_3 = Sí:\n")
## Diferencia de medias de p26 entre grupos p34_2 = Sí y p34_3 = Sí:
cat("Diferencia =", ic_diff_p26$estimate[1] - ic_diff_p26$estimate[2],"\n Intervalos de confianza96% = [", ic_diff_p26$conf.int[1], ",", ic_diff_p26$conf.int[2], "]\n")
## Diferencia = -2.250348
## Intervalos de confianza96% = [ -4.841348 , 0.3406526 ]
medias = c(ic_p4_g2$estimate, ic_p4_g3$estimate, ic_p26_g2$estimate, ic_p26_g3$estimate, ic_diff_p4$estimate[1] - ic_diff_p4$estimate[2], ic_diff_p26$estimate[1] - ic_diff_p26$estimate[2])
lim_inf = c(ic_p4_g2$conf.int[1], ic_p4_g3$conf.int[1], ic_p26_g2$conf.int[1], ic_p26_g3$conf.int[1], ic_diff_p4$conf.int[1], ic_diff_p26$conf.int[1])
lim_sup = c(ic_p4_g2$conf.int[2], ic_p4_g3$conf.int[2], ic_p26_g2$conf.int[2], ic_p26_g3$conf.int[2], ic_diff_p4$conf.int[2], ic_diff_p26$conf.int[2])
nombres = c("p4 (p34_2=Sí)", "p4 (p34_3=Sí)", "p26 (p34_2=Sí)", "p26 (p34_3=Sí)", "Dif p4", "Dif p26")
par(mfrow = c(3, 2)) # 3
#p4 (p34_2=Sí)
plot(0, xlim = c(ic_p4_g2$conf.int[1], ic_p4_g2$conf.int[2]), ylim = c(0, 1), yaxt = "n", ylab = "", xlab = "Libros leídos (p4)", main = "Intervalos de Confianza96%: p4 (p34_2 = Sí)")
arrows(ic_p4_g2$conf.int[1], 0.5, ic_p4_g2$conf.int[2], 0.5, angle = 90, code = 3, length = 0.1, lwd = 2)
points(ic_p4_g2$estimate, 0.5, pch = 19, cex = 1.2, col = "blue")
#p4 (p34_3=Sí)
plot(0, xlim = c(ic_p4_g3$conf.int[1], ic_p4_g3$conf.int[2]), ylim = c(0, 1), yaxt = "n", ylab = "", xlab = "Libros leídos (p4)", main = "Intervalos de Confianza96%: p4 (p34_3 = Sí)")
arrows(ic_p4_g3$conf.int[1], 0.5, ic_p4_g3$conf.int[2], 0.5, angle = 90, code = 3, length = 0.1, lwd = 2)
points(ic_p4_g3$estimate, 0.5, pch = 19, cex = 1.2, col = "blue")
#p26 (p34_2=Sí)
plot(0, xlim = c(ic_p26_g2$conf.int[1], ic_p26_g2$conf.int[2]), ylim = c(0, 1), yaxt = "n", ylab = "", xlab = "Minutos de lectura (p26)", main = "Intervalos de Confianza96%: p26 (p34_2 = Sí)")
arrows(ic_p26_g2$conf.int[1], 0.5, ic_p26_g2$conf.int[2], 0.5, angle = 90, code = 3, length = 0.1, lwd = 2)
points(ic_p26_g2$estimate, 0.5, pch = 19, cex = 1.2, col = "green")
#p26 (p34_3=Sí)
plot(0, xlim = c(ic_p26_g3$conf.int[1], ic_p26_g3$conf.int[2]), ylim = c(0, 1), yaxt = "n", ylab = "", xlab = "Minutos de lectura (p26)", main = "Intervalos de Confianza96%: p26 (p34_3 = Sí)")
arrows(ic_p26_g3$conf.int[1], 0.5, ic_p26_g3$conf.int[2], 0.5, angle = 90, code = 3, length = 0.1, lwd = 2)
points(ic_p26_g3$estimate, 0.5, pch = 19, cex = 1.2, col = "green")
#Diferencia p4
plot(0, xlim = c(ic_diff_p4$conf.int[1], ic_diff_p4$conf.int[2]), ylim = c(0, 1), yaxt = "n", ylab = "", xlab = "Diferencia en libros (p4)", main = "Intervalos de Confianza96%: Diferencia p4 (p34_2 vs p34_3)")
arrows(ic_diff_p4$conf.int[1], 0.5, ic_diff_p4$conf.int[2], 0.5, angle = 90, code = 3, length = 0.1, lwd = 2)
points(ic_diff_p4$estimate[1] - ic_diff_p4$estimate[2], 0.5, pch = 19, cex = 1.2, col = "red")
#Diferencia p26
plot(0, xlim = c(ic_diff_p26$conf.int[1], ic_diff_p26$conf.int[2]), ylim = c(0, 1), yaxt = "n", ylab = "", xlab = "Diferencia en minutos (p26)", main = "Intervalos de Confianza96%: Diferencia p26 (p34_2 vs p34_3)")
arrows(ic_diff_p26$conf.int[1], 0.5, ic_diff_p26$conf.int[2], 0.5, angle = 90, code = 3, length = 0.1, lwd = 2)
points(ic_diff_p26$estimate[1] - ic_diff_p26$estimate[2], 0.5, pch = 19, cex = 1.2, col = "red")
