Cooperativa Agroindutrial Nova Aliança - Safra 2018/19
Equipe: Débora Lang, Leonardo Reffatti, Mônica Rodrigues, Paulo Dullius e Waldir Schu. Autor do documento: Leonardo Reffatti
Março de 2019
Relatório agronômico referente a Safra 2018/19
O objetivo deste relatório é detalhar em seções as informações dos vinhedos dos cooperados da Cooperativa Agroindustrial Nova Aliança. Com a exploração de dados das áreas dos vinhedos (parcelas), variedades de uva, produção, produtividade e qualidade das uvas (grau glucométrico), rentabilidade, cadernos de campo e monitoramento do uso de agrotóxicos. Neste sentido, integramos informações componentes do processo de rastreabilidade agrícola da cooperativa Nova Aliança.
Dividimos este relatório em três seções listadas abaixo:
Seção 1: Mapeamento de vinhedos 2018
Seção 2: Cadernos de Campo 2018 e Monitoramento de aplicações de Agrotóxicos - Safra 2019
Seção 3: Análise de Rentabilidade em função de Produtividade e Grau Glucométrico - Safra 2018
Seção 1: Mapeamento de vinhedos 2018
Nesta seção apresentamos os dados coletados em campo registrados em parcelas de vinhedos organizados pela equipe de Agronomia da Coop. Nova Aliança. Para fazer a estruturação dos dados de campo foram utilizados unidades amostrais denominadas de parcelas.
Estas parcelas contém informações dos Grupos Familiares associados a Nova Aliança, das variedades/cultivares de uva, dos números dos cadastros vitícolas, das coordenadas geográficas, das altitudes, dos números dos cadernos de campo e das áreas em hectares. Com essas informações dos pontos é possível buscar correlações entre as variáveis e com outras bases de dados para gerar mais resultados e interpretações.
Neste relatório apresentamos algumas análises de estatísticas aplicadas nessas unidades amostrais. As análises aplicadas foram estatística descritiva básica, histogramas, gráficos de médias, teste T, Anova One-Way, Anova Fatorial, Regressão linear simples.
Partimos do ponto inicial de uma visualização geral da área de vinhedos da Cooperativa Nova Aliança no Mapa abaixo. Neste mapeamento (heatmap) quanto maior for a presença de vinhedos da cooperativa mais tons avermelhados são percebidos no mapa. A escala de intensidade vai do vermelho (maior quantidade) para o azul (menor quantidade).
Seguimos para detalhar as variáveis coletadas em cada uma das unidades amostrais (UA). Os resumos através das funções abaixo nos trazem alguns parâmetros estatísticos calculados para cada uma das variáveis nas UA chamadas de parcelas.
## 'data.frame': 3440 obs. of 13 variables:
## $ UA : Factor w/ 3440 levels "ADAIR_PASTORELLO|BORDO_(IVES)|P4|22530|0.53ha|n_cc1",..: 1 2 3 4 5 6 7 8 9 10 ...
## $ GRUPO_FAMILIAR: Factor w/ 444 levels "ADAIR_PASTORELLO",..: 1 1 1 1 1 1 1 1 2 2 ...
## $ VARIEDADE : Factor w/ 61 levels "ALICANTE_BOUCHET",..: 4 5 6 6 9 22 36 42 4 4 ...
## $ PARCELA : Factor w/ 135 levels "P1","P1.1","P1.2",..: 77 87 80 83 63 1 91 41 77 80 ...
## $ AREA : num 0.53 0.91 1.19 0.98 0.33 2.68 1.58 0.11 0.32 0.65 ...
## $ CV : int 22530 22530 22530 22530 22530 22530 22530 22530 4528 4528 ...
## $ ALTITUDE : int 510 450 531 474 486 470 599 512 580 621 ...
## $ LON : num -51.5 -51.5 -51.5 -51.5 -51.5 ...
## $ LAT : num -29 -29 -29 -29 -29 ...
## $ N_CC : int 1 2 3 3 4 5 6 7 9 9 ...
## $ NUCLEO : Factor w/ 6 levels "FC","JC","NP",..: 4 4 4 4 4 4 4 4 3 3 ...
## $ NUCLEO_N : int 1 1 1 1 1 1 1 1 5 5 ...
## $ COOP : int 1 1 1 1 1 1 1 1 1 1 ...## UA
## ADAIR_PASTORELLO|BORDO_(IVES)|P4|22530|0.53ha|n_cc1: 1
## ADAIR_PASTORELLO|BRS_CARMEM|P7|22530|0.91ha|n_cc2 : 1
## ADAIR_PASTORELLO|BRS_CORA|P5|22530|1.19ha|n_cc3 : 1
## ADAIR_PASTORELLO|BRS_CORA|P6|22530|0.98ha|n_cc3 : 1
## ADAIR_PASTORELLO|BRS_VIOLETA|P3|22530|0.33ha|n_cc4 : 1
## ADAIR_PASTORELLO|ISABEL|P1|22530|2.68ha|n_cc5 : 1
## (Other) :3434
## GRUPO_FAMILIAR VARIEDADE
## COOP._AGROINDUSTRIAL_NOVA_ALIANCA_LTDA: 36 ISABEL : 901
## FABIO_BRUSTOLIN : 34 BORDO_(IVES) : 574
## CELIO_FERRARI : 30 NIAGARA_BRANCA: 333
## SANTO_JOSE_BUSETTI : 25 NIAGARA_ROSADA: 222
## DEONILDA_MANARA_DETONI : 21 BRS_CORA : 186
## ADELAR_PARISOTTO : 20 BRS_CARMEM : 132
## (Other) :3274 (Other) :1092
## PARCELA AREA CV ALTITUDE
## P1 : 436 Min. :0.0100 Min. : 194 Min. : 22.0
## P2 : 414 1st Qu.:0.2200 1st Qu.: 3282 1st Qu.:566.0
## P3 : 385 Median :0.4050 Median : 6779 Median :623.0
## P4 : 356 Mean :0.5354 Mean : 8305 Mean :608.8
## P5 : 312 3rd Qu.:0.6900 3rd Qu.:13035 3rd Qu.:672.0
## P6 : 273 Max. :5.6500 Max. :27353 Max. :982.0
## (Other):1264
## LON LAT N_CC NUCLEO
## Min. :-55.42 Min. :-31.25 Min. : 0.0 FC:555
## 1st Qu.:-51.44 1st Qu.:-29.12 1st Qu.: 504.8 JC:681
## Median :-51.34 Median :-29.08 Median :1131.0 NP:652
## Mean :-51.39 Mean :-29.12 Mean :1161.6 PB:862
## 3rd Qu.:-51.25 3rd Qu.:-29.03 3rd Qu.:1774.2 SG:496
## Max. :-50.62 Max. :-28.78 Max. :2844.0 SV:194
##
## NUCLEO_N COOP
## Min. :1.00 Min. :1
## 1st Qu.:1.00 1st Qu.:1
## Median :3.00 Median :1
## Mean :3.22 Mean :1
## 3rd Qu.:5.00 3rd Qu.:1
## Max. :6.00 Max. :1
## ## [1] 1841.77## [1] 444## [1] 4.148131Conforme os “outputs”" das funções acima, os dados apresentam 1841.77 hectares de vinhedos distribuidos nos 444 grupos familiares, onde a área média por grupos familiares é de 4.14 hectares. São 4,5 variedades de uva a média no grupos familiares.
No próximo MAPA abaixo é possível visualizar a localização geográfica de todas as parcelas da Cooperativa Nova Aliança.
Cada ponto azul no mapa se refere ao ponto central do vinhedo (parcela). A marcação dos pontos no mapa é feita através do quadro “checkbox” posicionado na margem superior direita. Neste quadro criamos da opção de marcação OSM que é a definição padrão e Satélite que coloca a visualização em modo paisagem. Conforme nos aproximamos ou nos afastamos da visualização das parcelas, o mapa faz a separação ou união dos pontos, o valor informado internamente indica a soma de parcelas dentro daquele grupo “cluster”.
Ao clicar no ponto, cada uma das parcelas está identificada com o grupo familiar, cadastro vitícola, variedade/cultivar, número da parcela, área em hectares e número do caderno de campo correspondente. Por exemplo, DELVIO_WAIS|BORDO_(IVES)|P9|14779|1.23ha|n_cc582.
## [1] 3440## [1] 61São 3440 parcelas com mais de 60 variedades diferentes presentes. As variedades com maior número de parcelas são Isabel (901), Bordô (574) e Niágara Branca (333).
A variedade Isabel é a mais frentente nas parcelas, apresenta um somatório de área de 572.97 hectares, sendo a variedade com maior área plantada da Nova Aliança. Seguido de Bordô com 351.21 hectares e Niagara Branca 120.80 hectares.
## Variedade x
## 1 ALICANTE_BOUCHET 2.79
## 2 ANCELOTA 7.60
## 3 BAILEY 8.03
## 4 BORDO_(IVES) 351.21
## 5 BRS_CARMEM 63.66
## 6 BRS_CORA 96.18
## 7 BRS_LORENA 32.08
## 8 BRS_MAGNA 18.77
## 9 BRS_VIOLETA 44.28
## 10 CABERNET_FRANC 1.32
## 11 CABERNET_SAUVIGNON 29.99
## 12 CHARDONNAY 25.42
## 13 CONCORD_(FRANCESA) 15.94
## 14 CONCORD_CLONE_30 3.83
## 15 COUDERC_TINTA_(SEIBEL_1077) 29.69
## 16 COUDERC_TREZE_(BRANCA) 16.28
## 17 CYNTHIANA_(ZEPERINA_/_SANTIAGO) 8.74
## 18 DONA_ZILA 0.10
## 19 FREISA 1.15
## 20 GEWURZTRAMINER 1.66
## 21 HERBEMONT_(BORGONHA_/_CHAMPANHA) 10.97
## 22 ISABEL 572.97
## 23 ISABEL_PRECOCE 60.69
## 24 ITALIA 1.53
## 25 JACQUEZ_(SEIBEL_/_PICA_LONGA) 13.90
## 26 LONGANESI 1.34
## 27 MALBEC 2.50
## 28 MALVASIA_CHIANTI 1.05
## 29 MALVASIA_DE_CANDIA 24.03
## 30 MALVASIA_VERDE 0.45
## 31 MARSELAN 3.31
## 32 MERLOT 28.49
## 33 MOSCATO_BIANCO 0.68
## 34 MOSCATO_BRANCO 38.23
## 35 MOSCATO_DE_ALEXANDRIA 0.77
## 36 MOSCATO_EMBRAPA 59.92
## 37 MOSCATO_GIALLO 3.31
## 38 MOSCATO_HAMBURGO 0.59
## 39 MOSCATO_NAZARENO 0.67
## 40 MUDAS 0.19
## 41 NIAGARA_BRANCA 120.80
## 42 NIAGARA_ROSADA 65.68
## 43 PEVERELA 0.26
## 44 PINOT_GRIS 1.43
## 45 PINOT_NOIR 12.60
## 46 PORTA_ENXERTO 1.85
## 47 PROCECCO 0.34
## 48 RIESLING 1.66
## 49 RIESLING_ITALICO 0.55
## 50 RUBEA 10.29
## 51 RUBI 0.13
## 52 RUBI;_ITALIA 0.24
## 53 SAUVIGNON_BLANC 0.69
## 54 SEIBEL_2_(SEIBELETO_/_SEIBEL_CURTA) 11.65
## 55 SEYVE_VILLARD_5276_(BRANCA) 1.03
## 56 SEYVE_VILLARD_TINTA 0.60
## 57 TANNAT 12.11
## 58 TARDIA_DE_CAXIAS 0.60
## 59 TREBIANO 13.44
## 60 VENUS 0.50
## 61 VIOGNIER 1.01
## Variedade.1 x.1
## 1 ALICANTE_BOUCHET 3
## 2 ANCELOTA 4
## 3 BAILEY 25
## 4 BORDO_(IVES) 574
## 5 BRS_CARMEM 132
## 6 BRS_CORA 186
## 7 BRS_LORENA 65
## 8 BRS_MAGNA 44
## 9 BRS_VIOLETA 91
## 10 CABERNET_FRANC 5
## 11 CABERNET_SAUVIGNON 22
## 12 CHARDONNAY 12
## 13 CONCORD_(FRANCESA) 52
## 14 CONCORD_CLONE_30 8
## 15 COUDERC_TINTA_(SEIBEL_1077) 88
## 16 COUDERC_TREZE_(BRANCA) 35
## 17 CYNTHIANA_(ZEPERINA_/_SANTIAGO) 11
## 18 DONA_ZILA 1
## 19 FREISA 1
## 20 GEWURZTRAMINER 1
## 21 HERBEMONT_(BORGONHA_/_CHAMPANHA) 38
## 22 ISABEL 901
## 23 ISABEL_PRECOCE 115
## 24 ITALIA 6
## 25 JACQUEZ_(SEIBEL_/_PICA_LONGA) 37
## 26 LONGANESI 1
## 27 MALBEC 4
## 28 MALVASIA_CHIANTI 1
## 29 MALVASIA_DE_CANDIA 52
## 30 MALVASIA_VERDE 2
## 31 MARSELAN 1
## 32 MERLOT 31
## 33 MOSCATO_BIANCO 1
## 34 MOSCATO_BRANCO 60
## 35 MOSCATO_DE_ALEXANDRIA 2
## 36 MOSCATO_EMBRAPA 119
## 37 MOSCATO_GIALLO 8
## 38 MOSCATO_HAMBURGO 1
## 39 MOSCATO_NAZARENO 1
## 40 MUDAS 1
## 41 NIAGARA_BRANCA 333
## 42 NIAGARA_ROSADA 222
## 43 PEVERELA 3
## 44 PINOT_GRIS 1
## 45 PINOT_NOIR 10
## 46 PORTA_ENXERTO 4
## 47 PROCECCO 1
## 48 RIESLING 2
## 49 RIESLING_ITALICO 1
## 50 RUBEA 30
## 51 RUBI 1
## 52 RUBI;_ITALIA 2
## 53 SAUVIGNON_BLANC 1
## 54 SEIBEL_2_(SEIBELETO_/_SEIBEL_CURTA) 39
## 55 SEYVE_VILLARD_5276_(BRANCA) 3
## 56 SEYVE_VILLARD_TINTA 4
## 57 TANNAT 13
## 58 TARDIA_DE_CAXIAS 3
## 59 TREBIANO 21
## 60 VENUS 2
## 61 VIOGNIER 2As parcelas dos vinhedos apresentam área média de 0,53 hectare, a parcela com menor área registrada apresenta 0,01 hectare, e a parcela com área máxima apresenta 5,65 hectare.
Seguindo nas análises e resultados sobre as parcelas, as a próximas duas figuras representam um histograma e um “boxplot”, respectivamente. O histograma demonstra a frequência de parcelas com seus respectivos valores de área em hectares, a linha tracejada vermelha informa a média de área das parcelas em 0,53 hectare. Na sequencia o “boxplot” mostra os valores das áreas das parcelas entre os Núcleos da Cooperativa.
## [1] 0.01## [1] 0.5353983## [1] 5.65
Este “boxplot” nos mostra uma distribuição de valores de área muito semelhante e simétricos entre os Núcleos. Porém, essa semelhança de valores de área por parcelas dentro dos núcleos apresenta diferença significativa nas médias de áreas das parcelas nos Núcleos através da aplicação do teste ANOVA (p<0.001). Estes resultados podem ser verificados no Gráfico de médias de valores de áreas nos Núcleos. Portanto há diferença entre as áreas médias das parcelas quando comparamos os Núcleos da Cooperativa Nova Aliança. Lembrando que os dados tiveram resultados de não apresentar normalidade e homogeineidade nos testes de Shapiro-Wilk e Levene. Deduzimos que essa diferença acontece pois as parcelas dos grupos de Santana do Livramento estão localizadas dentro do Núcleo de Flores da Cunha, elevando a média das áreas das parcelas dos vinhedos deste núcleo, diferenciando estatisticamente FC dos demais Núcleos.
## FC JC NP PB SG SV
## 0.7053514 0.5096623 0.5312423 0.4837355 0.4722581 0.5444845## [1] "assign" "call" "coefficients" "contrasts"
## [5] "df.residual" "effects" "fitted.values" "model"
## [9] "qr" "rank" "residuals" "terms"
## [13] "xlevels"
##
## Shapiro-Wilk normality test
##
## data: resultado.anova.areaparcela.nucleo$residuals
## W = 0.77865, p-value < 2.2e-16## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 5 19.686 < 2.2e-16 ***
## 3434
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Df Sum Sq Mean Sq F value Pr(>F)
## NUCLEO 5 20.8 4.157 17.37 <2e-16 ***
## Residuals 3434 822.0 0.239
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1## Tukey multiple comparisons of means
## 95% family-wise confidence level
## factor levels have been ordered
##
## Fit: aov(formula = AREA ~ NUCLEO)
##
## $NUCLEO
## diff lwr upr p adj
## PB-SG 0.01147743 -0.06714391 0.09009878 0.9984085
## JC-SG 0.03740420 -0.04494491 0.11975330 0.7878807
## NP-SG 0.05898427 -0.02413302 0.14210155 0.3290480
## SV-SG 0.07222647 -0.04590562 0.19035856 0.5030388
## FC-SG 0.23309329 0.14689491 0.31929166 0.0000000
## JC-PB 0.02592676 -0.04559543 0.09744895 0.9066918
## NP-PB 0.04750683 -0.02489850 0.11991216 0.4204446
## SV-PB 0.06074904 -0.05010780 0.17160588 0.6234625
## FC-PB 0.22161585 0.14569347 0.29753824 0.0000000
## NP-JC 0.02158007 -0.05485678 0.09801692 0.9666868
## SV-JC 0.03482227 -0.07870877 0.14835332 0.9526041
## FC-JC 0.19568909 0.11591273 0.27546545 0.0000000
## SV-NP 0.01324220 -0.10084726 0.12733167 0.9994763
## FC-NP 0.17410902 0.09353994 0.25467810 0.0000000
## FC-SV 0.16086682 0.04451355 0.27722009 0.0011582
Outro resultado que compreendemos sobre as parcelas de vinhedos, são que estas apresentam variações de tamanho em área em função da altitude. Existe evidência significativa geral de todas as parcelas (p<0,001) de que a altitude influencia negativamente a área das parcelas dos vinhedos ( eq: Altitude = 0.9903-0.00075 Área ). Isso nos mostra que para cada 1m de aumento na altitude, temos uma diminuição de 7,5m2 na área da parcela de vinhedo. Este resultado pode ser visualizado através do Gráfico de dispersão abaixo.
## [1] 22## [1] 982## FC JC NP PB SG SV
## 649.0685 585.4097 633.0429 558.9304 609.8206 713.8299## FC JC NP PB SG SV
## 0.7053514 0.5096623 0.5312423 0.4837355 0.4722581 0.5444845
##
## Call:
## lm(formula = AREA ~ ALTITUDE)
##
## Residuals:
## Min 1Q Median 3Q Max
## -0.8246 -0.3156 -0.1231 0.1651 4.8973
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 1.017e+00 4.818e-02 21.11 <2e-16 ***
## ALTITUDE -7.908e-04 7.794e-05 -10.15 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 0.4879 on 3438 degrees of freedom
## Multiple R-squared: 0.02907, Adjusted R-squared: 0.02879
## F-statistic: 102.9 on 1 and 3438 DF, p-value: < 2.2e-16