##Ejercicio 1:
En un estudio para analizar la evolución de tubérculos almacenados, se deseaba comparar dos épocas de cosecha: Abril y Agosto, las que determinan diferentes periodos de almacenamiento. La variable en estudio fue la pérdida de peso por deshidratación (en gr). El objetivo es determinar si la época de cosecha de los tubérculos tiene efecto sobre la pérdida de peso promedio por deshidratación.
unidad Experimental: Tubérculos, n=20
Variable Respuesta: Perdida de peso por deshidratación en gramos.
Factor: Época de cosecha
Tratamientos: 2, se evaluaron los meses de Abril y Agosto
Réplicas: 10 tubérculos evaluados en cada época.
Construcción del diseño:
| Epoca | ||||||||||
|---|---|---|---|---|---|---|---|---|---|---|
| Agosto | 33.25 | 27.75 | 32.15 | 21.16 | 25.19 | 23.42 | 26.87 | 22.36 | 24.13 | 30.22 |
| Abril | 35.56 | 36.89 | 47.05 | 44.36 | 42.05 | 43.58 | 37.63 | 40.21 | 39.98 | 41.54 |
https://rpubs.com/Juve_Campos/tablasRMarkdown Método:
Para la comparación de la media de los tratamientos se usó la prueba t-Student. Para la verificación del supuesto de homogeneidad de varianzas se realizó la prueba F de razón de varianzas. Los tamaños de muestras fueron iguales. Se usó un nivel de significancia del 5%
Hipótesis de interés para la media
\[H_o:\mu_1=\mu_2\\ H_1:\mu_1\ne\mu_2\]
Hipótesis de interés para las varianzas
\[H_o: \frac{\sigma^2_1}{\sigma^2_2}=1\\ H_1:\frac{\sigma^2_1}{\sigma^2_2}\ne1\]
library(readxl)
setwd("C:/Users/UT/Desktop")
datos=read_xlsx("Epoca.xlsx")
str(datos)
## tibble [20 × 2] (S3: tbl_df/tbl/data.frame)
## $ Epoca: chr [1:20] "Abril" "Abril" "Abril" "Abril" ...
## $ Peso : num [1:20] 35.6 36.9 47 44.4 42 ...
Analsis descriptivo
`
aggregate(datos$Peso, by=list(datos$Epoca), mean)
## Group.1 x
## 1 Abril 40.885
## 2 Agosto 26.650
aggregate(datos$Peso, by=list(datos$Epoca), var)
## Group.1 x
## 1 Abril 12.81305
## 2 Agosto 17.25249
aggregate(datos$Peso, by=list(datos$Epoca), sd)
## Group.1 x
## 1 Abril 3.579532
## 2 Agosto 4.153612
boxplot(datos$Peso~datos$Epoca,
main="Box Plot para Peso según Época",
xlab = "Época",
ylab = "Peso en gramos")
Analsis inferencial:
# Verifico si se cumple el supuesto de homogeneidad de varianzas
var.test(Peso~Epoca, datos)
##
## F test to compare two variances
##
## data: Peso by Epoca
## F = 0.74268, num df = 9, denom df = 9, p-value = 0.6648
## alternative hypothesis: true ratio of variances is not equal to 1
## 95 percent confidence interval:
## 0.1844708 2.9900187
## sample estimates:
## ratio of variances
## 0.7426784