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"
