Descrição do crescimento de frutos de amora-preta por modelos de regressão não linear

rbf Revista Brasileira de Fruticultura Rev. Bras. Frutic. 0100-2945 1806-9967 Sociedade Brasileira de Fruticultura Resumo A amora-preta é um pequeno fruto com várias propriedades benéficas à saúde humana, e seu cultivo é uma alternativa para os pequenos produtores. Estudar o crescimento de frutos, ao longo do tempo, é extremamente importante para entender seu desenvolvimento, auxiliando no manejo mais adequado da cultura, evitando, por exemplo, perda após a colheita, que é um dos fatores agravantes da amora-preta, por ser um fruto de curto período de prateleira. Sendo assim, as curvas de crescimento ganham destaque, nesse tipo de estudo, e a modelagem, por meio de modelos estatísticos, ajuda a entender como tal crescimento acontece. Conduziu-se este trabalho com o objetivo de ajustar modelos não lineares, duplo Logístico e duplo Gompertz, para descrever o crescimento do diâmetro de quatro cultivares da amora-preta (Brazos, Choctaw, Guarani e Tupy). As estimativas dos parâmetros foram obtidas por meio do método de mínimos quadrados, utilizando o algoritmo de Gauss-Newton, por meio do software R. A comparação dos ajustes foi feitas pelo critério de informação Akaike (AICc), desvio padrão residual (DPR) e coeficiente de determinação ajustado (R2 aj). Os modelos descreveram de forma satisfatória os dados, sendo escolhido o modelo duplo Logístico para as cultivares Brazos e Guarani, e o duplo Gompertz para as cultivares Tupy e Choctaw. Introduction Small fruits, such as blackberry, strawberry, raspberry, currant, blueberry, among others, are growing in importance, quality and production volume. The cultivation of these fruits is characterized by the low implantation and production cost and good economic return. The ideal consumption of these fruits is fresh, maintaining their physicochemical, nutritional and biological properties. However, most of the production is destined to the consumption of the processed fruit, in the form of jellies, juices, ice cream, among others. The processing of these fruits is justified by the fact that they are highly perishable, that is, they have short shelf life (TULLIO, 2013). According to Antunes et al. (2014), the fresh fruit consists of 85% water, 10% carbohydrates, minerals, vitamins A and B and calcium, with some benefits to human health, such as bleeding control, blood pressure, sedative effect, action oxidant, among others. Blackberry belongs to the genus Rubus; it is a small fruit, with good adaptation to the mild climate. Its cultivation for commercialization began in Europe, in the 16th century, later taken to the USA and arriving in Brazil in the 1970s (RASEIRA; FRANZON, 2012). This is a culture of quick financial return for producers; however, it is still little explored in Brazil, with plantations located in the states of Rio Grande do Sul, Espírito Santo, Rio de Janeiro and Southern Minas Gerais (FACHINELLO et al., 2011; ANTUNES et al., 2014). Fruit harvesting is subjectively performed by color and size (CAVALINI et al., 2006). The preservation time of fresh blackberry fruits is short and knowing the best harvesting time is very important, increasing the consumption time of fresh fruits (GONÇALVES et al., 2012; SOUZA et al., 2018). With rapid loss of post-harvest quality, there is a very large consumption limitation in its natural form; therefore, it is very important to use additional post-harvest techniques to prolong this time. Therefore, choosing the right time to objectively harvest these fruits is very important, and having knowledge about the growth pattern of fruits can be useful in this process, reducing the rapid loss of post-harvest quality, which is the main problem found for the consumption of fresh blackberry fruits (GONÇALVES et al., 2012). Some authors have studied fruit development, observing the sigmoid growth pattern and using nonlinear regression models, as, in these cases, they are more suitable than linear ones, synthesizing data information in few parameters and with biological interpretations (MARO et al. , 2012; PRADO et al., 2013a; PRADO et al., 2013b; FERNANDES et al., 2014). Nonlinear models have been used in a fruit growth study with satisfactory results, providing good estimates of parameters and practical interpretation, with green dwarf coconut (PRADO et al., 2013b), cashew (MUIANGA et al., 2016), cacao (MUNIZ et al., 2017), Asian pear (RIBEIRO et al., 2018a), pequi (RIBEIRO et al., 2018b) and strawberry fruits (DIEL, 2019). Most fruit growth studies have used simple sigmoid models. However, in the case of stone fruits, growth occurs in two development stages, and in the first, growth is slow, in the initial stage, passing to fast and stabilizing at the end, with a second stage of the same form. This behavior makes the biometric measurements of fruits to present the shape of a double “S” over time, with characteristics of a curve with two sigmoids (COOMBE, 1976; FERNANDES et al., 2017). Double sigmoid models have been satisfactorily used in growth studies (ALVAREZ and BOCHE, 1999; GOMES et al., 2008; CUNHA JUNIOR et al., 2007; BEBBINGTON el al., 2009; FAMIANE et al., 2012; BRUCKNER et al., 2012; MARTINEZ et al., 2017; FERNANDES et al., 2017; SILVA et al., 2019). The estimation of parameters in nonlinear models, in general, is performed by minimizing the sum of squares of errors, leading to a system of normal nonlinear equations, requiring the use of iterative methods to obtain estimates. Among iterative methods, the most used is Gauss-Newton (ZEVIANI et al., 2012; SOUSA et al., 2014; CARNEIRO et al., 2014; FERNANDES et al., 2015; FURTADO et al., 2019; SILVA et al., 2019). Some authors do not take into account, when adjusting nonlinear models, the interpretation of the parameter associated with the inflection point, arguing that it is a scale parameter (SILVEIRA et al., 2011; CASSIANO and SAFADI, 2015). However, Fernandes et al (2015) suggest a parameterization for Logistic and Gompertz models that have a direct interpretation for this parameter. In addition, SARI et al. (2019) pointed out that the inflection point is the point at which the growth rate reaches its maximum, which can be estimated, therefore, by equating the first derivative in relation to the model time to zero. The aim of this study was to characterize the growth of blackberry fruits of four cultivars (‘Choctaw’, ‘Tupy’, ‘Brazos’ and ‘Guarani’), in relation to their diameter using double Logistic and double Gompertz nonlinear regression models. Material and methods Data were taken from Tadeu (2015), which experiment was carried out at the Federal University of Lavras from January 2012 to January 2014. The climate is altitude tropical, with dry winters and rainy summers (ALVARES et al., 2013). The variable under study was the diameter in millimeters (mm) of blackberry of ‘Brazos’, ‘Choctaw’, ‘Guarani’ and ‘Tupy’ cultivars. Seedlings were planted in 2009 with a T-wire cordon, 60 cm apart and 80 cm in height, with 3.0 m x 0.5 m spacing. From this study, the diameters of 12 fruits of each cultivar were measured. To obtain these measurements, digital caliper (King Tools 150 mm model, Cia, São Paulo, SP) was used. The choice of fruits, for the study, was made by the development stage, that is, fruits that lost petals and anthers (called anthesis) at the same time. The Logistic double (1) and double Gompertz (2) growth models were used, described by the following expressions, with parameters taken from Fernandes et. al. (2017): (1) Y i = A 1 / { 1 + exp ⁡ B 1 - K 1 t i + ( A 2 - A 1 ) / { 1 + exp ⁡ B 2 - K 2 t i + u i (2) Y i = A 1 exp ⁡ - exp ⁡ B 1 - k 1 t i + A 2 - A 1 exp ⁡ - exp ⁡ B 2 - k 2 t i + u i in which, ui = f1ui-1 +...+ fpui-p+ ei with i = 1, 2, ..., n and n = 13 the number of times in which blackberry fruits were measured; ui is the fit residue in the i-th time; f1 is the autoregressive parameter of order 1; ui-1 is the residue of the time fit immediately before the i-th measurement; fp is the autoregressive parameter of order p; ui-pis the fit residue in p times before the i-th measurement; ei is the model’s random error, which by assumption has normal distribution. When residues are independent, parameters f1 will be null and, consequently, ui = ei (FURTADO, 2019). In equations (1) and (2), Yi defines the average values of diameters (mm) of blackberry fruits in times ti in days; A1 indicates the maximum expected value of the first phase, A2 is the asymptotic value of the second growth phase, that is, the diameter value when the fruit is fully developed; B1 and B2 represent the abscissa of the inflection point of the respective curves; K1 and K2 represent the fruit maturity indices in the respective curves; ti refers to the time of the i-th measurement, expressed in days after anthesis (FERNANDES et al., 2015, FERNANDES et al., 2017). To adjust the models, the R software (R DEVELOPMENT CORE TEAM, 2018) was used to estimate parameters using least squares and the nlme, car, lmtest and qpcr packages. The initial values for the iterative method to start the iteration was a graphical analysis with the manipulate function of the R software (ALLAIRE, 2014; DIEL et al., 2019). Studying other points on the growth curve can be extremely important, as it can help the researcher to know the precocity, acceleration and growth rate of the phenomenon under study and to know these points, the derivatives of models are used in relation to the time variable (SARI et al., 2018). In this work, the first derivative of the double Logistic and double Gompertz models was made in order to know how the growth rate behaves in the two phases and in a complementary way to the Student t test to verify if the two phases are statistically equal. The Breusch-Pagan (BREUSCH; PAGAN,1979) and Shapiro-Wilk (SHAPIRO; WILK, 1965) tests were used to verify the homogeneity of variances and normality of residues, respectively. To assess the existence of residual autocorrelation, the Durbin-Watson test (DURBIN; WATSON, 1951) was used. In addition to hypothesis tests, graphical analysis was used to verify these assumptions. Models were compared using the following criteria: i) Residual standard deviation (RSD), calculated by the expression ,DPR= QME where QME is the residual mean square; ii) Adjusted determination coefficient (R2aj), obtained by R aj 2=1-(1-R2) ( n-1)n-p , where R2 is the determination coefficient, n the number of times in which measurements were taken and p the number of model parameters. iii) Corrected Akaike information criterion (AKAIKE, 1974), given by AICc=AIC+ 2p(p+1)n-p-1,with, AIC=nln⁡ SQRn+2p where SQR is the sum of the squares of residues, p is the number of model parameters, is the sample size and ln is the natural logarithmic operator. For all analyses, the significance level adopted was 5%. Results and discussion Table 1 shows the results obtained for the Shapiro- Wilk (SW), Breusch-Pagan (BP) and Durbin-Watson (DW) tests and Figures 1, 2, 3 and 4 show the graphical analysis of the residual vector of the double Logistic and double Gompertz models to diameter data of the four blackberry cultivars. The normality assumption was met in all analyses, except for the double Gompertz model of ‘Brazos’ cultivar; thus, the modeling of this cultivar was not presented and discussed in this study. Regarding variance homogeneity, for all models and all cultivars, this assumption was not violated. Table 1 P-value for Shapiro-Wilk (SW), Breusch-Pagan (BP) and Durbin-Watson (DW) tests in the residue analysis after adjusting the double Logistic and double Gompertz models to diameter data of ‘Tupy’, ‘Brazos’, ‘Guarani’ and ‘Choctaw’ cultivars. Cultivar Model Test SW BP DW P-value P-value P-value ‘Tupy’ double Logistic 0.9900 0.1002 0.0221 double Gompertz 0.9512 0.1301 0.0178 ‘Brazos’ double Logistic 0.0970 0.6996 0.7492 double Gompertz 0.0066 0.7004 0.9980 ‘Guarani’ double Logistic 0.9674 0.0536 0.8480 double Gompertz 0.6917 0.9800 0.0591 ‘Choctaw’ double Logistic 0.5120 0.4120 0.0001 double Gompertz 0.3974 0.4976 0.0025 Figure 1 Residual analysis for the diameter (mm) of ‘Tupy’ cultivar, where (a) and (c) represent the adjusted values in relation to residues, and (b) and (d) represent the residual values in relation to the theoretical quantiles for the double Logistic and double Gompertz models, respectively. Figure 2 Residual analysis for the diameter (mm) of ‘Brazos’ cultivar, where (a) and (c) represent the adjusted values in relation to residues, and (b) and (d) represent the residual values in relation to the theoretical quantiles for the double Logistic and double Gompertz models, respectively. Figure 3 Residual analysis for the diameter (mm) of ‘Guarani’ cultivar, where (a) and (c) represent the adjusted values in relation to residues, and (b) and (d) represent the residual values in relation to the theoretical quantiles for the double Logistic and double Gompertz models, respectively. Figure 4 Residual analysis for the diameter (mm) of ‘Guarani’ cultivar, where (a) and (c) represent the adjusted values in relation to residues, and in (b) and (d) represent the residual values in relation to the theoretical quantiles for the double Logistic and double Gompertz models, respectively. Analyzing the results presented in Table 1, for the independence of the residual vector, this assumption was violated for the double Logistic model and the double Gompertz model for ‘Tupy’ and ‘Choctaw’ cultivars. Therefore, these models were again adjusted and the first-order autoregressive term was added to model this dependency. This result was similar to that found by Lúcio et al. (2015), Muniz et al. (2017) and Ribeiro et al. (2018b) in tomato, cocoa and pear growth studies, respectively. Other studies with measures taken from the same individuals over time had results similar to those found for ‘Brazos’ and ‘Guarani’ cultivars, not violating the independence assumption, such as the works of Muianga et al. (2016) and Fernandes et al. (2014), for cashew and coffee, respectively. Table 3 presents the estimates of parameters and their respective 95% confidence intervals of double Logistic and double Gompertz models for blueberry cultivars. For ‘Tupy’ cultivar, both models had good adjustments, with R2aj values greater than 97%, and the relationship between diameter growth and time after anthesis was well explained. The comparison between models was made by fit quality evaluators, which values are shown in Table 2, with the double Gompertz model taking a slight advantage in relation to the double Logistic model, with higher R2aj value and lower AICc and RSD values; however any these models could be used to explain this relationship. Table 2 Fit quality evaluators, adjusted determination coefficient ( ), corrected Akaike information criterion ( and residual standard deviation (RSD) for adjusting the double Logistic and double Gompertz models to diameter data (mm) of blackberry fruits of ‘Tupy’, ‘Brazos’, ‘Guarani’ and ‘Choctaw’ cultivars. Cultivar Model Parameters RSD ‘Tupy’ double Logistic 0.9746 32.8388 0.7317 double Gompertz (AR1) 0.9751 32.1388 0.7304 ‘Brazos’ double Logistic 0.9871 20.5857 0.4248 ‘Guarani’ double Logistic 0.9888 18.7370 0.3194 ‘Choctaw’ double Logistic 0.9587 23.6090 0.6954 double Gompertz (AR1) 0.9635 23.3321 0.6382 Table 3 Estimates of parameters and their Confidence Intervals (LI; LS) for adjusting the double Logistic and double Gompertz models to diameter data (mm) of blackberry fruits of ‘Tupy’, ‘Brazos’, ‘Guarani’ and ‘Choctaw’ cultivars. Cultivar Model Parameters phi ‘Tupy’ Double Gompertz (AR1) 15.71 3.66 0.08 21.15 29.99 0.27 0,56 IC [11.29;20.13] [0.88;6.44] [0.03;0.13] [19.35;22.94] [28.09;31.89] [0.17;0.47] --------- ‘Brazos’ Double Logistic 11.71 2.76 0.18 19.31 32.08 0.28 --------- IC [9.75;13.68] [0.86;4.67] [0.07;0.28] [17.56;21.07] [29.78;34.37] [0.09;0.46] --------- ‘Guarani’ Double Logistic 12.74 3.58 0.11 16.52 32.94 0.43 --------- IC [10.75;14.73] [1.04;6.12] [0.06;0.16] [15.71;17.33] [30.86;35.02] [0.05;0.80] --------- ‘Choctaw’ Double Gompertz (AR1) 12.97 2.10 0.09 17.23 29.40 0.36 0,70 IC [10.73;15.22] [0.62;3.58] [0.04;0.14] [15.69;18.77] [27.85;30.95] [0.16;0.56] --------- Regarding the asymptotic value parameter (A2), a conclusion made concerns the overestimation of this value, as can be seen in Figure 5, and this estimate (A2) was close to the maximum growth found in data, with expected value of 21.15 mm Another parameter of the double Gompertz model is the growth rate, and for the first and second growth phases, values were, respectively, 0.0841mm / day and 0.2714mm / day and the Student’s t-test result was significant for these values, that is,K1 and K2 differ statistically. According to this result, it could be concluded that the second growth phase is more accelerated than the first, as observed by Silva et al (2013) for peach data. Figure 5 Growth rate graph (left) and adjustment graph of the double Gompertz model (right) to fruit diameter data for ‘Tupy’ cultivar. For ‘Brazos’ cultivar, normality was violated and the abscissa parameter of the inflection point for the first growth phase (B1) of the double Gompertz model was not significant, since the confidence interval contained zero; in addition, as can be seen in Table 2 and Table 3, respectively, only the double Logistic model was used. The model chosen for this cultivar obtained determination coefficient greater than 98%, indicating good fit. Checking the information of parameters of this model, the horizontal asymptote was 19.31 mm, being a good estimate of the diameter in the final stage of fruit maturation, as shown in Figure 6. Values are related to the growth rate of both phases, which are presented in Table 3, with values of 0.1829 mm / day and 0.2814 mm / day for K1 and K2, respectively, these two phases being statistically equal by the Student’s t test (p-value <0.05). This result is inconsistent with those found by Silva et al. (2013) and Fernandes et al (2017), because in both works, phases were different. Figure 6 Growth rate graph (left) and adjustment graph of the Logistic double model (right) to fruit diameter data for ‘Brazos’ cultivar. Analyzing results in Table 3, for ‘Guarani’ cultivar, the confidence interval of the B1 parameter of the double Gompertz model contained zero, this parameter being unreliable to draw conclusions about its result. Therefore, as verified for the double Gompertz model for ‘Brazos’ cultivar, the results of this model were not discussed below. As a result, it was not necessary to make a comparison with the double Logistic model by AICc and RSD evaluators, the double Logistic model being the only model tested and adequate to describe the diameter growth of the ‘Guarani’ cultivar. Even though there is no way to compare the models proposed in this work, Table 2 shows that the determination coefficient of the double Logistic model, the only one suitable for these data, was greater than 98.5%, suggesting that this model has excellent fit quality in the description of the diameter growth of this cultivar. Analyzing the results of the ‘Guarani’ cultivar, the asymptotic value (A2) had good estimate, with value close to the maximum growth found in data, with estimate of 16.52 mm, as can be seen in Figure 7. The values of estimates for both growth phases were 0.11 for the first phase and 0.43 for the second, with the Student t test being significant (p-value> 0.05), that is, phases are statistically different. Figure 7 Growth rate graph (left) and adjustment graph of the Logistic double model (right) to fruit diameter data for ‘Guarani’ cultivar. For ‘Choctaw’ cultivar, estimates and the 95% confidence interval are in Table 3, showing that all confidence intervals for parameters did not contain zero. Both models had good adjustments, with R2 aj showing values greater than 95%. As observed for ‘Tupy’ cultivar, any of the models could be chosen to characterize this growth, but the double Gompertz model presented lower AICc and RSD values and higher R2aj value, as can be seen in Table 2, which was chosen to comment about these results. The horizontal asymptote value (A2)did not overestimate this value, with estimate 17.2366 mm, this verification being made in Figure 8. Growth rate estimates were 0.0904 mm/ day and 0.3609 mm / day for the first and second growth phases, respectively, and Student’s t test was significant (p-value> 0.05), as in the works by Silva et al. (2013) and Fernandes et al. (2017), the second growth phase is faster than the first one. Figure 8 Growth rate graph (left) and adjustment graph of the double Gompertz model (right) to fruit diameter data for ‘Choctaw’ cultivar. The inflection point of the first growth phase occurred before the fifth day after anthesis, as can be seen in Figures 1, 2, 3 and 4 (left) and the second point occurred, approximately, on the thirtieth day after anthesis and other information was in relation to the growth pattern, as Figures 1, 2, 3 and 4 (left) corroborate pattern found in Figures 1, 2, 3 and 4 (right), that is, the growth of cultivars under study follows a characteristic of double sigmoid. Martins et al. (2019) evaluated the phenology of the ‘Tupy’ cultivar and concluded that it has the same double sigmoid growth pattern. For other fruits, such as coffee, Fernandes et al. (2017) and peach, Silva et al. (2013) and Silva et al. (2019) also observed the same curve shape in their studies. The model chosen to represent the nonlinear relationship between time and days after anthesis of the ‘Tupy’ cultivar presented final diameter growth estimate of 21.1522 mm; Martins et al. (2019) and Campagnolo and Pio (2012) observed value within the confidence interval of estimates found, and in this work, these values were, respectively, 21.40 mm and 21.00 mm. However, Ferreira et al. (2018) found value of 24.00 mm, which can be explained by the fact that these authors used castor bean cake as fertilizer. Conclusions The fruit diameter growth pattern of blackberry cultivars under study was the double sigmoidal, adequately described by nonlinear double Logistic and double Gompertz models. For ‘Tupy’ and ‘Choctaw’ cultivars, the most appropriate model was the Logistic double and for ‘Brazos’ and ‘Guarani’ cultivars, the double Gompertz was the most appropriate. ‘Tupy’ and ‘Choctaw’ cultivars presented the second growth phase more quickly than the first, while for ‘Brazos’ and ‘Guarani’ cultivars, these phases were statistically equal. At harvest time, fruits of the different cultivars presented diameters between 17.2 mm and 21.1 mm. AKAIKE, H. A new look at the statistical model identification. IEEE Transactions on Automatic Control, Boston, v.19, n.6, p. 716-723, 1974. (1) AKAIKE H. A new look at the statistical model identification. IEEE Transactions on Automatic Control Boston 19 6 1974 716 723 ALLAIRE, J.J. Manipulate: interactive plots for RStudio. R package version 1.0.1. Boston, 2014. Disponível em: https://rdrr.io/cran/manipulate/. ALLAIRE J.J. Manipulate: interactive plots for RStudio. 2014 Disponível em: https://rdrr.io/cran/manipulate/ 10 janeiro 2015 ALVARES, C.A.; STAPE, J.L.; SENTELHAS, P.C.; GONÇALVES, J.L.M.; SPAROVEK, G. Köppen's climate classification map for Brazil. Meteorologische Zeitschrift, Stuttgart, v.22, n.6, p.711-728, 2013. ALVARES C.A. STAPE J.L. SENTELHAS P.C. Köppen's climate classification map for Brazil. Meteorologische Zeitschrift Stuttgart 22 6 2013 711 728 ALVAREZ, A.; BOCHE, S. Modelos matemáticos para describir crescimientos doble sigmoideos en frutos de un nectarin tardio. Agrosur Journal, Valdivia, v.27, n.1, p.1-8, 1999. ALVAREZ A. BOCHE S. Modelos matemáticos para describir crescimientos doble sigmoideos en frutos de un nectarin tardio. Agrosur Journal Valdivia 27 1 1999 1 8 ANTUNES, L.E.C.; PEREIRA, I.S.; VIGNOLO, G.K.; GONÇALVES, M.A. Produção de amoreira-preta no Brasil. Revista Brasileira de Fruticultura, Jaboticabal, v.36, n.1, p.100-111, 2014. ANTUNES L.E.C. PEREIRA I.S. VIGNOLO G.K. Produção de amoreira-preta no Brasil. Revista Brasileira de Fruticultura Jaboticabal 36 1 2014 100 111 BEBBINGTON, M.; HALL, A.J.; LAI, C.D.; ZITIKIS, R. Dynamics and phases of kiwifruit (Actinidia deliciosa) growth curves. New Zealand Journal of Crop and Horticultural Science, Abingdon, v.37, n.3, p.179-188, 2009. BEBBINGTON M. HALL A.J. LAI C.D. Dynamics and phases of kiwifruit (Actinidia deliciosa) growth curves. New Zealand Journal of Crop and Horticultural Science Abingdon 37 3 2009 179 188 BREUSCH, T.S.; PAGAN, A.R. A simple test for heterocedasticity and random coeficiente variation. Econometrica, New York, v.47, n.5, 1287-1294, 1979. BREUSCH T.S. PAGAN A.R. A simple test for heterocedasticity and random coeficiente variation. Econometrica New York 47 5 1979 1287 1294 BRUCKNER, C.H.; SILVA, D.F.P.; SILVA, J.O.C.; MATIAS, R.G.P.; RIBEIRO, M.R. Curva de crescimento e padrão respiratório de frutos de genótipos de pessegueiro em região de clima subtropical. Revista Brasileira de Fruticultura, Jaboticabal, v.35, n.2, p.235–241, 2012. BRUCKNER C.H. SILVA D.F.P. SILVA J.O.C. Curva de crescimento e padrão respiratório de frutos de genótipos de pessegueiro em região de clima subtropical. Revista Brasileira de Fruticultura Jaboticabal 35 2 2012 235 241 CAMPAGNOLO, M.A.; PIO, R. Produção da amoreira-preta'Tupy'sob diferentes épocas de poda. Ciência Rural, Santa Maria, v.42, n .2, p.225-231, 2012. CAMPAGNOLO M.A. PIO R. Produção da amoreira-preta'Tupy'sob diferentes épocas de poda. Ciência Rural Santa Maria 42 2 2012 225 231 CARNEIRO, A.P.S.; MUNIZ, J.A.; CARNEIRO, P.L.S.; MALHADO, C.H.M.; FILHO, R.M.; SILVA, F.F. Identidade de modelos não lineares para comparar curvas de crescimento de bovinos da raça Tabapuã. Pesquisa Agropecuária Brasileira, Brasília, DF, v.49, n.1, p.57-62, 2014. CARNEIRO A.P.S. MUNIZ J.A. CARNEIRO P.L.S. Identidade de modelos não lineares para comparar curvas de crescimento de bovinos da raça Tabapuã. Pesquisa Agropecuária Brasileira J.A. 49 1 2014 57 62 CASSIANO, F.R.; SÁFADI, T. Modelos de crescimento animal para tempos irregulares. Pesquisa Agropecuária Brasileira, Brasília, DF, v.50, n.11, p.1114-1119, 2015. CASSIANO F.R. SÁFADI T. Modelos de crescimento animal para tempos irregulares. Pesquisa Agropecuária Brasileira Brasília 50 11 2015 1114 1119 CAVALINI, F.C.; JACOMINO, A.P.; LOCHOSKI, M.A.; KLUGE, R.A.; ORTEGA, E.M.M. Maturity indexes for’kumagai’and’paluma’guavas. Revista Brasileira de Fruticultura, Jaboticabal, v.28, n.2, p.176–179, 2006. CAVALINI F.C. JACOMINO A.P. LOCHOSKI M.A. Maturity indexes for’kumagai’and’paluma’guavas. Revista Brasileira de Fruticultura Jaboticabal 28 2 2006 176 179 COOMBE, B.G. The development of fleshy fruits. Annual Review of Plant Physiology, Palo Alto, v.27, p.507-528, 1976. COOMBE B.G. The development of fleshy fruits. Annual Review of Plant Physiology Palo Alto 27 1976 507 528 CUNHA JUNIOR, L.C.; MATTIUZ, B.; DURIGAN, M.F.B.; MARTINS, R.N. Caracterização da curva de maturação de pêssegos’ aurora-1’, na região de jaboticabal-sp. Revista Brasileira de Fruticultura, Jaboticabal, v.29, n.3, p.661–665, 2007. CUNHA L.C. JUNIOR MATTIUZ B. DURIGAN M.F.B. Caracterização da curva de maturação de pêssegos’ aurora-1’, na região de jaboticabal-sp. Revista Brasileira de Fruticultura Jaboticabal 29 3 2007 661 665 DIEL, M.I.; SARI, B.G.; KRYSCZUN, D.K.; OLIVOTO, T.; PINHEIRO, M.V.M.; MEIRA, D.; SHIMIDT, D.; LÚCIO, A.D.C. Nonlinear regression for description of strawberry (Fragaria x ananassa) production. The Journal of Horticultural Science and Biotechnology, Londres, v.94, n.2, p.259-273, 2019. DIEL M.I. SARI B.G. KRYSCZUN D.K. Nonlinear regression for description of strawberry (Fragaria x ananassa) production. The Journal of Horticultural Science and Biotechnology Londres 94 2 2019 259 273 DURBIN, J.; WATSON, G.S. Testing for serial correlation in last squares regression. Biometrika, London, v.38, n.1/2, p.159-177, 1951. DURBIN J. WATSON G.S. Testing for serial correlation in last squares regression. Biometrika London 38 1/2 1951 159 177 FACHINELLO, J.C.; PASA, M.S.; SHIMITZ, J.D.; BETEMPS, D.L. Situação e perspectivas da fruticultura de clima temperado no brasil. Revista Brasileira de Fruticultura, Jaboticabal, v.33, n.1, p.109–120, 2011. FACHINELLO J.C. PASA M.S. SHIMITZ J.D. Situação e perspectivas da fruticultura de clima temperado no brasil. Revista Brasileira de Fruticultura Jaboticabal 33 1 2011 109 120 FAMIANI, F.; CASULLI, V.; BALDICCHI, A.; BATTISTELLI, A.; MOSCATELLO, S.; WALKER, R.P. Development and metabolism of the fruit seed of the Japanese plum Ozark Premier. Journal of Plant Physiology, Freiburg, v.169, p.551-560, 2012. FAMIANI F. CASULLI V. BALDICCHI A. Development and metabolism of the fruit seed of the Japanese plum Ozark Premier. Journal of Plant Physiology Freiburg 169 2012 551 560 FERNANDES, T.J.; MUNIZ, J.A.; PEREIRA, A.A.; MUNIZ, F.R.; MUIANGA, C.A. Parameterization effects in nonlinear models to describe growth curves. Acta Scientiarum, Maringá, v.37, n.4, p.397, 2015. FERNANDES T.J. MUNIZ J.A. PEREIRA A.A. Parameterization effects in nonlinear models to describe growth curves. Acta Scientiarum Maringá 37 4 2015 397 FERNANDES, T.J.; PEREIRA, A.A.; MUNIZ, J.A. Double sigmoidal models describing the growth of coffee berries. Ciência Rural, Santa Maria, v.47, n.8, 2017. FERNANDES T.J. PEREIRA A.A. MUNIZ J.A. Double sigmoidal models describing the growth of coffee berries. Ciência Rural Santa Maria 47 8 2017 FERNANDES, T.J.; PEREIRA, A.A.; MUNIZ, J.A.; SAVIAN, T.V. Seleção de modelos não lineares para a descrição das curvas de crescimento do fruto do cafeeiro. Coffee Science, Lavras, v.9, n.2, p.207-2015, 2014. FERNANDES T.J. PEREIRA A.A. MUNIZ J.A. Seleção de modelos não lineares para a descrição das curvas de crescimento do fruto do cafeeiro. Coffee Science Lavras 9 2 2014 207 2015 FERREIRA, L.V.; COCCO, C.; FINKENAUER, D.; PICOLOTTO, L.; ANTUNES, L.E.C. Adubação com torta de mamona sobre o crescimento e produção da amoreira-preta. Cultura Agronômica: Revista de Ciências Agronômicas, Ilha Solteira, v.27, n.1, p.34-43, 2018. FERREIRA L.V. COCCO C. FINKENAUER D. Adubação com torta de mamona sobre o crescimento e produção da amoreira-preta. Cultura Agronômica: Revista de Ciências Agronômicas Ilha Solteira 27 1 2018 34 43 FURTADO, T.D.R.; MUNIZ, J.A.; SILVA, E.M.; FERNANDES, J.G. Drying kinetics of jabuticaba pulp by regression models. Revista Brasileira de Fruticultura, Jaboticabal, v.41, n.1, 2019. FURTADO T.D.R. MUNIZ J.A. SILVA E.M. Drying kinetics of jabuticaba pulp by regression models. Revista Brasileira de Fruticultura Jaboticabal 41 1 2019 GOMES, R.; PINHEIRO, M.C.B.; LIMA, H.A. Fenologia reprodutiva de quatro espécies de Sapotaceae da restinga de Marica, RJ. Revista Brasileira de Botânica, São Paulo.v.31, n.4, p.679-687, 2008. GOMES R. PINHEIRO M.C.B. LIMA H.A. Fenologia reprodutiva de quatro espécies de Sapotaceae da restinga de Marica, RJ. Revista Brasileira de Botânica São Paulo 31 4 2008 679 687 GONÇALVES, E.D.; PIMENTEL, R.M.A.; LIMA, L.C.O.; CASTRICINI, A.; ZAMBON, C.R.; ANTUNES, L.E.C.; TREVISAN, R.Manutenção da qualidade pós-colheita das pequenas frutas.Informe Agropecuário, Belo Horizonte, v.33, n.268, p.89-95, 2012. GONÇALVES E.D. PIMENTEL R.M.A. LIMA L.C.O. Manutenção da qualidade pós-colheita das pequenas frutas. Informe Agropecuário Belo Horizonte 33 268 2012 89 95 LUCIO, A.D.; SARI, B.G.; RODRIGUES, M.; BEVILAQUA, L.M.; VOSS, H.M.G.; COPETTI, D.; FAÉ, M. Modelos não-lineares para a estimativa da produção de tomate do tipo cereja. Ciência Rural, Santa Maria, v.46, n.2, p.233-241, 2015. LUCIO A.D. SARI B.G. RODRIGUES M. Modelos não-lineares para a estimativa da produção de tomate do tipo cereja. Ciência Rural Santa Maria 46 2 2015 233 241 MARO, L.A.C.; PIO, R.; SILVA, T.C.; PATTO, L.S. Ciclo de produção de cultivares de framboeseiras (Rubusidaeus) submetidas a poda drástica nas condições do sul de Minas Gerais. Revista Brasileira de Fruticultura, Jaboticabal, v.34, n.2, p.435-441, 2012. MARO L.A.C. PIO R. SILVA T.C. Ciclo de produção de cultivares de framboeseiras (Rubusidaeus) submetidas a poda drástica nas condições do sul de Minas Gerais. Revista Brasileira de Fruticultura Jaboticabal 34 2 2012 435 441 MARTÍNEZ, J.D.; HERNANDEZ, I.C.; REYES, N.G.; NAJERA, J.B.P.; GUERECA M.C.G.; PEREZ, E.C. Growth models of peach fruit prunus persica (l) in three handling systems. Interciencia, Caracas, v.42, n.9, p.596 – 602, 2017. MARTÍNEZ J.D. HERNANDEZ I.C. REYES N.G. Growth models of peach fruit prunus persica (l) in three handling systems. Interciencia Caracas 42 9 2017 596 602 MARTINS, W.A., SANTOS, S.C., JARA, R.S., SOUZA, J.L.A.C., GALVÃO, J.R.; BISCARO, G.A. Fenologia e demanda térmica de amoreira-preta cv.Tupy. Revista de Ciências Agrárias, Lisboa, v.42, n.3, p.720-730, 2019. ( MARTINS W.A. SANTOS S.C. JARA R.S. Fenologia e demanda térmica de amoreira-preta cv.Tupy. Revista de Ciências Agrárias Lisboa 42 3 2019 720 730 MUIANGA, C.A.; MUNIZ, J.A.; NASCIMENTO, M.S.; FERNANDES, T.J.; SAVIAN, T.V. Descrição da curva de crescimento de frutos do cajueiro por modelos não lineares. Revista Brasileira de Fruticultura, Jaboticabal, v.38, n.1, p.22-32, 2016. MUIANGA C.A. MUNIZ J.A. NASCIMENTO M.S. Descrição da curva de crescimento de frutos do cajueiro por modelos não lineares. Revista Brasileira de Fruticultura Jaboticabal 38 1 2016 22 32 MUNIZ, J.A.; NASCIMENTO, M.S.; FERNANDES, T.J. Nonlinear models for description of cacao fruit growth assumpition violations. Revista Caatinga, Mossoró, v.30, n.1, p.250-257, 2017. MUNIZ J.A. NASCIMENTO M.S. FERNANDES T.J. Nonlinear models for description of cacao fruit growth assumpition violations. Revista Caatinga Mossoró 30 1 2017 250 257 PRADO, T.K.L.; MUNIZ, J.A.; SAVIAN, T.V.; SÁFADI, T. Ajuste do modelo logístico na descrição do crescimento de frutos de coqueiro anão por meio de algoritmos iterativos MCMC. Revista Brasileira de Biometria, Lavras, v.31, n.2, p.216-232, 2013b. PRADO T.K.L. MUNIZ J.A. SAVIAN T.V. Ajuste do modelo logístico na descrição do crescimento de frutos de coqueiro anão por meio de algoritmos iterativos MCMC. Revista Brasileira de Biometria Lavras 31 2 2013b 216 232 PRADO, T.K.L.; SAVIAN, T.V.; MUNIZ, J.A. Ajuste dos modelos Gompertz e Logístico aos dados de crescimento de frutos de coqueiro anão verde. Ciência Rural, Santa Maria,v.43, n.5, p.803-809, 2013a. PRADO T.K.L. SAVIAN T.V. MUNIZ J.A. Ajuste dos modelos Gompertz e Logístico aos dados de crescimento de frutos de coqueiro anão verde. Ciência Rural Santa Maria 43 5 2013a 803 809 R Core Team. R: a language and environment for statistical computing.Vienna, 2017. Disponível em: https://www.R-project.org/. R Core Team R: a language and environment for statistical computing. Vienna 2017 RASEIRA, M.C.B.; FRANZON, R.C. Melhoramento genético e cultivares de amora-preta e mirtilo. Informe Agropecuário, Belo Horizonte, v.33, n.268, p.11–20, 2012. RASEIRA M.C.B. FRANZON R.C. Melhoramento genético e cultivares de amora-preta e mirtilo. Informe Agropecuário Belo Horizonte 33 268 2012 11 20 RIBEIRO, T.D.; MATTOS, R.V.P.; MORAIS, A.R.; MUNIZ, J.A. Description of the growth of pequi fruits by nonlinear models. Revista Brasileira de Fruticultura, Jaboticabal, v.40, n.4, 2018b. RIBEIRO T.D. MATTOS R.V.P. MORAIS A.R. Description of the growth of pequi fruits by nonlinear models. Revista Brasileira de Fruticultura Jaboticabal 40 4 2018b RIBEIRO, T.D.; SAVIAN, T.V.; FERNANDES, T.J.; MUNIZ, J.A. The use of the nonlinear models in the growth of pears of ‘Shinseiki’cultivar. Ciência Rural, Santa Maria, v.48, n.1, 2018a. RIBEIRO T.D. SAVIAN T.V. FERNANDES T.J. The use of the nonlinear models in the growth of pears of ‘Shinseiki’cultivar. Ciência Rural Santa Maria 48 1 2018a SARI, B.G.; LÚCIO, A.D.C.; SAMANTA, C.S.; SAVIAN, T.V. Describing tomato plant production using growth models. Scientia Horticulturae, Amsterdam, v.246, p.146-154, 2019. SARI B.G. LÚCIO A.D.C. SAMANTA C.S. Describing tomato plant production using growth models. Scientia Horticulturae Amsterdam 246 2019 146 154 SHAPIRO, S.S.; WILK, M.B. An analysis of variance test for normality (complete samples). Biometrika, Cambridge, v.52, n.3-4, p.591-611, 1965 SHAPIRO S.S. WILK M.B. An analysis of variance test for normality (complete samples). Biometrika Cambridge 52 3-4 1965 591 611 SILVA, D.F.P.; SILVA, J.O.C.; MATIAS, R.G.P.; RIBEIRO, M.R.; BRUCKNER, C.H. Growth curve and respiratory pattern of genotypes of peach fruit in subtropical region. Revista Brasileira de Fruticultura , Jaboticabal, v.35, n.2, p.642-649, 2013. SILVA D.F.P. SILVA J.O.C. MATIAS R.G.P. Growth curve and respiratory pattern of genotypes of peach fruit in subtropical region. Revista Brasileira de Fruticultura Jaboticabal 35 2 2013 642 649 SILVA, É.M.; SILVA, V.F.; FERNANDES, F.A.; MUNIZ, J.A.; FERNANDES, T.J. O crescimento de frutos de pêssegos caracterizados por modelos de regressão não lineares. Sigmae, Alfenas, v.8, n.2, p.290-294, 2019. SILVA É.M. SILVA V.F. FERNANDES F.A. O crescimento de frutos de pêssegos caracterizados por modelos de regressão não lineares. Sigmae Alfenas 8 2 2019 290 294 SILVEIRA, F.G.SILVA, F.F.; CARNEIRO P.L.S.; MALHADO, C.H.M.; MUNIZ, J.A. Análise de agrupamento na seleção de modelos de regressão não-lineares para curvas de crescimento de ovinos cruzados. Ciência Rural, Santa Maria, v.41, n.4, p.692-698, 2011. SILVEIRA F.G. SILVA F.F. CARNEIRO P.L.S. Análise de agrupamento na seleção de modelos de regressão não-lineares para curvas de crescimento de ovinos cruzados. Ciência Rural Santa Maria 41 4 2011 692 698 SOUSA, I.A, KUNZLE NETO, J.E.; MUNIZ, J.A.; GUIMARÃES, R.M.; SAVIAN, T.V.; MUNIZ, F.R. Fitting nonlinear autoregressive models to describe coffee seed germination. Ciência Rural, Santa Maria, v.44, p.2016-2021, 2014. SOUSA I.A. KUNZLE J.E. NETO MUNIZ J.A. Fitting nonlinear autoregressive models to describe coffee seed germination. Ciência Rural Santa Maria 44 2014 2016 2021 SOUZA, A.V.; SILVA, V.; RIBEIRO, M.; VIEITES, R.L. Evolução da coloração de frutos e geleias de amora-preta ao longo do período de armazenamento. Revista Iberoamericana de Tecnología Postcosecha, Cidade do México, v.19, n.2, 2018. SOUZA A.V. SILVA V. RIBEIRO M. Evolução da coloração de frutos e geleias de amora-preta ao longo do período de armazenamento. Revista Iberoamericana de Tecnología Postcosecha Cidade do México 19 2 2018 TADEU, M.H. Poda drástica de verão na produção e avaliação de doenças em amoreiras-pretas em região subtropical. 2014. Dissertação (Mestrado em Agronomia) – Universidade Federal de Lavras, Lavras, 2014. TADEU M.H. Poda drástica de verão na produção e avaliação de doenças em amoreiras-pretas em região subtropical. Dissertação (Mestrado em Agronomia) Universidade Federal de Lavras Lavras 2014 TULLIO, L.; AYUB, R.A. Produção da amora-preta cv tupy, em função da intensidade da poda. Semina: Ciências Agrárias, Londrina, v.34, n.3, p.1147-1152, 2013. TULLIO L. AYUB R.A. Produção da amora-preta cv tupy, em função da intensidade da poda. Semina: Ciências Agrárias Londrina 34 3 2013 1147 1152 ZEVIANI, W.M.; SILVA, C.A; CARNEIRO, W.J.O.; MUNIZ, J.A. Modelos não lineares para a liberação de potássio de estercos animais em latossolos. Ciência Rural, Santa Maria, v.42, n.10, p.1789-1796, 2012. ZEVIANI W.M. SILVA C.A. CARNEIRO W.J.O. Modelos não lineares para a liberação de potássio de estercos animais em latossolos. Ciência Rural Santa Maria 42 10 2012 1789 1796