Actividad 3

# Enunciado del estudio
cat("En un estudio geológico, se quiere investigar si tres tipos diferentes de rocas sedimentarias (arenisca, lutita y caliza) presentan diferencias significativas en su porosidad. La porosidad es una medida importante porque afecta la capacidad de las rocas para almacenar y transmitir agua y petróleo. Se seleccionaron aleatoriamente cinco muestras de cada tipo de roca de diferentes yacimientos geológicos para estudiar su porosidad.\n\n")

## En un estudio geológico, se quiere investigar si tres tipos diferentes de rocas sedimentarias (arenisca, lutita y caliza) presentan diferencias significativas en su porosidad. La porosidad es una medida importante porque afecta la capacidad de las rocas para almacenar y transmitir agua y petróleo. Se seleccionaron aleatoriamente cinco muestras de cada tipo de roca de diferentes yacimientos geológicos para estudiar su porosidad.

# Datos
porosidad <- data.frame(
  Tipo_Roca = rep(c("Arenisca", "Lutita", "Caliza"), each = 5),
  Porosidad = c(14, 15, 13, 16, 14, 6, 7, 8, 5, 6, 12, 11, 12, 10, 13)
)

# Análisis de varianza (ANOVA) para comparar los tres tipos de rocas
anova_result <- aov(Porosidad ~ Tipo_Roca, data = porosidad)
summary(anova_result)

##             Df Sum Sq Mean Sq F value   Pr(>F)    
## Tipo_Roca    2  164.8    82.4   63.38 4.18e-07 ***
## Residuals   12   15.6     1.3                     
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

# Test de Tukey para comparaciones múltiples
tukey_result <- TukeyHSD(anova_result)
print(tukey_result)

##   Tukey multiple comparisons of means
##     95% family-wise confidence level
## 
## Fit: aov(formula = Porosidad ~ Tipo_Roca, data = porosidad)
## 
## $Tipo_Roca
##                 diff       lwr        upr     p adj
## Caliza-Arenisca -2.8 -4.723824 -0.8761762 0.0057015
## Lutita-Arenisca -8.0 -9.923824 -6.0761762 0.0000003
## Lutita-Caliza   -5.2 -7.123824 -3.2761762 0.0000296

# Gráfico de barras para visualizar la porosidad promedio por tipo de roca
library(ggplot2)
ggplot(porosidad, aes(x = Tipo_Roca, y = Porosidad)) +
  geom_bar(stat = "summary", fun = "mean", fill = "skyblue") +
  labs(title = "Porosidad Promedio por Tipo de Roca",
       x = "Tipo de Roca",
       y = "Porosidad Promedio")

# Respuestas a las preguntas
cat("\nPreguntas:\n")

## 
## Preguntas:

cat("• ¿Hay diferencias significativas en la porosidad entre los tres tipos de rocas?\n")

## • ¿Hay diferencias significativas en la porosidad entre los tres tipos de rocas?

if (summary(anova_result)[[1]][[5]][[1]] < 0.05) {
  cat("   Sí, hay diferencias significativas en la porosidad entre los tres tipos de rocas.\n")
} else {
  cat("   No se encontraron diferencias significativas en la porosidad entre los tres tipos de rocas.\n")
}

##    Sí, hay diferencias significativas en la porosidad entre los tres tipos de rocas.

cat("• ¿Cuál de los tipos de rocas presenta mayor porosidad en promedio?\n")

## • ¿Cuál de los tipos de rocas presenta mayor porosidad en promedio?

means <- aggregate(Porosidad ~ Tipo_Roca, data = porosidad, FUN = mean)
max_mean <- max(means$Porosidad)
roca_max_mean <- means$Tipo_Roca[which.max(means$Porosidad)]
cat("   En promedio, la", roca_max_mean, "presenta la mayor porosidad.\n")

##    En promedio, la Arenisca presenta la mayor porosidad.