#Punto 2.1

#a
round (((24/-8+7)/(9-6*(2))),1)
## [1] -1.3
#b
round ((2^3*-7+4-((1)/(3)+((1)/(2)))),1)
## [1] -52.8
#c
round (((sqrt(16)+5*4-3^-2)/(3*(4-8)+1)),1)
## [1] -2.2
#d
round ((1)/(2)*(((5)/(4)-2^-1)*((13)/(4))-sqrt(7)+8/(2-((1)/(5)))-3^2),1)
## [1] -2.4
#e
round ((2*sin((pi)/(3))+5*cos((pi)/(4))-((tan((pi)/(6))/(4)))),1)
## [1] 5.1
#f
round ((4*log2(5)-3*log(7,base=3)+((1)/(3))*log(8)^3),1)
## [1] 7
#Punto 2.3
GC = c(14.2, 8.0, 9.2, 12.1, 8.7, 11.6, 11.0, 12.5, 10.0, 9.0, 8.5, 13.1, 12.9, 8.4, 11.2, 9.8, 12.4, 11.0, 13.0, 8.6)
NP = c(11, 10, 6, 9, 10, 12, 14, 10, 8, 9, 11, 12, 11, 11, 15, 13, 12, 12, 8, 13)

#media
prom_GC = mean(GC)
print(prom_GC)
## [1] 10.76
prom_NP = mean(NP)
print(prom_NP)
## [1] 10.85
#moda 
moda_GC = mode(GC)
print(moda_GC)
## [1] "numeric"
moda_NP = mode(NP)
print(moda_NP)
## [1] "numeric"
#mediana
acom_GC = sort(GC)
mediana_GC = median(acom_GC)
print(mediana_GC)
## [1] 11
acom_NP = sort(NP)
mediana_NP = median(acom_NP)
print(mediana_NP)
## [1] 11
#desviación estandar
de_GC = sd(GC)
d_1 = round(de_GC, 2)
print(d_1)
## [1] 1.91
de_NP = sd(NP)
d_2 = round(de_NP, 2)
print(d_2)
## [1] 2.18
#varianza
varia_GC = var(GC)
v_1 = round(varia_GC, 2)
print(v_1)
## [1] 3.67
varia_NP = var(NP)
v_2 = round(varia_NP, 2)
print(v_2)
## [1] 4.77
#cuartiles
cuartiles_GC = quantile(GC)
print(cuartiles_GC, 2)
##   0%  25%  50%  75% 100% 
##  8.0  8.9 11.0 12.4 14.2
cuartiles_NP = quantile(NP)
print(cuartiles_NP, 2)
##   0%  25%  50%  75% 100% 
##  6.0  9.8 11.0 12.0 15.0
#gráficas 
ID = c(1:20)
barplot(GC, names.arg = ID, col = "red", xlab = "ID", ylab = "GC", main = "ID vs GC")
mtext("Gráfica", side = 3, line = 3, cex = 1.2, font = 2)

barplot(NP, names.arg = ID, col = "orange", xlab = "ID", ylab = "NP", main = "ID vs NP")
mtext("Gráfica", side = 3, line = 3, cex = 1.2, font = 2)

#Punto 2.4
PHP = c(14.0, 24.7, 16.4, 26.0, 25.7, 24.6)
TMI =c(8.8, 10.2, 8.0, 9.1, 8.2, 9.4)

#media
prom_PHP = mean(PHP)
print(prom_PHP)
## [1] 21.9
prom_TMI = mean(TMI)
print(prom_TMI)
## [1] 8.95
#moda 
moda_PHP = mode(PHP)
print(moda_PHP)
## [1] "numeric"
acom_TMI = sort(TMI)
mediana_TMI = median(acom_TMI)
print(mediana_TMI)
## [1] 8.95
#desviación estandar
de_PHP = sd(PHP)
de_1 = round(de_PHP, 2)
print(de_1)
## [1] 5.27
de_TMI = sd(TMI)
de_2 = round(de_TMI, 2)
print(de_2)
## [1] 0.81
#Gráficas
INDICA = c("Central", "Norte", "Sur", "Oriente", "Oeste", "Rural")
barplot(PHP, names.arg = INDICA, col = "Red", xlab = "Indicador", ylab = "PHP", main = "Zona vs PHP")
mtext("Gráfica", side = 3, line = 3, cex = 1.2, font = 2)

barplot(TMI, names.arg = INDICA, col = "Orange", xlab = "Indicador", ylab = "TMI", main = "Zona vs TMI")
mtext("Gráfica", side = 3, line = 3, cex = 1.2, font = 2)

#Punto 2.5
datos = matrix(c(8, 12, 14, 9, 5, 10), nrow = 2, byrow = FALSE)
rownames(datos) = c("Masculino", "Femenino")
colnames(datos) = c("Casado", "Soltero", "Otro")

tabla = as.data.frame(datos)

rownames(tabla) = c("Masculino", "Femenino")
colnames(tabla) = c("Casado", "Soltero", "Otro")

print(tabla)
##           Casado Soltero Otro
## Masculino      8      14    5
## Femenino      12       9   10
#Matris
datos = matrix(c(8, 12, 14, 9, 5, 10), nrow = 2, byrow = FALSE)

rownames(datos) <- c("Masculino", "Femenino")
colnames(datos) <- c("Casado", "Soltero", "Otro")

#Crear data frame
df <- as.data.frame(datos)

print(df)
##           Casado Soltero Otro
## Masculino      8      14    5
## Femenino      12       9   10
#Gráfico de barras apiladas
barplot(as.matrix(df), beside = TRUE, legend.text = FALSE, col = c("cyan", "brown"),
        main = "Estado Conyugal por Género (Barras Apiladas)", xlab = "Género", ylab = "Frecuencia")

#Gráfico de barras agrupadas
barplot(as.matrix(df), beside = FALSE, col = c("purple", "yellow"),
        main = "Estado Conyugal por Género (Barras Agrupadas)", xlab = "Género", ylab = "Frecuencia")