JEISON

library(readxl)

## Warning: package 'readxl' was built under R version 3.4.2

DATOS_BIOSTA <- read_excel("C:/Users/usuario/Desktop/DATOS BIOSTA.xlsx", 
    sheet = "Hoja1")
View(DATOS_BIOSTA)

attach(DATOS_BIOSTA)

tapply(Promedio_global,Edad,mean)

##        1        2        3 
## 7.599600 7.297273 7.103333

tapply(Promedio_global,Edad,sd)

##         1         2         3 
## 0.9266089 0.6880073 1.3202399

tapply(Promedio_global,Edad,length)

##  1  2  3 
## 25 22  3

ANOVA1<-aov(Promedio_global~Edad)
summary(ANOVA1)

##             Df Sum Sq Mean Sq F value Pr(>F)
## Edad         1   1.40  1.4033   1.978  0.166
## Residuals   48  34.05  0.7095

shapiro.test(residuals(ANOVA1))

## 
##  Shapiro-Wilk normality test
## 
## data:  residuals(ANOVA1)
## W = 0.98376, p-value = 0.7174

library(car)

## Warning: package 'car' was built under R version 3.4.2

levene.test(Promedio_global~as.factor(Edad))

## Warning: 'levene.test' is deprecated.
## Use 'leveneTest' instead.
## See help("Deprecated") and help("car-deprecated").

## Levene's Test for Homogeneity of Variance (center = median)
##       Df F value Pr(>F)
## group  2  0.8227 0.4455
##       47

drop1(ANOVA1,test = "F")

## Single term deletions
## 
## Model:
## Promedio_global ~ Edad
##        Df Sum of Sq    RSS     AIC F value Pr(>F)
## <none>              34.054 -15.204               
## Edad    1    1.4033 35.457 -15.184   1.978  0.166

Se encontro que el promedio global medio de los estudiantes de biologia y ensenanza de las ciencias, fue mayor en la primera categoria de edad (17-20 anos) (7.6±0.93, n=25). Seguido de la segunda categoria de edad (21-24 anos) (7.3±0.69, n=22). Por ultimo fue la tercera categoria de edad (mas de 25 anos)(7.1±1.32, n=3).