DADOS
dados=read.csv("C:/Users/Samsung/Documents/Doutorado CENA_USP/Aulas e disciplinas/Metodos agronomicos/Pratica 2/data.csv")
str(dados)
## 'data.frame': 6 obs. of 10 variables:
## $ amostra : chr "A1" "A2" "A3" "A1" ...
## $ grupo : int 1 1 1 2 2 2
## $ MF.fresca : num 206 175 232 142 131 ...
## $ MF.ms : num 82 71.5 96.9 60.1 54.5 ...
## $ FRACIO.fresca: num 60.7 53.1 60 43.2 35.9 ...
## $ RSA : num 19.3 19.8 19.4 15.1 12.5 ...
## $ LFSA : num 3.47 2.6 3.43 2.58 2.27 2.1
## $ RA.fresca : num 144.4 121.1 170.7 98.9 95.1 ...
## $ RA.ms : num 57.4 49.5 71.3 41.7 39.6 32.4
## $ IAF.LI : num 394 305 376 242 210 ...
dados
C - Área específica de folhas (AEF, cm2 g-1)
dados$AEF=dados$IAF.LI/dados$LFSA
dados
D - Proporção de folhas na sub-amostra (FSA)
dados$FSA=dados$LFSA/(dados$LFSA+dados$RSA)
dados
E - Peso seco da fração de folha na amostra (FA, g
MS)
dados$FA=(dados$LFSA+dados$RSA+dados$RA.ms)*dados$FSA
dados
F - Área total de folhas na amostra (AFA, cm2)
dados$AFA=(dados$FA*dados$AEF)
dados
F - Área total de folhas na amostra (AFA, cm2)
dados$IAF=(dados$AFA/2500)
dados
Média
mean1=mean(dados[-c(4:6),]$IAF)
mean1
## [1] 0.5223502
mean2=mean(dados[-c(1:3),]$IAF)
mean2
## [1] 0.2993721
meangeral=mean(dados$IAF)
meangeral
## [1] 0.4108611
Desvio padrão
dp1=sd(dados[-c(4:6),]$IAF)
dp1
## [1] 0.1176421
dp2=sd(dados[-c(1:3),]$IAF)
dp2
## [1] 0.03208892
dpgeral=sd(dados$IAF)
dpgeral
## [1] 0.1444421
Coeficiente de variação
cv1=sd(dados[-c(4:6),]$IAF)/mean(dados[-c(4:6),]$IAF)*100
cv1
## [1] 22.52169
cv2=sd(dados[-c(1:3),]$IAF)/mean(dados[-c(1:3),]$IAF)*100
cv2
## [1] 10.71874
cvgeral=sd(dados$IAF)/mean(dados$IAF)*100
cvgeral
## [1] 35.15594
Anova e Tukey
library(emmeans)
mod = lm(IAF~grupo, data = dados)
anova(mod)
medias=emmeans(mod,~grupo)
pairs(medias)
## contrast estimate SE df t.ratio p.value
## 1 - 2 0.223 0.0704 4 3.167 0.0339
pwpp(medias)
multcomp::cld(medias)