# ABSTRACT:

Assessing sugarcane (Saccharum spp.) stalk growth helps to adequately manage the phenological stages of the crop. The aim of this study was to describe the height-growth curve of four sugarcane varieties (RB92579, RB93509, RB931530 and SP79-1011), in irrigated plant-cane and ratoon cane plantations, using the Logistic and Gompertz nonlinear models, while considering all deviations from assumptions. The model parameters were estimated based on the least squares method using the Gauss-Newton algorithm. To select the most suitable model, nonlinear measures, adjusted coefficient of determination (R2 adj), residual standard deviation (RSD), and corrected Akaike information criterion (AICc) were used. Based on the best models, stalk height growth rates and crop phenological stages were determined using critical points. All tests were performed in the free software environment for statistical computing and graphics, R. In general, the Logistic and Gompertz models without AR(1) better described the plant-cane and ratoon cane stalk height, respectively. All varieties showed early growth, and the RB92579 variety presented higher rates in both cycles.

Key words:
logistic model; gompertz model; growth rates

# RESUMO:

Palavras-chave:
modelo logístico; modelo gompertz; taxas de crescimento

# INTRODUCTION:

Sugarcane (Saccharum spp.) is a crop grown in more than 100 countries with tropical and subtropical climates and is particularly important for worldwide agribusinesses. Brazil is the largest sugarcane producer, with 8.59 million hectares and a total production of 620.44 million tons, which is ahead of India, China, and Thailand. Brazilian production is equivalent to 39% of the total sugarcane produced in the world, and the Central South is the largest Brazilian sugarcane-producing region, with over 92% of Brazilian sugarcane production, of which the state of São Paulo accounts for 58% (CONAB, 2019CONAB (COMPANHIA NACIONAL DE ABASTECIMENTO). Acompanhamento da safra brasileira: cana-de-açúcar, segundo levantamento - safra 2018-2019. 2019. Available from: <Available from: https://www.conab.gov.br/info-agro/safras/cana/boletim-da-safra-de-cana-de-acucar >. Acessed: Mar. 9, 2019.
https://www.conab.gov.br/info-agro/safra...
).

This crop generates jobs/income and is the main raw material used in sugarcane, ethanol, and cachaça production. The low yield of sugarcane varieties affects the amount of sugarcane derivatives. Thus, studying crop growth and adapting variety management to phenological stages becomes important, particularly for maximizing plant growth, which coincides with periods of increased water availability. Water deficiency, in addition to limiting crop yields, reduces the genetic potential of the crop.

The yield of sugarcane varieties is determined based on the morphophysiological variables of the plant. TOPPA et al. (2010TOPPA, E. V. B. et al. Physiology development in the vegetative stage of sugarcane. Pesquisa Aplicada & Agrotecnologia, v.3, n.2, p.169-186. 2010. Available from: <Available from: https://www.researchgate.net/publication/277820032_Physiology_development_in_the_vegetative_stage_of_sugarcane >. Accessed: May, 27, 2019. doi: 10.5777/paet.v3i2.1064.
https://www.researchgate.net/publication...
) considered stalk height and diameter as the main variables, thereby conducting studies to understand the dynamics of their growth. According to ARCHONTOULIS & MIGUEZ (2015ARCHONTOULIS, S. V; MIGUEZ, F. E. Nonlinear regression models and applications in agricultural research. Agronomy Journal, v.107, n.2, p.786-798. 2015. Available from: <Available from: https://dl.sciencesocieties.org/publications/aj/abstracts/107/2/786 >. Accessed: Feb. 10, 2019. doi: 10.2134/agronj2012.0506.
https://dl.sciencesocieties.org/publicat...
), the most widely used standard method is quantitative growth analysis, often performed using multiple linear and polynomial regression models. Plant growth shows a sigmoid pattern: slow at first, undergoing a faster phase, and eventually slowing down again. Therefore, growth should be described using nonlinear models, which have better horizontal asymptotes and fewer parameters. Additionally, nonlinear models provide a better fit than linear models and their estimates may have biological and practical interpretations.

Several authors have used these models to describe plant and fruit growth (MUIANGA et al., 2016MUIANGA, C. A. et al.Description of the growth curve of cashew fruits in nonlinear models. Revista Brasileira de Fruticultura, v.38, n.1, p.22-32. 2016. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S0100-29452016000100022&script=sci_abstract&tlng=pt >. Accessed: May, 25, 2019. doi: 10.1590/0100-2945-295/14.
http://www.scielo.br/scielo.php?pid=S010...
; PEREIRA et al., 2016PEREIRA, A. A. et al. Modeling of the diameter of the canopy pruned coffee tree cultivated in different densities and water regimes.Lavras: Coffee Science, v.11, n.4, p.495-501. 2016. Available from: <Available from: http://www.coffeescience.ufla.br/index.php/Coffeescience/article/view/1145/pdf_1145 >. Accessed: May, 25, 2019. doi: 10.25186/cs.v11i4.1145.
http://www.coffeescience.ufla.br/index.p...
; RIBEIRO et al., 2018aRIBEIRO, T. D. et al.Description of the growth of pequi fruits by nonlinear models. Jaboticabal: Revista Brasileira de Fruticultura, v.40, n.4, e-949. 2018a. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-29452018000400705 >. Accessed: May, 13, 2019. doi: 10.1590/0100-29452018949.
http://www.scielo.br/scielo.php?script=s...
; RIBEIRO et al., 2018bRIBEIRO, T. D. et al. The use of the nonlinear models in the growth of pears of ‘Shinseiki’ cultivar. Santa Maria: Ciência Rural, v.48, n.1, e20161097. 2018b. Availale from: < Availale from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-84782018000100202 >. Accessed: May, 25, 2019. doi: 10.1590/0103-8478cr20161097.
http://www.scielo.br/scielo.php?script=s...
). The most frequently used models for sugarcane are the Logistic and Gompertz models, commonly adjusted for stalk growth (SILVA et al., 2012SILVA, T. G. F. da et al. Biometrics of the sugarcane shoot during irrigated ratoon cycle in the Submedio of the Vale do São Francisco. Revista Ciência Agronômica, v.43, n.3, p.500-509. 2012. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1806-66902012000300012&script=sci_arttext >. Accessed: May, 27, 2019. doi: 10.1590/S1806-66902012000300012.
http://www.scielo.br/scielo.php?pid=S180...
; BATISTA et al., 2013BATISTA, E. L. da S. et al.Modeling of growth of sugarcane cultivars duringthe crop formation period. Revista Brasileira de Engenharia Agrícola e Ambiental, v.17, n.10, p.1080-1087. 2013. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1415-43662013001000009&script=sci_arttext >. Acessed: May, 26, 2019. doi: 10.1590/S1415-43662013001000009.
http://www.scielo.br/scielo.php?pid=S141...
). Based on these models, plant growth rates can be determined (GAVA et al., 2001GAVA, G. J. de C. et al. Growth and accumulation of nitrogen by sugarcane cultivated in soil covered with cane trash.Pesquisa Agropecuária Brasileira, v.36, n.11, p.1347-1354. 2001. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S0100-204X2001001100004&script=sci_arttext&tlng=pt >. Accessed: May, 23, 2019. doi: 10.1590/S0100-204X2001001100004.
http://www.scielo.br/scielo.php?pid=S010...
; SANTOS et al., 2009SANTOS, V. R. dos et al. Growth and yield of sugarcane under different phosphorus sources. Revista Brasileira de Engenharia Agrícola e Ambiental, v.13, n.4, p.389-396. 2009. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-43662009000400004 >. Accessed: May, 27, 2019. doi: 10.1590/S1415-43662009000400004.
http://www.scielo.br/scielo.php?script=s...
; PEDULA et al., 2016PEDULA, R. O. et al.Growth analysis of sugarcane inoculated with diazotrophic bacteria and nitrogen fertilization. African Journal of Agricultural Research, v.11, n.30, p.2786-2795. 2016. Available from: <Available from: https://academicjournals.org/journal/AJAR/article-abstract/4E1ABBB59770 >. Accessed: May, 28, 2019. doi: 10.5897/AJAR2016.11141.
), indicating increase in growth over time. These rates help determined the duration of the phenological stages, model curve inflection point (IP), points of maximum acceleration (PMA), minimum deceleration point (PMD), and asymptotic deceleration (PDA). According to OLIVEIRA et al. (2010OLIVEIRA, E. C. A. de et al. Growth and dry matter production in sugarcane varieties grown under full irrigation.Revista Brasileira de Engenharia Agrícola e Ambiental, v.14, n.9, p.951-960. 2010. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-43662010000900007 >. Accessed: Apr. 24, 2019. doi: 10.1590/S1415-43662010000900007.
http://www.scielo.br/scielo.php?script=s...
), the IP is highly important for sugarcane because it indicates the date when the greatest plant growth has been achieved and from when growth starts to increase at decreasing rates until nearing zero, indicating that the plant is already in the maturation phase. According to MISCHAN et al. (2011MISCHAN, M. M. et al. Determination of a point sufficiently close to the asymptote in nonlinear growth functions. Scientia Agricola, v.68, n.1, p.109-114. 2011. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S0103-90162011000100016&script=sci_arttext&tlng=pt >. Accessed: May, 10, 2019. doi: 10.1590/S0103-90162011000100016.
http://www.scielo.br/scielo.php?pid=S010...
), the PDA and PMD indicate the end of the slow and exponential growth phases; respectively, whereas the PDA indicated growth stabilization.

According to MUIANGA et al. (2016MUIANGA, C. A. et al.Description of the growth curve of cashew fruits in nonlinear models. Revista Brasileira de Fruticultura, v.38, n.1, p.22-32. 2016. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S0100-29452016000100022&script=sci_abstract&tlng=pt >. Accessed: May, 25, 2019. doi: 10.1590/0100-2945-295/14.
http://www.scielo.br/scielo.php?pid=S010...
) and RIBEIRO et al. (2018bRIBEIRO, T. D. et al. The use of the nonlinear models in the growth of pears of ‘Shinseiki’ cultivar. Santa Maria: Ciência Rural, v.48, n.1, e20161097. 2018b. Availale from: < Availale from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-84782018000100202 >. Accessed: May, 25, 2019. doi: 10.1590/0103-8478cr20161097.
http://www.scielo.br/scielo.php?script=s...
), the nonlinear model parameters are not estimated based on the least squares method. However, because the system of equations has no explicit solution, iterative methods are used to obtain the estimates, most commonly the Gauss-Newton algorithm. According to these authors, once the model parameters are estimated, the residuals are analyzed, usually applying the Shapiro-Wilk, Durbin-Watson, and Breush-Pagan tests, which assessed the assumptions of normality, independence of residuals, and variance homogeneity, respectively. When assessing the failure to fulfill an assumption, weighted or generalized least squares methods are used to estimate the parameters (MUNIZ et al., 2017MUNIZ, J. A. et. al. Nonlinear models for description of cacao fruit growth with assumption violations. Revista Caatinga, v.30, n.1, p.250-257. 2017. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1983-21252017000100250 >. Accessed: May, 13, 2019. doi: 10.1590/1983-21252017v30n128rc.
http://www.scielo.br/scielo.php?script=s...
). Through nonlinear measures, the nonlinear behavior of the models are evaluated and their adjustment was evaluated based on quality evaluators, with the determination coefficient, residual standard deviation, and Akaike information criterion being the most cited.

Thus, this study aimed to describe the sugarcane stalk height growth curve, for the RB92579, RB93509, RB931530, and SP79-1011 varieties, in plant-cane and ratoon cane cycles, using the Logistic and Gompertz nonlinear models, considering any deviations from assumptions.

# MATERIALS AND METHODS:

The analyzed data was derived from an experiment conducted in the experimental Agrometeorology area of the Center for Agricultural Sciences (Centro de Ciências Agrárias), Federal University of Alagoas (Universidade Federal de Alagoas - UFAL), Delza Gitaí Campus, Rio Largo, in the Coastal Tablelands of Alagoas (ALMEIDA et al., 2008ALMEIDA, A. C. dos S. et al. Vegetative development and production of sugarcane varieties as a function of water availability and thermic units. Ciência e Agrotecnologia, v.32, n.5, p.1441-1448. 2008. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1413-70542008000500013&script=sci_arttext&tlng=es >. Acessed: May, 26, 2019. doi: 10.1590/S1413-70542008000500013.
http://www.scielo.br/scielo.php?pid=S141...
). The experiment was conducted using the randomized block design, with five replicates and four treatments, including the RB92579, RB93509, RB931530, and SP79-1011 varieties, in plant-cane and ratoon cane crop cycles. The soil of the experimental area was classified as cohesive argisolic Yellow Latosol of a medium/clayey texture, with mild topography and a slope lower than 2%. The varieties were planted in October 2003 within row spacings of 1.25 m, in 4 m plots with 5 furrows, corresponding to 25 m² per plot and a total area of 500 m². Two t ha-1 of dolomitic limestone (PRNT=85%) was added, with 70, 60, and 120 kg ha-1 of N, P2O5, and K2O base fertilization, respectively.

The plants were irrigated during the first plant growth stages (from sprouting to tillering), based on crop evapotranspiration (ETc) data. In this study, 17 plant-cane and ratoon cane stalk height measurements were taken 30 days after planting (DAP) and after cutting (DAC) for 360 days. The interval between observations was 15 days until the plants completed half of the growth cycle and 30 days in the final 180 DAP and DAC. All plants in the central row of each plot were evaluated, measuring the stalks from the soil to the insertion of the first fully developed leaf (leaf +1) with a visible auricle.

The Logistic and Gompertz nonlinear models were evaluated according to the parameterization suggested by FERNANDES et al. (2015FERNANDES, T. J. et al.Parameterization effects in nonlinear models to describe growth curves. Acta Scientiarum. Technology, v.37, n.4, p.397-402. 2015. Available from: <Available from: http://periodicos.uem.br/ojs/index.php/ActaSciTechnol/article/view/27855 >. Accessed: May, 13, 2019. doi: 10.4025/actascitechnol.v37i4.27855.
http://periodicos.uem.br/ojs/index.php/A...
). These models are represented by the following equations:

${Y}_{i}=A/\left(1+exp\left(k\left(B-{x}_{i}\right)\right)\right)+{\epsilon }_{i}$ and ${Y}_{i}=Aexp\left(-exp\left(k\left(B-{x}_{i}\right)\right)\right)+{\epsilon }_{i}$

respectively. For both models, Y i represents the observed height (m) values in observations i = 1,...,17; A represents the maximum horizontal asymptote, that is, the maximum stalk height to be reached; k represents the growth rate (the higher k is, the faster the plants will reach A); and B is interpreted as the abscissa of the inflection point of the Logistic model, from which growth slows down. The inflection point occurs at 50% of the horizontal asymptote (0.5A) in the Logistic model and 37% in the Gompertz model (0.37A). corresponds to the random error, which is assumed to be independently and identically distributed following a normal distribution with a mean of zero and constant variance, that is, $\epsilon ~N\left(0,{\sigma }^{2}\right)$ .

The graphical analysis of the residuals as a function of time (DAP and DAC) was performed, indicating the order (p) of the autocorrelation parameters (AR) through the partial autocorrelation function. When assessing the correlation between observations, the parameters were estimated using the generalized least squares method, incorporating AR(p) into the residuals. Thus, the residual value is provided using ${\epsilon }_{t}={\varnothing }_{1}{\epsilon }_{t-1}+\dots +{\varnothing }_{p}{\epsilon }_{t-p}+{e}_{t}$ . . Wherein, (t and (t-1 are the residuals at times t=1,2...p and t-p, respectively; is the autocorrelation parameter of order 1; while e t is the white noise, that is, ${e}_{t}~N\left(0;{\sigma }_{e}^{2}\right)$

being, $E\left({e}_{t}\right)=0$, $E\left({e}_{e}^{2}\right)={\sigma }_{e}^{2}$ , and , $E\left({e}_{t};{e}_{t-h}\right)=0$

for h≠0 (MORETTIN & TOLOI, 2006MORETTIN, P. A.; TOLOI, C. Análise de séries temporais. 2. ed. São Paulo: Blucher, 2006.). Subsequently, Shapiro-Wilk, Breusch-Pagan, and Durbin-Watson tests were used to validate the non-violation of the assumptions of normality, variance homogeneity, and independence of residuals.

The best model was selected using intrinsic nonlinear measures (IN) and the parametric effect (PE), assessed with the rms.curv( ) function of the R software (R CORE TEAM, 2018R CORE TEAM. R: A Language and Environment for Statistical Computing. 2018. Available from: <Available from: http://www.R-project.org/ >. Accessed: Oct. 20, 2018.
http://www.R-project.org/...
). Nonlinearity is used to indicate how distant a nonlinear model is from a linear model (FERNANDES et al., 2015FERNANDES, T. J. et al.Parameterization effects in nonlinear models to describe growth curves. Acta Scientiarum. Technology, v.37, n.4, p.397-402. 2015. Available from: <Available from: http://periodicos.uem.br/ojs/index.php/ActaSciTechnol/article/view/27855 >. Accessed: May, 13, 2019. doi: 10.4025/actascitechnol.v37i4.27855.
http://periodicos.uem.br/ojs/index.php/A...
). The quality of the model fit was assessed based on the adjusted coefficient of determination ${R}_{adj}^{2}=1-\left[\left(1-{R}^{2}\right)\left(n-i\right)\right]/\left(n-p\right)$ , residual standard deviation: $RSD=\sqrt{SSR/\left(n-p\right)}$ , and corrected Akaike Information Criterion: AICc = nIn (SSR/n) + [2 p (p + 1)] / (n - p - 1). Wherein, SSR and SST are the sum of the squared residuals and total square sum, respectively; n is the number of observations; p is the number of adjusted model parameters; i is related to the fit of the intercept to the curve, which is 1 if there is an intercept and 0 otherwise. R 2 = 1 - (SSR / SST) is the coefficient of determination, which indicates the variation of the data explained by the model, whereas ${R}_{adj}^{2}$ is only used to selected the models. The most suitable model is that with the highest ${R}_{adj}^{2}$ and lowest AICc values, evaluated using the Rsq.ad( ) and AICc( ) functions of the R software, respectively.

The model parameters were estimated based on the least squares method, using the Gauss-Newton algorithm, based on the Taylor series expansion, and implemented in the nls ( ) function of the R software. Stalk growth rates (SGR) were determined through the first derivative of the best growth over time (DAP or DAC) for each variety. The first derivatives of the Logistic and Gompertz models are provided by $y=Akexp\left(k\left(B-{x}_{i}\right)\right)/{\left(\left(1+exp\left(B-{x}_{i}\right)\right)\right)}^{2}$ and, $y=Akexp\left(-exp\left(k\left(B-x\right)\right)\right)exp\left(k\left(B-x\right)\right)$ respectively. The coordinates of the points of maximum acceleration (PMA), minimum deceleration (PMD), and asymptotic deceleration (PDA) of stalk growth were also calculated using the second, third, and fourth derivatives in relation to the independent variable (time), respectively. Based on the Logistic model, the PMA, PMD, and PDA coordinates are obtained by: ${x}_{PMA}=\left(-\mathrm{ln}\left(2+\sqrt{3}\right)-kB\right)/k$ , ${y}_{PMA}=A/\left(3+\sqrt{3}\right)$,

${x}_{PMD}=\left(-\mathrm{ln}\left(2+\sqrt{3}\right)-kB\right)/k$ , ${y}_{PMD}=A/\left(3+\sqrt{3}\right)$ ,

; and ${x}_{PDA}=\left(-\mathrm{ln}\left(5+2\sqrt{6}\right)-kB\right)/k$ ,

${y}_{PDA}=A/\left(2\left(3-\sqrt{6}\right)\right)$ respectively. Conversely, the Gompertz model is used to determine: ${x}_{PMA}=\left(-\mathrm{ln}\left(\left(3+\sqrt{5}\right)/2\right)-kB\right)/k$ , ${y}_{PMA}=Aexp\left(-\left(\left(3+\sqrt{5}\right)/2\right)\right)$ ; ${x}_{PMD}=\left(-\mathrm{ln}\left(\left(3+\sqrt{5}\right)/2\right)-kB\right)/k$ , ${y}_{PMD}=Aexp\left(-\left(3+\sqrt{5}\right)/2\right)$ and ${x}_{PDA}=\left(-\mathrm{ln}\left(0,165758\right)-kB\right)/k$ , ${y}_{PMA}=Aexp\left(-0,165758\right)$

, respectively (MISCHAN et al., 2011MISCHAN, M. M. et al. Determination of a point sufficiently close to the asymptote in nonlinear growth functions. Scientia Agricola, v.68, n.1, p.109-114. 2011. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S0103-90162011000100016&script=sci_arttext&tlng=pt >. Accessed: May, 10, 2019. doi: 10.1590/S0103-90162011000100016.
http://www.scielo.br/scielo.php?pid=S010...
; KLEINPAUL et al., 2018KLEINPAUL, J. A. et al. Nonlinear Models for Plant Height of Rye Cultivars at Sowing Times. Journal of Agricultural Science, v.10, n.12, p.157-170. 2018. Available from: <Available from: http://www.ccsenet.org/journal/index.php/jas/article/download/0/0/37368/37661 >. Accessed: Jun. 6, 2019. doi: 10.5539/jas.v10n12p157.
http://www.ccsenet.org/journal/index.php...
). The significance of the model parameters was assessed using the lack-of-fit (F) test at the 1% significance level.

# RESULTS AND DISCUSSION:

Figure 1 shows the graphs concerning the analysis of the residuals of the Logistic and Gompertz models adjusted to the average stalk height of plant-cane and ratoon cane. The trend in the Gompertz model, as shown in figure 1 (e, f, g, and h) in the plant-cane cycle, indicated residual autocorrelation. Several authors have also observed residual autocorrelation in growth studies, including PEREIRA et al. (2014PEREIRA, A. A. et al. Description vegetative growth of coffee tree farming ruby MG 1192 using regression models. Coffee Science, v.9, n.2, p.266-274. 2014. Available from: <Available from: http://www.sbicafe.ufv.br:80/handle/123456789/8051 >. Accessed: Apr. 23, 2019.
http://www.sbicafe.ufv.br:80/handle/1234...
), MUIANGA et al. (2016MUIANGA, C. A. et al.Description of the growth curve of cashew fruits in nonlinear models. Revista Brasileira de Fruticultura, v.38, n.1, p.22-32. 2016. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S0100-29452016000100022&script=sci_abstract&tlng=pt >. Accessed: May, 25, 2019. doi: 10.1590/0100-2945-295/14.
http://www.scielo.br/scielo.php?pid=S010...
), and MUNIZ et al. (2017MUNIZ, J. A. et. al. Nonlinear models for description of cacao fruit growth with assumption violations. Revista Caatinga, v.30, n.1, p.250-257. 2017. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1983-21252017000100250 >. Accessed: May, 13, 2019. doi: 10.1590/1983-21252017v30n128rc.
http://www.scielo.br/scielo.php?script=s...
) when evaluating the height of coffee plants, cashew nuts, and cocoa beans, respectively. The failure to fulfill the assumption of the independence of the residuals (autocorrelation) primarily occurs when the measurements are taken from the same individual over time.

Figure 1
Graphs of the residuals of the fit of the Logistic (a, b, c, and d) and Gompertz (e, f, g, and h) models and of the Logistic (i, j, k, and l) and Gompertz (m, n, o, and p) models adjusted to the average stalk height of four sugarcane varieties (RB92579, RB93509, RB931530, and SP79-1011), produced in the plant-cane and ratoon cane cycles, respectively.

The autocorrelation function of the residuals also indicated residual dependence for the Gompertz model when adjusted to plant-cane stalk height data (Figure 2; e, f, g, and h). As lags 1 and 2 were significant, second-order autocorrelation was modeled. However, the F test indicated the non-significance of the parameter, which showed a confidence interval with the value zero. Thus, to ensure a greater accuracy of estimates, the first-order autoregressive parameter AR(1) was added. Table 1 outlines the results from the analysis of the assumptions of the normality of the residuals, variance homogeneity, and independence of the residuals assessed using the Shapiro-Wilk (SW), Breusch-Pagan (BP), and Durbin-Watson (DW) tests, respectively, at the 1% significance level. All assumptions were fulfilled, indicating that the estimates obtained are reliable and the inferences are valid.

Figure 2
Graphs of the partial autocorrelation function of the residuals of the fit of the Logistic (a, b, c, and d) and Gompertz (e, f, g, and h) models and of the Logistic (i, j, k, and l) and Gompertz (m, n, o, and p) models adjusted to the average stalk height of four sugarcane varieties (RB92579, RB931530, and SP79-1011), produced in the plant-cane and ratoon cane cycles, respectively.

Table 1
P-values for the Shapiro-Wilk (SW), Breusch-Pagan (BP), and Durbin-Watson (DW) tests, intrinsic nonlinearity (IN), parametric effect (PE), adjusted coefficient of determination (${R}_{adj}^{2}$), residual standard deviation (RSD), and corrected Akaike information criterion (AICc) of the Logistic and Gompertz models adjusted to the average stalk height of four sugarcane varieties (RB92579, RB93509, RB931530, and SP79-1011), produced in the plant-cane and ratoon cane cycles.

Table 1 also outlines the values for the assessment of the intrinsic nonlinear measures (IN) and parametric effects (PE). Neither the intrinsic nonlinearity nor parametric effect were significant in either model (Logistic and Gompertz), in the four varieties and both cycles of cultivation. According to FERNANDES et al. (2015FERNANDES, T. J. et al.Parameterization effects in nonlinear models to describe growth curves. Acta Scientiarum. Technology, v.37, n.4, p.397-402. 2015. Available from: <Available from: http://periodicos.uem.br/ojs/index.php/ActaSciTechnol/article/view/27855 >. Accessed: May, 13, 2019. doi: 10.4025/actascitechnol.v37i4.27855.
http://periodicos.uem.br/ojs/index.php/A...
), the values of these measures are significant when higher than 0.5, thus indicating deviation from the linearity of the adjusted models. In this study, the Logistic model presented the lower values for the two measures in plant-cane, whereas the Gompertz model, concerning ratoon cane, presented lower values in three varieties (RB92579, RB931530, and SP79-1011). Thus, these models are close to the linear model, which ensures a greater reliability of their estimates.

The fit quality evaluators, outlined in Table 1, also indicated higher values and lower RSD and AICc values of the Logistic model in plant-cane. The sugarcane crop requires proper management during the first growth stage, especially in the plant-cane cycle. According to OLIVEIRA et al. (2010OLIVEIRA, E. C. A. de et al. Growth and dry matter production in sugarcane varieties grown under full irrigation.Revista Brasileira de Engenharia Agrícola e Ambiental, v.14, n.9, p.951-960. 2010. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-43662010000900007 >. Accessed: Apr. 24, 2019. doi: 10.1590/S1415-43662010000900007.
http://www.scielo.br/scielo.php?script=s...
), at first, plant-cane tillers grow slowly because they rely on the adaption of the seed piece reserves. Therefore, the Logistic model is perfectly suited to describe its growth, corroborating the results from this study. The Logistic model can capture longer sprouting and tillering stages because its inflection point is reported at the midpoint of the asymptote (0.5A). MUIANGA et al. (2016MUIANGA, C. A. et al.Description of the growth curve of cashew fruits in nonlinear models. Revista Brasileira de Fruticultura, v.38, n.1, p.22-32. 2016. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S0100-29452016000100022&script=sci_abstract&tlng=pt >. Accessed: May, 25, 2019. doi: 10.1590/0100-2945-295/14.
http://www.scielo.br/scielo.php?pid=S010...
), PEREIRA et al. (2016PEREIRA, A. A. et al. Modeling of the diameter of the canopy pruned coffee tree cultivated in different densities and water regimes.Lavras: Coffee Science, v.11, n.4, p.495-501. 2016. Available from: <Available from: http://www.coffeescience.ufla.br/index.php/Coffeescience/article/view/1145/pdf_1145 >. Accessed: May, 25, 2019. doi: 10.25186/cs.v11i4.1145.
http://www.coffeescience.ufla.br/index.p...
), MUNIZ et al. (2017MUNIZ, J. A. et. al. Nonlinear models for description of cacao fruit growth with assumption violations. Revista Caatinga, v.30, n.1, p.250-257. 2017. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1983-21252017000100250 >. Accessed: May, 13, 2019. doi: 10.1590/1983-21252017v30n128rc.
http://www.scielo.br/scielo.php?script=s...
), and RIBEIRO et al. (2018bRIBEIRO, T. D. et al. The use of the nonlinear models in the growth of pears of ‘Shinseiki’ cultivar. Santa Maria: Ciência Rural, v.48, n.1, e20161097. 2018b. Availale from: < Availale from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0103-84782018000100202 >. Accessed: May, 25, 2019. doi: 10.1590/0103-8478cr20161097.
http://www.scielo.br/scielo.php?script=s...
) also reported better fits in the Logistic model than in the Gompertz model when studying cashew nut and cacao bean growth, coffee canopy diameter, and pears, respectively.

For the ratoon cane cycle, the Gompertz model showed higher values and lower RSD and AICc values in three varieties (RB92579, RB931530, and SP79-1011), corroborating the results of the nonlinear measures. Ratoon cane usually grows fast by using the seed piece reserve after the stalks have been cut, thus shortening phenological stages. Therefore, the Gompertz model describes ratoon cane growth better than the Logistic model, because the inflection point of the former precedes that of the latter. Curiously, the RB93509 variety presented better nonlinear measure and fit evaluator results in the Logistic model. This indicated that the variety grew slowly, perhaps due to water stress and soil compaction as a result of mechanized harvesting. These factors strongly affect ratoon sprouting and tiller growth. In the literature, some authors have also reported good fits of the Gompertz model when evaluating ratoon cane growth, such as BATISTA et al. (2013BATISTA, E. L. da S. et al.Modeling of growth of sugarcane cultivars duringthe crop formation period. Revista Brasileira de Engenharia Agrícola e Ambiental, v.17, n.10, p.1080-1087. 2013. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1415-43662013001000009&script=sci_arttext >. Acessed: May, 26, 2019. doi: 10.1590/S1415-43662013001000009.
http://www.scielo.br/scielo.php?pid=S141...
), who studied the varieties RB92579, RB867515, RB928064, and RB855453 during the crop development period.

Thus, the Logistic model was selected as the best model to describe plant-cane growth. Whereas, the Gompertz model was selected for the ratoon cane varieties RB92579, RB931530, and SP79-1011, while the Logistic model was selected for the ratoon cane variety RB93509. The analysis of the parameter estimates of the selected models shows that in both plant-cane and ratoon cane, the RB92579 variety had the highest A and k values (Table 2), followed by SP79-1011. The results also showed that the higher the k value is, the lower the B estimate will be, indicating that the variety grows fast and it takes less time to reach the maximum growth rate. The inflection points of the growth curves occurred approximately from 129 to 146 DAP in plant-cane and from 79 to 90 DAC in ratoon cane. All estimated parameters were significant according to the t test at the 1% significance level, thus rejecting the hypothesis that model parameter estimates are zero.

Table 2
Estimates of the parameters (A, B, and k), standard error of the estimate in parentheses, and autoregressive parameters (ϕ 1 ) of the Logistic and Gompertz models adjusted to the average stalk height (m) of four sugarcane varieties (RB92579, RB93509, RB931530, and SP79-1011), produced in the plant-cane and ratoon cane cycles.

The results of the plant-cane estimates are similar to those observed by FREITAS et al. (2012FREITAS, C. A. S. de et al.Vegetative growth of sugar cane (Saccharum officinarum L.) irrigated with water from treated domestic sewage. Conexões-Ciência e Tecnologia, v.6, n.1, p.27-43. 2012. Available from: <Available from: http://conexoes.ifce.edu.br/index.php/conexoes/article/view/456 >. Accessed: May, 24, 2019. doi: 10.21439/conexoes.v6i1.456.
http://conexoes.ifce.edu.br/index.php/co...
) when analyzing irrigated sugarcane growth. These authors assessed average stalk height values ranging from 2.435 to 2.743 m. SANTOS et al. (2009SANTOS, V. R. dos et al. Growth and yield of sugarcane under different phosphorus sources. Revista Brasileira de Engenharia Agrícola e Ambiental, v.13, n.4, p.389-396. 2009. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-43662009000400004 >. Accessed: May, 27, 2019. doi: 10.1590/S1415-43662009000400004.
http://www.scielo.br/scielo.php?script=s...
) reported results close to the estimates of the RB93509 and SP79-1011 varieties when assessing the growth of the RB75126 variety, with the average stalk height ranging from 2.312 to 2.420 m. COSTA et al. (2011COSTA, C. T. S. et al. Growth and yield of four varieties of sugarcane in the third ratoon. Revista Caatinga, v.24, n. 3, p.56-63. 2011. Available from: <Available from: http://www.redalyc.org/articulo.oa?id=237119874008 >. Acessed: May, 9, 2019.
http://www.redalyc.org/articulo.oa?id=23...
) conducted the same experiment in the fourth crop cycle (3rd ratoon), and reported results similar to the findings of the present study (2nd ratoon) using the Gompertz model. The authors obtained estimates of 2.750 and 2.381 m for the asymptotic parameter of the RB92579 and SP79-1011 varieties, respectively. According to SILVA et al. (2012SILVA, T. G. F. da et al. Biometrics of the sugarcane shoot during irrigated ratoon cycle in the Submedio of the Vale do São Francisco. Revista Ciência Agronômica, v.43, n.3, p.500-509. 2012. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1806-66902012000300012&script=sci_arttext >. Accessed: May, 27, 2019. doi: 10.1590/S1806-66902012000300012.
http://www.scielo.br/scielo.php?pid=S180...
), the RB92579 variety showed higher shoot growth and is one of the varieties with the highest dry mass production, especially in ratoon cane, using soil nutrient reserves after the stalks of the previous cycle are cut.

The graphs concerning the fit of the best models are shown in figures 3A and 3C for plant-cane and ratoon cane, respectively. The growth of each variety is explained by figures 3B and 3D, which showed the stalk growth rates (SGR) of the four varieties. In these figures, the highest SGR indicated the inflection point (IP) of the model curve, which coincides with the estimate of parameter B, according to the variety and respective model. In plant-cane, the IP of the RB93509, RB92579, SP79-1011 and RB931530 varieties occurred when they showed an SGR of 2.096, 2.201, 2.618, and 1.678 cm day-1 (Figure 3B), respectively. According to the parameter B estimates (129 to 146 DAP), the IP occurred in February for the first three varieties, and in March for the last, thus coinciding with the rainy season in the Northeast region.

Figure 3
Graphs of the fit of the best models adjusted to the average sugarcane stalk height. Logistic model (1A) for plant-cane. Gompertz (varieties RB92579, RB931530. and SP79-1011) and Logistic (RB93509) models for ratoon cane (1C). and respective stalk growth rates (1B and 1D).

Conversely, in ratoon cane, the IP (from 79 to 83 DAC) occurred in December when the RB92579, RB931530, and SP79-1011 varieties reached an SGR of 2.009, 1.510, and 1.546 cm day-1 (Figure 3D), respectively. COSTA et al. (2011COSTA, C. T. S. et al. Growth and yield of four varieties of sugarcane in the third ratoon. Revista Caatinga, v.24, n. 3, p.56-63. 2011. Available from: <Available from: http://www.redalyc.org/articulo.oa?id=237119874008 >. Acessed: May, 9, 2019.
http://www.redalyc.org/articulo.oa?id=23...
) also observed low SGR values in relation to the IP: 1.210, 0.900, and 0.800 cm day-1 in the RB92579, RB931530, and SP79-1011 varieties, respectively. SILVA et al. (2012SILVA, T. G. F. da et al. Biometrics of the sugarcane shoot during irrigated ratoon cycle in the Submedio of the Vale do São Francisco. Revista Ciência Agronômica, v.43, n.3, p.500-509. 2012. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1806-66902012000300012&script=sci_arttext >. Accessed: May, 27, 2019. doi: 10.1590/S1806-66902012000300012.
http://www.scielo.br/scielo.php?pid=S180...
) observed low SGR in the IP of the height increase of the irrigated RB92579 variety, approximately 1.837 cm day-1, at 218 DAC in the first ratoon cane cycle in the Sub-middle São Francisco Valley, Bahia State. The low rates may have been caused by the fact that with each sugarcane cut the regrowth capacity of the sugarcane decreases, due to factors such as soil compaction, variety characteristics, and edaphoclimatic conditions.

According to critical point estimates, three stalk growth stages were identified. Plant-cane PMA estimates indicated that the first stage (from sprouting to tillering) ended from approximately 93 to 103 DAP, when the stalks showed 21% of the asymptote (A) in the four study varieties. The vegetative growth phase ended approximately between 165 and 188 DAP, as determined by the PMD, with 79% of the asymptotic parameter estimate. Conversely, growth stabilization (PDA) started from 190 to 220 DAP, with 91% of the A estimate, indicating the asymptotic deceleration point (PDA). When stalk growth reaches the PDA, crop management should be restricted to promote sucrose accumulation over vegetative growth, as reported by OLIVEIRA et al. (2010OLIVEIRA, E. C. A. de et al. Growth and dry matter production in sugarcane varieties grown under full irrigation.Revista Brasileira de Engenharia Agrícola e Ambiental, v.14, n.9, p.951-960. 2010. Available from: <Available from: http://www.scielo.br/scielo.php?script=sci_arttext&pid=S1415-43662010000900007 >. Accessed: Apr. 24, 2019. doi: 10.1590/S1415-43662010000900007.
http://www.scielo.br/scielo.php?script=s...
), SILVA et al. (2012SILVA, T. G. F. da et al. Biometrics of the sugarcane shoot during irrigated ratoon cycle in the Submedio of the Vale do São Francisco. Revista Ciência Agronômica, v.43, n.3, p.500-509. 2012. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1806-66902012000300012&script=sci_arttext >. Accessed: May, 27, 2019. doi: 10.1590/S1806-66902012000300012.
http://www.scielo.br/scielo.php?pid=S180...
), and BATISTA et al. (2013BATISTA, E. L. da S. et al.Modeling of growth of sugarcane cultivars duringthe crop formation period. Revista Brasileira de Engenharia Agrícola e Ambiental, v.17, n.10, p.1080-1087. 2013. Available from: <Available from: http://www.scielo.br/scielo.php?pid=S1415-43662013001000009&script=sci_arttext >. Acessed: May, 26, 2019. doi: 10.1590/S1415-43662013001000009.
http://www.scielo.br/scielo.php?pid=S141...
).

The critical points assessed using the Gompertz model in three varieties (RB92579, RB931530, and SP79-1011) indicated that the first stalk growth stage ended on the first month (from 29 to 30 DAC), whereas the second stage ended from 130 to 136 DAC, with a stalk height of 7 and 68%, respectively. Conversely, at 85% height, the growth of these varieties began to stabilize, from 175 to 183 DAC. These results indicated small differences between the PMA, PDM, and PDA estimates in the ratoon cane cycle, most likely because ratoon cane is more easily managed than plant-cane.

At the maturation stage, the vegetative growth of sugarcane stalks decreases. The results of the critical point estimates indicated that sugarcane maturation accounted for approximately 21 and 32% of plant-cane and ratoon cane stalk height growth, respectively. Ratoon cane grows very fast, closing the inter rows very early on. Thus, weeds are suffocated, thereby reducing moisture losses during the first two stalk growth stages.

# CONCLUSION:

The four study varieties showed growth with a sigmoid pattern in both crops. In general, the Logistic and Gompertz models without AR(1) better described the data of the mean sugarcane stalk height. The Logistic model was better suited to describe the plant-cane growth curve, whereas the Gompertz model better described the growth of the three varieties in the ratoon cane cycle.

The maximum stalk height growth rates occurred from 129 to 146 DAP and from 79 to 90 DAC, when the stalks reached from 1.678 to 2.618 cm day-1 and from 1.510 to 2.009 cm day-1 for plant-cane and ratoon cane, respectively.

# ACKNOWLEDGMENTS

To the Mozambique Scholarship Institute (Instituto de Bolsas de Estudo de Moçambique - IBE), Coordination for the Improvement of Higher Education Personnel (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES), and National Council for Scientific and Technological Development (Conselho Nacional de Desenvolvimento Científico e Tecnológico - CNPq) for the scholarships granted to the first, second, and third authors, respectively.

# REFERENCES

• CR-2019-0408.R1

# Publication Dates

• Publication in this collection
06 Apr 2020
• Date of issue
2020