Teste de normalidade:

knitr::kable(shapiro_geral)
variavel n estatistica p_valor
Al 54 0.9245 0.0022
Areia 54 0.7817 0.0000
Argila 54 0.7581 0.0000
CBM 54 0.7311 0.0000
CTC.a.pH7.cmolc 54 0.9808 0.5330
CTC.efetiva 54 0.8848 0.0001
Ca 54 0.8375 0.0000
Corg 54 0.9571 0.0510
Densidade 54 0.9372 0.0071
EC 54 0.9271 0.0028
H.Al 54 0.9450 0.0152
K 54 0.8264 0.0000
Macroporosidade 54 0.7193 0.0000
Mat.Org 54 0.9583 0.0576
Mg 54 0.8113 0.0000
Microporosidade 54 0.6347 0.0000
P.Mehlich 54 0.9443 0.0141
PT 54 0.7755 0.0000
RBS 54 0.9310 0.0040
SB 54 0.8421 0.0000
Sat.Al.M 54 0.9172 0.0012
Sat.Base 54 0.8899 0.0001
Silte 54 0.6309 0.0000
Tar 54 0.8790 0.0001
Ts 45 0.5884 0.0000
UR 54 0.8165 0.0000
Us 45 0.7892 0.0000
pH.CaCl2 54 0.8504 0.0000

Obs: Grande parte das variáveis não se adequaram à premissa de normalidade.

Analise Descritiva geral

Estatísticas descritivas dos atributos do solo por uso do solo.
Uso.do.solo variavel media desvio_padrao coef_variacao minimo q1 mediana q3 maximo
Abandonada Al 0.66 0.21 32.77 0.40 0.50 0.60 0.70 1.20
Abandonada Areia 787.88 0.00 0.00 787.88 787.88 787.88 787.88 787.88
Abandonada Argila 152.00 0.00 0.00 152.00 152.00 152.00 152.00 152.00
Abandonada CBM 111.06 74.91 67.45 16.40 68.22 111.80 130.90 283.60
Abandonada CTC.a.pH7.cmolc 3.77 0.70 18.43 2.80 3.23 3.60 4.47 4.80
Abandonada CTC.efetiva 1.35 0.27 19.72 1.00 1.20 1.30 1.50 2.00
Abandonada Ca 0.43 0.10 23.75 0.30 0.40 0.40 0.50 0.60
Abandonada Corg 9.41 3.32 35.23 6.20 7.00 7.80 10.93 16.50
Abandonada Densidade 1.47 0.13 8.66 1.20 1.40 1.50 1.60 1.60
Abandonada EC 13.58 4.43 32.59 7.40 11.00 12.25 14.65 24.40
Abandonada H.Al 3.08 0.63 20.54 2.20 2.52 2.90 3.77 4.00
Abandonada K 0.09 0.03 30.24 0.06 0.07 0.08 0.09 0.16
Abandonada Macroporosidade 0.17 0.05 29.10 0.10 0.10 0.20 0.20 0.20
Abandonada Mat.Org 16.21 5.72 35.30 10.70 12.00 13.50 18.85 28.50
Abandonada Mg 0.16 0.05 32.87 0.10 0.10 0.20 0.20 0.20
Abandonada Microporosidade 0.23 0.05 20.79 0.20 0.20 0.20 0.30 0.30
Abandonada P.Mehlich 4.04 1.61 39.69 1.80 2.25 4.00 5.00 7.00
Abandonada PT 0.41 0.08 19.78 0.30 0.32 0.40 0.50 0.50
Abandonada RBS 81.11 45.21 55.74 17.10 45.97 75.75 114.85 155.10
Abandonada SB 0.69 0.14 20.85 0.48 0.59 0.70 0.79 0.92
Abandonada Sat.Al.M 48.57 8.52 17.54 36.60 42.80 46.45 53.65 69.80
Abandonada Sat.Base 18.29 3.50 19.13 10.80 15.67 19.65 20.80 23.10
Abandonada Silte 60.12 0.00 0.00 60.12 60.12 60.12 60.12 60.12
Abandonada Tar 28.53 1.58 5.55 26.99 26.99 28.53 30.07 30.07
Abandonada Ts 32.91 0.00 0.00 32.91 32.91 32.91 32.91 32.91
Abandonada UR 78.97 5.68 7.19 73.45 73.45 78.97 84.49 84.49
Abandonada Us 0.15 0.00 0.00 0.15 0.15 0.15 0.15 0.15
Abandonada pH.CaCl2 4.04 0.11 2.83 3.90 4.00 4.00 4.10 4.40
Floresta Al 0.56 0.11 20.68 0.30 0.50 0.55 0.60 0.80
Floresta Areia 727.88 0.00 0.00 727.88 727.88 727.88 727.88 727.88
Floresta Argila 209.00 0.00 0.00 209.00 209.00 209.00 209.00 209.00
Floresta CBM 131.83 166.32 126.17 5.50 34.08 54.55 188.18 556.40
Floresta CTC.a.pH7.cmolc 4.17 0.55 13.30 3.30 4.00 4.10 4.30 5.50
Floresta CTC.efetiva 1.49 0.37 24.91 0.90 1.25 1.50 1.60 2.60
Floresta Ca 0.63 0.23 36.73 0.30 0.50 0.60 0.70 1.30
Floresta Corg 12.50 1.99 15.92 9.30 10.93 12.50 13.97 16.20
Floresta Densidade 1.38 0.15 11.16 1.10 1.30 1.40 1.40 1.80
Floresta EC 17.47 4.37 25.00 12.20 14.55 16.85 18.67 29.70
Floresta H.Al 3.23 0.42 13.12 2.70 2.90 3.25 3.40 4.40
Floresta K 0.08 0.03 38.04 0.05 0.06 0.07 0.10 0.17
Floresta Macroporosidade 0.15 0.05 34.30 0.10 0.10 0.15 0.20 0.20
Floresta Mat.Org 21.56 3.43 15.89 16.10 18.88 21.50 24.12 28.00
Floresta Mg 0.19 0.12 62.30 0.10 0.10 0.20 0.20 0.60
Floresta Microporosidade 0.27 0.05 16.93 0.20 0.23 0.30 0.30 0.30
Floresta P.Mehlich 3.12 1.76 56.28 0.80 1.40 3.00 4.00 7.00
Floresta PT 0.42 0.05 12.35 0.30 0.40 0.40 0.40 0.50
Floresta RBS 117.65 60.27 51.23 26.30 95.18 105.50 130.52 304.70
Floresta SB 0.93 0.36 38.57 0.46 0.70 0.93 1.08 2.07
Floresta Sat.Al.M 39.19 9.07 23.15 19.50 33.97 38.60 46.82 56.60
Floresta Sat.Base 21.67 6.28 28.99 13.80 16.92 20.95 24.60 37.60
Floresta Silte 63.12 0.00 0.00 63.12 63.12 63.12 63.12 63.12
Floresta Tar 26.95 1.22 4.54 25.76 25.76 26.95 28.14 28.14
Floresta Ts 26.78 0.38 1.40 26.42 26.42 26.78 27.15 27.15
Floresta UR 87.00 7.78 8.95 79.44 79.44 87.00 94.57 94.57
Floresta Us 0.20 0.04 20.58 0.16 0.16 0.20 0.24 0.24
Floresta pH.CaCl2 4.10 0.15 3.74 3.80 4.03 4.10 4.10 4.60
Mangueira Al 0.17 0.07 41.16 0.10 0.10 0.20 0.20 0.30
Mangueira Areia 764.24 0.00 0.00 764.24 764.24 764.24 764.24 764.24
Mangueira Argila 126.00 0.00 0.00 126.00 126.00 126.00 126.00 126.00
Mangueira CBM 314.21 314.14 99.98 5.50 61.38 215.25 514.07 976.40
Mangueira CTC.a.pH7.cmolc 5.03 0.66 13.19 4.10 4.53 5.05 5.40 6.70
Mangueira CTC.efetiva 2.26 0.58 25.75 1.40 1.83 2.15 2.75 3.20
Mangueira Ca 1.52 0.48 31.59 0.80 1.12 1.40 1.80 2.40
Mangueira Corg 9.40 1.24 13.21 7.00 8.55 9.30 10.35 11.20
Mangueira Densidade 1.44 0.12 8.64 1.30 1.40 1.40 1.50 1.80
Mangueira EC 13.52 2.10 15.50 9.70 12.45 13.80 14.90 17.70
Mangueira H.Al 2.93 0.49 16.86 2.40 2.62 2.75 3.05 4.20
Mangueira K 0.08 0.01 19.53 0.06 0.06 0.07 0.09 0.11
Mangueira Macroporosidade 0.17 0.07 38.85 0.10 0.10 0.20 0.20 0.30
Mangueira Mat.Org 16.20 2.15 13.27 12.00 14.82 16.00 17.92 19.40
Mangueira Mg 0.47 0.17 36.90 0.20 0.32 0.40 0.60 0.80
Mangueira Microporosidade 0.26 0.05 20.01 0.20 0.20 0.30 0.30 0.30
Mangueira P.Mehlich 3.98 2.16 54.40 1.00 2.00 4.00 5.75 8.00
Mangueira PT 0.42 0.05 12.99 0.30 0.40 0.40 0.48 0.50
Mangueira RBS 155.46 83.58 53.76 34.90 92.83 167.75 212.48 354.20
Mangueira SB 2.09 0.63 30.14 1.20 1.61 2.00 2.63 3.11
Mangueira Sat.Al.M 8.53 5.17 60.55 3.10 3.68 8.10 11.68 20.50
Mangueira Sat.Base 40.86 9.30 22.76 21.70 35.55 40.80 46.08 56.20
Mangueira Silte 109.76 0.00 0.00 109.76 109.76 109.76 109.76 109.76
Mangueira Tar 27.65 0.55 1.97 27.12 27.12 27.65 28.18 28.18
Mangueira Ts 26.78 0.38 1.40 26.42 26.42 26.78 27.15 27.15
Mangueira UR 76.31 3.06 4.01 73.33 73.33 76.31 79.28 79.28
Mangueira Us 0.17 0.07 38.22 0.11 0.11 0.17 0.24 0.24
Mangueira pH.CaCl2 4.71 0.20 4.18 4.40 4.53 4.70 4.80 5.00

Descritiva: “Uso de Solo” como fator

Média ± desvio padrão dos atributos do solo por uso do solo.
Uso.do.solo n Al Areia Argila CBM CTC.a.pH7.cmolc CTC.efetiva Ca Corg Densidade EC H.Al K Macroporosidade Mat.Org Mg Microporosidade P.Mehlich PT RBS SB Sat.Al.M Sat.Base Silte Tar Ts UR Us pH.CaCl2
Abandonada 18 0,66 ± 0,21 787,88 ± 0,00 152,00 ± 0,00 111,06 ± 74,91 3,77 ± 0,70 1,35 ± 0,27 0,43 ± 0,10 9,41 ± 3,32 1,47 ± 0,13 13,58 ± 4,43 3,08 ± 0,63 0,09 ± 0,03 0,17 ± 0,05 16,21 ± 5,72 0,16 ± 0,05 0,23 ± 0,05 4,04 ± 1,61 0,41 ± 0,08 81,11 ± 45,21 0,69 ± 0,14 48,57 ± 8,52 18,29 ± 3,50 60,12 ± 0,00 28,53 ± 1,58 32,91 ± 0,00 78,97 ± 5,68 0,15 ± 0,00 4,04 ± 0,11
Floresta 18 0,56 ± 0,11 727,88 ± 0,00 209,00 ± 0,00 131,83 ± 166,32 4,17 ± 0,55 1,49 ± 0,37 0,63 ± 0,23 12,50 ± 1,99 1,38 ± 0,15 17,47 ± 4,37 3,23 ± 0,42 0,08 ± 0,03 0,15 ± 0,05 21,56 ± 3,43 0,19 ± 0,12 0,27 ± 0,05 3,12 ± 1,76 0,42 ± 0,05 117,65 ± 60,27 0,93 ± 0,36 39,19 ± 9,07 21,67 ± 6,28 63,12 ± 0,00 26,95 ± 1,22 26,79 ± 0,38 87,00 ± 7,78 0,20 ± 0,04 4,10 ± 0,15
Mangueira 18 0,17 ± 0,07 764,24 ± 0,00 126,00 ± 0,00 314,21 ± 314,14 5,03 ± 0,66 2,26 ± 0,58 1,52 ± 0,48 9,40 ± 1,24 1,44 ± 0,12 13,52 ± 2,10 2,93 ± 0,49 0,08 ± 0,01 0,17 ± 0,07 16,20 ± 2,15 0,47 ± 0,17 0,26 ± 0,05 3,98 ± 2,16 0,42 ± 0,05 155,46 ± 83,58 2,09 ± 0,63 8,53 ± 5,17 40,86 ± 9,30 109,76 ± 0,00 27,65 ± 0,55 26,79 ± 0,38 76,31 ± 3,06 0,17 ± 0,07 4,71 ± 0,20

Teste estatístico

Aplicação do teste de Kruskal-Wallis com comparações múltiplas de Dunn

Difrenças entre Periodos

Comparação dos atributos do solo entre períodos (Kruskal-Wallis + Dunn).
Periodo Al Areia Argila CBM CTC.a.pH7.cmolc CTC.efetiva Ca Corg Densidade EC H.Al K Macroporosidade Mat.Org Mg Microporosidade P.Mehlich PT RBS SB Sat.Al.M Sat.Base Silte Tar Ts UR Us pH.CaCl2
Inverno 0,50 a 764,24 a 152,00 a 218,20 a 4,10 a 1,50 a 0,70 a 9,30 a 1,50 a 13,60 a 2,90 a 0,06 a 0,20 a 16,00 a 0,20 a 0,20 a 2,00 a 0,40 a 125,20 a 1,00 a 36,30 a 23,60 a 63,12 a 26,99 a 26,42 a 84,49 a 0,24 a 4,10 a
verão 0,60 a 764,24 a 152,00 a 65,50 a 4,60 a 1,60 a 0,60 a 10,80 a 1,30 a 14,60 a 2,90 a 0,09 a 0,20 a 18,70 a 0,20 a 0,30 a 5,00 a 0,40 a 75,40 a 0,91 a 39,50 a 21,10 a 63,12 a 28,18 a 27,15 a 73,45 a 0,15 a 4,10 a

Difrenças entre Periodos e Uso do Solo

Comparação dos atributos do solo por uso do solo (Kruskal-Wallis + Dunn) em cada período.
Periodo Uso.do.solo Al Areia Argila CBM CTC.a.pH7.cmolc CTC.efetiva Ca Corg Densidade EC H.Al K Macroporosidade Mat.Org Mg Microporosidade P.Mehlich PT RBS SB Sat.Al.M Sat.Base Silte Tar UR pH.CaCl2 Ts Us
Inverno Abandonada 0.51 a 787.88 a 152 a 123.33 a 3.42 a 1.19 a 0.41 a 7.29 a 1.57 a 11.27 a 2.74 a 0.07 a 0.17 a 12.56 a 0.17 a 0.2 a 2.76 a 0.4 a 106.41 a 0.68 a 44.28 a 19.26 a 60.12 a 26.99 a 84.49 a 3.99 a NA NA
Inverno Floresta 0.53 a 727.88 b 209 b 193.94 a 4.12 a 1.49 ab 0.68 b 12.88 b 1.48 b 19.16 b 3.17 b 0.07 a 0.12 a 22.22 b 0.18 a 0.26 b 1.8 a 0.4 a 145.42 ab 0.96 a 37.16 a 22.61 a 63.12 b 25.76 b 94.57 b 4.06 a NA NA
Inverno Mangueira 0.19 b 764.24 c 126 c 566.59 b 4.72 b 1.88 b 1.21 c 9.17 c 1.51 ab 13.79 c 3.03 ab 0.06 a 0.19 a 15.78 c 0.38 b 0.22 ab 2.29 a 0.41 a 194.02 b 1.69 b 10.88 b 35.4 b 109.76 c 27.12 c 79.28 c 4.59 b NA NA
verão Abandonada 0.8 a 787.88 a 152 a 98.79 a 4.12 a 1.51 a 0.46 a 11.53 a 1.38 a 15.9 a 3.41 a 0.11 a 0.17 a 19.87 a 0.14 a 0.27 a 5.33 a 0.41 a 55.81 a 0.71 a 52.87 a 17.33 a 60.12 a 30.07 a 73.45 a 4.09 a 32.91 a 0.15 a
verão Floresta 0.58 a 727.88 b 209 b 69.71 a 4.21 a 1.49 a 0.59 a 12.12 a 1.29 a 15.78 a 3.3 a 0.1 a 0.18 a 20.9 a 0.21 a 0.29 a 4.44 a 0.43 a 89.88 a 0.9 a 41.22 a 20.73 a 63.12 b 28.14 b 79.44 b 4.14 a 27.15 b 0.16 b
verão Mangueira 0.14 b 764.24 c 126 c 61.82 a 5.33 b 2.64 b 1.83 b 9.63 a 1.37 a 13.26 a 2.83 a 0.09 a 0.16 a 16.62 a 0.57 b 0.29 a 5.67 a 0.43 a 116.9 a 2.49 b 6.19 b 46.32 b 109.76 c 28.18 c 73.33 c 4.83 b 27.15 b 0.11 c

PCA por Uso do Solo, avaliando diferenças entre os periodos:

Tabelas

Resumo da PCA para o uso do solo Abandonada.
Atributo CP1 CP2 CP3 CP4
Autovalores 9.328 4.949 2.492 1.871
Proporção 0.406 0.215 0.108 0.081
Acumulado 0.406 0.621 0.729 0.810
UR 0.290 0.078 -0.197 0.150
Tar -0.290 -0.078 0.197 -0.150
pH.CaCl2 -0.060 -0.145 0.488 -0.128
Ca -0.088 0.368 0.272 -0.020
Mg 0.009 0.360 0.108 0.264
SB -0.060 0.403 0.239 0.003
Al -0.290 -0.052 -0.087 0.145
H.Al -0.262 0.141 -0.214 -0.143
CTC.a.pH7.cmolc -0.253 0.208 -0.151 -0.129
CTC.efetiva -0.272 0.168 0.043 0.118
P.Mehlich -0.265 -0.169 0.041 -0.109
K -0.157 0.058 0.102 -0.536
Mat.Org -0.286 0.097 -0.099 0.206
Sat.Al.M -0.200 -0.298 -0.189 0.110
Sat.Base 0.148 0.237 0.392 0.094
Corg -0.286 0.099 -0.097 0.204
EC -0.251 0.181 -0.109 0.215
RBS 0.166 0.104 -0.221 -0.241
CBM 0.052 -0.082 -0.046 0.364
PT 0.005 -0.249 0.279 0.351
Macroporosidade 0.067 -0.219 0.260 0.056
Microporosidade -0.229 -0.147 0.156 0.145
Densidade 0.235 0.258 -0.098 0.092
Resumo da PCA para o uso do solo Floresta.
Atributo CP1 CP2 CP3 CP4
Autovalores 8.315 6.361 3.167 2.057
Proporção 0.333 0.254 0.127 0.082
Acumulado 0.333 0.587 0.714 0.796
Us 0.235 0.261 -0.115 0.060
Ts -0.235 -0.261 0.115 -0.060
UR 0.235 0.261 -0.115 0.060
Tar -0.235 -0.261 0.115 -0.060
pH.CaCl2 -0.246 0.132 0.195 -0.207
Ca -0.173 0.333 -0.020 0.009
Mg -0.245 0.235 0.161 -0.051
SB -0.210 0.311 0.042 0.004
Al -0.093 -0.111 -0.281 0.292
H.Al -0.121 -0.063 -0.251 0.375
CTC.a.pH7.cmolc -0.229 0.151 -0.163 0.278
CTC.efetiva -0.233 0.263 -0.043 0.076
P.Mehlich -0.268 -0.113 0.214 0.078
K -0.295 0.026 0.130 0.193
Mat.Org 0.194 -0.038 0.402 0.210
Sat.Al.M 0.103 -0.320 -0.192 0.100
Sat.Base -0.156 0.319 0.110 -0.077
Corg 0.194 -0.041 0.400 0.210
EC 0.216 0.056 0.369 0.239
RBS 0.059 0.193 -0.201 0.247
CBM 0.125 0.143 -0.076 0.250
PT -0.122 -0.104 0.113 0.406
Macroporosidade -0.187 -0.126 -0.053 0.283
Microporosidade -0.154 -0.047 0.176 0.155
Densidade 0.213 0.172 0.219 0.184
Resumo da PCA para o uso do solo Mangueira.
Atributo CP1 CP2 CP3 CP4
Autovalores 12.129 3.561 2.889 1.911
Proporção 0.485 0.142 0.116 0.076
Acumulado 0.485 0.628 0.743 0.820
Us 0.255 0.204 0.075 -0.078
Ts -0.255 -0.204 -0.075 0.078
UR 0.255 0.204 0.075 -0.078
Tar -0.255 -0.204 -0.075 0.078
pH.CaCl2 -0.198 0.018 0.064 -0.017
Ca -0.260 0.133 0.138 0.009
Mg -0.238 0.195 0.035 -0.174
SB -0.263 0.160 0.116 -0.043
Al 0.174 -0.135 -0.317 0.274
H.Al 0.085 -0.114 -0.267 -0.201
CTC.a.pH7.cmolc -0.188 0.063 -0.090 -0.189
CTC.efetiva -0.266 0.152 0.087 -0.012
P.Mehlich -0.195 -0.266 -0.199 -0.030
K -0.244 -0.075 -0.073 0.047
Mat.Org -0.091 0.276 -0.438 0.020
Sat.Al.M 0.214 -0.155 -0.303 0.149
Sat.Base -0.242 0.174 0.202 0.049
Corg -0.090 0.281 -0.432 0.016
EC -0.014 0.350 -0.375 -0.204
RBS 0.189 -0.112 0.019 -0.050
CBM 0.207 0.211 0.184 0.107
PT -0.025 -0.287 0.014 -0.569
Macroporosidade 0.127 -0.224 0.028 -0.439
Microporosidade -0.189 -0.180 -0.139 -0.113
Densidade 0.117 0.252 -0.022 -0.422

PCA com apenas as 10 variáveis que mais contribíram em cada período

Analise Geral (Periodos e Usos do Solo)

Tabela

Resumo da PCA geral dos atributos do solo.
Atributo CP1 CP2 CP3 CP4
Autovalores 9.475 5.481 4.124 1.955
Proporção (%) 33.839 19.577 14.728 6.981
Acumulado (%) 33.839 53.416 68.144 75.125
Us 0.068 -0.333 0.100 -0.065
Ts 0.105 0.322 -0.091 -0.241
UR 0.146 -0.324 0.122 0.075
Tar -0.002 0.391 -0.086 -0.114
pH.CaCl2 -0.304 0.025 0.026 0.038
Ca -0.310 0.006 0.106 -0.033
Mg -0.299 0.032 0.101 -0.011
SB -0.312 0.014 0.106 -0.030
Al 0.268 0.165 0.031 -0.056
H.Al 0.088 0.130 0.198 -0.089
CTC.a.pH7.cmolc -0.226 0.098 0.224 -0.083
CTC.efetiva -0.281 0.092 0.150 -0.061
P.Mehlich -0.068 0.359 0.029 0.083
K 0.003 0.285 0.048 0.006
Mat.Org 0.107 0.056 0.417 -0.111
Sat.Al.M 0.303 0.108 -0.060 -0.013
Sat.Base -0.310 -0.029 0.060 -0.012
Corg 0.107 0.055 0.417 -0.113
EC 0.097 -0.022 0.385 -0.224
RBS -0.082 -0.234 -0.028 0.119
CBM -0.079 -0.231 -0.045 -0.203
PT -0.030 0.079 -0.017 0.518
Macroporosidade -0.018 0.038 -0.192 0.386
Microporosidade -0.038 0.209 0.247 0.257
Densidade 0.004 -0.216 -0.086 -0.323
Argila 0.219 -0.063 0.258 0.249
Areia -0.034 0.144 -0.364 -0.325
Silte -0.297 -0.060 0.002 -0.026

Gráficos

Versao com apenas as 10 variáveis mais importantes:

Abordagem hierárquica multivariada por blocos funcionais (Como sugerida no texto, via cluters).

Realizei inicialmente apenas para os Periodos “inverno” e “verão”:

Obs: Como demonstrado no código abaixo, não consegui encontrar na planilha algumas das variaveis sugeridas no texto enviado, por isso não estão nos clusters. São elas: “V”,“m”, “MO”, “CO”,“Ds”. Aguardo retorno para ajuste

clusters_variaveis <- list(
  
  acidez = c("pH.CaCl2", "Al", "H.Al"),
  
  fertilidade = c("Ca", "Mg", "K",  "CTC.a.pH7.cmolc", "P.Mehlich"),
  
  biologico = c( "RBS", "CBM", "EC"),
  
  fisico = c("PT", "Macroporosidade", "Microporosidade",  "Areia", "Argila", "Silte")
)

Clusters (Distância de Ward.D2):

gráfico PCA por cluster

Tabela da Permanova

Resultados da PERMANOVA para cada grupo de variáveis em função do período.
Cluster F_modelo R2 p_valor
acidez 2.1165 0.0391 0.112
fertilidade 14.6686 0.2200 0.001
biologico 8.6736 0.1430 0.001
fisico 3.1903 0.0578 0.014

Indice de Qualidade do Solo

Essa etapa, considero embrionária, pois preciso das informações sobre os pesos. Dessa forma, calculei índices teóricos a título de exemplo, para sugestões posteiores.

SQI-1 - baseado na contribuição das variáveis na PCA;

SQI-2 - baseado na redução de redundância via correlação;

SQI-3 - baseado em indicadores-chave de fácil mensuração (chutei alguns).

SQI-1 (PCA)

  • Variáveis selecionadas por PCA:
## [1] "Tar"       "P.Mehlich" "UR"        "Us"        "Ts"

SQI-2 (Correlação)

  • Variáveis selecionadas por correlação:
## [1] "Mg"    "K"     "PT"    "Areia"

SQI-3 (Campo)

  • Variáveis selecionadas por conveniência (mensuração):
## [1] "Ts"              "Corg"            "Densidade"       "Macroporosidade"
## [5] "CBM"

Definir comportamento das variáveis. O código diz o seguinte, por exemplo, para a variável “Tar”: “mais_melhor”, quanto maior o valor dessa variável maior será o valor do índice. (aqui, é onde é preciso conhecimento técnico sobre os atributos, por isso coloquei atribuições arbitrárias que precisam ser ajustadas).

tipos_sqi1 <- list(
  Tar = "mais_melhor",        
  P.Mehlich = "mais_melhor",
  UR = "mais_melhor",
  Us = "menos_melhor",
  Ts = "mais_melhor"

)

tipos_sqi2 <- list(
  Mg = "mais_melhor",
  K  = "mais_melhor",
  PT  = "mais_melhor",
  Areia = "menos_melhor"
)

tipos_sqi3 <- list(
  Ts = "mais_melhor",        
  Corg = "mais_melhor",
  Densidade = "menos_melhor",
  Macroporosidade = "mais_melhor",
  CBM = "mais_melhor"
)

Gráfico

Média, erro padrão e letras de comparação para os índices de qualidade do solo.
indice grupo n media sd se letra
SQI-1 Abandonada 18 0.4230909 0.1114820 0.0262766 a
SQI-1 Floresta 18 0.3878429 0.0483683 0.0114005 a
SQI-1 Mangueira 18 0.2609476 0.0669928 0.0157904 b
SQI-2 Abandonada 18 0.1111277 0.0652106 0.0153703 a
SQI-2 Floresta 18 0.2531200 0.1468130 0.0346042 b
SQI-2 Mangueira 18 0.5230600 0.2067442 0.0487301 c
SQI-3 Abandonada 18 0.4134699 0.2920667 0.0688408 ab
SQI-3 Floresta 18 0.3938526 0.0727594 0.0171496 a
SQI-3 Mangueira 18 0.2798983 0.0641016 0.0151089 b