Graficas de muestreos
Andres Felipe Cortes Caucali
2022-06-18
#Librerias
library(readr)
library(dplyr)##
## Attaching package: 'dplyr'## The following objects are masked from 'package:stats':
##
## filter, lag## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, unionlibrary(ggplot2)
library(lattice)
library(tigerstats)## Loading required package: abd## Loading required package: nlme##
## Attaching package: 'nlme'## The following object is masked from 'package:dplyr':
##
## collapse## Loading required package: grid## Loading required package: mosaic## Registered S3 method overwritten by 'mosaic':
## method from
## fortify.SpatialPolygonsDataFrame ggplot2##
## The 'mosaic' package masks several functions from core packages in order to add
## additional features. The original behavior of these functions should not be affected by this.##
## Attaching package: 'mosaic'## The following object is masked from 'package:Matrix':
##
## mean## The following object is masked from 'package:ggplot2':
##
## stat## The following objects are masked from 'package:dplyr':
##
## count, do, tally## The following objects are masked from 'package:stats':
##
## binom.test, cor, cor.test, cov, fivenum, IQR, median, prop.test,
## quantile, sd, t.test, var## The following objects are masked from 'package:base':
##
## max, mean, min, prod, range, sample, sum## Welcome to tigerstats!
## To learn more about this package, consult its website:
## http://homerhanumat.github.io/tigerstatslibrary(multcompView)
library(readxl)Potasio <- read_excel("C:\\Users\\USER\\Documents\\UNIVERSIDAD F\\2022-1\\Fisiologia vegetal\\Analisis-de-datos_files\\Cuaderno de campo virtual, datos.xlsx",
sheet = "Datos totales")
df1=data.frame(Potasio)#Analisis de Varianza de dos factores
#Organizacion de datos
head(df1)str(df1)## 'data.frame': 64 obs. of 19 variables:
## $ Muestreo : chr "Muestreo 1" "Muestreo 1" "Muestreo 1" "Muestreo 1" ...
## $ Tratamientos : chr "Control" "Control" "Control" "Control" ...
## $ Repeticion : num 1 2 3 4 1 2 3 4 1 2 ...
## $ Clorofila : num 38.3 35.7 40 38.9 36.6 36.6 32.8 41.4 46.2 40.7 ...
## $ Temperatura : num 18.4 19.6 18 17.8 15.8 20.4 20.6 20.6 18.4 22.8 ...
## $ AreaF : num 253 292 256 173 190 ...
## $ NumH : num 6 7 8 6 6 6 5 7 9 8 ...
## $ LongitudAerea: num 7.2 8.3 7.1 7.5 18.9 20.8 17 22.7 16.8 19 ...
## $ EstomasA : num 4 6 11 5 3 5 1 2 6 15 ...
## $ EstomasC : num 14 3 7 4 1 4 6 3 3 1 ...
## $ EstomasT : num 39 39 41 40 31 35 34 25 14 34 ...
## $ PesoFH : num 10 15 14 13 17 18 19 15 20 21 ...
## $ PesoSH : num 0.8 1.3 1.2 0.6 0.9 1.8 0.8 1.7 1.4 1.2 ...
## $ DiametroR : num 31.7 30.3 33.1 33.7 37.6 30.1 36 40 44.1 27.6 ...
## $ PesoFR : num 18.2 15.9 17.8 20 34.5 22.9 27 27.1 27.4 13.9 ...
## $ PesoSR : num 1.6 1.2 1.3 1.4 2.4 2.7 2.6 2.8 2.8 1.9 ...
## $ CRA : num 65 68 70 70 75 70 80 70 65 60 ...
## $ RAF : num 8.97 9.44 8.05 5.24 3.68 ...
## $ AFE : num 25.3 19.4 18.3 13.3 11.2 ...Los Factores Muestreo y tratamientos que estan como caracteres , se transforman en Factores para el analisis de varianzas
# creating a variable as factor for the ANOVA
df1$Muestreo <- as.factor(df1$Muestreo)
df1$Tratamientos <- as.factor(df1$Tratamientos)
str(df1)## 'data.frame': 64 obs. of 19 variables:
## $ Muestreo : Factor w/ 4 levels "Muestreo 1","Muestreo 2",..: 1 1 1 1 2 2 2 2 3 3 ...
## $ Tratamientos : Factor w/ 4 levels "Control","NaCl",..: 1 1 1 1 1 1 1 1 1 1 ...
## $ Repeticion : num 1 2 3 4 1 2 3 4 1 2 ...
## $ Clorofila : num 38.3 35.7 40 38.9 36.6 36.6 32.8 41.4 46.2 40.7 ...
## $ Temperatura : num 18.4 19.6 18 17.8 15.8 20.4 20.6 20.6 18.4 22.8 ...
## $ AreaF : num 253 292 256 173 190 ...
## $ NumH : num 6 7 8 6 6 6 5 7 9 8 ...
## $ LongitudAerea: num 7.2 8.3 7.1 7.5 18.9 20.8 17 22.7 16.8 19 ...
## $ EstomasA : num 4 6 11 5 3 5 1 2 6 15 ...
## $ EstomasC : num 14 3 7 4 1 4 6 3 3 1 ...
## $ EstomasT : num 39 39 41 40 31 35 34 25 14 34 ...
## $ PesoFH : num 10 15 14 13 17 18 19 15 20 21 ...
## $ PesoSH : num 0.8 1.3 1.2 0.6 0.9 1.8 0.8 1.7 1.4 1.2 ...
## $ DiametroR : num 31.7 30.3 33.1 33.7 37.6 30.1 36 40 44.1 27.6 ...
## $ PesoFR : num 18.2 15.9 17.8 20 34.5 22.9 27 27.1 27.4 13.9 ...
## $ PesoSR : num 1.6 1.2 1.3 1.4 2.4 2.7 2.6 2.8 2.8 1.9 ...
## $ CRA : num 65 68 70 70 75 70 80 70 65 60 ...
## $ RAF : num 8.97 9.44 8.05 5.24 3.68 ...
## $ AFE : num 25.3 19.4 18.3 13.3 11.2 ...#Anova 1 CRA 2 factores
# analysis of variance
anova1 <- aov(df1$AFE ~ df1$Tratamientos*df1$Muestreo, data = df1)
summary(anova1)## Df Sum Sq Mean Sq F value Pr(>F)
## df1$Tratamientos 3 201.9 67.3 3.487 0.0227 *
## df1$Muestreo 3 972.6 324.2 16.802 1.35e-07 ***
## df1$Tratamientos:df1$Muestreo 9 33.2 3.7 0.191 0.9941
## Residuals 48 926.2 19.3
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1#Crear una tabla con factores, medias y desviación estándar
# table with factors, means and standard deviation
data_summary <- group_by(df1, Muestreo, Tratamientos) %>%
summarise(mean=mean(AFE), sd=sd(AFE)) %>%
arrange(desc(mean))## `summarise()` has grouped output by 'Muestreo'. You can override using the
## `.groups` argument.print(data_summary)## # A tibble: 16 × 4
## # Groups: Muestreo [4]
## Muestreo Tratamientos mean sd
## <fct> <fct> <dbl> <dbl>
## 1 Muestreo 1 NaCl 24.7 4.84
## 2 Muestreo 1 NaCl+ KCL 2 24.5 3.01
## 3 Muestreo 2 NaCl 20.7 4.87
## 4 Muestreo 1 NaCl+ KCL 1 20.3 5.26
## 5 Muestreo 1 Control 19.1 4.93
## 6 Muestreo 2 NaCl+ KCL 2 18.4 1.92
## 7 Muestreo 3 NaCl+ KCL 2 17.8 3.63
## 8 Muestreo 2 NaCl+ KCL 1 17.8 6.02
## 9 Muestreo 3 NaCl 16.1 3.68
## 10 Muestreo 2 Control 15.8 4.99
## 11 Muestreo 3 NaCl+ KCL 1 14.1 6.21
## 12 Muestreo 4 NaCl+ KCL 2 13.3 6.99
## 13 Muestreo 3 Control 13.2 4.03
## 14 Muestreo 4 NaCl 12.4 1.08
## 15 Muestreo 4 NaCl+ KCL 1 11.0 1.73
## 16 Muestreo 4 Control 9.31 1.09#Prueba de Tukey, para comparacion de medias
# Tukey's test
tukey<- TukeyHSD(anova1)
print(tukey)## Tukey multiple comparisons of means
## 95% family-wise confidence level
##
## Fit: aov(formula = df1$AFE ~ df1$Tratamientos * df1$Muestreo, data = df1)
##
## $`df1$Tratamientos`
## diff lwr upr p adj
## NaCl-Control 4.124957451 -0.0082484820 8.258163 0.0506348
## NaCl+ KCL 1-Control 1.462762577 -2.6704433564 5.595969 0.7825210
## NaCl+ KCL 2-Control 4.133761163 0.0005552301 8.266967 0.0499575
## NaCl+ KCL 1-NaCl -2.662194874 -6.7954008076 1.471011 0.3276916
## NaCl+ KCL 2-NaCl 0.008803712 -4.1244022212 4.142010 0.9999999
## NaCl+ KCL 2-NaCl+ KCL 1 2.670998586 -1.4622073468 6.804205 0.3248282
##
## $`df1$Muestreo`
## diff lwr upr p adj
## Muestreo 2-Muestreo 1 -4.001618 -8.134823 0.1315884 0.0609949
## Muestreo 3-Muestreo 1 -6.832029 -10.965235 -2.6988232 0.0003417
## Muestreo 4-Muestreo 1 -10.655787 -14.788993 -6.5225809 0.0000001
## Muestreo 3-Muestreo 2 -2.830412 -6.963618 1.3027943 0.2755416
## Muestreo 4-Muestreo 2 -6.654169 -10.787375 -2.5209634 0.0004939
## Muestreo 4-Muestreo 3 -3.823758 -7.956964 0.3094483 0.0791072
##
## $`df1$Tratamientos:df1$Muestreo`
## diff lwr
## NaCl:Muestreo 1-Control:Muestreo 1 5.67119436 -5.550083
## NaCl+ KCL 1:Muestreo 1-Control:Muestreo 1 1.21873127 -10.002546
## NaCl+ KCL 2:Muestreo 1-Control:Muestreo 1 5.40326507 -5.818012
## Control:Muestreo 2-Control:Muestreo 1 -3.29237062 -14.513648
## NaCl:Muestreo 2-Control:Muestreo 1 1.60405678 -9.617221
## NaCl+ KCL 1:Muestreo 2-Control:Muestreo 1 -1.30798750 -12.529265
## NaCl+ KCL 2:Muestreo 2-Control:Muestreo 1 -0.71697802 -11.938255
## Control:Muestreo 3-Control:Muestreo 1 -5.86038895 -17.081666
## NaCl:Muestreo 3-Control:Muestreo 1 -2.98596539 -14.207243
## NaCl+ KCL 1:Muestreo 3-Control:Muestreo 1 -4.92785136 -16.149129
## NaCl+ KCL 2:Muestreo 3-Control:Muestreo 1 -1.26072030 -12.481998
## Control:Muestreo 4-Control:Muestreo 1 -9.76496441 -20.986242
## NaCl:Muestreo 4-Control:Muestreo 1 -6.70717993 -17.928457
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 1 -8.04956609 -19.270844
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 1 -5.80824608 -17.029524
## NaCl+ KCL 1:Muestreo 1-NaCl:Muestreo 1 -4.45246309 -15.673741
## NaCl+ KCL 2:Muestreo 1-NaCl:Muestreo 1 -0.26792929 -11.489207
## Control:Muestreo 2-NaCl:Muestreo 1 -8.96356498 -20.184842
## NaCl:Muestreo 2-NaCl:Muestreo 1 -4.06713758 -15.288415
## NaCl+ KCL 1:Muestreo 2-NaCl:Muestreo 1 -6.97918186 -18.200459
## NaCl+ KCL 2:Muestreo 2-NaCl:Muestreo 1 -6.38817238 -17.609450
## Control:Muestreo 3-NaCl:Muestreo 1 -11.53158331 -22.752861
## NaCl:Muestreo 3-NaCl:Muestreo 1 -8.65715976 -19.878437
## NaCl+ KCL 1:Muestreo 3-NaCl:Muestreo 1 -10.59904572 -21.820323
## NaCl+ KCL 2:Muestreo 3-NaCl:Muestreo 1 -6.93191466 -18.153192
## Control:Muestreo 4-NaCl:Muestreo 1 -15.43615877 -26.657436
## NaCl:Muestreo 4-NaCl:Muestreo 1 -12.37837429 -23.599652
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 1 -13.72076045 -24.942038
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 1 -11.47944044 -22.700718
## NaCl+ KCL 2:Muestreo 1-NaCl+ KCL 1:Muestreo 1 4.18453380 -7.036744
## Control:Muestreo 2-NaCl+ KCL 1:Muestreo 1 -4.51110189 -15.732379
## NaCl:Muestreo 2-NaCl+ KCL 1:Muestreo 1 0.38532551 -10.835952
## NaCl+ KCL 1:Muestreo 2-NaCl+ KCL 1:Muestreo 1 -2.52671877 -13.747996
## NaCl+ KCL 2:Muestreo 2-NaCl+ KCL 1:Muestreo 1 -1.93570929 -13.156987
## Control:Muestreo 3-NaCl+ KCL 1:Muestreo 1 -7.07912022 -18.300398
## NaCl:Muestreo 3-NaCl+ KCL 1:Muestreo 1 -4.20469666 -15.425974
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 1:Muestreo 1 -6.14658263 -17.367860
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 1:Muestreo 1 -2.47945156 -13.700729
## Control:Muestreo 4-NaCl+ KCL 1:Muestreo 1 -10.98369568 -22.204973
## NaCl:Muestreo 4-NaCl+ KCL 1:Muestreo 1 -7.92591119 -19.147189
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 1:Muestreo 1 -9.26829735 -20.489575
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 1 -7.02697735 -18.248255
## Control:Muestreo 2-NaCl+ KCL 2:Muestreo 1 -8.69563569 -19.916913
## NaCl:Muestreo 2-NaCl+ KCL 2:Muestreo 1 -3.79920829 -15.020486
## NaCl+ KCL 1:Muestreo 2-NaCl+ KCL 2:Muestreo 1 -6.71125257 -17.932530
## NaCl+ KCL 2:Muestreo 2-NaCl+ KCL 2:Muestreo 1 -6.12024309 -17.341521
## Control:Muestreo 3-NaCl+ KCL 2:Muestreo 1 -11.26365402 -22.484931
## NaCl:Muestreo 3-NaCl+ KCL 2:Muestreo 1 -8.38923046 -19.610508
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 2:Muestreo 1 -10.33111643 -21.552394
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 2:Muestreo 1 -6.66398536 -17.885263
## Control:Muestreo 4-NaCl+ KCL 2:Muestreo 1 -15.16822948 -26.389507
## NaCl:Muestreo 4-NaCl+ KCL 2:Muestreo 1 -12.11044499 -23.331722
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 2:Muestreo 1 -13.45283115 -24.674109
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 2:Muestreo 1 -11.21151115 -22.432789
## NaCl:Muestreo 2-Control:Muestreo 2 4.89642740 -6.324850
## NaCl+ KCL 1:Muestreo 2-Control:Muestreo 2 1.98438312 -9.236894
## NaCl+ KCL 2:Muestreo 2-Control:Muestreo 2 2.57539260 -8.645885
## Control:Muestreo 3-Control:Muestreo 2 -2.56801833 -13.789296
## NaCl:Muestreo 3-Control:Muestreo 2 0.30640523 -10.914872
## NaCl+ KCL 1:Muestreo 3-Control:Muestreo 2 -1.63548074 -12.856758
## NaCl+ KCL 2:Muestreo 3-Control:Muestreo 2 2.03165033 -9.189627
## Control:Muestreo 4-Control:Muestreo 2 -6.47259379 -17.693871
## NaCl:Muestreo 4-Control:Muestreo 2 -3.41480930 -14.636087
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 2 -4.75719546 -15.978473
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 2 -2.51587546 -13.737153
## NaCl+ KCL 1:Muestreo 2-NaCl:Muestreo 2 -2.91204428 -14.133322
## NaCl+ KCL 2:Muestreo 2-NaCl:Muestreo 2 -2.32103480 -13.542312
## Control:Muestreo 3-NaCl:Muestreo 2 -7.46444573 -18.685723
## NaCl:Muestreo 3-NaCl:Muestreo 2 -4.59002217 -15.811300
## NaCl+ KCL 1:Muestreo 3-NaCl:Muestreo 2 -6.53190814 -17.753186
## NaCl+ KCL 2:Muestreo 3-NaCl:Muestreo 2 -2.86477707 -14.086055
## Control:Muestreo 4-NaCl:Muestreo 2 -11.36902119 -22.590299
## NaCl:Muestreo 4-NaCl:Muestreo 2 -8.31123670 -19.532514
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 2 -9.65362286 -20.874900
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 2 -7.41230286 -18.633580
## NaCl+ KCL 2:Muestreo 2-NaCl+ KCL 1:Muestreo 2 0.59100948 -10.630268
## Control:Muestreo 3-NaCl+ KCL 1:Muestreo 2 -4.55240145 -15.773679
## NaCl:Muestreo 3-NaCl+ KCL 1:Muestreo 2 -1.67797789 -12.899255
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 1:Muestreo 2 -3.61986386 -14.841141
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 1:Muestreo 2 0.04726721 -11.174010
## Control:Muestreo 4-NaCl+ KCL 1:Muestreo 2 -8.45697691 -19.678254
## NaCl:Muestreo 4-NaCl+ KCL 1:Muestreo 2 -5.39919242 -16.620470
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 1:Muestreo 2 -6.74157858 -17.962856
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 2 -4.50025858 -15.721536
## Control:Muestreo 3-NaCl+ KCL 2:Muestreo 2 -5.14341093 -16.364688
## NaCl:Muestreo 3-NaCl+ KCL 2:Muestreo 2 -2.26898737 -13.490265
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 2:Muestreo 2 -4.21087334 -15.432151
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 2:Muestreo 2 -0.54374227 -11.765020
## Control:Muestreo 4-NaCl+ KCL 2:Muestreo 2 -9.04798639 -20.269264
## NaCl:Muestreo 4-NaCl+ KCL 2:Muestreo 2 -5.99020190 -17.211479
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 2:Muestreo 2 -7.33258806 -18.553866
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 2:Muestreo 2 -5.09126806 -16.312546
## NaCl:Muestreo 3-Control:Muestreo 3 2.87442356 -8.346854
## NaCl+ KCL 1:Muestreo 3-Control:Muestreo 3 0.93253759 -10.288740
## NaCl+ KCL 2:Muestreo 3-Control:Muestreo 3 4.59966866 -6.621609
## Control:Muestreo 4-Control:Muestreo 3 -3.90457546 -15.125853
## NaCl:Muestreo 4-Control:Muestreo 3 -0.84679097 -12.068068
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 3 -2.18917713 -13.410455
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 3 0.05214287 -11.169135
## NaCl+ KCL 1:Muestreo 3-NaCl:Muestreo 3 -1.94188596 -13.163163
## NaCl+ KCL 2:Muestreo 3-NaCl:Muestreo 3 1.72524510 -9.496032
## Control:Muestreo 4-NaCl:Muestreo 3 -6.77899902 -18.000276
## NaCl:Muestreo 4-NaCl:Muestreo 3 -3.72121453 -14.942492
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 3 -5.06360069 -16.284878
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 3 -2.82228069 -14.043558
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 1:Muestreo 3 3.66713106 -7.554146
## Control:Muestreo 4-NaCl+ KCL 1:Muestreo 3 -4.83711305 -16.058391
## NaCl:Muestreo 4-NaCl+ KCL 1:Muestreo 3 -1.77932857 -13.000606
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 1:Muestreo 3 -3.12171473 -14.342992
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 3 -0.88039472 -12.101672
## Control:Muestreo 4-NaCl+ KCL 2:Muestreo 3 -8.50424411 -19.725522
## NaCl:Muestreo 4-NaCl+ KCL 2:Muestreo 3 -5.44645963 -16.667737
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 2:Muestreo 3 -6.78884579 -18.010123
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 2:Muestreo 3 -4.54752579 -15.768803
## NaCl:Muestreo 4-Control:Muestreo 4 3.05778449 -8.163493
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 4 1.71539832 -9.505879
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 4 3.95671833 -7.264559
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 4 -1.34238616 -12.563664
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 4 0.89893384 -10.322344
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 4 2.24132000 -8.979957
## upr p adj
## NaCl:Muestreo 1-Control:Muestreo 1 16.892471836 0.8945413
## NaCl+ KCL 1:Muestreo 1-Control:Muestreo 1 12.440008743 1.0000000
## NaCl+ KCL 2:Muestreo 1-Control:Muestreo 1 16.624542543 0.9249395
## Control:Muestreo 2-Control:Muestreo 1 7.928906852 0.9993361
## NaCl:Muestreo 2-Control:Muestreo 1 12.825334251 0.9999999
## NaCl+ KCL 1:Muestreo 2-Control:Muestreo 1 9.913289972 1.0000000
## NaCl+ KCL 2:Muestreo 2-Control:Muestreo 1 10.504299453 1.0000000
## Control:Muestreo 3-Control:Muestreo 1 5.360888522 0.8690852
## NaCl:Muestreo 3-Control:Muestreo 1 8.235312080 0.9997887
## NaCl+ KCL 1:Muestreo 3-Control:Muestreo 1 6.293426116 0.9632551
## NaCl+ KCL 2:Muestreo 3-Control:Muestreo 1 9.960557179 1.0000000
## Control:Muestreo 4-Control:Muestreo 1 1.456313064 0.1534731
## NaCl:Muestreo 4-Control:Muestreo 1 4.514097550 0.7189487
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 1 3.171711389 0.4277295
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 1 5.413031391 0.8764285
## NaCl+ KCL 1:Muestreo 1-NaCl:Muestreo 1 6.768814382 0.9847997
## NaCl+ KCL 2:Muestreo 1-NaCl:Muestreo 1 10.953348182 1.0000000
## Control:Muestreo 2-NaCl:Muestreo 1 2.257712491 0.2592755
## NaCl:Muestreo 2-NaCl:Muestreo 1 7.154139890 0.9936020
## NaCl+ KCL 1:Muestreo 2-NaCl:Muestreo 1 4.242095611 0.6614535
## NaCl+ KCL 2:Muestreo 2-NaCl:Muestreo 1 4.833105092 0.7816333
## Control:Muestreo 3-NaCl:Muestreo 1 -0.310305840 0.0384530
## NaCl:Muestreo 3-NaCl:Muestreo 1 2.564117719 0.3104395
## NaCl+ KCL 1:Muestreo 3-NaCl:Muestreo 1 0.622231754 0.0826840
## NaCl+ KCL 2:Muestreo 3-NaCl:Muestreo 1 4.289362817 0.6716440
## Control:Muestreo 4-NaCl:Muestreo 1 -4.214881297 0.0008760
## NaCl:Muestreo 4-NaCl:Muestreo 1 -1.157096812 0.0181351
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 1 -2.499482973 0.0050465
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 1 -0.258162970 0.0402083
## NaCl+ KCL 2:Muestreo 1-NaCl+ KCL 1:Muestreo 1 15.405811274 0.9915389
## Control:Muestreo 2-NaCl+ KCL 1:Muestreo 1 6.710175583 0.9828755
## NaCl:Muestreo 2-NaCl+ KCL 1:Muestreo 1 11.606602983 1.0000000
## NaCl+ KCL 1:Muestreo 2-NaCl+ KCL 1:Muestreo 1 8.694558703 0.9999734
## NaCl+ KCL 2:Muestreo 2-NaCl+ KCL 1:Muestreo 1 9.285568184 0.9999992
## Control:Muestreo 3-NaCl+ KCL 1:Muestreo 1 4.142157253 0.6396992
## NaCl:Muestreo 3-NaCl+ KCL 1:Muestreo 1 7.016580812 0.9911357
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 1:Muestreo 1 5.074694847 0.8244965
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 1:Muestreo 1 8.741825910 0.9999791
## Control:Muestreo 4-NaCl+ KCL 1:Muestreo 1 0.237581795 0.0608208
## NaCl:Muestreo 4-NaCl+ KCL 1:Muestreo 1 3.295366281 0.4537085
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 1:Muestreo 1 1.952980120 0.2142764
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 1 4.194300122 0.6510819
## Control:Muestreo 2-NaCl+ KCL 2:Muestreo 1 2.525641783 0.3036972
## NaCl:Muestreo 2-NaCl+ KCL 2:Muestreo 1 7.422069183 0.9968084
## NaCl+ KCL 1:Muestreo 2-NaCl+ KCL 2:Muestreo 1 4.510024904 0.7181112
## NaCl+ KCL 2:Muestreo 2-NaCl+ KCL 2:Muestreo 1 5.101034384 0.8288928
## Control:Muestreo 3-NaCl+ KCL 2:Muestreo 1 -0.042376547 0.0482600
## NaCl:Muestreo 3-NaCl+ KCL 2:Muestreo 1 2.832047012 0.3597806
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 2:Muestreo 1 0.890161047 0.1016052
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 2:Muestreo 1 4.557292110 0.7277795
## Control:Muestreo 4-NaCl+ KCL 2:Muestreo 1 -3.946952004 0.0011594
## NaCl:Muestreo 4-NaCl+ KCL 2:Muestreo 1 -0.889167519 0.0231235
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 2:Muestreo 1 -2.231553680 0.0065628
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 2:Muestreo 1 0.009766323 0.0504088
## NaCl:Muestreo 2-Control:Muestreo 2 16.117704874 0.9651440
## NaCl+ KCL 1:Muestreo 2-Control:Muestreo 2 13.205660595 0.9999989
## NaCl+ KCL 2:Muestreo 2-Control:Muestreo 2 13.796670076 0.9999660
## Control:Muestreo 3-Control:Muestreo 2 8.653259145 0.9999672
## NaCl:Muestreo 3-Control:Muestreo 2 11.527682704 1.0000000
## NaCl+ KCL 1:Muestreo 3-Control:Muestreo 2 9.585796739 0.9999999
## NaCl+ KCL 2:Muestreo 3-Control:Muestreo 2 13.252927802 0.9999985
## Control:Muestreo 4-Control:Muestreo 2 4.748683687 0.7656557
## NaCl:Muestreo 4-Control:Muestreo 2 7.806468173 0.9989963
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 2 6.464082012 0.9726609
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 2 8.705402014 0.9999748
## NaCl+ KCL 1:Muestreo 2-NaCl:Muestreo 2 8.309233195 0.9998438
## NaCl+ KCL 2:Muestreo 2-NaCl:Muestreo 2 8.900242676 0.9999911
## Control:Muestreo 3-NaCl:Muestreo 2 3.756831745 0.5544150
## NaCl:Muestreo 3-NaCl:Muestreo 2 6.631255304 0.9799921
## NaCl+ KCL 1:Muestreo 3-NaCl:Muestreo 2 4.689369339 0.7541511
## NaCl+ KCL 2:Muestreo 3-NaCl:Muestreo 2 8.356500402 0.9998720
## Control:Muestreo 4-NaCl:Muestreo 2 -0.147743712 0.0441645
## NaCl:Muestreo 4-NaCl:Muestreo 2 2.910040773 0.3748871
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 2 1.567654612 0.1657903
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 2 3.808974615 0.5659942
## NaCl+ KCL 2:Muestreo 2-NaCl+ KCL 1:Muestreo 2 11.812286955 1.0000000
## Control:Muestreo 3-NaCl+ KCL 1:Muestreo 2 6.668876024 0.9814100
## NaCl:Muestreo 3-NaCl+ KCL 1:Muestreo 2 9.543299583 0.9999999
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 1:Muestreo 2 7.601413618 0.9980923
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 1:Muestreo 2 11.268544681 1.0000000
## Control:Muestreo 4-NaCl+ KCL 1:Muestreo 2 2.764300567 0.3469219
## NaCl:Muestreo 4-NaCl+ KCL 1:Muestreo 2 5.822085052 0.9253507
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 1:Muestreo 2 4.479698891 0.7118498
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 2 6.721018894 0.9832449
## Control:Muestreo 3-NaCl+ KCL 2:Muestreo 2 6.077866544 0.9482325
## NaCl:Muestreo 3-NaCl+ KCL 2:Muestreo 2 8.952290102 0.9999934
## NaCl+ KCL 1:Muestreo 3-NaCl+ KCL 2:Muestreo 2 7.010404138 0.9910092
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 2:Muestreo 2 10.677535201 1.0000000
## Control:Muestreo 4-NaCl+ KCL 2:Muestreo 2 2.173291086 0.2462151
## NaCl:Muestreo 4-NaCl+ KCL 2:Muestreo 2 5.231075572 0.8497424
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 2:Muestreo 2 3.888689411 0.5837051
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 2:Muestreo 2 6.130009413 0.9522085
## NaCl:Muestreo 3-Control:Muestreo 3 14.095701034 0.9998667
## NaCl+ KCL 1:Muestreo 3-Control:Muestreo 3 12.153815069 1.0000000
## NaCl+ KCL 2:Muestreo 3-Control:Muestreo 3 15.820946132 0.9796155
## Control:Muestreo 4-Control:Muestreo 3 7.316702017 0.9957623
## NaCl:Muestreo 4-Control:Muestreo 3 10.374486503 1.0000000
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 3 9.032100342 0.9999959
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 3 11.273420344 1.0000000
## NaCl+ KCL 1:Muestreo 3-NaCl:Muestreo 3 9.279391510 0.9999992
## NaCl+ KCL 2:Muestreo 3-NaCl:Muestreo 3 12.946522573 0.9999998
## Control:Muestreo 4-NaCl:Muestreo 3 4.442278458 0.7040642
## NaCl:Muestreo 4-NaCl:Muestreo 3 7.500062944 0.9974353
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 3 6.157676783 0.9542275
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 3 8.398996785 0.9998934
## NaCl+ KCL 2:Muestreo 3-NaCl+ KCL 1:Muestreo 3 14.888408538 0.9978062
## Control:Muestreo 4-NaCl+ KCL 1:Muestreo 3 6.384164423 0.9685129
## NaCl:Muestreo 4-NaCl+ KCL 1:Muestreo 3 9.441948909 0.9999998
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 1:Muestreo 3 8.099562748 0.9996418
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 3 10.340882750 1.0000000
## Control:Muestreo 4-NaCl+ KCL 2:Muestreo 3 2.717033360 0.3381002
## NaCl:Muestreo 4-NaCl+ KCL 2:Muestreo 3 5.774817846 0.9204850
## NaCl+ KCL 1:Muestreo 4-NaCl+ KCL 2:Muestreo 3 4.432431685 0.7020051
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 2:Muestreo 3 6.673751687 0.9815879
## NaCl:Muestreo 4-Control:Muestreo 4 14.279061960 0.9997195
## NaCl+ KCL 1:Muestreo 4-Control:Muestreo 4 12.936675799 0.9999999
## NaCl+ KCL 2:Muestreo 4-Control:Muestreo 4 15.177995802 0.9951478
## NaCl+ KCL 1:Muestreo 4-NaCl:Muestreo 4 9.878891314 1.0000000
## NaCl+ KCL 2:Muestreo 4-NaCl:Muestreo 4 12.120211316 1.0000000
## NaCl+ KCL 2:Muestreo 4-NaCl+ KCL 1:Muestreo 4 13.462597477 0.9999944##Pantalla de letras compacta para indicar diferencias significativas
# creating the compact letter display
tukey.cld <- multcompLetters4(anova1, tukey)
print(tukey.cld)## $`df1$Tratamientos`
## NaCl+ KCL 2 NaCl NaCl+ KCL 1 Control
## "a" "ab" "ab" "b"
##
## $`df1$Muestreo`
## Muestreo 1 Muestreo 2 Muestreo 3 Muestreo 4
## "a" "ab" "bc" "c"
##
## $`df1$Tratamientos:df1$Muestreo`
## NaCl:Muestreo 1 NaCl+ KCL 2:Muestreo 1 NaCl:Muestreo 2
## "a" "ab" "abc"
## NaCl+ KCL 1:Muestreo 1 Control:Muestreo 1 NaCl+ KCL 2:Muestreo 2
## "abcd" "abcd" "abcd"
## NaCl+ KCL 2:Muestreo 3 NaCl+ KCL 1:Muestreo 2 NaCl:Muestreo 3
## "abcd" "abcd" "abcd"
## Control:Muestreo 2 NaCl+ KCL 1:Muestreo 3 NaCl+ KCL 2:Muestreo 4
## "abcd" "abcd" "bcd"
## Control:Muestreo 3 NaCl:Muestreo 4 NaCl+ KCL 1:Muestreo 4
## "cd" "cd" "cd"
## Control:Muestreo 4
## "d"##Codigo para añadir la letra compacta a la tabla de datos
# adding the compact letter display to the table with means and sd
cld <- as.data.frame.list(tukey.cld$`df1$Tratamientos:df1$Muestreo`)
data_summary$Tukey <- cld$Letters
print(data_summary)## # A tibble: 16 × 5
## # Groups: Muestreo [4]
## Muestreo Tratamientos mean sd Tukey
## <fct> <fct> <dbl> <dbl> <chr>
## 1 Muestreo 1 NaCl 24.7 4.84 a
## 2 Muestreo 1 NaCl+ KCL 2 24.5 3.01 ab
## 3 Muestreo 2 NaCl 20.7 4.87 abc
## 4 Muestreo 1 NaCl+ KCL 1 20.3 5.26 abcd
## 5 Muestreo 1 Control 19.1 4.93 abcd
## 6 Muestreo 2 NaCl+ KCL 2 18.4 1.92 abcd
## 7 Muestreo 3 NaCl+ KCL 2 17.8 3.63 abcd
## 8 Muestreo 2 NaCl+ KCL 1 17.8 6.02 abcd
## 9 Muestreo 3 NaCl 16.1 3.68 abcd
## 10 Muestreo 2 Control 15.8 4.99 abcd
## 11 Muestreo 3 NaCl+ KCL 1 14.1 6.21 abcd
## 12 Muestreo 4 NaCl+ KCL 2 13.3 6.99 bcd
## 13 Muestreo 3 Control 13.2 4.03 cd
## 14 Muestreo 4 NaCl 12.4 1.08 cd
## 15 Muestreo 4 NaCl+ KCL 1 11.0 1.73 cd
## 16 Muestreo 4 Control 9.31 1.09 dGrafico presentacion articulo
# gray-scale barplot
ggplot(data_summary, aes(x = Muestreo, y = mean, fill =Tratamientos)) +
geom_bar(stat = "identity", position = "dodge", alpha = 0.5, colour = "gray25") +
geom_errorbar(aes(ymin=mean-sd, ymax=mean+sd), position = position_dodge(0.9), width = 0.25,
show.legend = FALSE, colour = "gray25") +
labs(x="Muestreos", y="Estomas Abiertos") +
theme_bw() +
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank()) +
theme(legend.position = c(0.11, 0.75)) +
geom_text(aes(label=c("a")), position = position_dodge(0.90), size = 3,
vjust=-0.8, hjust=-0.5, colour = "gray25") +
ylim(0, 100) +
scale_fill_grey()
###grafico a color
# coloured barplot
ggplot(data_summary, aes(x = Muestreo, y = mean, fill =Tratamientos, colour = Tratamientos)) +
geom_bar(stat = "identity", position = "dodge", alpha = 0.5) +
geom_errorbar(aes(ymin=mean-sd, ymax=mean+sd), position = position_dodge(0.9), width = 0.25) +
labs(x="Tratamientos", y="CRA") +
theme_bw() +
theme(panel.grid.major = element_blank(), panel.grid.minor = element_blank()) +
theme(legend.position = c(0.1, 0.75)) +
geom_text(aes(label=Tukey), position = position_dodge(0.90), size = 3,
vjust=-0.8, hjust=-0.5, colour = "gray25") +
ylim(0, 100) +
geom_text(aes(label=Tratamientos, y = 5), position = position_dodge(1))