No se encontrÓ diferencia significativa entre el nÚmero de individuos de la zona alterada ZA(3.17±3.04, n=30), y el nÚmero de individuos del bosque secundario BS(2.33±1.58, n=30) (chi-squared=0.5,gl=1, p>0.05). De igual forma se realizó una comparación no planeada Bonferroni y no se encontró diferencias significativas entre ZA-BS(P>0.05).
read.table("clipboard",header = T)->bichos
tapply(bichos$abun,bichos$parc,mean)
## BS ZA
## 2.333333 3.166667
tapply(bichos$abun,bichos$parc,sd)
## BS ZA
## 1.582955 3.040909
##normalidad de los residuos
reg<-lm(bichos$abun~bichos$parc)
residuos<-residuals(reg)
shapiro.test(residuos)
##
## Shapiro-Wilk normality test
##
## data: residuos
## W = 0.86719, p-value = 1.018e-05
##homogeneidad de los residuos
fligner.test(bichos$abun~bichos$parc)
##
## Fligner-Killeen test of homogeneity of variances
##
## data: bichos$abun by bichos$parc
## Fligner-Killeen:med chi-squared = 1.9075, df = 1, p-value = 0.1672
kruskal.test(bichos$abun~bichos$parc)
##
## Kruskal-Wallis rank sum test
##
## data: bichos$abun by bichos$parc
## Kruskal-Wallis chi-squared = 0.50943, df = 1, p-value = 0.4754
pairwise.wilcox.test(bichos$abun,bichos$parc, p.adj = "b", exact=F)
##
## Pairwise comparisons using Wilcoxon rank sum test
##
## data: bichos$abun and bichos$parc
##
## BS
## ZA 0.48
##
## P value adjustment method: bonferroni
boxplot(bichos$abun,bichos$parc,xlab="Zonas de muestreo",names=c("Zona alterada","Bosque secundario"),ylab="número de individuos")