library(readxl)
Estrellas <- read_excel("C:/Users/una/Downloads/Estrellas.xlsx",
col_types = c("numeric", "numeric", "numeric",
"text", "text", "numeric"))
attach(Estrellas)
tapply(Talla_cm,Especie,length)
## Acanthaster planci Pharia pyramidatus Phataria unifascialis
## 93 969 1477
anova1<-aov(Talla_cm~Especie)
residuos<-residuals(anova1)
shapiro.test(residuos)
##
## Shapiro-Wilk normality test
##
## data: residuos
## W = 0.96641, p-value < 2.2e-16
summary(anova1)
## Df Sum Sq Mean Sq F value Pr(>F)
## Especie 2 7001 3500 223.5 <2e-16 ***
## Residuals 2536 39713 16
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
drop1(anova1, test = "F")
## Single term deletions
##
## Model:
## Talla_cm ~ Especie
## Df Sum of Sq RSS AIC F value Pr(>F)
## <none> 39713 6988.0
## Especie 2 7000.8 46714 7396.3 223.53 < 2.2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
anova2<-kruskal.test(Talla_cm,as.factor(Especie))
pairwise.wilcox.test(Talla_cm,as.factor(Especie), p.adj = "none", exact=F)
##
## Pairwise comparisons using Wilcoxon rank sum test
##
## data: Talla_cm and as.factor(Especie)
##
## Acanthaster planci Pharia pyramidatus
## Pharia pyramidatus 3.6e-15 -
## Phataria unifascialis 0.26 < 2e-16
##
## P value adjustment method: none
tapply(Talla_cm,Especie, median)
## Acanthaster planci Pharia pyramidatus Phataria unifascialis
## 14 16 14
tapply(Talla_cm,Especie, sd)
## Acanthaster planci Pharia pyramidatus Phataria unifascialis
## 4.509405 4.245992 3.716870
tapply(Talla_cm,Especie, length)
## Acanthaster planci Pharia pyramidatus Phataria unifascialis
## 93 969 1477
anova2
##
## Kruskal-Wallis rank sum test
##
## data: Talla_cm and as.factor(Especie)
## Kruskal-Wallis chi-squared = 405.14, df = 2, p-value < 2.2e-16
Se encontró que el tamaño medio de las estrellas fue mayor en la especie Pharia pyramidatus (16±4.25,n= 969), seguido de las especies Phataria unifascialis (14±3.72, n= 1477) y Acanthaster planci(14± 4.51, n= 93 ) de manera significativa (chi =400, gl= 2 p< 0.05). Así mismo se realizaron pruebas de pairwise.wilcox y se encontrarón diferencias significativas entre Pharia pyramidatus y Phataria unifascialis (p< 0.05) y entre Pharia pyramidatus y Acanthaster planci (p<0.05).