###Prueba t-student comparar 2 medias independientes
set.seed(123)
rto_criolla <- rnorm(60,2.8,0.25) #rendimiento
rto_past <- rnorm(60,2.72,0.27)
ge_criolla <- rnorm(60,1.08,0.05)  #gravedad especifica
ge_past <- rnorm(60,1.03,0.03)

df_var <- data.frame(rto_criolla,rto_past,ge_criolla,ge_past)
boxplot(rto_criolla,rto_past, names = c("criolla", "pastusa"))
boxplot(ge_criolla,ge_past, names = c("criolla", "pastusa"))

par(mfrow = c(1,2))
hist(rto_criolla)
hist(rto_past)

par(mfrow = c(1,2))
hist(ge_criolla)
hist(ge_past)

par(mfrow = c(1,2))
plot(rto_criolla,ge_criolla)
plot(rto_past,ge_past)

cor(ge_criolla,rto_criolla)
cor(ge_past,rto_past)

library(corrplot)
M = cor(df_var)
corrplot(M, method = 'number')

# Pruebas de hipotesis
# H0: mu_rtoc = mu_rtop, Ha: mu_rtoc != mu_rtop
# H0: mu_gec = mu_gep, Ha: mu_gec != mu_gep

t.test(df_var$rto_criolla,df_var$rto_past)
t.test(df_var$rto_criolla,df_var$rto_past,alternative = "t")
t.test(df_var$rto_criolla,df_var$rto_past,alternative = "l")
t.test(df_var$rto_criolla,df_var$rto_past,alternative = "g")

p1 <- t.test(df_var$rto_criolla,df_var$rto_past,alternative = "t");p1
p2 <- t.test(sort(df_var$rto_criolla),sort(df_var$rto_past),alternative = "t");p2

# Analisis de dependencia espacial
library(spdep)

n <- 60
distancia_paso <- 0.30
observaciones <- sort(df_var$rto_criolla)

coordenadas <- matrix(c(seq(0, by = distancia_paso, length.out = n), rep(0, n)), ncol = 2)

vecinos <- dnearneigh(coordenadas, 0, 0.31)
pesos <- nb2listw(vecinos, style = "W", zero.policy = TRUE)

resultado_moran <- moran.test(observaciones, pesos)
print(resultado_moran)

plot(sp.correlogram(vecinos, observaciones, order = 5, method = "I", style = "W"))

# Prueba T^2 de Hotelling para muestras independientes
library(Hotelling)

set.seed(123)
grupo1 <- data.frame(rto_criolla, ge_criolla)
grupo2 <- data.frame(rto_past, ge_past)

resultado_t2 <- hotelling.test(grupo1, grupo2)
print(resultado_t2)