set.seed(123)
cra = rnorm(80,2.8,0.2)
hist(cra, xlim = c(2.2, 3.4), ylim = c(0, 20), col = "green")
abline(v = mean(cra),col='pink', lwd= 5)
library(readxl)
## Warning: package 'readxl' was built under R version 4.2.2
Libro1 <- read_excel("Libro1.xlsx")#
## New names:
## • `cra60` -> `cra60...1`
## • `cra80` -> `cra80...2`
## • `cra60` -> `cra60...3`
## • `cra80` -> `cra80...5`
#View(Libro1)
plot(Libro1$cra60...3,Libro1$cra80...5,
pch = 16, cex = 1.5,
xlab = "cra 60",ylab = "cra 80")
# coeficiente de correlacion
cor(Libro1$cra60...3,Libro1$cra80...5)
## [1] 0.9770091
Ejemplos en algunas areas
\[H_0: \mu_{20}=\mu_{40} \]
Frutales:se quiere medir el diametro de duraznos a dos fechas, seguramenta estaran correlacionadas ya que a mayor numero de dias se espera tener mayor diametro. Microbiologia: al medir el diametro de colonias creciendo en agar, interezado en comparar su diametro inical con su diametro final, se encuentra patron similar. * ecologia:al medir el oxigeno disuelto en agua,el cual depnede de la temperatura de la laguna y la profundida de la medida, por lo cual es probable que se encuentre una concentracion diferente a una mayor profundida. * cafe : indice *la cra son igual en 60 grados y en 80 grados \[H_0: \mu_{cra_{60}} = \mu_{cra{80}}\]
#prueba t-student de muestras relacionadas.
prueba3 = t.test(Libro1$cra60...3,
Libro1$cra80...5,
alternative = "t",
mu = 0,
paired = TRUE)
ifelse(prueba3$p.value<0.05,"rechazo Ho", "no rechazoHo")
## [1] "rechazo Ho"
boxplot(Libro1$cra60...3,Libro1$cra80...5)