Analisis de datos - Experimento 1

##    Genotype              CC                 Rep             ma       
##  Length:48          Length:48          Min.   :1.00   Min.   :23.40  
##  Class :character   Class :character   1st Qu.:1.75   1st Qu.:29.68  
##  Mode  :character   Mode  :character   Median :2.50   Median :31.73  
##                                        Mean   :2.50   Mean   :31.77  
##                                        3rd Qu.:3.25   3rd Qu.:34.42  
##                                        Max.   :4.00   Max.   :37.63  
##        mb           mclorofila         ma2             mb2        
##  Min.   : 6.030   Min.   :33.40   Min.   :19.60   Min.   : 6.300  
##  1st Qu.: 8.957   1st Qu.:37.62   1st Qu.:28.92   1st Qu.: 9.335  
##  Median :10.415   Median :42.42   Median :30.90   Median :10.885  
##  Mean   :10.602   Mean   :42.81   Mean   :30.97   Mean   :11.020  
##  3rd Qu.:12.230   3rd Qu.:46.55   3rd Qu.:32.69   3rd Qu.:11.848  
##  Max.   :15.600   Max.   :76.90   Max.   :49.30   Max.   :34.000  
##   mclorofila2         Afol            nfolha             Ra        
##  Min.   :27.30   Min.   : 20.88   Min.   : 5.000   Min.   :0.0800  
##  1st Qu.:38.27   1st Qu.: 62.84   1st Qu.: 7.000   1st Qu.:0.1700  
##  Median :41.12   Median : 90.72   Median : 8.000   Median :0.2350  
##  Mean   :40.91   Mean   : 91.01   Mean   : 8.979   Mean   :0.2405  
##  3rd Qu.:43.58   3rd Qu.:119.79   3rd Qu.:11.250   3rd Qu.:0.3000  
##  Max.   :49.90   Max.   :169.57   Max.   :16.000   Max.   :0.5500  
##       Folh             Tal              Flo               Mpv       
##  Min.   :0.1400   Min.   :0.0700   Min.   :0.00000   Min.   :504.0  
##  1st Qu.:0.2400   1st Qu.:0.1875   1st Qu.:0.04200   1st Qu.:520.8  
##  Median :0.3400   Median :0.2500   Median :0.07000   Median :534.5  
##  Mean   :0.3473   Mean   :0.2669   Mean   :0.07187   Mean   :541.4  
##  3rd Qu.:0.4400   3rd Qu.:0.3600   3rd Qu.:0.09000   3rd Qu.:555.2  
##  Max.   :0.6500   Max.   :0.4800   Max.   :0.36000   Max.   :621.0  
##      Mpinv         viabpolen     
##  Min.   : 0.00   Min.   : 95.45  
##  1st Qu.:12.00   1st Qu.: 96.49  
##  Median :14.50   Median : 97.31  
##  Mean   :14.73   Mean   : 97.28  
##  3rd Qu.:20.00   3rd Qu.: 97.70  
##  Max.   :24.00   Max.   :100.00

Dados

#install.packages("agricolae")
#install.packages("stats")
#install.packages("tidyverse")
#library(tidyverse)
#library(stats)
library(agricolae)
#install.packages(c("PerformanceAnalytics","plotly","reshape2","ggrepel","rgl",
                #  "dplyr", "car","sf","tmap","rgdal","esquisse"))
#install.packages("dplyr")
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, union
library(PerformanceAnalytics)
## Loading required package: xts
## Loading required package: zoo
## 
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric
## 
## Attaching package: 'xts'
## The following objects are masked from 'package:dplyr':
## 
##     first, last
## 
## Attaching package: 'PerformanceAnalytics'
## The following objects are masked from 'package:agricolae':
## 
##     kurtosis, skewness
## The following object is masked from 'package:graphics':
## 
##     legend
library(plotly)
## Loading required package: ggplot2
## 
## Attaching package: 'plotly'
## The following object is masked from 'package:ggplot2':
## 
##     last_plot
## The following object is masked from 'package:stats':
## 
##     filter
## The following object is masked from 'package:graphics':
## 
##     layout
library(tidyverse)
## -- Attaching packages --------------------------------- tidyverse 1.3.0 --
## v tibble  3.0.3     v purrr   0.3.4
## v tidyr   1.1.1     v stringr 1.4.0
## v readr   1.3.1     v forcats 0.5.0
## -- Conflicts ------------------------------------ tidyverse_conflicts() --
## x plotly::filter() masks dplyr::filter(), stats::filter()
## x xts::first()     masks dplyr::first()
## x dplyr::lag()     masks stats::lag()
## x xts::last()      masks dplyr::last()
library(reshape2)
## 
## Attaching package: 'reshape2'
## The following object is masked from 'package:tidyr':
## 
##     smiths
library(ggrepel)
library(rgl)
library(car)
## Loading required package: carData
## 
## Attaching package: 'car'
## The following object is masked from 'package:purrr':
## 
##     some
## The following object is masked from 'package:dplyr':
## 
##     recode
library(sf)
## Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1
library(tmap)
library(rgdal)
## Loading required package: sp
## rgdal: version: 1.5-16, (SVN revision 1050)
## Geospatial Data Abstraction Library extensions to R successfully loaded
## Loaded GDAL runtime: GDAL 3.0.4, released 2020/01/28
## Path to GDAL shared files: C:/Users/Alex/Documents/R/win-library/4.0/rgdal/gdal
## GDAL binary built with GEOS: TRUE 
## Loaded PROJ runtime: Rel. 6.3.1, February 10th, 2020, [PJ_VERSION: 631]
## Path to PROJ shared files: C:/Users/Alex/Documents/R/win-library/4.0/rgdal/proj
## Linking to sp version:1.4-2
## To mute warnings of possible GDAL/OSR exportToProj4() degradation,
## use options("rgdal_show_exportToProj4_warnings"="none") before loading rgdal.
library(esquisse)
library(ggplot2)
library(gridExtra)
## 
## Attaching package: 'gridExtra'
## The following object is masked from 'package:dplyr':
## 
##     combine
dt<-data.frame(lapply(dt, type.convert), stringsAsFactors=FALSE)
summary(aov(dt[,17]~ Genotype*CC, data = dt))
##             Df Sum Sq Mean Sq F value Pr(>F)  
## Genotype     3  140.1   46.69   1.465 0.2403  
## CC           2   20.3   10.15   0.318 0.7294  
## Genotype:CC  6  393.9   65.65   2.060 0.0827 .
## Residuals   36 1147.2   31.87                 
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
Genotype<-as.factor(Genotype)
CC<-as.factor(CC)
Inter<- as.factor(Genotype:CC)

Loops para ANOVA e normalidade

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

## Warning in res[i] <- sort(unlist(aov(dt[, i + 3] ~ Genotype + CC + Inter, :
## number of items to replace is not a multiple of replacement length

visualizar datos

print(cbind.data.frame(names(dt)[4:18], results<0.05))
##    names(dt)[4:18]     1     2     3
## 1               ma FALSE FALSE FALSE
## 2               mb FALSE FALSE FALSE
## 3       mclorofila FALSE FALSE FALSE
## 4              ma2 FALSE FALSE FALSE
## 5              mb2 FALSE FALSE  TRUE
## 6      mclorofila2 FALSE  TRUE  TRUE
## 7             Afol  TRUE  TRUE  TRUE
## 8           nfolha  TRUE  TRUE FALSE
## 9               Ra FALSE  TRUE FALSE
## 10            Folh FALSE  TRUE  TRUE
## 11             Tal FALSE  TRUE  TRUE
## 12             Flo FALSE  TRUE FALSE
## 13             Mpv FALSE FALSE  TRUE
## 14           Mpinv FALSE FALSE FALSE
## 15       viabpolen FALSE FALSE FALSE

Significancia e proba de normalidade

test_norm=cbind.data.frame(names(dt)[4:18], results, shapt, t_kruskal)
names(test_norm)=c("Variavel","Vetor","Riego","Vetor*Riego","Shapiro", "Kolmogorov")
test_norm
##       Variavel Vetor Riego Vetor*Riego Shapiro Kolmogorov
## 1           ma 0.996 0.300       0.359    0.04       0.44
## 2           mb 0.596 0.325       0.146    0.04       0.17
## 3   mclorofila 0.446 0.183       0.933    0.00       0.62
## 4          ma2 0.292 0.739       0.050    0.00       0.04
## 5          mb2 0.146 0.490       0.036    0.00       0.01
## 6  mclorofila2 0.121 0.026       0.025    0.07       0.03
## 7         Afol 0.014 0.000       0.001    0.73       0.00
## 8       nfolha 0.000 0.000       0.537    0.13       0.00
## 9           Ra 0.877 0.001       0.477    0.00       0.02
## 10        Folh 0.333 0.000       0.001    0.48       0.00
## 11         Tal 0.688 0.000       0.000    0.73       0.00
## 12         Flo 0.179 0.002       0.519    0.00       0.00
## 13         Mpv 0.407 0.988       0.045    0.79       0.08
## 14       Mpinv 0.240 0.729       0.083    0.38       0.30
## 15   viabpolen 0.264 0.742       0.078    0.68       0.29

Comparación de medias

t_Inter
## [[1]]
##              dt[, i + 3] groups
## WT:100%CC        34.9075      a
## WT-V:30%CC       32.9250     ab
## PR1-F:100%CC     32.9175     ab
## PR1-G:100%CC     32.1675     ab
## PR1-G:30%CC      32.0325     ab
## WT-V:60%CC       31.8875     ab
## PR1-G:60%CC      31.5925     ab
## PR1-F:30%CC      31.3175     ab
## PR1-F:60%CC      31.0150     ab
## WT:30%CC         30.8875     ab
## WT-V:100%CC      30.5650     ab
## WT:60%CC         29.0675      b
## 
## [[2]]
##              dt[, i + 3] groups
## WT:100%CC        12.5225      a
## PR1-F:30%CC      12.1325      a
## WT-V:30%CC       11.8225     ab
## PR1-G:30%CC      11.5150    abc
## PR1-F:100%CC     11.1075    abc
## WT-V:60%CC       10.5175    abc
## PR1-G:60%CC      10.3675    abc
## PR1-G:100%CC     10.3150    abc
## PR1-F:60%CC      10.2925    abc
## WT-V:100%CC       9.3500    abc
## WT:30%CC          8.7750     bc
## WT:60%CC          8.5025      c
## 
## [[3]]
##              dt[, i + 3] groups
## WT-V:100%CC      49.1025      a
## WT:100%CC        46.2675     ab
## WT-V:30%CC       44.7525     ab
## PR1-F:100%CC     44.0250     ab
## PR1-G:30%CC      43.5475     ab
## PR1-G:100%CC     43.5250     ab
## WT-V:60%CC       42.8975     ab
## PR1-F:60%CC      41.9750     ab
## PR1-G:60%CC      41.9600     ab
## PR1-F:30%CC      40.1475     ab
## WT:30%CC         38.0000      b
## WT:60%CC         37.5650      b
## 
## [[4]]
##              dt[, i + 3] groups
## PR1-F:60%CC      36.4250      a
## WT:100%CC        33.5675     ab
## WT-V:30%CC       32.9250     ab
## PR1-G:100%CC     32.0075    abc
## WT-V:60%CC       30.7400     bc
## WT:30%CC         30.5400     bc
## PR1-F:100%CC     30.3925     bc
## WT-V:100%CC      30.3500     bc
## PR1-F:30%CC      30.2075     bc
## PR1-G:30%CC      28.8500     bc
## WT:60%CC         28.5175     bc
## PR1-G:60%CC      27.1425      c
## 
## [[5]]
##              dt[, i + 3] groups
## PR1-F:60%CC      17.6750      a
## WT:100%CC        13.3900     ab
## WT-V:60%CC       11.1650     bc
## WT-V:100%CC      11.1250     bc
## WT:30%CC         10.6650     bc
## PR1-F:100%CC     10.6475     bc
## PR1-G:100%CC     10.4925     bc
## WT-V:30%CC       10.3450     bc
## PR1-F:30%CC       9.8575     bc
## PR1-G:30%CC       9.8175     bc
## WT:60%CC          9.2750     bc
## PR1-G:60%CC       7.7900      c
## 
## [[6]]
##              dt[, i + 3] groups
## WT:100%CC        46.9575      a
## PR1-F:60%CC      43.9600     ab
## PR1-G:100%CC     42.5000    abc
## WT-V:60%CC       41.9075    abc
## WT-V:100%CC      41.4725     bc
## WT:30%CC         41.2100     bc
## PR1-F:100%CC     41.0400     bc
## WT-V:30%CC       40.4600     bc
## PR1-F:30%CC      40.0675    bcd
## PR1-G:30%CC      38.6675     cd
## WT:60%CC         37.7900     cd
## PR1-G:60%CC      34.9350      d
## 
## [[7]]
##              dt[, i + 3] groups
## WT-V:100%CC     145.3975      a
## PR1-F:100%CC    132.6700     ab
## PR1-G:100%CC    132.4575     ab
## PR1-F:60%CC     112.1725     bc
## WT:100%CC       105.3200     cd
## WT:60%CC         90.7325     cd
## PR1-G:60%CC      88.6750      d
## PR1-G:30%CC      64.9425      e
## WT-V:30%CC       61.2075      e
## PR1-F:30%CC      58.6350      e
## WT-V:60%CC       55.1250      e
## WT:30%CC         44.7675      e
## 
## [[8]]
##              dt[, i + 3] groups
## PR1-F:100%CC       13.50      a
## WT-V:100%CC        12.50     ab
## PR1-F:60%CC        11.00    abc
## PR1-G:100%CC       10.50     bc
## PR1-F:30%CC         9.25     cd
## WT-V:30%CC          8.50    cde
## WT-V:60%CC          8.50    cde
## PR1-G:60%CC         7.25     de
## WT:100%CC           7.25     de
## PR1-G:30%CC         6.75     de
## WT:60%CC            6.50     de
## WT:30%CC            6.25      e
## 
## [[9]]
##              dt[, i + 3] groups
## WT-V:100%CC      0.33750      a
## WT:100%CC        0.32000      a
## PR1-G:100%CC     0.28500     ab
## PR1-F:60%CC      0.27075    abc
## PR1-F:100%CC     0.26500   abcd
## WT:60%CC         0.26250   abcd
## PR1-G:60%CC      0.25225   abcd
## PR1-F:30%CC      0.22825   abcd
## PR1-G:30%CC      0.18700    bcd
## WT-V:60%CC       0.18450    bcd
## WT-V:30%CC       0.15250     cd
## WT:30%CC         0.14125      d
## 
## [[10]]
##              dt[, i + 3] groups
## WT-V:100%CC       0.5550      a
## PR1-G:100%CC      0.4925     ab
## PR1-F:100%CC      0.4525    abc
## PR1-F:60%CC       0.4175    bcd
## WT:60%CC          0.3875    bcd
## PR1-G:60%CC       0.3625    cde
## WT:100%CC         0.3325    def
## PR1-G:30%CC       0.2600    efg
## WT-V:30%CC        0.2400     fg
## WT-V:60%CC        0.2325     fg
## WT:30%CC          0.2300     fg
## PR1-F:30%CC       0.2050      g
## 
## [[11]]
##              dt[, i + 3] groups
## WT-V:100%CC       0.4175      a
## PR1-G:100%CC      0.3925      a
## PR1-F:60%CC       0.3475     ab
## WT:60%CC          0.3325     ab
## PR1-F:100%CC      0.2950      b
## PR1-G:60%CC       0.2750     bc
## WT:100%CC         0.2725     bc
## WT-V:30%CC        0.1975     cd
## PR1-G:30%CC       0.1900     cd
## WT:30%CC          0.1850     cd
## WT-V:60%CC        0.1725      d
## PR1-F:30%CC       0.1250      d
## 
## [[12]]
##              dt[, i + 3] groups
## PR1-F:100%CC     0.14275      a
## PR1-G:100%CC     0.12500     ab
## PR1-F:60%CC      0.09675    abc
## WT-V:100%CC      0.08450    abc
## WT:60%CC         0.07750    abc
## PR1-G:60%CC      0.07050     bc
## WT:100%CC        0.06500     bc
## PR1-G:30%CC      0.05000      c
## WT:30%CC         0.03875      c
## WT-V:30%CC       0.03850      c
## WT-V:60%CC       0.03800      c
## PR1-F:30%CC      0.03525      c
## 
## [[13]]
##              dt[, i + 3] groups
## WT-V:30%CC        566.50      a
## PR1-F:60%CC       561.50     ab
## WT:100%CC         556.50    abc
## WT-V:60%CC        551.25    abc
## PR1-G:30%CC       550.75    abc
## WT:60%CC          539.25   abcd
## PR1-F:100%CC      538.50   abcd
## WT-V:100%CC       537.00   abcd
## PR1-G:100%CC      536.50   abcd
## WT:30%CC          525.75    bcd
## PR1-F:30%CC       519.75     cd
## PR1-G:60%CC       513.00      d
## 
## [[14]]
##              dt[, i + 3] groups
## PR1-G:100%CC       18.75      a
## WT-V:30%CC         18.25      a
## WT-V:60%CC         17.75      a
## WT:60%CC           17.50      a
## PR1-G:30%CC        17.00      a
## PR1-F:100%CC       16.50      a
## WT:30%CC           13.75     ab
## PR1-G:60%CC        13.50     ab
## WT-V:100%CC        13.25     ab
## PR1-F:60%CC        12.25     ab
## PR1-F:30%CC        11.50     ab
## WT:100%CC           6.75      b
## 
## [[15]]
##              dt[, i + 3] groups
## WT:100%CC       98.79175      a
## PR1-F:60%CC     97.82937     ab
## PR1-F:30%CC     97.80828     ab
## WT-V:100%CC     97.56636     ab
## WT:30%CC        97.38299     ab
## PR1-G:60%CC     97.37673     ab
## PR1-F:100%CC    96.92815      b
## PR1-G:30%CC     96.91871      b
## WT-V:60%CC      96.82122      b
## WT:60%CC        96.75317      b
## WT-V:30%CC      96.73132      b
## PR1-G:100%CC    96.50449      b
t_CC
## [[1]]
##        dt[, i + 3] groups
## 100%CC    32.63938      a
## 30%CC     31.79062      a
## 60%CC     30.89062      a
## 
## [[2]]
##        dt[, i + 3] groups
## 30%CC     11.06125      a
## 100%CC    10.82375      a
## 60%CC      9.92000      a
## 
## [[3]]
##        dt[, i + 3] groups
## 100%CC    45.73000      a
## 30%CC     41.61187      a
## 60%CC     41.09938      a
## 
## [[4]]
##        dt[, i + 3] groups
## 100%CC    31.57937      a
## 60%CC     30.70625      a
## 30%CC     30.63062      a
## 
## [[5]]
##        dt[, i + 3] groups
## 60%CC     11.47625      a
## 100%CC    11.41375      a
## 30%CC     10.17125      a
## 
## [[6]]
##        dt[, i + 3] groups
## 100%CC    42.99250      a
## 30%CC     40.10125      b
## 60%CC     39.64813      b
## 
## [[7]]
##        dt[, i + 3] groups
## 100%CC   128.96125      a
## 60%CC     86.67625      b
## 30%CC     57.38813      c
## 
## [[8]]
##        dt[, i + 3] groups
## 100%CC     10.9375      a
## 60%CC       8.3125      b
## 30%CC       7.6875      b
## 
## [[9]]
##        dt[, i + 3] groups
## 100%CC    0.301875      a
## 60%CC     0.242500      a
## 30%CC     0.177250      b
## 
## [[10]]
##        dt[, i + 3] groups
## 100%CC    0.458125      a
## 60%CC     0.350000      b
## 30%CC     0.233750      c
## 
## [[11]]
##        dt[, i + 3] groups
## 100%CC    0.344375      a
## 60%CC     0.281875      b
## 30%CC     0.174375      c
## 
## [[12]]
##        dt[, i + 3] groups
## 100%CC   0.1043125      a
## 60%CC    0.0706875     ab
## 30%CC    0.0406250      b
## 
## [[13]]
##        dt[, i + 3] groups
## 100%CC    542.1250      a
## 60%CC     541.2500      a
## 30%CC     540.6875      a
## 
## [[14]]
##        dt[, i + 3] groups
## 60%CC      15.2500      a
## 30%CC      15.1250      a
## 100%CC     13.8125      a
## 
## [[15]]
##        dt[, i + 3] groups
## 100%CC    97.44769      a
## 30%CC     97.21032      a
## 60%CC     97.19512      a
t_Genotype
## [[1]]
##       dt[, i + 3] groups
## PR1-G    31.93083      a
## WT-V     31.79250      a
## PR1-F    31.75000      a
## WT       31.62083      a
## 
## [[2]]
##       dt[, i + 3] groups
## PR1-F   11.177500      a
## PR1-G   10.732500      a
## WT-V    10.563333      a
## WT       9.933333      a
## 
## [[3]]
##       dt[, i + 3] groups
## WT-V     45.58417      a
## PR1-G    43.01083      a
## PR1-F    42.04917      a
## WT       40.61083      a
## 
## [[4]]
##       dt[, i + 3] groups
## PR1-F    32.34167      a
## WT-V     31.33833      a
## WT       30.87500      a
## PR1-G    29.33333      a
## 
## [[5]]
##       dt[, i + 3] groups
## PR1-F   12.726667      a
## WT      11.110000     ab
## WT-V    10.878333     ab
## PR1-G    9.366667      b
## 
## [[6]]
##       dt[, i + 3] groups
## WT       41.98583      a
## PR1-F    41.68917     ab
## WT-V     41.28000     ab
## PR1-G    38.70083      b
## 
## [[7]]
##       dt[, i + 3] groups
## PR1-F   101.15917      a
## PR1-G    95.35833     ab
## WT-V     87.24333     bc
## WT       80.27333      c
## 
## [[8]]
##       dt[, i + 3] groups
## PR1-F   11.250000      a
## WT-V     9.833333      a
## PR1-G    8.166667      b
## WT       6.666667      b
## 
## [[9]]
##       dt[, i + 3] groups
## PR1-F   0.2546667      a
## PR1-G   0.2414167      a
## WT      0.2412500      a
## WT-V    0.2248333      a
## 
## [[10]]
##       dt[, i + 3] groups
## PR1-G   0.3716667      a
## PR1-F   0.3583333      a
## WT-V    0.3425000      a
## WT      0.3166667      a
## 
## [[11]]
##       dt[, i + 3] groups
## PR1-G   0.2858333      a
## WT      0.2633333      a
## WT-V    0.2625000      a
## PR1-F   0.2558333      a
## 
## [[12]]
##       dt[, i + 3] groups
## PR1-F  0.09158333      a
## PR1-G  0.08183333      a
## WT     0.06041667      a
## WT-V   0.05366667      a
## 
## [[13]]
##       dt[, i + 3] groups
## WT-V     551.5833      a
## WT       540.5000      a
## PR1-F    539.9167      a
## PR1-G    533.4167      a
## 
## [[14]]
##       dt[, i + 3] groups
## PR1-G    16.41667      a
## WT-V     16.41667      a
## PR1-F    13.41667      a
## WT       12.66667      a
## 
## [[15]]
##       dt[, i + 3] groups
## WT       97.64264      a
## PR1-F    97.52193      a
## WT-V     97.03963      a
## PR1-G    96.93331      a

Teste de Kruskal-Wallis

t_krus=cbind.data.frame(names(dt)[4:18], t_kruskal)
names(test_norm)=c("Variavel","Test Kruskal-Wallis ")
t_krus
##    names(dt)[4:18] t_kruskal
## 1               ma      0.44
## 2               mb      0.17
## 3       mclorofila      0.62
## 4              ma2      0.04
## 5              mb2      0.01
## 6      mclorofila2      0.03
## 7             Afol      0.00
## 8           nfolha      0.00
## 9               Ra      0.02
## 10            Folh      0.00
## 11             Tal      0.00
## 12             Flo      0.00
## 13             Mpv      0.08
## 14           Mpinv      0.30
## 15       viabpolen      0.29
med_kruskal
## [[1]]
##              dt[, i + 3] groups
## WT:100%CC         37.875      a
## WT-V:30%CC        30.625      a
## PR1-F:100%CC      29.375     ab
## PR1-G:100%CC      27.250     ab
## PR1-F:30%CC       26.500     ab
## PR1-G:30%CC       25.750     ab
## WT-V:60%CC        25.250     ab
## PR1-G:60%CC       22.875     ab
## PR1-F:60%CC       20.250     ab
## WT:30%CC          19.500     ab
## WT-V:100%CC       18.750     ab
## WT:60%CC          10.000      b
## 
## [[2]]
##              dt[, i + 3] groups
## PR1-F:30%CC       36.000      a
## WT:100%CC         35.250      a
## WT-V:30%CC        31.750      a
## PR1-G:30%CC       30.750     ab
## PR1-F:100%CC      27.125    abc
## PR1-G:60%CC       24.250    abc
## WT-V:60%CC        23.250    abc
## PR1-G:100%CC      23.125    abc
## PR1-F:60%CC       22.000    abc
## WT-V:100%CC       17.750    abc
## WT:30%CC          12.750     bc
## WT:60%CC          10.000      c
## 
## [[3]]
##              dt[, i + 3] groups
## WT:100%CC         33.000      a
## WT-V:30%CC        31.500      a
## PR1-G:100%CC      27.875      a
## PR1-F:100%CC      27.500      a
## PR1-G:30%CC       27.250      a
## WT-V:60%CC        26.500      a
## WT-V:100%CC       26.000      a
## PR1-F:60%CC       24.000      a
## PR1-G:60%CC       24.000      a
## PR1-F:30%CC       20.625      a
## WT:30%CC          13.250      a
## WT:60%CC          12.500      a
## 
## [[4]]
##              dt[, i + 3] groups
## WT:100%CC         39.125      a
## PR1-F:60%CC       37.000      a
## PR1-G:100%CC      31.250      a
## WT-V:30%CC        28.375     ab
## WT-V:60%CC        25.875     ab
## PR1-G:30%CC       24.625     ab
## WT:30%CC          23.625     ab
## WT-V:100%CC       23.125     ab
## PR1-F:100%CC      22.750     ab
## PR1-F:30%CC       22.000     ab
## WT:60%CC          11.625     bc
## PR1-G:60%CC        4.625      c
## 
## [[5]]
##              dt[, i + 3] groups
## WT:100%CC         43.500      a
## PR1-F:60%CC       41.500     ab
## WT-V:60%CC        29.250    abc
## WT-V:100%CC       28.000    bcd
## WT:30%CC          25.125     cd
## PR1-F:100%CC      24.750     cd
## PR1-G:100%CC      24.500     cd
## WT-V:30%CC        22.625     cd
## PR1-G:30%CC       18.500    cde
## PR1-F:30%CC       17.250    cde
## WT:60%CC          13.750     de
## PR1-G:60%CC        5.250      e
## 
## [[6]]
##              dt[, i + 3] groups
## WT:100%CC          43.00      a
## PR1-F:60%CC        35.25     ab
## PR1-G:100%CC       29.00    abc
## WT-V:60%CC         28.00    abc
## WT-V:100%CC        26.75    abc
## WT:30%CC           25.00     bc
## PR1-F:100%CC       24.50     bc
## WT-V:30%CC         22.75     bc
## PR1-G:30%CC        22.25     bc
## PR1-F:30%CC        20.25    bcd
## WT:60%CC           13.25     cd
## PR1-G:60%CC         4.00      d
## 
## [[7]]
##              dt[, i + 3] groups
## WT-V:100%CC        43.00      a
## PR1-G:100%CC       41.00     ab
## PR1-F:100%CC       40.50     ab
## PR1-F:60%CC        33.75     bc
## WT:100%CC          31.75     cd
## WT:60%CC           25.25      d
## PR1-G:60%CC        25.00      d
## PR1-G:30%CC        14.25      e
## WT-V:30%CC         12.25      e
## PR1-F:30%CC        12.00      e
## WT-V:60%CC          9.25      e
## WT:30%CC            6.00      e
## 
## [[8]]
##              dt[, i + 3] groups
## PR1-F:100%CC      42.875      a
## WT-V:100%CC       40.000     ab
## PR1-F:60%CC       35.625    abc
## PR1-G:100%CC      33.500   abcd
## PR1-F:30%CC       28.000   bcde
## WT-V:60%CC        26.500    cde
## WT-V:30%CC        22.250    def
## PR1-G:60%CC       16.750    efg
## WT:100%CC         16.250    efg
## PR1-G:30%CC       13.125     fg
## WT:60%CC          10.750     fg
## WT:30%CC           8.375      g
## 
## [[9]]
##              dt[, i + 3] groups
## WT-V:100%CC       39.625      a
## WT:100%CC         34.125      a
## PR1-G:100%CC      33.125     ab
## WT:60%CC          30.375    abc
## PR1-F:100%CC      30.250    abc
## PR1-F:60%CC       29.750    abc
## PR1-G:60%CC       27.625    abc
## PR1-G:30%CC       17.250    bcd
## PR1-F:30%CC       16.875    bcd
## WT-V:60%CC        15.625     cd
## WT-V:30%CC        10.625      d
## WT:30%CC           8.750      d
## 
## [[10]]
##              dt[, i + 3] groups
## WT-V:100%CC       43.750      a
## PR1-G:100%CC      39.875     ab
## PR1-F:100%CC      36.250    abc
## PR1-F:60%CC       33.750   abcd
## WT:60%CC          31.250    bcd
## PR1-G:60%CC       27.500     cd
## WT:100%CC         23.375     de
## PR1-G:30%CC       15.250     ef
## WT-V:30%CC        11.750      f
## WT-V:60%CC        11.625      f
## WT:30%CC          11.000      f
## PR1-F:30%CC        8.625      f
## 
## [[11]]
##              dt[, i + 3] groups
## WT-V:100%CC       44.000      a
## PR1-G:100%CC      40.875     ab
## PR1-F:60%CC       33.750    abc
## WT:60%CC          32.625    abc
## PR1-F:100%CC      28.250    bcd
## WT:100%CC         25.875    cde
## PR1-G:60%CC       25.750    cde
## WT-V:30%CC        15.875    def
## PR1-G:30%CC       14.250     ef
## WT:30%CC          13.625     ef
## WT-V:60%CC        12.375      f
## PR1-F:30%CC        6.750      f
## 
## [[12]]
##              dt[, i + 3] groups
## PR1-G:100%CC      43.500      a
## PR1-F:60%CC       37.375     ab
## WT-V:100%CC       34.125     ab
## PR1-F:100%CC      33.250     ab
## WT:60%CC          31.375     ab
## PR1-G:60%CC       26.000     bc
## WT:100%CC         25.125     bc
## PR1-G:30%CC       17.375     cd
## WT:30%CC          14.000     cd
## WT-V:60%CC        11.000      d
## WT-V:30%CC        10.875      d
## PR1-F:30%CC       10.000      d
## 
## [[13]]
##              dt[, i + 3] groups
## WT-V:30%CC        37.125      a
## PR1-F:60%CC       33.500     ab
## WT:100%CC         30.750     ab
## PR1-G:30%CC       30.000    abc
## WT-V:60%CC        28.375    abc
## PR1-F:100%CC      27.000    abc
## WT:60%CC          26.000    abc
## WT-V:100%CC       25.250    abc
## PR1-G:100%CC      20.500   abcd
## WT:30%CC          17.000    bcd
## PR1-F:30%CC       12.375     cd
## PR1-G:60%CC        6.125      d
## 
## [[14]]
##              dt[, i + 3] groups
## PR1-G:100%CC      34.000      a
## WT-V:30%CC        32.125      a
## WT-V:60%CC        31.500      a
## WT:60%CC          31.250      a
## PR1-F:100%CC      28.625      a
## PR1-G:30%CC       28.500      a
## WT-V:100%CC       20.875     ab
## PR1-F:60%CC       19.875     ab
## WT:30%CC          19.750     ab
## PR1-G:60%CC       19.375     ab
## PR1-F:30%CC       19.125     ab
## WT:100%CC          9.000      b
## 
## [[15]]
##              dt[, i + 3] groups
## WT:100%CC          40.00      a
## PR1-F:60%CC        31.25     ab
## PR1-F:30%CC        29.25     ab
## WT-V:100%CC        28.50     ab
## PR1-G:60%CC        28.25     ab
## WT:30%CC           27.75     ab
## PR1-G:30%CC        22.50     ab
## PR1-F:100%CC       20.75     ab
## WT-V:60%CC         18.25      b
## WT-V:30%CC         17.50      b
## WT:60%CC           16.25      b
## PR1-G:100%CC       13.75      b

Correlaciones

##  [1] "Genotype"    "CC"          "Rep"         "ma"          "mb"         
##  [6] "mclorofila"  "ma2"         "mb2"         "mclorofila2" "Afol"       
## [11] "nfolha"      "Ra"          "Folh"        "Tal"         "Flo"        
## [16] "Mpv"         "Mpinv"       "viabpolen"

##           Var1 Var2 value
## 1         SPAD SPAD  1.00
## 2    Leaf area SPAD  0.19
## 3        Leafs SPAD  0.15
## 4 Root biomass SPAD  0.21
## 5 Leaf biomass SPAD  0.09
## 6 Stem biomass SPAD  0.10

Area folhar

##  [1] WT-V:100%CC  WT-V:100%CC  WT-V:100%CC  WT-V:100%CC  WT-V:60%CC  
##  [6] WT-V:60%CC   WT-V:60%CC   WT-V:60%CC   WT-V:30%CC   WT-V:30%CC  
## [11] WT-V:30%CC   WT-V:30%CC   WT:100%CC    WT:100%CC    WT:100%CC   
## [16] WT:100%CC    WT:60%CC     WT:60%CC     WT:60%CC     WT:60%CC    
## [21] WT:30%CC     WT:30%CC     WT:30%CC     WT:30%CC     PR1-F:100%CC
## [26] PR1-F:100%CC PR1-F:100%CC PR1-F:100%CC PR1-F:60%CC  PR1-F:60%CC 
## [31] PR1-F:60%CC  PR1-F:60%CC  PR1-F:30%CC  PR1-F:30%CC  PR1-F:30%CC 
## [36] PR1-F:30%CC  PR1-G:100%CC PR1-G:100%CC PR1-G:100%CC PR1-G:100%CC
## [41] PR1-G:60%CC  PR1-G:60%CC  PR1-G:60%CC  PR1-G:60%CC  PR1-G:30%CC 
## [46] PR1-G:30%CC  PR1-G:30%CC  PR1-G:30%CC 
## 12 Levels: WT:100%CC WT:60%CC WT:30%CC WT-V:100%CC WT-V:60%CC ... PR1-G:30%CC
## The following objects are masked _by_ .GlobalEnv:
## 
##     CC, Genotype, viabpolen
## The following objects are masked from dt (pos = 28):
## 
##     Afol, CC, Flo, Folh, Genotype, ma, ma2, mb, mb2, mclorofila,
##     mclorofila2, Mpinv, Mpv, nfolha, Ra, Rep, Tal
##  [1] "Genotype"    "CC"          "Rep"         "ma"          "mb"         
##  [6] "mclorofila"  "ma2"         "mb2"         "mclorofila2" "Afol"       
## [11] "nfolha"      "Ra"          "Folh"        "Tal"         "Flo"        
## [16] "Mpv"         "Mpinv"       "viabpolen"
##  [1] WT-V:100%CC  WT-V:100%CC  WT-V:100%CC  WT-V:100%CC  WT-V:60%CC  
##  [6] WT-V:60%CC   WT-V:60%CC   WT-V:60%CC   WT-V:30%CC   WT-V:30%CC  
## [11] WT-V:30%CC   WT-V:30%CC   WT:100%CC    WT:100%CC    WT:100%CC   
## [16] WT:100%CC    WT:60%CC     WT:60%CC     WT:60%CC     WT:60%CC    
## [21] WT:30%CC     WT:30%CC     WT:30%CC     WT:30%CC     PR1-F:100%CC
## [26] PR1-F:100%CC PR1-F:100%CC PR1-F:100%CC PR1-F:60%CC  PR1-F:60%CC 
## [31] PR1-F:60%CC  PR1-F:60%CC  PR1-F:30%CC  PR1-F:30%CC  PR1-F:30%CC 
## [36] PR1-F:30%CC  PR1-G:100%CC PR1-G:100%CC PR1-G:100%CC PR1-G:100%CC
## [41] PR1-G:60%CC  PR1-G:60%CC  PR1-G:60%CC  PR1-G:60%CC  PR1-G:30%CC 
## [46] PR1-G:30%CC  PR1-G:30%CC  PR1-G:30%CC 
## 12 Levels: PR1-F:100%CC PR1-F:30%CC PR1-F:60%CC PR1-G:100%CC ... WT-V:60%CC

Biomassa de tallo

##  [1] "Genotype"    "CC"          "Rep"         "ma"          "mb"         
##  [6] "mclorofila"  "ma2"         "mb2"         "mclorofila2" "Afol"       
## [11] "nfolha"      "Ra"          "Folh"        "Tal"         "Flo"        
## [16] "Mpv"         "Mpinv"       "viabpolen"

Numero de folhas

## The following objects are masked _by_ .GlobalEnv:
## 
##     CC, Genotype, viabpolen
## The following objects are masked from dt (pos = 3):
## 
##     Afol, CC, Flo, Folh, Genotype, ma, ma2, mb, mb2, mclorofila,
##     mclorofila2, Mpinv, Mpv, nfolha, Ra, Rep, Tal, viabpolen
## The following objects are masked from dt (pos = 29):
## 
##     Afol, CC, Flo, Folh, Genotype, ma, ma2, mb, mb2, mclorofila,
##     mclorofila2, Mpinv, Mpv, nfolha, Ra, Rep, Tal
##  [1] "Genotyp"     "CC"          "Rep"         "ma"          "mb"         
##  [6] "mclorofila"  "ma2"         "mb2"         "mclorofila2" "Afol"       
## [11] "nfolha"      "Ra"          "Folh"        "Tal"         "Flo"        
## [16] "Mpv"         "Mpinv"       "viabpolen"

Biomassa de raiz

Clororfila

##   [1] "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V" 
##  [10] "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT"    "WT"    "WT"   
##  [19] "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"   
##  [28] "WT"    "WT"    "WT"    "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F"
##  [37] "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F"
##  [46] "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G"
##  [55] "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "WT-V"  "WT-V"  "WT-V"  "WT-V" 
##  [64] "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V" 
##  [73] "WT-V"  "WT-V"  "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"   
##  [82] "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "PR1-F"
##  [91] "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F"
## [100] "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-G" "PR1-G" "PR1-G" "PR1-G"
## [109] "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G"
## [118] "PR1-G" "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V" 
## [127] "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT"    "WT"   
## [136] "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"   
## [145] "WT"    "WT"    "WT"    "WT"    "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F"
## [154] "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F"
## [163] "PR1-F" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G"
## [172] "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "WT-V"  "WT-V"  "WT-V" 
## [181] "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V" 
## [190] "WT-V"  "WT-V"  "WT-V"  "WT"    "WT"    "WT"    "WT"    "WT"    "WT"   
## [199] "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"   
## [208] "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F"
## [217] "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-G" "PR1-G" "PR1-G"
## [226] "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G"
## [235] "PR1-G" "PR1-G" "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V" 
## [244] "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT-V"  "WT"   
## [253] "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"    "WT"   
## [262] "WT"    "WT"    "WT"    "WT"    "WT"    "PR1-F" "PR1-F" "PR1-F" "PR1-F"
## [271] "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F" "PR1-F"
## [280] "PR1-F" "PR1-F" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G"
## [289] "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G" "PR1-G"
## Warning: Using size for a discrete variable is not advised.
## Warning: Removed 1 rows containing missing values (geom_point).

IRGA

## The following objects are masked from dt (pos = 3):
## 
##     ma2, mb2
## The following objects are masked from dt (pos = 29):
## 
##     ma2, mb2
##  [1] "Vaso"           "Genotypei"      "Cci"            "Rept biológica"
##  [5] "E1"             "GS1"            "A1"             "E2"            
##  [9] "GS2"            "A2"             "E3"             "GS3"           
## [13] "A3"             "total E"        "TotalGS"        "TotalA"        
## [17] "ma1"            "mb1"            "mcl1"           "ma2"           
## [21] "mb2"            "mcl2"
##  [1] WT    WT    WT    WT    WT    WT    WT    WT    WT    PR1-F PR1-F PR1-F
## [13] PR1-F PR1-F PR1-F PR1-F PR1-F PR1-F PR1-G PR1-G PR1-G PR1-G PR1-G PR1-G
## [25] PR1-G PR1-G PR1-G
## Levels: WT PR1-F PR1-G
##  [1] "Vaso"           "Genotypei"      "Cci"            "Rept biológica"
##  [5] "E1"             "GS1"            "A1"             "E2"            
##  [9] "GS2"            "A2"             "E3"             "GS3"           
## [13] "A3"             "total E"        "TotalGS"        "TotalA"        
## [17] "ma1"            "mb1"            "mcl1"           "ma2"           
## [21] "mb2"            "mcl2"
## The following objects are masked from dt (pos = 3):
## 
##     ma2, mb2
## 
## The following objects are masked from dt (pos = 29):
## 
##     ma2, mb2
##  [1] "Vaso"           "Genotypei"      "Cci"            "Rept biológica"
##  [5] "E1"             "GS1"            "A1"             "E2"            
##  [9] "GS2"            "A2"             "E3"             "GS3"           
## [13] "A3"             "total E"        "TotalGS"        "TotalA"        
## [17] "ma1"            "mb1"            "mcl1"           "ma2"           
## [21] "mb2"            "mcl2"
##  [1] "Vaso"           "Genotype"       "Soil moinsture" "Rept biológica"
##  [5] "E1"             "GS1"            "A1"             "E2"            
##  [9] "GS2"            "A2"             "E3"             "GS3"           
## [13] "A3"             "total E"        "TotalGS"        "TotalA"        
## [17] "ma1"            "mb1"            "mcl1"           "ma2"           
## [21] "mb2"            "mcl2"