library(RcmdrMisc)
## Warning: package 'RcmdrMisc' was built under R version 4.0.5
## Loading required package: car
## Loading required package: carData
## Loading required package: sandwich
CADIZ <-c(182,170,175,167,171,174,181,169,
174,174,170,176,168,178,180)
MALAGA <-c(181,173,177,170,170,175,169,169,
171,173,177,182,179,165,174)
DB_ESTATURA <- data.frame(CADIZ =CADIZ, MALAGA =MALAGA)
DB_EST_APILADA <- stack(DB_ESTATURA[, c("CADIZ","MALAGA")])
names(DB_EST_APILADA) <- c("ESTATURA", "CIUDAD")
DB_EST_APILADA
## ESTATURA CIUDAD
## 1 182 CADIZ
## 2 170 CADIZ
## 3 175 CADIZ
## 4 167 CADIZ
## 5 171 CADIZ
## 6 174 CADIZ
## 7 181 CADIZ
## 8 169 CADIZ
## 9 174 CADIZ
## 10 174 CADIZ
## 11 170 CADIZ
## 12 176 CADIZ
## 13 168 CADIZ
## 14 178 CADIZ
## 15 180 CADIZ
## 16 181 MALAGA
## 17 173 MALAGA
## 18 177 MALAGA
## 19 170 MALAGA
## 20 170 MALAGA
## 21 175 MALAGA
## 22 169 MALAGA
## 23 169 MALAGA
## 24 171 MALAGA
## 25 173 MALAGA
## 26 177 MALAGA
## 27 182 MALAGA
## 28 179 MALAGA
## 29 165 MALAGA
## 30 174 MALAGA
Boxplot(ESTATURA~CIUDAD, data=DB_EST_APILADA, id=list(method="y"))
numSummary(DB_EST_APILADA[,"ESTATURA", drop=FALSE], groups=DB_EST_APILADA$CIUDAD,
statistics=c("mean", "sd", "IQR", "quantiles"), quantiles=c(0,.25,.5,.75,1))
## mean sd IQR 0% 25% 50% 75% 100% ESTATURA:n
## CADIZ 173.9333 4.787882 7 167 170 174 177 182 15
## MALAGA 173.6667 4.850135 7 165 170 173 177 182 15
t.test(ESTATURA~CIUDAD, alternative='two.sided', conf.level=.95, var.equal=FALSE, data=DB_EST_APILADA)
##
## Welch Two Sample t-test
##
## data: ESTATURA by CIUDAD
## t = 0.15154, df = 27.995, p-value = 0.8806
## alternative hypothesis: true difference in means is not equal to 0
## 95 percent confidence interval:
## -3.337921 3.871255
## sample estimates:
## mean in group CADIZ mean in group MALAGA
## 173.9333 173.6667
#95 percent confidence interval: # -3.337921 3.871255
βFin