library(readxl)
AFF <- read_excel("I:/Personal/biologia/Bioestadistica prac/datos2.xlsx", sheet = "AF-2Fac")
attach(AFF)
cond=as.factor(Condic)
tapply(AF,cond,length)
## Enfermo con Hongo Enfermo sin Hongo Sano con Hongo Sano sin Hongo
## 46 14 47 13
tapply(AF,cond,mean)
## Enfermo con Hongo Enfermo sin Hongo Sano con Hongo Sano sin Hongo
## 39.01974 40.17000 61.14213 68.25000
tapply(AF,cond,sd)
## Enfermo con Hongo Enfermo sin Hongo Sano con Hongo Sano sin Hongo
## 9.380503 10.695835 24.375706 25.888197
shapiro.test(aov(AF~cond)$residuals)
##
## Shapiro-Wilk normality test
##
## data: aov(AF ~ cond)$residuals
## W = 0.97411, p-value = 0.02041
fligner.test(AF~cond)
##
## Fligner-Killeen test of homogeneity of variances
##
## data: AF by cond
## Fligner-Killeen:med chi-squared = 27.994, df = 3, p-value =
## 3.643e-06
kruskal.test(AF~cond)
##
## Kruskal-Wallis rank sum test
##
## data: AF by cond
## Kruskal-Wallis chi-squared = 34.103, df = 3, p-value = 1.884e-07
pairwise.wilcox.test(AF,cond,p.adj= "none",exact=FALSE)
##
## Pairwise comparisons using Wilcoxon rank sum test
##
## data: AF and cond
##
## Enfermo con Hongo Enfermo sin Hongo Sano con Hongo
## Enfermo sin Hongo 0.85439 - -
## Sano con Hongo 2.7e-06 0.00385 -
## Sano sin Hongo 1.1e-05 0.00074 0.38908
##
## P value adjustment method: none
library(ggplot2)
b = ggplot(AFF,aes(factor(Condic),AF))
b +scale_y_continuous(name = "Área Foliar") + scale_x_discrete(name="Condicion")+ geom_boxplot()+geom_point()
Se encontró que el área foliar en presencia del Matapalo y el Hongo(39.020±9.3805,n= 46) y el Matapalo sin Hongo(40.170±10.6958,n= 14) fue menor que las hojas del árbol sin dicho parasito y con Hongo(61.1420±24.3757,n= 47) y sin Hongo (68.250±25.8882,n= 13). Así mismo se encontraron diferencias significativas para todos los casos exepto para las condiciones Enfermo con Hongo-Enfermo sin Hongo (P= 0.854), Sano sin Hongo-Sano con Hongo(P= 0.389).