#packages
library("tidyverse")
#cargo el df
#los descriptivos (como en articulo original), por genero
df3 %>%
group_by(Género) %>%
summarise(n=n(), Porcentaje = `n`/89*100)
`summarise()` ungrouping output (override with `.groups` argument)
#los descriptivos (como en articulo original), por grupo de edad
df3 %>%
group_by(Grupo.de.edad) %>%
summarise(n=n(), Porcentaje = `n`/89*100)
`summarise()` ungrouping output (override with `.groups` argument)
#agrupo los datos por genero, saco promedios y sd
df3 %>%
group_by(Género) %>%
summarise(n=n(), Promedioananguloderecho = mean(derecho), Desvanguloderecho = sd(derecho), Promedioanguloizquierdo= mean(izquierdo), Desvanguloizquierdo= sd(izquierdo), promedioconcavidadderecha= mean(concavidad_derecha), desvconcavidadderecha= sd(concavidad_derecha), promedioconcavidadizquierda= mean(concavidad_izquierda), desvconcavidadizquierda= sd(concavidad_izquierda)) %>%
ungroup()
`summarise()` ungrouping output (override with `.groups` argument)
#agrupo los datos por grupo de edad, saco promedios y sd
df3 %>%
group_by(Grupo.de.edad) %>%
summarise(n=n(), Promedioananguloderecho = mean(derecho), Desvanguloderecho = sd(derecho), Promedioanguloizquierdo= mean(izquierdo), Desvanguloizquierdo= sd(izquierdo), promedioconcavidadderecha= mean(concavidad_derecha), desvconcavidadderecha= sd(concavidad_derecha), promedioconcavidadizquierda= mean(concavidad_izquierda), desvconcavidadizquierda= sd(concavidad_izquierda)) %>%
ungroup()
`summarise()` ungrouping output (override with `.groups` argument)
#grafica para angulo derecho

#grafica para angulo izquierdo

#grafica para concavidad derecha

#grafica para concavidad izquierda

#graficas de boxplot para angulo derecho

#graficas de boxplot para angulo izquierdo

#graficas de boxplot para concavidad derecha

#graficas de boxplot para concavidad izquierda

#existe diferencia entre angulo derecho y genero ???
t.test(df3$derecho~df3$Género)#existen
Welch Two Sample t-test
data: df3$derecho by df3$Género
t = 3.077, df = 87.97, p-value = 0.00279
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
1.62398 7.54842
sample estimates:
mean in group Femenino mean in group Masculino
125.062 120.475
#existe diferencia entre angulo izquierdo y genero ???
t.test(df3$izquierdo~df3$Género)#no existen
Welch Two Sample t-test
data: df3$izquierdo by df3$Género
t = 2.01, df = 87.46, p-value = 0.0475
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
0.034948 6.131585
sample estimates:
mean in group Femenino mean in group Masculino
124.957 121.874
#existe diferencia entre concavidad derecha y genero ???
t.test(df3$concavidad_derecha~df3$Género) #no existen
Welch Two Sample t-test
data: df3$concavidad_derecha by df3$Género
t = -1.786, df = 76.02, p-value = 0.0781
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-5.539904 0.301771
sample estimates:
mean in group Femenino mean in group Masculino
14.0772 16.6963
#existe diferencia entre concavidad izquierda y genero ???
t.test(df3$concavidad_izquierda~df3$Género)# existen
Welch Two Sample t-test
data: df3$concavidad_izquierda by df3$Género
t = -2.829, df = 73.33, p-value = 0.00602
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-6.18350 -1.07219
sample estimates:
mean in group Femenino mean in group Masculino
12.9645 16.5924
LS0tCnRpdGxlOiAiMjAyMF90ZXNpc19nZW5lcm8geSBwYW5vcmFtaWNhcyIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQojcGFja2FnZXMKYGBge3J9CmxpYnJhcnkoInRpZHl2ZXJzZSIpCmBgYAojY2FyZ28gZWwgZGYKYGBge3J9CmRmMyA8LSByZWFkLmNzdigiaHR0cHM6Ly9kb2NzLmdvb2dsZS5jb20vc3ByZWFkc2hlZXRzL2QvZS8yUEFDWC0xdlNwUDd4bEI5YURLYmpPNWtkdVgxSE5ZSzl1UFlnbW5raGhoVXJVNTEycW5PN1dUeEJKQ3FYRmNBVTFZS25yXzQyaTdLRkFtbUtpQzNMay9wdWI/Z2lkPTE1Njk4NjE3NTgmc2luZ2xlPXRydWUmb3V0cHV0PWNzdiIpCmBgYAojbG9zIGRlc2NyaXB0aXZvcyAoY29tbyBlbiBhcnRpY3VsbyBvcmlnaW5hbCksIHBvciBnZW5lcm8gCmBgYHtyfQpvcHRpb25zKGRpZ2l0cyA9IDQpCmRmMyAlPiUgCiAgZ3JvdXBfYnkoR8OpbmVybykgJT4lIAogIHN1bW1hcmlzZShuPW4oKSwgUG9yY2VudGFqZSA9ICBgbmAvODkqMTAwKQpgYGAKI2xvcyBkZXNjcmlwdGl2b3MgKGNvbW8gZW4gYXJ0aWN1bG8gb3JpZ2luYWwpLCBwb3IgZ3J1cG8gZGUgZWRhZApgYGB7cn0Kb3B0aW9ucyhkaWdpdHMgPSA0KQpkZjMgJT4lIAogIGdyb3VwX2J5KEdydXBvLmRlLmVkYWQpICU+JSAKICBzdW1tYXJpc2Uobj1uKCksIFBvcmNlbnRhamUgPSAgYG5gLzg5KjEwMCkKYGBgCiNhZ3J1cG8gbG9zIGRhdG9zIHBvciBnZW5lcm8sIHNhY28gcHJvbWVkaW9zIHkgc2QKYGBge3J9Cm9wdGlvbnMoZGlnaXRzID0gNikKZGYzICU+JSAKICBncm91cF9ieShHw6luZXJvKSAlPiUgCiAgc3VtbWFyaXNlKG49bigpLCBQcm9tZWRpb2FuYW5ndWxvZGVyZWNobyA9IG1lYW4oZGVyZWNobyksIERlc3Zhbmd1bG9kZXJlY2hvID0gc2QoZGVyZWNobyksIFByb21lZGlvYW5ndWxvaXpxdWllcmRvPSBtZWFuKGl6cXVpZXJkbyksIERlc3Zhbmd1bG9penF1aWVyZG89IHNkKGl6cXVpZXJkbyksIHByb21lZGlvY29uY2F2aWRhZGRlcmVjaGE9IG1lYW4oY29uY2F2aWRhZF9kZXJlY2hhKSwgZGVzdmNvbmNhdmlkYWRkZXJlY2hhPSBzZChjb25jYXZpZGFkX2RlcmVjaGEpLCBwcm9tZWRpb2NvbmNhdmlkYWRpenF1aWVyZGE9IG1lYW4oY29uY2F2aWRhZF9penF1aWVyZGEpLCBkZXN2Y29uY2F2aWRhZGl6cXVpZXJkYT0gc2QoY29uY2F2aWRhZF9penF1aWVyZGEpKSAlPiUgCiAgdW5ncm91cCgpCmBgYAojYWdydXBvIGxvcyBkYXRvcyBwb3IgZ3J1cG8gZGUgZWRhZCwgc2FjbyBwcm9tZWRpb3MgeSBzZApgYGB7cn0Kb3B0aW9ucyhkaWdpdHMgPSA2KQpkZjMgJT4lIAogIGdyb3VwX2J5KEdydXBvLmRlLmVkYWQpICU+JSAKICBzdW1tYXJpc2Uobj1uKCksIFByb21lZGlvYW5hbmd1bG9kZXJlY2hvID0gbWVhbihkZXJlY2hvKSwgRGVzdmFuZ3Vsb2RlcmVjaG8gPSBzZChkZXJlY2hvKSwgUHJvbWVkaW9hbmd1bG9penF1aWVyZG89IG1lYW4oaXpxdWllcmRvKSwgRGVzdmFuZ3Vsb2l6cXVpZXJkbz0gc2QoaXpxdWllcmRvKSwgcHJvbWVkaW9jb25jYXZpZGFkZGVyZWNoYT0gbWVhbihjb25jYXZpZGFkX2RlcmVjaGEpLCBkZXN2Y29uY2F2aWRhZGRlcmVjaGE9IHNkKGNvbmNhdmlkYWRfZGVyZWNoYSksIHByb21lZGlvY29uY2F2aWRhZGl6cXVpZXJkYT0gbWVhbihjb25jYXZpZGFkX2l6cXVpZXJkYSksIGRlc3Zjb25jYXZpZGFkaXpxdWllcmRhPSBzZChjb25jYXZpZGFkX2l6cXVpZXJkYSkpICU+JSAKICB1bmdyb3VwKCkKYGBgCiNncmFmaWNhIHBhcmEgYW5ndWxvIGRlcmVjaG8gCmBgYHtyfQpkZjMgJT4lIAogIGdncGxvdChhZXMoeD1FZGFkLCB5PWRlcmVjaG8sIGNvbG9yPUfDqW5lcm8pKSArCiAgZ2VvbV9wb2ludCgpKwogIGdlb21fc21vb3RoKCkrCiAgdGhlbWVfbWluaW1hbCgpCmBgYAojZ3JhZmljYSBwYXJhIGFuZ3VsbyBpenF1aWVyZG8KYGBge3J9CmRmMyAlPiUgCiAgZ2dwbG90KGFlcyh4PUVkYWQsIHk9aXpxdWllcmRvLCBjb2xvcj1Hw6luZXJvKSkgKwogIGdlb21fcG9pbnQoKSsKICBnZW9tX3Ntb290aCgpKwogIHRoZW1lX21pbmltYWwoKQpgYGAKI2dyYWZpY2EgcGFyYSBjb25jYXZpZGFkIGRlcmVjaGEKYGBge3J9CmRmMyAlPiUgCiAgZ2dwbG90KGFlcyh4PUVkYWQsIHk9Y29uY2F2aWRhZF9kZXJlY2hhLCBjb2xvcj1Hw6luZXJvKSkgKwogIGdlb21fcG9pbnQoKSsKICBnZW9tX3Ntb290aCgpKwogIHRoZW1lX21pbmltYWwoKQpgYGAKI2dyYWZpY2EgcGFyYSBjb25jYXZpZGFkIGl6cXVpZXJkYQpgYGB7cn0KZGYzICU+JSAKICBnZ3Bsb3QoYWVzKHg9RWRhZCwgeT1jb25jYXZpZGFkX2l6cXVpZXJkYSwgY29sb3I9R8OpbmVybykpICsKICBnZW9tX3BvaW50KCkrCiAgZ2VvbV9zbW9vdGgoKSsKICB0aGVtZV9taW5pbWFsKCkKYGBgCiNncmFmaWNhcyBkZSBib3hwbG90IHBhcmEgYW5ndWxvIGRlcmVjaG8gCmBgYHtyfQpkZjMgJT4lIAogIGdncGxvdChhZXMoeD1Hw6luZXJvLCB5PWRlcmVjaG8pKSArCiAgZ2VvbV9ib3hwbG90KCkrCiAgdGhlbWVfbWluaW1hbCgpCmBgYAojZ3JhZmljYXMgZGUgYm94cGxvdCBwYXJhIGFuZ3VsbyBpenF1aWVyZG8KYGBge3J9CmRmMyAlPiUgCiAgZ2dwbG90KGFlcyh4PUfDqW5lcm8sIHk9aXpxdWllcmRvKSkgKwogIGdlb21fYm94cGxvdCgpKwogIHRoZW1lX21pbmltYWwoKQpgYGAKI2dyYWZpY2FzIGRlIGJveHBsb3QgcGFyYSBjb25jYXZpZGFkIGRlcmVjaGEKYGBge3J9CmRmMyAlPiUgCiAgZ2dwbG90KGFlcyh4PUfDqW5lcm8sIHk9Y29uY2F2aWRhZF9kZXJlY2hhKSkgKwogIGdlb21fYm94cGxvdCgpKwogIHRoZW1lX21pbmltYWwoKQpgYGAKI2dyYWZpY2FzIGRlIGJveHBsb3QgcGFyYSBjb25jYXZpZGFkIGl6cXVpZXJkYQpgYGB7cn0KZGYzICU+JSAKICBnZ3Bsb3QoYWVzKHg9R8OpbmVybywgeT1jb25jYXZpZGFkX2l6cXVpZXJkYSkpICsKICBnZW9tX2JveHBsb3QoKSsKICB0aGVtZV9taW5pbWFsKCkKYGBgCiNleGlzdGUgZGlmZXJlbmNpYSBlbnRyZSBhbmd1bG8gZGVyZWNobyB5IGdlbmVybyA/Pz8gCmBgYHtyfQp0LnRlc3QoZGYzJGRlcmVjaG9+ZGYzJEfDqW5lcm8pI2V4aXN0ZW4gCmBgYAojZXhpc3RlIGRpZmVyZW5jaWEgZW50cmUgYW5ndWxvIGl6cXVpZXJkbyB5IGdlbmVybyA/Pz8gCmBgYHtyfQp0LnRlc3QoZGYzJGl6cXVpZXJkb35kZjMkR8OpbmVybykjbm8gZXhpc3RlbiAKYGBgCiNleGlzdGUgZGlmZXJlbmNpYSBlbnRyZSBjb25jYXZpZGFkIGRlcmVjaGEgeSBnZW5lcm8gPz8/IApgYGB7cn0KdC50ZXN0KGRmMyRjb25jYXZpZGFkX2RlcmVjaGF+ZGYzJEfDqW5lcm8pICNubyBleGlzdGVuIApgYGAKI2V4aXN0ZSBkaWZlcmVuY2lhIGVudHJlIGNvbmNhdmlkYWQgaXpxdWllcmRhIHkgZ2VuZXJvID8/PyAKYGBge3J9CnQudGVzdChkZjMkY29uY2F2aWRhZF9penF1aWVyZGF+ZGYzJEfDqW5lcm8pIyBleGlzdGVuIApgYGAKCgoKCgoKCgoKCgoKCg==