##Caso 2: Pruebas de hipotesis para dos muestras pareadas
set.seed(123)
cra = rnorm(n = 80, mean = 2.8, sd= 0.2)
hist(cra, xlim = c(2.2, 3.4), ylim = c(0,20))
abline(v = mean(cra), col= 'orange',lwd=4)
library(readxl)
## Warning: package 'readxl' was built under R version 4.2.2
libro1 <- read_excel("C:/Users/fcecursos/Desktop/libro1.xlsx")
## New names:
## • `` -> `...3`
View(libro1)
head(libro1)
## # A tibble: 6 × 4
## cra60 cra80 ...3 ruido
## <dbl> <dbl> <lgl> <dbl>
## 1 1.92 2.62 NA 0.200
## 2 1.98 2.67 NA 0.190
## 3 2.06 2.81 NA 0.245
## 4 1.91 2.61 NA 0.202
## 5 1.81 2.47 NA 0.165
## 6 2.07 2.77 NA 0.207
plot(libro1$cra60,
libro1$cra80,
pch = 16, cex=0.9, xlab = 'CRA60', ylab = 'CRA80')
#coeficiente de correlñacion
cor(libro1$cra60, libro1$cra80)
## [1] 0.9739356
Ejemplos en algunas areas * Suelos: si se mide la MO dwe los primeros 20cm y luego a 40cm posiblemente se encuentre correlacion entre ambas medidas, a mayor profundidad habra menor MO.
*Frutales: Se quiere medir el diametro de durasnos en dois fechas, seguraqmente estaran correlacionadas, ya que a mayor numero de dias se espera tenewr mayor diametro.
*Microbiologia: Al medir el diametro medio a colonias creciendo en agar, se puede estar intyerezado en comparar un diametyro inicial y uno finmal, al relizar una medicion posterior y se encuentre un patron similar.
*Ecologia: Al medir el oxogeno disuelto en agua, el cual depende de la temperatura de la laguna y la profundidad de la medida, por lo cual es probable a una profumdidadf se encuentre una concentraciuon diferente a una profumndidad mayor, seguramente estara correlacionadas. \[H-0: \mu_{oxigeno_{50cm}} = \mu{oxigeno_{100cm}}\]
Cafe: En la maduracion de grano de cafe, cuando se m9ide color exitente en diferentes esopacios nde color, entre la cualñes puyedan estar LAB o RGB con los que a su vez se puden calcular diferentes indices (por ejemplo $GLI=/frac{2G -(R+B)})
\(H-0: \mu_{oxigeno_{50cm}} = \mu{oxigeno_{100cm}}\)
\[H_0: \mu_{CRA_{60}} = \mu{CRA{80}}\]
#Prueba 3
#Prube t-student para dos muestras relacionadas o pareadas
prueba3 = t.test(libro1$cra60,
libro1$cra80,
alternative = 't', # t = two.sided
mu = 0,
paired = TRUE,)
ifelse(prueba3$p.value<0.05, 'Rechazo H0', 'NO rechazo H0')
## [1] "Rechazo H0"
boxplot(libro1)