Acessibilidade / Reportar erro

Classification of the coefficient of variation for experiments with eucalyptus seedlings in greenhouse1 1 Parte da Dissertação da primeira autora, apresentada ao Programa de Pós-Graduação em Agronomia da Universidade Estadual Paulista “Julio de Mesquita Filho”, Faculdade de Engenharia de Ilha Solteira, com apoio da CAPES (Código de Financiamento 001)

Classificação do coeficiente de variação para experimentos com mudas de eucalipto em casa de vegetação

ABSTRACT

In experimental statistics, the variability, which allows for conclusions related to experimental accuracy, is usually demonstrated by means of experimental coefficient of variation (CV). The aims of this study were to evaluate the distribution of CV in several experiments of eucalyptus growth, compare the existent classification ranges proposed in the literature (methods proposed by Garcia; Costa, Seraphin and Zimmermann, and Pimentel-Gomes), and define the better one to be used as reference for eucalyptus seedlings cultures in greenhouses or protected cultivation. For that, a blank test with 100 repetitions was conducted and evaluations were performed in 30, 60, 90 and 120 days. Moreover, previous data from literature were collected, for comparison proposes. The obtained results were submitted to Shapiro-Wilk and Kolmogorov-Smirnov tests, and the classification tables were submitted to adherence and heterogeneity tests. As a result, it was possible to conclude that data obtained with the experiments with eucalyptus in the present research did not fit the classification of Pimentel-Gomes and tables obtained by Costa, Seraphin and Zimmermann, while the tables related to the method proposed by Garcia would be preferably used.

Key words:
Experimental precision; Experimental variation; Parametric tests; Non-parametric tests; E. grandis x E. urophylla

RESUMO

Na estatística experimental, a variabilidade, que permite conclusões relacionadas à precisão experimental, é geralmente demonstrada por meio do coeficiente de variação experimental (CV). O objetivo deste estudo foi avaliar a distribuição do CV em diversos experimentos de crescimento de eucalipto, comparar as faixas de classificação existentes propostas na literatura (métodos propostos por Garcia; Costa, Seraphin e Zimmermann e Pimentel-Gomes), e definir o melhor a ser utilizado como referência para cultivo de mudas de eucalipto em estufas ou cultivo protegido. Para tanto, foi realizado um teste em branco com 100 repetições e avaliações realizadas em 30, 60, 90 e 120 dias. Além disso, foram coletados dados anteriores da literatura, para fins de comparação. Os resultados obtidos foram submetidos aos testes de Shapiro-Wilk e Kolmogorov-Smirnov, e as tabelas de classificação submetidas aos testes de aderência e heterogeneidade. Como resultado, foi possível concluir que os dados obtidos com os experimentos com eucalipto na presente pesquisa não se enquadram na classificação de Pimentel-Gomes e nas tabelas obtidas por Costa, Seraphin e Zimmermann enquanto as tabelas por meio do método proposto por Garcia seriam preferencialmente utilizadas.

Palavras-chave:
Precisão experimental; Variação experimental; Testes paramétricos; Testes não paramétricos; E. grandis x E. urophylla

INTRODUCTION

The genus Eucalyptus is native of Australia and presents more than 700 cataloged species (SANTAROSA; PENTEADO; GOULART, 2014SANTAROSA, E.; PENTEADO JUNIOR, J. F.; GOULART, I. C. G. R. Cultivo de eucalipto em propriedades rurais: diversificação da produção e renda. Colombo: Embrapa Florestas, 2014.). Due to its economic importance, it is crucial to understand the possible variables related to the growing of eucalyptus crops.

In the experimental statistics, the variability allows the conclusion regarding the experimental precision, and it is frequently demonstrated through the coefficient of experimental variation (CV), which is an estimate of the experimental error in relation to the overall average of the experiment (CRUZ et al., 2012CRUZ, E. A. et al. Coeficiente de variação como medida de precisão em experimentos com tomate em ambiente protegido. Enciclopédia Biosfera, v. 8, n. 14, p. 220-233, 2012.). Statistic tests may help the researcher to understand the behavior and make inferences related to the results obtained in their projects (CARGNELUTTI FILHO; LÚCIO; LOPES, 2009CARGNELUTTI FILHO, A.; LÚCIO, A. D.; LOPES, S. J. Experimentação agrícola e florestal. Santa Maria: UFSM, 2009. 204 p.).

Lower estimated CVs correspond to higher experimental precision, as well as the higher precision, corresponds to greater quality of the experiments, leading to small differences between estimated averages (CRUZ et al., 2012CRUZ, E. A. et al. Coeficiente de variação como medida de precisão em experimentos com tomate em ambiente protegido. Enciclopédia Biosfera, v. 8, n. 14, p. 220-233, 2012.). The CV is a measure of dispersion and its main advantage is the possibility of comparison of variables of different natures, as well as the comparison of different studies involving the same variable, which allows the quantification of the experimental accuracy.

According to Garcia (1989)GARCIA, C. H. Tabelas para classificação do coeficiente de variação. Piracicaba: IPEF, 1989. 11 p., it is important to analyze the CV because when obtaining measures of dispersion or variability, the analysis becomes more informative. The author emphasizes the need of evaluating not only the variable in question, but also the nature of the data and the type of the experiment implemented, in addition to the number of repetitions used in the experiment. The importance of CV classification tables are also highlighted, because they are used as reference to the researcher to verify if the results obtained are in a range of expected values, which may reflect in a good accuracy of the analysis.

The importance of normal distribution is evident, being an essential condition for analysis by several methods, such as parametric statistics. In those cases, when the assumption of normality is invalidated, corresponding non-parametric tests should be use, even if they are not so robust, because if the parametric technique is used, the interpretation and inferences may not be safe or adequate.

Pimentel-Gomes (2009) correlates the ranges of experimental accuracy according to his CV classification table, being high for low CV, good for medium CV, low for high CV and very low accuracy for very high CV. Cargnelutti Filho et al. (2018)CARGNELUTTI FILHO, A. et al. Dimensionamento amostral para avaliação de altura e diâmetro de plantas de timbaúva. Floresta e Ambiente, v. 25, n. 1, 2018. uses this classification and states that CV lower than 10% can guarantee high accuracy of the experiment.

The CV classification was used in several cultures, such as tomato (CRUZ et al., 2012CRUZ, E. A. et al. Coeficiente de variação como medida de precisão em experimentos com tomate em ambiente protegido. Enciclopédia Biosfera, v. 8, n. 14, p. 220-233, 2012.), sugarcane (COUTO; PETERNELLI; BARBOSA, 2013COUTO, M. F.; PETERNELLI, L. A.; BARBOSA, M. H. P. Classification of the coefficients of variation for sugarcane crops. Ciência rural, v. 43, n. 6, p. 957-961, 2013.), eucalyptus (MORA; ARRIAGADA, 2016MORA, F.; ARRIAGADA, O. A classification proposal for coefficients of variation in Eucalyptus experiments involving survival, growth and wood quality variables. Bragantia, v. 75, n. 3, p. 263-267, 2016.), papaya (FERREIRA et al., 2016FERREIRA, J. P. et al. Comparison of methods for classification of the coefficient of variation in papaya. Revista Ceres, v. 63, n. 2, p. 138-144, 2016.) and lettuce (SCHMILDT et al., 2017SCHMILDT, E. R. et al. Coeficiente de variação como medida da precisão em experimentos de alface. Revista Agro@mbiente On-line, v. 11, n. 4, p. 290-295, 2017.).

The aims of this study were to evaluate the distribution of CV in several experiments of eucalyptus seedlings growth, compare the existing classification ranges proposed in the literature, and define the better one to be used as reference for eucalyptus cultures in greenhouses or controlled environments.

MATERIAL AND METHODS

Plant material

The Coefficients of Variation (CV) were obtained from simulations of a blank test and from literature data. To obtain the data related to the blank test, 100 basic units (bu) of clones of the hybrid VR3 (Eucalyptus grandis x Eucalyptus urophylla) were used. The experiment was set up at UNESP/Campus Ilha Solteira (SP, Brazil), in protected cultivation conditions, from February to July 2014. The soil used was collected at a depth of 0.0 to 0.20 m, under eucalyptus plantations, located in the municipality of Três Lagoas (MS, Brazil) (Latitude: 20º 59'S and Longitude 51º 48' W).

The soil was sieved (4 mm) and packed in seedling bags 0.40 m high (8 kg of soil). Limestone and NPK were incorporated in the soil, according to the recommendation of use for that region, for the cultivation of eucalyptus (1.5 mg ha-1 of limestone, 0.5 mg ha-1 NPK - 12-20-16). The eucalyptus seedlings (Eucalyptus spp.), used as indicators, were donated by FIBRIA Ltda. The variables evaluated were: height (H), stem diameter at the base (DB), chlorophyll content in old leaves (ChlorO) and new leaves (ChlorN) at 30, 60, 90 and 120 days, dry mass of the aerial part (DMAP) and the root system (DMRS) at 120 days.

Data from literature were collected from 104 articles in national and international journals from the last 30 years with eucalyptus seedlings under protected cultivation. The variables evaluated were: height (H), stem at base diameter (DB), dry mass of leaves (DML), dry mass of the aerial part (DMAP), dry mass of the root system (DMRS), total dry mass (TDM), number of leaves (NL), leaf area (LA) and root length (RL).

Statistical analysis

The average CV values were obtained for each variable in the study, as well as their respective standard deviations (SD). The CVs were classified as: low, medium, high and very high according to the definition of limits proposed by Garcia (1989)GARCIA, C. H. Tabelas para classificação do coeficiente de variação. Piracicaba: IPEF, 1989. 11 p., Pimentel-Gomes (2009)PIMENTEL-GOMES, F. Curso de estatística experimental. 15. ed. Piracicaba: FEALQ, 2009. 451 p. and Costa, Seraphin and Zimmermann (2002)COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002. (Table 1).

Table 1
Classification ranges of coefficients of variation (CV), according to Garcia (1989), Pimentel-Gomes (2009) and Costa, Seraphin and Zimmermann (2002)

To analyze the normality of the data, two parametric adherence tests were used, which are commonly applied regardless of the sample size, being: Kolmogorov-Smirnov test (1933KOLMOGOROV, A. Sulla determinazione empirica di una lgge di distribuzione. Giornale dell'Istituto Italiano degli Attuari, v. 4, p. 83-91, 1933.; 1936)SMIRNOV, N. V. Sui la distribution de w2 (Criterium de MRv Mises). Comptes Rendus , v. 202, p. 449-452, 1936. and Shapiro-Wilk test (1965)SHAPIRO, S. S.; WILK, M. B. An analysis of variance test for normality (complete samples). Biometrika, v. 52, n. 3/4, p. 591-611, 1965..

The CVs of the experiments were classified into ranges according to each method (Table 1), regardless the distribution presented by the normality test. Since the variables may present equality between the frequencies observed and expected by different methods (CRUZ et al., 2012CRUZ, E. A. et al. Coeficiente de variação como medida de precisão em experimentos com tomate em ambiente protegido. Enciclopédia Biosfera, v. 8, n. 14, p. 220-233, 2012.), the CV classification ranges were calculated by the methodologies of Garcia (1989)GARCIA, C. H. Tabelas para classificação do coeficiente de variação. Piracicaba: IPEF, 1989. 11 p. and Costa, Seraphin and Zimmermann (2002)COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002..

Based on this assumption, classification tables were created for all methods regardless of the result of the normality test and, for greater robustness of the analysis, from the formulation, the tables for each method were compared in relation to the frequency of values of the original data by the adherence tests, to compare with the normal distribution and the heterogeneity test, to verify the differences between methods.

The adherence test was performed under the hypothesis of normality, with approximately 68% of the observations included in an arithmetic average (X) range, plus or minus the standard deviation (SD), 95.4% of the observations in the interval xand 99.7% of those in the interval between 3SD ± x of the total area, corresponds to 1 (FONTELLES, 2012FONTELLES, M. J. Bioestatística aplicada à pesquisa experimental. São Paulo: Livraria da Física, 2012. v. 2.). For this test it was used the chi-square test (X2) (PEARSON, 1900PEARSON, K. X. On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, v. 50, n. 302, p. 157-175, 1900.), which is a non-parametric test that evaluates if the observed distribution of frequencies in the CV classification ranges in each method adjusted to the expected distribution (normal distribution). The CV values obtained were correlated with experimental accuracy from the corresponding expected frequencies, being: low - 15.9%, medium - 68.3%, high - 13.6% and very high - 2.3% according to De Mesquita Lopes, Branco and Soares (2013)DE MESQUITA LOPES, M.; BRANCO, V. T. F. C.; SOARES, J. B. Utilização dos testes estatísticos de Kolmogorov-Smirnov e Shapiro-Wilk para verificação da normalidade para materiais de pavimentação. Transportes, v. 21, n. 1, p. 59-66, 2013..

To compare the tests used, the heterogeneity test based on x2 and the Fisher's exact test (FISHER, 1922FISHER, R. A. On the mathematical foundations of theoretical statistics. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, v. 222, n. 594/604, p. 309-368, 1922.) was calculated for each CV classification range. The x2 was used when the cells had values greater than or equal to 6 and the number of samples was greater than 40. When the cells had values less than 5 and the number of samples less than 40, the Fisher's exact test was used (FONTELLES, 2012FONTELLES, M. J. Bioestatística aplicada à pesquisa experimental. São Paulo: Livraria da Física, 2012. v. 2.).

After performing those tests, a table was formulated for recommendation from the best method, verifying if it corroborated with the criteria for using the methods as recommended by the authors from the classification of the normality tests. In order to verify if there were differences between the recommendations of the classification tables from the blank test and from the literature data, the heterogeneity test was performed on the variables common to these two tables. Statistical analyses were performed using Microsoft Excel and R version 3.6.1 (R CORE TEAM, 2019R CORE TEAM. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2019.).

RESULTS AND DISCUSSION

Kolmogorov-Smirnov and Shapiro-Wilk statistic tests for the blank test

According to the application of the Kolmogorov-Smirnov normality test (KS), all variables presented normal distribution. For the Shapiro-Wilk test (SW), only the variables height - H (at 30 days), stem diameter at the base - DB (at 30, 60 and 90 days), chlorophyll content in old leaves - ChlorO (at 30 days), chlorophyll content in new leaves ChlorN (at 30, 90 and 120 days), dry mass of the aerial part - DMAP (at 120 days) and dry mass of the root system - DMRS (at 120 days) presented a normal distribution; the others had a non-normal distribution, at 5% significance level, suggesting that the variables were analyzed by appropriate methodologies according to the distribution (normal and non-normal), since the observed and expected frequencies may or may not present equality in each method used (Table 2).

Table 2
Kolmogorov-Smirnov and Shapiro-Wilk normality tests for the variables studied

Regarding the classification of the distribution based on the data, the variables chosen were those presenting normality in both methods (SW and KS). Consequently, the classification was recommended following the SW (Table 2), the only one that did not show non-normality in the results, probably because it is more effective for all sample sizes. It is worth to mention that, despite robustness of the SW test being considered low for small sample sizes (TORMAN; COSTER; RIBOLDI, 2012TORMAN, V. B. L.; COSTER, R.; RIBOLDI, J. Normalidade de variáveis: métodos de verificação e comparação de alguns testes não-paramétricos por simulação. Revista HCPA, v. 32, n. 2, p. 227-234, 2012.), the test is still recommended for samples larger than 10 to continue the normality test (TORMAN; COSTER; RIBOLDI, 2012TORMAN, V. B. L.; COSTER, R.; RIBOLDI, J. Normalidade de variáveis: métodos de verificação e comparação de alguns testes não-paramétricos por simulação. Revista HCPA, v. 32, n. 2, p. 227-234, 2012.).

Tables of classification of the coefficient of variation by the methods of Garcia (1989)GARCIA, C. H. Tabelas para classificação do coeficiente de variação. Piracicaba: IPEF, 1989. 11 p., Costa, Seraphin and Zimmermann (2002)COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002. and Pimentel-Gomes (2009)PIMENTEL-GOMES, F. Curso de estatística experimental. 15. ed. Piracicaba: FEALQ, 2009. 451 p. for the blank test

According to the methodology 1, during all period of evaluation, for the variable H, the coefficients of variation presented values ranging from 5.67% to 19.93%. For DB, the experimental values ranged from 3.37% to 9.98%, while for ChlorO and ChlorN, values varied from 1.19% to 7.96% and from 3.71% to 9.00% respectively. At 120 days, the CV values for the variable DMAP and DMRS varied from 10.43% to 32.84% and from 6.44% to 21.59%, respectively (Table 3).

Table 3
CV classification ranges calculated by Garcia's method (1989) at 30, 60, 90 and 120 days for the variables analyzed

It was found that values obtained in each variation coefficient classification range are very particular for each variable, indicating that the variation classification proposed by Pimentel-Gomes (2009), commonly used in several types of research, presents results quite distinct to those obtained in this study, since values for the very high range, for example, are, in general, not greater than 30% (Tables 3 and 4). Thus, it is necessary to classify the CV considering each variable under study (COSTA; SERAPHIN; ZIMMERMANN, 2002COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002.; GARCIA, 1989GARCIA, C. H. Tabelas para classificação do coeficiente de variação. Piracicaba: IPEF, 1989. 11 p.).

However, the classification ranges proposed by Pimentel-Gomes (2009) were based not only on agricultural characteristics, but also on experimental designs, and should not be extended to any experiments.

According to methodology 3, during the evaluation period, the coefficients of variation for H range from 6.53% to 16.35%; for DB from 3.78% to 8.68%; for ChlorO 1.75% to 6.30%; for ChlorN from 4.2% to 7.77%; for DMAP from 10.84% to 29.86% and for DMRS from 8.09% to 18.68% (Table 4). Compared to the methodology of Pimentel-Gomes (2009), it was possible to state that the classification ranges are distinct since the low to very high ranges for the method of Costa, Seraphin, and Zimmermann (2002)COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002. would correspond to the low to average for Pimentel-Gomes (2009).

Table 4
CV classification ranges calculated by the Costa, Seraphin and Zimmermann (2002) method at 30, 60, 90 and 120 days for the variables analyzed, and CV classification ranges recommended by Pimentel-Gomes (2009)

The amplitudes of the CV classifications between methods 1 and 3 presented 28 CV values below 1% and the highest value was 4.16, while for 1 and 2 only a unique CV value was below 1% and 28 were above 10%. For 2 and 3 no CV values below 1% were registered, while 32 CV values above 10% were found. This indicated the similarities between the methodologies calculated in this study. However, when compared to method 2, it can be suggest that the CV classifications are discrepant.

Only the variable DMAP calculated in this study fits the method variable for both methods calculated in this study, the ranges found were equal to method 2; for the other variables the values were close between the calculated methods, but distinct from method 2.

Adherence and heterogeneity test for the blank test

To increase the experimental accuracy, the adherence and heterogeneity test was performed using Fisher's Exact Test for all variables in the blank test (Tables 5, 6, 7 and 8).

Table 5
Comparison of the frequency of results of coefficients of variation (CV) for each variable in the classifications of Garcia (1989), Pimentel-Gomes (2009) and Costa, Seraphin and Zimmermann (2002) and the adherence and heterogeneity tests applied to variables at 30 days
Table 6
Comparison of the frequency of results of coefficients of variation (CV) for each variable within the classifications of Garcia (1989), Pimentel-Gomes (2009) and Costa, Seraphin and Zimmermann (2002) and the adherence and heterogeneity tests applied to variables at 60 days
Table 7
Comparison of the frequency of results of coefficients of variation (CV) for each variable within the classifications of Garcia (1989), Pimentel-Gomes (2009) and Costa, Seraphin and Zimmermann (2002) and the adherence and heterogeneity tests applied to variables at 90 days
Table 8
Comparison of the frequency of results of coefficients of variation (CV) for each variable within the classifications of Garcia (1989), Pimentel-Gomes (2009) and Costa, Seraphin and Zimmermann (2002) and the adherence and heterogeneity tests applied to variables at 120 days

By the adherence test, at 30 days, all variables accepted the null hypothesis (frequencies observed equal to expected frequencies), i.e., CV values following a normal distribution (CRUZ et al., 2012CRUZ, E. A. et al. Coeficiente de variação como medida de precisão em experimentos com tomate em ambiente protegido. Enciclopédia Biosfera, v. 8, n. 14, p. 220-233, 2012.), by methods 1 and 3, corroborating results obtained by the KS test; therefore, it is stated that these variables, by methods 1 and 3, did not depend on the KS test to obtain the CV classification table. For method 2, none of the variables followed a normal distribution (Table 5), which means that the observed frequencies were statistically different from the frequencies expected for a normal distribution. Moreover, for this method, normality tests are necessary to obtain the CV classification table.

At 30 days, for the heterogeneity test, all the variables analyzed were not significant between methods 1 and 3, indicating similarity between the proportions of CVs in the classification ranges of the compared methods. Therefore, it was concluded that there was no statistical difference in relation to the low, medium, high and very high limits for the methods used (CRUZ et al., 2012CRUZ, E. A. et al. Coeficiente de variação como medida de precisão em experimentos com tomate em ambiente protegido. Enciclopédia Biosfera, v. 8, n. 14, p. 220-233, 2012.). Thus, the classification ranges are similar. If no statistical differences exist, any of the methods for the classification of coefficient of variation ranges can be used for the data analyzed (Table 5).

For the other method comparison (2 and 3, 1 and 2) all the variables analyzed were significant different among each other (Table 5), showing that the classification ranges between the methods are different. Thus, it was possible to assert that these methods could be recommended, considering the nature of the data.

At the end of the 60 days, the adherence test shows the same behavior as in the 30 days i.e, the CV values calculated followed a normal distribution using methods 1 and 3. Regarding the method 2, none of the variables followed a normal distribution (Tabel 6). Thus, it is possible to affirm that even at 60 days, it is not necessary to perform normality tests by methods 1 and 3. However, the same pattern was not recorded for method 2.

By the heterogeneity test, at 60 days, the same pattern recorded in 30 days was observed. Between methods 1 and 3, it was possible to conclude that none were different and any of the methods could be recommended. For comparisons between 2 and 3 and 1 and 2, the proportions between the methods were similar, suggesting a specific approach in relation to the nature of the data (Table 6).

At 90 days, for the adherence test, by methods 1 and 3, the calculated CV values followed a normal distribution. Only for method 2 the variable did not follow a normal distribution (Table 7). This assessment followed the same behavior as previously recorded. Thus, it is possible to affirm that, considering methods 1 and 3, it is not necessary to perform normality tests. However, the same does not happen with method 2.

By the heterogeneity test, it was possible to concluded that, between methods 1 and 3, any of the methods can be recommended, while for the comparisons between methods 2 and 3 and 1 and 2, a more specific approach is recommended regarding the nature of the data (Table 7).

At 120 days, for the adherence test, using methods 1 and 3, all variables presented similarity with the normal distribution, while for method 2 only the variables DMAP and DMRS followed a normal distribution (Table 8).

From the heterogeneity test, it was possible to state that, between methods 1 and 3, a more specific approach is recommended in relation to the nature of the data, in all variables analyzed. For methods 2 and 3, any of the proposed methods could be used for the variable DMAP. Thus, it was concluded that between methods 1 and 2, a more specific approach is recommended in relation to variables H and DB. For the other variables, it was possible to concluded that any of the suggested methods can be used. Therefore, could be determine the method that best fits the data set, as well as the level of dependence on the classifications based on the normality tests (Table 8).

Since methods 1 and 3 were the most similar to normal (adherence test) and there was no difference among them or among evaluations (heterogeneity test), the recommendation of the CV classification table was based on the highest calculated values of adherence (Table 9). Therefore, the method 1 was recommended for the data calculated in the blank test.

Table 9
Recommendation of the coefficient of variation classification ranges for the blank test with Eucalyptus in a greenhouse

Kolmogorov-Smirnov and Shapiro-Wilk for literature data

Regarding the data obtained in the literature, all variables presented normal distribution by the KS test. For the SW test, only the variables DML, LA and RL had a normal distribution, the others presented a non-normal distribution, at a significance level of 5% (Table 10).

Since for these data there were variables presenting different values in the two methods, following the criteria adopted, it was recommended the classification from the SW (Table 10), because the tests with the blank test were the only ones where it was found non-normality for some variables.

Table 10
Kolmogorov-Smirnov and Shapiro-Wilk normality tests for the variables found in the literature

Tables of classification of the coefficient of variation by the methods of Garcia (1989)GARCIA, C. H. Tabelas para classificação do coeficiente de variação. Piracicaba: IPEF, 1989. 11 p., Costa, Seraphin and Zimmermann (2002)COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002. and Pimentel-Gomes (2009)PIMENTEL-GOMES, F. Curso de estatística experimental. 15. ed. Piracicaba: FEALQ, 2009. 451 p. for literature data

The CVs, by method 1, ranged from 9.08% to 41.19% for DMAP, 13.25% to 39.65%, 9.28% to 39.06% for DML, 9.25% to 31.12% for TDM, 9.11% to 42.15% for LA and 8.12% to 41.08% for RL. The smallest amplitudes are in variables H (19.92), DB (19.82) and TDM (21.87) (Table 11). As well as in the literature data related to the experiment in a blank test, it can also be seen that the values obtained in each range of classification of the variation coefficient are very peculiar for each variable, which indicates that the proposed variation of Pimentel-Gomes classification (2009) presents very different results when comparing to those obtained in this study (Tables 11 and 12). It can be highlighted the variables LA, RL, DMAP, and DMRS presented values greater than 30% in the average range, whereas for Pimentel-Gomes (2009)PIMENTEL-GOMES, F. Curso de estatística experimental. 15. ed. Piracicaba: FEALQ, 2009. 451 p. values greater than or equal to 30% should be considered very high.

Table 11
CV classification ranges calculated by Garcia's method (1989) for the variables found in the literature

The variations found by method 3 were: 3.81% to 24.54% for H, 6.06% to 22.11% for DB, 6.39% to 45.61% for DML, 8.71% to 35.12% DMAP, 13.39% to 38.42% for DMRS, 8.43% to 31.74% for TDM, 9.12% to 44.46% for LA and 12.31% to 31.67% for RL. The smallest amplitudes are in variables H (20.73), DB (16.05) and RL (19.36) (Table 12). Compared to the methodology of Pimentel-Gomes (2009), it was possible to assert that the classification ranges are distinct, with emphasis on the variables DMRS and RL, which for the low, medium, and high ranges present higher values than those proposed by Pimentel-Gomes (2009).

Table 12
CV classifi cation ranges calculated by the Costa, Seraphin and Zimmermann methods (2002) for the variables found in the literature, and CV classifi cation ranges recommended by Pimentel-Gomes (2009)

The amplitudes of the CV classifications between methods 1 and 3 presented 8 CV values below 1%, with the highest value registered of 9.41, while methods 1, 2 and 3 presented 6 CV values below 1% and 10 above 5%. In method 2 only one CV value was below 1%, while 12 CV values above 5% were registered. Thus, for most of the variables evaluated, there were similarities between the calculated classifications. Despite that, it was possible to affirm that those two were distinct from 2, because only for the variable TDM, in both tables, the ranges found were equal to 2, while for the other variables the values were different between the calculated methodologies and in relation to method 2.

Adherence and heterogeneity test for literature data

Using the adherence test, the probabilities for all variables were similar to the normal distribution by method 1, corroborating the KS test. For method 3, the variables DML, DMAP, DMRS, TDM, LA and RL were similar to normal distribution, while for method 2, only the variables DML, DMAP, TDM and LA. Variables H and DB were not similar to the normal distribution by methods 2 and 3, corroborating SW test (Table 13). Therefore, it is suggested that for the variables showing a normal distribution by the adherence test, the obtained data classification table does not depend on the normality of the test. Different behavior was recorded for other variables not preseting normality, whose observed frequencies were statistically different from those expected for a normal distribution, which means that only the adherence test is not enough to come to a conclusion and that normality tests are needed.

Table 13
Comparison of the frequency of results of coefficients of variation (CV) for each variable within the classifications of Garcia (1989), Pimentel-Gomes (2009) and Costa, Seraphin and Zimmermann (2002) and the adherence and heterogeneity tests applied to variables in the literature

Although considering the heterogeneity test, only variables H, DB and RL were significant in the comparison between the tables of methods 2 and 3. For the comparison between methods 1 and 2, the same variables (H, DB, and RL) were significant. For comparison between methods 1 and 3, variables H, DB and RL were significant, meaning that the proportions between classification ranges among the methods tested presented different behavior for those variables, in all comparisons, for the highlighted variables. Therefore, it is necessary to consider the characteristics of the analyzed variables to determine the appropriate classification range. For the remaining variables, no significant differences were registered between the methods tested. Thus, it was possible to suggest that the classification ranges were similar (Table 13).

The method that best fits the data obtained in the literature was method 1, since it was the only method that resembled the normal distribution in all variables by the adherence test (Table 14).

Table 14
Recommendation of the classifi cation ranges of the coeffi cient of variation with Eucalyptus data obtained in the literature

Comparison between data obtained with blank test and literature data

In the heterogeneity test comparing the classification tables from the blank test (Table 9) and from the literature data (Table 14) no significant differences between the tables generated for the variables tested were registered. Therefore, using only these variables as standard, it was possible to suggest that there was no difference in using either of the two recommendation tables for the variables H, DB, DMAP and DMRS.

For the variables not evaluated in this study, care should be taken with the inferences made, because there are many factors that may change them. Because the nature of the data and the place where the study was conducted should be taken into account (CARGNELUTTI FILHO et al., 2011CARGNELUTTI FILHO, A. et al. Métodos de estimativa do tamanho ótimo de parcelas experimentais de híbridos de milho simples, triplo e duplo. Ciência Rural, v. 41, n. 9, p. 1509-1516, 2011.), requiring, therefore, a specific classification of the CVs associated with each variable studied (SCHMILDT et al., 2017SCHMILDT, E. R. et al. Coeficiente de variação como medida da precisão em experimentos de alface. Revista Agro@mbiente On-line, v. 11, n. 4, p. 290-295, 2017.).

According to the classification recommended in this study, the following studies presented a high precision, for the following variables: Tatagiba et al. (2015)TATAGIBA, S. D. et al. Determinação da máxima capacidade de retenção de água no substrato para produção de mudas de eucalipto em viveiro. Floresta, v. 45, n. 4, p. 745-754, 2015. and Salles et al. (2012)SALLES, T. T. et al. Modelo de Clutter na modelagem de crescimento e produção de eucalipto em sistemas de integração lavoura-pecuária-floresta. Pesquisa Agropecuária Brasileira , v. 47, n. 2, p. 253-260, 2012. with coefficients of variation of 3.5% and 6.75%, respectively for the variable H; Santos Junior et al. (2015)SANTOS JUNIOR, A. et al. Glyphosate drift in eucalyptus plants. Planta Daninha, v. 33, n. 3, p. 615-621, 2015. and Medeiros et al. (2016)MEDEIROS, W. N. et al. Crescimento inicial e concentração de nutrientes em clones de Eucalyptus urophylla x Eucalyptus grandis sob interferência de plantas daninhas. Ciência Florestal , v. 26, n. 1, p. 147-157, 2016. with a CV of 3.79% and 7.95%, respectively for the variable DB; Stahl et al. (2013)STAHL, J. et al. Produção de massa seca e eficiência nutricional de clones de Eucalyptus dunnii e Eucalyptus benthamii em função da adição de doses de fósforo ao solo. Ciência Florestal , v. 23, n. 2, p. 287-295, 2013. and Dias et al. (2015)DIAS, L. P. R. et al. Eficiência relativa de fosfatos naturais na adubação de plantio de mudas de Eucalyptus dunnii Maiden e Eucalyptus benthamii Maiden et Cambage em solo sem e com calagem. Ciência Florestal, v. 25, n. 1, p. 37-48, 2015. with a CV of 6.7% and 10.5%, in sequence, for the variable DMAP and Xavier et al. (2013)XAVIER, T. M. T. et al. Influence of water deficit and season on biomass yield from clonal stands of eucalyptus. Cerne, v. 19, n. 3, p. 473-479, 2013. with a CV of 2.46% for the variable DMRS.

CONCLUSION

  1. Data obtained with the experiments with eucalyptus in the present research did not fit the classification of Pimentel-Gomes (2009)PIMENTEL-GOMES, F. Curso de estatística experimental. 15. ed. Piracicaba: FEALQ, 2009. 451 p. and tables obtained by Costa, Seraphin and Zimmermann (2002)COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002., whereas the tables obtained by means of the method proposed by Garcia (1989)GARCIA, C. H. Tabelas para classificação do coeficiente de variação. Piracicaba: IPEF, 1989. 11 p. would be preferably used;

  2. Results obtained with the present research may serve as a basis for assisting other researchers in choosing appropriate CV classification ranges to data sets with eucalyptus seedlings.

  • 1
    Parte da Dissertação da primeira autora, apresentada ao Programa de Pós-Graduação em Agronomia da Universidade Estadual Paulista “Julio de Mesquita Filho”, Faculdade de Engenharia de Ilha Solteira, com apoio da CAPES (Código de Financiamento 001)

ACKNOWLEDGEMENTS

The authors thank the Coordination for the Improvement of Higher Education Personnel (CAPES) and the State of São Paulo Research Foundation (FAPESP grant 2015/18225-4) for the project financing.

REFERENCES

  • CARGNELUTTI FILHO, A. et al Dimensionamento amostral para avaliação de altura e diâmetro de plantas de timbaúva. Floresta e Ambiente, v. 25, n. 1, 2018.
  • CARGNELUTTI FILHO, A. et al Métodos de estimativa do tamanho ótimo de parcelas experimentais de híbridos de milho simples, triplo e duplo. Ciência Rural, v. 41, n. 9, p. 1509-1516, 2011.
  • CARGNELUTTI FILHO, A.; LÚCIO, A. D.; LOPES, S. J. Experimentação agrícola e florestal Santa Maria: UFSM, 2009. 204 p.
  • COSTA, N. H. A. D.; SERAPHIN, J. C.; ZIMMERMANN, F. J. P. Novo método de classificação de coeficientes de variação para a cultura do arroz de terras altas. Pesquisa Agropecuária Brasileira, v. 37, n. 3, p. 243-249, 2002.
  • COUTO, M. F.; PETERNELLI, L. A.; BARBOSA, M. H. P. Classification of the coefficients of variation for sugarcane crops. Ciência rural, v. 43, n. 6, p. 957-961, 2013.
  • CRUZ, E. A. et al Coeficiente de variação como medida de precisão em experimentos com tomate em ambiente protegido. Enciclopédia Biosfera, v. 8, n. 14, p. 220-233, 2012.
  • DE MESQUITA LOPES, M.; BRANCO, V. T. F. C.; SOARES, J. B. Utilização dos testes estatísticos de Kolmogorov-Smirnov e Shapiro-Wilk para verificação da normalidade para materiais de pavimentação. Transportes, v. 21, n. 1, p. 59-66, 2013.
  • DIAS, L. P. R. et al Eficiência relativa de fosfatos naturais na adubação de plantio de mudas de Eucalyptus dunnii Maiden e Eucalyptus benthamii Maiden et Cambage em solo sem e com calagem. Ciência Florestal, v. 25, n. 1, p. 37-48, 2015.
  • FERREIRA, J. P. et al Comparison of methods for classification of the coefficient of variation in papaya. Revista Ceres, v. 63, n. 2, p. 138-144, 2016.
  • FISHER, R. A. On the mathematical foundations of theoretical statistics. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, v. 222, n. 594/604, p. 309-368, 1922.
  • FONTELLES, M. J. Bioestatística aplicada à pesquisa experimental São Paulo: Livraria da Física, 2012. v. 2.
  • GARCIA, C. H. Tabelas para classificação do coeficiente de variação Piracicaba: IPEF, 1989. 11 p.
  • KOLMOGOROV, A. Sulla determinazione empirica di una lgge di distribuzione. Giornale dell'Istituto Italiano degli Attuari, v. 4, p. 83-91, 1933.
  • MEDEIROS, W. N. et al Crescimento inicial e concentração de nutrientes em clones de Eucalyptus urophylla x Eucalyptus grandis sob interferência de plantas daninhas. Ciência Florestal , v. 26, n. 1, p. 147-157, 2016.
  • MORA, F.; ARRIAGADA, O. A classification proposal for coefficients of variation in Eucalyptus experiments involving survival, growth and wood quality variables. Bragantia, v. 75, n. 3, p. 263-267, 2016.
  • PEARSON, K. X. On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling. The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, v. 50, n. 302, p. 157-175, 1900.
  • PIMENTEL-GOMES, F. Curso de estatística experimental 15. ed. Piracicaba: FEALQ, 2009. 451 p.
  • R CORE TEAM. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing, 2019.
  • SALLES, T. T. et al Modelo de Clutter na modelagem de crescimento e produção de eucalipto em sistemas de integração lavoura-pecuária-floresta. Pesquisa Agropecuária Brasileira , v. 47, n. 2, p. 253-260, 2012.
  • SANTAROSA, E.; PENTEADO JUNIOR, J. F.; GOULART, I. C. G. R. Cultivo de eucalipto em propriedades rurais: diversificação da produção e renda. Colombo: Embrapa Florestas, 2014.
  • SANTOS JUNIOR, A. et al Glyphosate drift in eucalyptus plants. Planta Daninha, v. 33, n. 3, p. 615-621, 2015.
  • SCHMILDT, E. R. et al Coeficiente de variação como medida da precisão em experimentos de alface. Revista Agro@mbiente On-line, v. 11, n. 4, p. 290-295, 2017.
  • SHAPIRO, S. S.; WILK, M. B. An analysis of variance test for normality (complete samples). Biometrika, v. 52, n. 3/4, p. 591-611, 1965.
  • SMIRNOV, N. V. Sui la distribution de w2 (Criterium de MRv Mises). Comptes Rendus , v. 202, p. 449-452, 1936.
  • STAHL, J. et al Produção de massa seca e eficiência nutricional de clones de Eucalyptus dunnii e Eucalyptus benthamii em função da adição de doses de fósforo ao solo. Ciência Florestal , v. 23, n. 2, p. 287-295, 2013.
  • TATAGIBA, S. D. et al Determinação da máxima capacidade de retenção de água no substrato para produção de mudas de eucalipto em viveiro. Floresta, v. 45, n. 4, p. 745-754, 2015.
  • TORMAN, V. B. L.; COSTER, R.; RIBOLDI, J. Normalidade de variáveis: métodos de verificação e comparação de alguns testes não-paramétricos por simulação. Revista HCPA, v. 32, n. 2, p. 227-234, 2012.
  • XAVIER, T. M. T. et al Influence of water deficit and season on biomass yield from clonal stands of eucalyptus. Cerne, v. 19, n. 3, p. 473-479, 2013.

Edited by

Editor-in-Chief: Prof. Carlos Alexandre Gomes Costa - costacag@gmail.com

Publication Dates

  • Publication in this collection
    27 Sept 2021
  • Date of issue
    2021

History

  • Received
    27 Aug 2020
  • Accepted
    31 May 2021
Universidade Federal do Ceará Av. Mister Hull, 2977 - Bloco 487, Campus do Pici, 60356-000 - Fortaleza - CE - Brasil, Tel.: (55 85) 3366-9702 / 3366-9732, Fax: (55 85) 3366-9417 - Fortaleza - CE - Brazil
E-mail: ccarev@ufc.br