Ademas es muy importanyte resaltar que, en todos los casos, los
sesgos son relativamente pequeños en valor absoluto. Esto sugiere que
las estimaciones tienden a estar bastante cerca de la proporción real,
lo que es una indicación positiva de la validez del Teorema del Límite
Central. Sin embargo, los sesgos varían en dirección según el porcentaje
real de plantas enfermas en la población.
tamanos_muestra <- c(5, 10, 15, 20, 30, 50, 60, 100, 200, 500)
print(paste("Porcentaje de plantas enfermas: 10% \n"))
## [1] "Porcentaje de plantas enfermas: 10% \n"
for (n_muestra in tamanos_muestra) {
resultados_estimador1 <- replicate(n_simulaciones, muestra_estimador(poblacion1, n_muestra))
# Prueba de Shapiro-Wilk
shapiro_test <- shapiro.test(resultados_estimador1)
print(shapiro_test)
# Gráfico de normalidad
qqnorm(resultados_estimador1, main = paste("QQ-Plot, n =", n_muestra))
qqline(resultados_estimador1)
# Comentar resultados
print(paste("Tamaño de muestra:", n_muestra))
print(paste("Shapiro-Wilk p-value:", shapiro_test$p.value))
}
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.67709, p-value < 2.2e-16

## [1] "Tamaño de muestra: 5"
## [1] "Shapiro-Wilk p-value: 9.33521571718581e-30"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.84053, p-value < 2.2e-16

## [1] "Tamaño de muestra: 10"
## [1] "Shapiro-Wilk p-value: 5.19252371656721e-22"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.89264, p-value < 2.2e-16

## [1] "Tamaño de muestra: 15"
## [1] "Shapiro-Wilk p-value: 3.29458069470375e-18"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.91935, p-value = 1.06e-15

## [1] "Tamaño de muestra: 20"
## [1] "Shapiro-Wilk p-value: 1.05987147471932e-15"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.94875, p-value = 3.845e-12

## [1] "Tamaño de muestra: 30"
## [1] "Shapiro-Wilk p-value: 3.84458935706495e-12"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.96532, p-value = 1.757e-09

## [1] "Tamaño de muestra: 50"
## [1] "Shapiro-Wilk p-value: 1.7566701412616e-09"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.97437, p-value = 1.12e-07

## [1] "Tamaño de muestra: 60"
## [1] "Shapiro-Wilk p-value: 1.11996897807292e-07"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.98213, p-value = 8.088e-06

## [1] "Tamaño de muestra: 100"
## [1] "Shapiro-Wilk p-value: 8.08818320440799e-06"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.99158, p-value = 0.006184

## [1] "Tamaño de muestra: 200"
## [1] "Shapiro-Wilk p-value: 0.00618393390659231"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador1
## W = 0.99178, p-value = 0.007249

## [1] "Tamaño de muestra: 500"
## [1] "Shapiro-Wilk p-value: 0.00724919889773392"
print(paste("Porcentaje de plantas enfermas: 90% \n"))
## [1] "Porcentaje de plantas enfermas: 90% \n"
for (n_muestra in tamanos_muestra) {
resultados_estimador2 <- replicate(n_simulaciones, muestra_estimador(poblacion2, n_muestra))
# Prueba de Shapiro-Wilk
shapiro_test <- shapiro.test(resultados_estimador2)
print(shapiro_test)
# Gráfico de normalidad
qqnorm(resultados_estimador2, main = paste("QQ-Plot, n =", n_muestra))
qqline(resultados_estimador2)
# Comentar resultados
print(paste("Tamaño de muestra:", n_muestra))
print(paste("Shapiro-Wilk p-value:", shapiro_test$p.value))
}
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.71518, p-value < 2.2e-16

## [1] "Tamaño de muestra: 5"
## [1] "Shapiro-Wilk p-value: 2.74710109939113e-28"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.84354, p-value < 2.2e-16

## [1] "Tamaño de muestra: 10"
## [1] "Shapiro-Wilk p-value: 8.08398554433445e-22"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.88333, p-value < 2.2e-16

## [1] "Tamaño de muestra: 15"
## [1] "Shapiro-Wilk p-value: 5.63563487015288e-19"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.92573, p-value = 5.124e-15

## [1] "Tamaño de muestra: 20"
## [1] "Shapiro-Wilk p-value: 5.12426447231444e-15"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.9576, p-value = 8.367e-11

## [1] "Tamaño de muestra: 30"
## [1] "Shapiro-Wilk p-value: 8.36707513619835e-11"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.97377, p-value = 8.301e-08

## [1] "Tamaño de muestra: 50"
## [1] "Shapiro-Wilk p-value: 8.30131986767019e-08"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.97971, p-value = 1.951e-06

## [1] "Tamaño de muestra: 60"
## [1] "Shapiro-Wilk p-value: 1.95076256948616e-06"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.98543, p-value = 6.655e-05

## [1] "Tamaño de muestra: 100"
## [1] "Shapiro-Wilk p-value: 6.65540609182676e-05"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.99005, p-value = 0.001846

## [1] "Tamaño de muestra: 200"
## [1] "Shapiro-Wilk p-value: 0.00184621253529977"
##
## Shapiro-Wilk normality test
##
## data: resultados_estimador2
## W = 0.99435, p-value = 0.0611

## [1] "Tamaño de muestra: 500"
## [1] "Shapiro-Wilk p-value: 0.0611002537570506"