# ABSTRACT:

The objectives of this study were to fit the Gompertz and Logistic models for the fresh and dry matter of leaves and the fresh and dry matter of shoots of three lettuce cultivars and indicate the best model to describe their growth in autumn-winter. The lettuce cultivars Gloriosa, Pira Verde, and Stella were evaluated in the autumn-winter of 2016 and 2017, in soilless in a protected environment. After transplantation, the fresh and dry matter of leaves and shoots were weighed every seven days. These dependent variables were fit using the accumulated thermal sum. The parameters of the Gompertz and Logistic models were estimated, the assumptions of the models were verified, the indicators of fit quality and critical points were calculated and the parametric and intrinsic curvature measures quantified. The Logistic and Gompertz models presented a satisfactory adjustment for the fresh and dry matter of leaves and the fresh and dry matter of shoots, for the lettuce cultivars Gloriosa, Pira Verde and Stella, in autumn-winter. The Logistic model best describes the growth of the lettuce cultivars.

Key words:
dry matter; fresh matter; Gompertz; Lactuca sativa L.; Logistic.

# RESUMO:

Palavras-chave:
Gompertz; Lactuca sativa L.; Logístico; massa de matéria fresca; massa de matéria seca.

# INTRODUCTION:

Lettuce (Lactuca sativa L.) is a temperate leafy green vegetable (SALA & COSTA, 2012SALA, C. F.; COSTA, C. P. Retrospective and trends of Brazilian lettuce crop. Horticultura Brasileira, v.30, p.187-194, 2012 Available from: <Available from: http://dx.doi.org/10.1590/S0102-05362012000200002 >. Accessed: Mar, 12, 2018. doi: 10.1590/S0102-05362012000200002.
http://dx.doi.org/10.1590/S0102-05362012...
). Its leaves are consumed as raw salads, soups, and creams, and is a source of dietary fibers, vitamins, and minerals (NTSOANE et al., 2016NTSOANE, L. L. M. et al. Variety-specific responses of lettuce grown under the different coloured shade nets on phytochemical quality after postharvest storage. The Journal of Horticultural Science and Biotechnology, v.91, p.520-528, 2016. Available from: <Available from: https://doi.org/10.1080/14620316.2016.1178080 >. Accessed: Sep. 10, 2018. doi: 10.1080/14620316.2016.1178080.
https://doi.org/10.1080/14620316.2016.11...
). It is the main leafy green vegetable sold and consumed in Brazil, mainly because of its ease of production and acquisition.

Lettuce cultivars are classified as crisphead, iceberg, or butterhead and other types (‘mimosa’, romaine, baby, and purple), corresponding to 43.3%, 41.2%, 5.0% and 10.5%, respectively, of the lettuce traded at the General Warehousing Company of São Paulo (CEAGESP, 2017CEAGESP. Companhia de Entrepostos e Armazéns Gerais de São Paulo. 2017. Available from: <Available from: http://www.ceagesp.gov.br/guia-ceagesp/alface-crespa/ >. Accessed: Jan. 09, 2018.
http://www.ceagesp.gov.br/guia-ceagesp/a...
). Thus, different genetic materials exist for the leaf morphological characteristics or head shape and also the growing seasons. In Rio Grande do Sul, the most favorable season to the growth of the crop is winter, because in this period temperatures vary between -3 to 18 ºC (ALVARES et al., 2013ALVARES, C. A. et al. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, v.22, p.711-728, 2013. Available from: <Available from: https://doi.org/10.1127/0941-2948/2013/0507 >. Accessed: May, 7, 2018. doi: 10.1127/0941-2948/2013/0507.
https://doi.org/10.1127/0941-2948/2013/0...
).

One way to characterize plant growth is via modeling (STRECK et al., 2008STRECK, N. A. et al. Modeling leaf appearance in cultivated rice and red rice. Pesquisa Agropecuária Brasileira, v.43, p.559-567, 2008. Available from: <Available from: http://dx.doi.org/10.1590/S0100-204X2008000500002 >. Accessed: Oct. 20, 2018. doi: S0100-204X2008000500002.
http://dx.doi.org/10.1590/S0100-204X2008...
). Adjusting growth models to plant species helps in the evaluation of plant response to environmental conditions, as well as understanding its growth pattern (LYRA et al., 2003LYRA, G. B. et al. Fitting growth models to hydroponic lettuce (Lactuca sativa L.) grown under greenhouse conditions. Revista Brasileira de Agrometeorologia, v.11, p.69-77, 2003. Available from: <Available from: https://www.researchgate.net/publication/230996209_Modelos_de_crescimento_para_alface_Lactuca_sativa_L_cultivada_em_sistema_hidroponico_sob_condicoes_de_casa_de_vegetacao >. Accessed: Jan. 18, 2018.
https://www.researchgate.net/publication...
). Growth models using the accumulated thermal sum allow to make inferences on precocity, velocity and stabilization of the plant growth through the interpretation of parameters and critical points of the adjusted model curve (MISCHAN & PINHO, 2014MISCHAN, M. M.; PINHO, S. Z. Modelos não lineares: funções assintóticas de crescimento. Cultura Acadêmica: São Paulo, 2014.).

The accumulated thermal sum is a biological time measure in plants, being more accurate than days in the civil calendar or days after sowing/transplant (GILMORE & ROGERS, 1958GILMORE, E. C.; ROGERS, J. S. Heat units as a method of measuring maturity in corn. Agronomy Journal, v.50, p.611-615, 1958. Available from: <Available from: https://doi.org/10.2134/agronj1958.00021962005000100014x >. Accessed: Jan. 22, 2018. doi: 10.1080/14620316.2018.1472045.
https://doi.org/10.2134/agronj1958.00021...
; MCMASTER & SMIKA, 1988MCMASTER, G. S.; SMIKA, D. E. Estimation and evaluation of winter wheat phenology in the Central Great Plains. Agricultural and Forest Meteorology, v.43, p.1-18, 1988. Available from: <Available from: https://doi.org/10.1016/0168-1923(88)90002-0 >. Accessed: Sep. 10, 2018. doi: 10.1016/0168-1923(88)90002-0.
https://doi.org/10.1016/0168-1923(88)900...
). The use of accumulated thermal sum as elapsed time of the crop cycle assumes a linear relation between growth or plant development and temperature (BONHOMME, 2000BONHOMME, R. Bases and limits to using ‘degree.day’ units. European Journal of Agronomy, v.13, p.1-10, 2000. Available from: <Available from: https://doi.org/10.1016/S1161-0301(00)00058-7 >. Accessed: Oct, 02, 2018. doi: 10.1016/S1161-0301(00)00058-7.
https://doi.org/10.1016/S1161-0301(00)00...
). However, this would not be realistic from the biological point of view, since the plant growth in response to the thermal accumulation is nonlinear. Therefore, nonlinear models are more often used to describe the growth of plants, generally, faster in its initial phase, then decreasing its speed and, finally, tending to a stability in the adult phase (PAINE et al., 2012PAINE, C. E. T. et al. How to fit nonlinear plant growth models and calculate growth rates: an update for ecologists. Methods in Ecology and Evolution, v.3, p.245-256, 2012. Available from: <Available from: https://doi.org/10.1111/j.2041-210X.2011.00155.x >. Accessed: Oct. 21, 2018. doi: 10.1111/j.2041-210X.2011.00155.x.
https://doi.org/10.1111/j.2041-210X.2011...
; MISCHAN & PINHO, 2014MISCHAN, M. M.; PINHO, S. Z. Modelos não lineares: funções assintóticas de crescimento. Cultura Acadêmica: São Paulo, 2014.).

Mathematical models must be able to reproduce the plants behavior as closest as possible to the real. The adjustment of nonlinear models have been applied to describe the growth of Allium sativum L. (PUIATTI et al., 2013PUIATTI, A. G. et al. Cluster analysis applied to nonlinear regression models selection for the description of dry matter accumulation of garlic plants. Revista Brasileira de Biometria, v.31, p.337-351, 2013. Available from: <Available from: http://jaguar.fcav.unesp.br/RME/fasciculos/v31/v31_n3/A2_Guilherme_PauloCecon.pdf >. Accessed: Oct. 21, 2018.
http://jaguar.fcav.unesp.br/RME/fascicul...
; REIS et al., 2014REIS, R. M. et al. Nonlinear regression models applied to clusters of garlic accessions. Horticultura Brasileira, 32(2), 178-183, 2014. Available from: <Available from: http://dx.doi.org/10.1590/S0102-05362014000200010 >. Accessed: May, 7, 2018. doi: 10.1590/S0102-05362014000200010.
http://dx.doi.org/10.1590/S0102-05362014...
) and production of the Cucurbita pepo and Capisicum annuum (LÚCIO et al., 2015LÚCIO, A. D. et al. Nonlinear models to describe production of fruit in Cucurbita pepo and Capiscum annuum. Scientia Horticulturae, v.193, p.286-293, 2015. Available from: <Available from: https://doi.org/10.1016/j.scienta.2015.07.021 >. Accessed: May, 19, 2018. doi: 10.1016/j.scienta.2015.07.021.
https://doi.org/10.1016/j.scienta.2015.0...
), cherry tomatoes (LÚCIO et al., 2016LÚCIO, A. D. et al. Nonlinear models for estimating cherry tomato yield. Ciência Rural, v.46, p.233-241, 2016. Available from: <Available from: http://dx.doi.org/10.1590/0103-8478cr20150067 >. Accessed: May, 19, 2018. doi: 10.1590/0103-8478cr20150067.
http://dx.doi.org/10.1590/0103-8478cr201...
) and strawberry (DIEL et al., 2018DIEL, M. I. et al. Nonlinear regression for description of strawberry (Fragaria x ananassa) production. The Journal of Horticultural Science and Biotechnology, p.1-15, 2018. Available from: <Available from: https://doi.org/10.1080/14620316.2018.1472045 >. Accessed: Dec. 18, 2018. doi: 10.1080/14620316.2018.1472045.
https://doi.org/10.1080/14620316.2018.14...
). According to TERRA et al. (2010TERRA, M. F. et al. (2010). Fitting Logistic and Gompertz models to the growth data of dwarf date palm (Phoenix roebelenni) fruits. Magistra, v.22, p.1-7. Available from: <Available from: https://www.ufrb.edu.br/magistra/2000-atual/volume-22-ano-2010/996-numero-1-jan-mar >. Accessed: May, 9, 2018.
https://www.ufrb.edu.br/magistra/2000-at...
), models allow condensing information from a series of data, taken over time, into a small set of biologically interpretable parameters.

It has been shown that for lettuce, the models Gompertz, Logistic, and Expolinear fit well the cultivars Grand Rapids, Regina, and Great Lakes, in a hydroponic system during the summer (LYRA et al., 2003LYRA, G. B. et al. Fitting growth models to hydroponic lettuce (Lactuca sativa L.) grown under greenhouse conditions. Revista Brasileira de Agrometeorologia, v.11, p.69-77, 2003. Available from: <Available from: https://www.researchgate.net/publication/230996209_Modelos_de_crescimento_para_alface_Lactuca_sativa_L_cultivada_em_sistema_hidroponico_sob_condicoes_de_casa_de_vegetacao >. Accessed: Jan. 18, 2018.
https://www.researchgate.net/publication...
). However, no report has been reported describing growth using nonlinear models in other seasons or with different cultivars in protected environment.

Suppose that the Gompertz and Logistic models are suitable to describe the growth of three lettuce cultivars during the autumn-winter season and that it is possible to select the most appropriate model. The objectives of this research were to adjust the Gompertz and Logistic models for the fresh and dry matter of leaves and shoots of three lettuce cultivars (Gloriosa, Pira Verde, and Stella) and indicated the model that best describes the growth in autumn-winter.

# MATERIALS AND METHODS:

Two experiments were carried out with lettuce cultivars, one in the autumn-winter of 2016 (experiment 1) and other in the autumn-winter of 2017 (experiment 2). The plants were grown using a closed soilles system, in a protected environment of umbrella type, with 115 m² (5×23 m) environment covered with 150 μm anti-UV polyethylene. The location is between coordinates 29º42’S, 53º49’W and 95 m altitude. According to Köppen’s classification, the climate of the region is humid subtropical Cfa, with hot summers and no defined dry season (ALVARES et al., 2013ALVARES, C. A. et al. Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift, v.22, p.711-728, 2013. Available from: <Available from: https://doi.org/10.1127/0941-2948/2013/0507 >. Accessed: May, 7, 2018. doi: 10.1127/0941-2948/2013/0507.
https://doi.org/10.1127/0941-2948/2013/0...
).

The lettuce cultivars evaluated were: Gloriosa (iceberg - light green leaves, crisp, consistent, prominent ribs, compact head), Pira Verde (crisp green - consistent and loose leaves that do not form head), and Stella (butterhead - delicate and smooth leaves with loosely formed head). These cultivars were recommended by the seed companies for autumn-winter conditions. The seedlings were produced in the floating system in 200-cell expanded polystyrene trays filled with commercial Plantmax® substrate. Transplanting was carried out when the plants developed four to five leaves, on 30/Jun/2016 (experiment 1) and 04/Jun/2017 (experiment 2).

Plants were grown in six benches of corrugated fiber cement sheets, 3.66 m long, 1.10 m wide, 6 mm thick, with six troughs of 5 cm in depth. The troughs were covered with clear 100-μm-thick plastic film and filled with washed gravel number two. The benches were raised (0.85 m) on fixed masonry blocks at the two end portions, with slope of 2%. This slope allowed the nutrient solution to return to the 500 L plastic storage tank. The solution was pumped by a low-power submersible motor pump (with a timer) to a PVC pipe (25 mm diameter). From this pipe derived four drip hoses with pots placed under the drippers at a distance of 30 cm between the plants in the row, to a density of 11.11 m-2 plants. Each bench consisted of four rows, totaling 44 three-liter volume pots (11 pots per row), filled with washed sieved coarse sand, with 0 dS m-1 electrical conductivity.

The nutrient solution consisted of the following macronutrient composition (in mmol L-1): 10.36 NO3 -; 1.0 H2PO4 -; 3.36 NH4; 1.0 SO4; 4.0 de K+; 2.0 Ca2+; 1.0 Mg2+; and micronutrients (mg L-1): 1.0 Fe; 0.50 Mn; 0.22 Zn; 0.26 B; 0.06 Cu, and 0.03 de Mo, for lettuce, with 1.33 dS m-1 electrical conductivity (EC) and pH 5.5 to 6.5. The EC and pH were monitored throughout the growing cycle and corrected if there was variation of 20%, higher or lower than the standard.

All cultivars were grown in the same environment, according to the information mentioned above. Seven days after transplantation, seven to twelve plants per cultivar were evaluated (nine times) in experiment 1 (total of 172 plants) and in experiment 2 six plants per cultivar were used and 10 support points (total of 180 plants) until the beginning of flowering. The variables fresh leaf matter (FLM, as g plant-1), dry leaf matter (DLM, as g plant-1), sum fresh stem matter and fresh leaf matter = fresh shoot matter (FSM, as g plant-1), and sum dry stem matter and dry leaf matter = dry shoot matter (DSM, as g plant-1) were weighed with a digital scale. For dry matter, the samples were packed in paper bags and incubated in a forced air circulation oven (60 ± 5 °C) until obtaining constant mass.

The indoor air temperature data were recorded every three hours by a digital data logger (0.1 °C resolution and 0.5 °C accuracy) installed in a weather-proof shelter located inside the umbrella greenhouse. These data were used to calculate the daily thermal sum by the method of GILMORE & ROGERS (1958GILMORE, E. C.; ROGERS, J. S. Heat units as a method of measuring maturity in corn. Agronomy Journal, v.50, p.611-615, 1958. Available from: <Available from: https://doi.org/10.2134/agronj1958.00021962005000100014x >. Accessed: Jan. 22, 2018. doi: 10.1080/14620316.2018.1472045.
https://doi.org/10.2134/agronj1958.00021...
) and ARNOLD (1959ARNOLD, C. T. The determination and significance of the base temperature in a linear heat unit system. Proceedings of the American Society for Horticultural Science, v.74, p.430-455, 1959.), using equations 1 and 2:

$STd=\left(Tmax+Tmin\right)/2-Tb$(1)

Where:

Tmax: maximum daily temperature as ºC;

Tmin: minimum daily temperature as ºC;

Tb: lettuce base temperature = 10ºC (BRUNINI, 1976BRUNINI, O. Base temperature for lettuce in a heat-unit system. Bragantia, v.35, p.213-219, 1976. Available from: <Available from: https://doi.org/10.1590/S0006-87051976000100019 >. Accessed: Oct. 15, 2018. doi: 10.1590/S0006-87051976000100019.
https://doi.org/10.1590/S0006-8705197600...
)

$\mathrm{S}\mathrm{T}\mathrm{a}\mathrm{}=\sum \mathrm{}\mathrm{S}\mathrm{T}\mathrm{d}$(2)

Where:

STa: accumulated thermal sum;

∑ STd: sum of the daily thermal sum.

The fit of the Gompertz (WINDSOR, 1932WINDSOR, C. P. The Gompertz Curve as a Growth Curve. Proceedings of the National Academy of Sciences of the United States of America, 18, 1-8, 1932. Available from: <Available from: https://doi.org/10.1073/pnas.18.1.1 >. Accessed: May, 9, 2018. doi: 10.1073/pnas.18.1.1.
https://doi.org/10.1073/pnas.18.1.1...
) and Logistic (NELDER, 1961NELDER, J. A. The fitting of a generalization of the Logistic curve. Biometrics, v.17, p.89-110, 1961. Available from: <Available from: https://doi.org/10.2307/2527498 >. Accessed: May, 16, 2018. doi: 10.2307/2527498.
https://doi.org/10.2307/2527498...
) models for each character (dependent variable) was performed using the repetitions of each evaluation (for each cultivar x experiment separately), using the raw data, as a function of the accumulated thermal sum (independent variable). The equation used for the Gompertz model was: ${y}_{i}=aexp\left[-exp\left(b-c{x}_{i}\right)+{\epsilon }_{i}\right]$ , and for the Logistic: ${y}_{i}=a/\left[1+exp\left(-b-c{x}_{i}\right)+{\epsilon }_{i}\right]$y i is the i-th observation of the dependent variable with i = 1, 2, ..., n; x i is the i-th observation of the independent variable; a is the asymptotic value; b is a location parameter, important for maintaining the sigmoidal shape of the model; c is associated with growth, indicating the precocity index. The higher the value of c the less time will be required for the plant to reach the asymptotic value (a).

The assumptions of normality, independence, and homogeneity of the model residuals were tested using the Shapiro-Wilk (SHAPIRO & WILK, 1965SHAPIRO, S. S.; WILK, M. B. An analysis of variance test for normality. Biometrika, v.52, p.591-611, 1965. Available from: <Available from: http://dx.doi.org/10.2307/2333709 >. Accessed: Feb. 5, 2018. doi: 10.2307/2333709.
http://dx.doi.org/10.2307/2333709...
), Durbin-Watson (DURBIN & WATSON, 1950), and Breusch-Pagan tests (BREUSCH & PAGAN, 1979BREUSCH, T.; PAGAN, A. A Simple test for heteroscedasticity and random coefficient variation. Sociedade Econométrica, v.47, p.1287-1294, 1979. Available from: <Available from: http://dx.doi.org/10.2307/1911963 >. Accessed: Jul. 15, 2018. doi: 10.2307/1911963.
http://dx.doi.org/10.2307/1911963...
) respectively.

The estimates of the parameters (a, b e c) for each response were compared between the experiments for each cultivar, and between the cultivars in each experiment, by overlapping confidence intervals (CI) of the parameter estimates in each model. For this purpose, we calculated the lower and upper limits of the 95% confidence interval.

The coefficient of determination ${\mathrm{R}}^{2}=1-\frac{\mathrm{S}\mathrm{Q}\mathrm{R}}{\mathrm{S}\mathrm{Q}\mathrm{T}}$ was used to assess the quality of fit of the models, and the best fit was considered the model with the coefficient closest to 1 or 100%. The Akaike Information Criterion $\mathrm{A}\mathrm{I}\mathrm{C}=\mathrm{ln}\left(\mathrm{\sigma }²\right)+2\left(\mathrm{p}+1\right)/\mathrm{n}$ in which the lower its value the better the model (that is, the more suitable the model is to describe the study), and the Residual Standard Deviation $\mathrm{R}\mathrm{S}\mathrm{D}=\sqrt{\frac{\mathrm{S}\mathrm{Q}\mathrm{R}}{\mathrm{n}-\mathrm{p}}}$ , define the best fit of the model with values closer to zero. The intrinsic curvature measures (ICM) and curvature measures of the parameter effect (PE) were quantified using the geometric concept of curvature (BATES & WATTS, 1998BATES, D. M.; WATTS, D. G. Nonlinear regression analysis and its applications. New York: John Wiley & Sons, 1998.). The selection of the best model to describe plant growth, is based on the one that provides the lowest values of intrinsic or parametric curvature measures. Were calculated according to the equations described in MISCHAN & PINHO (2014MISCHAN, M. M.; PINHO, S. Z. Modelos não lineares: funções assintóticas de crescimento. Cultura Acadêmica: São Paulo, 2014.), the inflection point (IP) to Gompertz ($IPx=\frac{b}{c}$and $IPy=\frac{a}{e}\right)$) and to Logistc $\left(IPx=\frac{-b}{c}$ and $IPy=\frac{a}{2}$), the maximum acceleration point (MAP) to Gompertz $\left(MAPx=\frac{b-ln\left(\frac{3+\sqrt{5}}{2}\right)}{c}$ and $MAPy=aexp\left(-\frac{3+\sqrt{5}}{2}\right)$ ) to Logistc $MAPx=\frac{1}{c}\left[-b-ln\left(2+\sqrt{3}\right)\right]$ and $MAPy=\frac{a}{3+\sqrt{3}}$ the maximum deceleration point (MDP) to Gompertz $PDMx=\frac{b-ln\left(\frac{3-\sqrt{5}}{2}\right)}{c}$ and $PDMy=aexp\left(-\frac{3-\sqrt{5}}{2}\right)\right)$

and to Logistc $PDMx=\frac{1}{c}\left[-b-ln\left(2-\sqrt{3}\right)\right]$ and $PDMy=\frac{a}{3-\sqrt{3}}$ .

Inferences about plant growth were made from these critical points. The calculations were performed using the Microsoft Office Excel® applications and the R software, with the nls function (R DEVELOPMENT CORE TEAM, 2018R development Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria, 2018. Available from <Available from http://www.R-project.org/ >. Accessed: Oct. 21, 2018.
http://www.R-project.org/...
).

# RESULTS AND DISCUSSION:

The assumptions of normality, homogeneity, and independence of errors were met in both the Gompertz and Logistic models for fresh and dry matter of leaves and fresh and dry matter of shoots of lettuce cultivars in both experiments, as the Shapiro-Wilk, Durbin-Watson, and Breusch-Pagan tests had p-values equal to or greater than 0.05. Similar results were reported by RIBEIRO et al. (2018), in which the assumptions were taken to nonlinear models.

The estimates of a are the asymptotic values, that is, in the case of lettuce, represent the maximum mass achieved. For all the characters of the cultivars, the a values for the Gompertz model were higher than for the Logistic model (Tables 1 and 2). The estimation of b, in theory, provides a concept of the ratio between the initial values and the missing amount to reach the asymptote. The estimate of parameter c, represents the growth speed, which was higher in the Logistic model (Tables 1 and 2).

Table 1
Estimates of the parameters a, b, and c, lower limit (LL) and upper limit (UL) of the confidence interval (CI 95%) of the Gompertz model for the characters as a function of accumulated thermal sum (in °C) of lettuce cultivars (Gloriosa, Pira Verde, and Stella) in two experiments.
Table 2
Estimation of the parameters a, b, and c, lower limit (LL) and upper limit (UL) of the confidence interval (CI 95%) of the Logistic model for the characters as a function of accumulated thermal sum (in °C) of lettuce cultivars (Gloriosa, Pira Verde, and Stella) in two experiments.

The estimates of the parameters (a, b and c) of each character for the Gompertz and Logistic models were compared between each experiment (Tables 1 and 2) and between the each cultivars (Tables 3), by the criterion of overlapping confidence intervals. This comparison criterion was used by WHEELERN et al. (2006WHEELER, M. W. et al. Comparing median lethal concentration values using confidence interval overlap or ratio tests. Environmental Toxicology and Chemistry, v.25, p.1441-1444, 2006. Available from: <Available from: http://dx.doi.org/10.1897/05-320R.1 >. Accessed: Oct. 20, 2018. doi: 10.1897/05-320R.1.
http://dx.doi.org/10.1897/05-320R.1...
) and by BEM et al. (2018BEM, C. M. et al. Gompertz and Logistic models to the productive traits of sunn hemp. Journal of Agricultural Science, v.10, p.225-238, 2018. Available from: <Available from: https://doi.org/10.5539/jas.v10n1p225 >. Accessed: May, 18, 2018. doi: 10.5539/jas.v10n1p225.
https://doi.org/10.5539/jas.v10n1p225...
), to verify if the growth curves have differed according to the treatments.

Table 3
Comparison of estimates of parameters (a, b and c) in the Gompertz and Logistic models for characters as a function of cumulative thermal sum based on the confidence interval (CI 95%), between lettuce cultivars Gloriosa, Pira Verde, and Stella in two experiments.

To clarify the comparison by the criterion of overlapping 95% confidence intervals (CI), the FLM of cv. Pira Verde will be used as an example to compare the estimate of the parameter a of the Logistic model between experiments 1 and 2 (Table 2). The following results were reported: the estimate of parameter a (354.7561) in experiment 1 is within the confidence interval of the estimate of parameter a in experiment 2 (329.5156 to 385.3911). Also, the estimate of parameter a (357.4533) in experiment 2 is within the confidence interval of the estimate of parameter a of experiment 1 (225.9895 to 483.5227). Therefore, the estimates of the parameter a are not different between the experiments. When at least one of the estimates is within the CI of the other, it can be concluded that the effect is not significant. If the two parameter estimates are outside the CI of the other, it can be concluded that the effect is significant.

In the Gompertz model, the parameters b and c for FLM and FSM of cv. Gloriosa were not different between the experiments (Table 1). However, the parameter a differed in all the characters, with higher values of FLM and FSM in experiment 1, which indicates higher matter production in relation to experiment 2. Opposite behavior was observed for DLM and DSM. For cv. Pira Verde, the estimates were not different for the characters except for DLM in relation to parameter c. These results indicated that, for this cultivar, there was no difference of the Gompertz models between the experiments. However, for cv. Stella, no differences were observed for FLM and FSM between the experiments. The DLM and DSM were not different asymptotically.

The Logistic model of cv. Gloriosa showed that DLM and DSM differed between experiments for parameters a, b and c (Table 2). FLM and FSM showed differences only in the asymptotic values and were higher in the experiment 1 than in the experiment 2, which indicated that the plants had higher production of green matter in experiment 1. Characters differed of cultivar Pira Verde with respect to the parameters b and c, and did not differ for parameter a. Asymptotic values of cultivar Stella were not different between experiments for the all the characters, the estimate of b was similar for FLM and FSM. However, all characters differed for growth rate.

These results suggested that the growth models had different behaviors between experiments 1 and 2. Similar results were reported for genotype tomato in two years, in the Cordillera and Ellen genotypes were more premature in 2015/2016 crop, and the Gaucho genotype was more premature in the 2016/2017 crop (SARI et al. 2019SARI, B. G. et al. Describing tomato plant production using growth models. Scientia Horticulturae, v.246, p.146-154, 2019. Available from: <Available from: https://doi.org/10.1016/j.scienta.2018.10.044 >. Accessed: Jun. 3, 2019. doi: 10.1016/j.scienta.2018.10.044.
https://doi.org/10.1016/j.scienta.2018.1...
).

Comparing the cultivars in each experiment, we found that cvs. Gloriosa and Pira Verde, in the Gompertz model, experiment 1, showed no difference between the characters (Table 3). This means that the Gompertz model makes no difference between these cultivars. Conversely, Gompertz model differed for all the characters of cvs. Pira Verde and Stella, since at least one of the three parameters (a, b and c) was significant. This same behavior was observed between cvs. Gloriosa and Stella. In experiment 2, the cultivars Gloriosa and Pira Verde were not different for FLM, while Pira Verde and Stella did not differ for DLM and DSM.

The estimates of the Logistic model parameters, for FSM of cvs. Pira Verde and Stella in experiment 1, for FSM of cvs. Gloriosa and Pira Verde in experiment 2, and for DLM and DSM of cvs. Pira Verde and Stella in experiment 2, were not different (Table 3). All other comparisons differed in at least one of the three parameters of the Logistic model. For the Gompertz and Logistic models there was a predominance of differences, which indicated the need of specific models per character and cultivar. Different models were also required in groups of garlic accesses (REIS et al., 2014REIS, R. M. et al. Nonlinear regression models applied to clusters of garlic accessions. Horticultura Brasileira, 32(2), 178-183, 2014. Available from: <Available from: http://dx.doi.org/10.1590/S0102-05362014000200010 >. Accessed: May, 7, 2018. doi: 10.1590/S0102-05362014000200010.
http://dx.doi.org/10.1590/S0102-05362014...
).

Goodness-of-fit indicators are used to define the most suitable model. The Logistic and Gompertz models presented acceptable goodness-of-fit values (high R², low AIC and RSD) and close to each other (Tables 4 and 5). The R² indicator was used by LIRA et al. (2003) to study the growth curve of lettuce cultivars. However, it is recommended to use more than one fit quality indicator to increase the reliability of the model choice.

Table 4
Coefficient of determination (R2), Akaike information criterion (AIC), residual standard deviation (RSD), intrinsic curvature measures (ICM), curvature measures of the parameter effect (PE), inflection point (IP), maximum acceleration point (MAP), and maximum deceleration point (MDP) of the Gompertz model for characters (1) as a function of the accumulated thermal sum (in °C) of lettuce cultivars (Gloriosa, Pira Verde, and Stella) in two experiments.
Table 5
Coefficient of determination (R2), Akaike information criterion (AIC), residual standard deviation (RSD), intrinsic curvature measures (ICM), curvature measures of the parameter effect (PE), inflection point (IP), maximum acceleration point (MAP), and maximum deceleration point (MDP) of the Logistc model for characters (1) as a function of the accumulated thermal sum (in °C) of lettuce cultivars (Gloriosa, Pira Verde, and Stella) in two experiments.

The Gompertz and Logistic models satisfactorily described the growth curve of lettuce cultivars, with R² values equal to or higher than 0.913, except for cv. Stella, which showed lower Goodness-of-fit in experiment 2 (0.769 ≤ R² ≤ 0.826), of both models. LYRA et al. (2003LYRA, G. B. et al. Fitting growth models to hydroponic lettuce (Lactuca sativa L.) grown under greenhouse conditions. Revista Brasileira de Agrometeorologia, v.11, p.69-77, 2003. Available from: <Available from: https://www.researchgate.net/publication/230996209_Modelos_de_crescimento_para_alface_Lactuca_sativa_L_cultivada_em_sistema_hidroponico_sob_condicoes_de_casa_de_vegetacao >. Accessed: Jan. 18, 2018.
https://www.researchgate.net/publication...
) adjusted growth models for dry leaf matter in lettuce cultivars in the summer in Viçosa, Minas Gerais, Brazil, and obtained results partially similar to the present study, with a coefficient of determination equal to or greater than 0.98.

Although, the models presented satisfactory Goodness-of-fit for FLM and FSM of cv. Gloriosa in experiment 1, the Gompertz model overestimated the parameter a with 929.8712 for FLM and 984.0988 for FSM (Table 1), that is, these estimates were higher than the maximum values observed in the data set, which were 647.4 g plant-1 of FLM and 665.77 g plant-1 of FSM. A greater overestimation was found for FLM and FSM of cv. Pira Verde in experiment 1, where the Gompertz model estimated the parameter a as 1234.9167 for FLM and as 1334.4813 for FSM, while the maximum value observed was 276.16 g plant-1 for FLM and 288.02 g plant-1 for FSM. Cases of overestimation of parameters in the Gompertz model were also described for the dry matter of bulbs of garlic accesses (REIS et al., 2014REIS, R. M. et al. Nonlinear regression models applied to clusters of garlic accessions. Horticultura Brasileira, 32(2), 178-183, 2014. Available from: <Available from: http://dx.doi.org/10.1590/S0102-05362014000200010 >. Accessed: May, 7, 2018. doi: 10.1590/S0102-05362014000200010.
http://dx.doi.org/10.1590/S0102-05362014...
).

Intrinsic curvature measures (ICM) and curvature measures of the parameter effect (PE) help us choose the best model. We found that the Logistic model had lower ICM for most of the characters of the cultivars in the two experiments and the PE was always smaller than in the Gompertz model (Tables 4 and 5). The lower ICM and, especially, the lower PE indicate better suitability of the Logistic model, when compared to the Gompertz model.

Considering the five Goodness-of-fit indicators (R2, AIC, RSD, ICM, and PE), we can infer that the Logistic model had suitable behavior regardless of cultivar, character, and experiment and is the best indicated to describe the growth of the lettuce cultivars. To exemplify the growth curve shape of the Logistic model for each character, with the respective critical points, we selected cv. Gloriosa of experiment 2 (Figure 1). The other growth curves can be constructed with the respective estimates of the parameters (Table 2).

Figure 1
Logistic model plot for fresh leaf matter (FLM, as g plant-1), dry leaf matter (DLM as g plant-1), fresh shoot matter (FSM, as g plant-1), and dry shoot matter (DSM, as g plant-1) as a function of the accumulated thermal sum (STa, in °C), for the cultivar Gloriosa, in experiment 2.

Inflection points, maximum acceleration and maximum deceleration are used to infer the crop growth, having as a base the general behavior to the cultivars Gloriosa, Pira Verde and Stella (Tables 4 and 5). The maximum acceleration point occurred at the beginning of the curve, when the plants showed slow growth, which is related to smaller plants and still young leaves. For most cultivars, in both experiments, the inflection point (IP) coincided with the phase close to harvest point, with the appearance of senescent basal leaves, which in practice is one of the criteria to classify the produce commercially. In general, independent of the experiment, in the Gompertz model, the cultivars reached the IP with lower STa than in the Logistic model. Among the cultivars, Gloriosa required the greatest accumulation of thermal sum and had the highest dry and fresh matter compared with cvs. Stella and Pira Verde. The results showed that iceberg cultivars need higher thermal sum during the autumn-winter period, due to the process of formation of a compact commercial head (YURI et al., 2017YURI, J. E. et al. Agronomic performance of crisphead lettuce genotypes at Sub-Middle São Francisco Valley. Horticultura Brasileira, p.35, p.292-297, 2017. Available from: <Available from: http://dx.doi.org/10.1590/s0102-053620170222 >. Accessed: Oct. 14, 2018. doi: 10.1590/s0102-053620170222.
http://dx.doi.org/10.1590/s0102-05362017...
).

The characters fresh leaf matter and fresh shoot matter represent the edible part of the lettuce, that is, the part of major commercial interest. Thus, the fresh mass has greater relevance than the dry matter. Accordingly, the inflection points of the fresh leaf mass and aerial part can be used in a practical way, since the accumulated thermal sum (IPx) reflects the amount of mass accumulated (IPy) near the harvest point. Therefore, this information is useful for producers and researchers that work with this crop.

The results of this study showed that the Logistic nonlinear growth model and its critical points are relevant to help the selection of promising lettuce cultivars. The Logistic model was also used to describe the growth curve of dry matter of the aerial part, the bulb and the whole plant of the onion culture (PÔRTO et al., 2006PÔRTO, D. R. Q. et al. Macronutrients accumulation by onion ‘Optima’ established by direct sowing. Horticultura Brasileira, 24, 470-475, 2006. Available from: <Available from: http://dx.doi.org/10.1590/S0102-05362006000400015 >. Accessed: Oct. 21, 2018. doi: 10.1590/S0102-05362006000400015.
http://dx.doi.org/10.1590/S0102-05362006...
), the production of genotype tomato (SARI et al., 2019SARI, B. G. et al. Describing tomato plant production using growth models. Scientia Horticulturae, v.246, p.146-154, 2019. Available from: <Available from: https://doi.org/10.1016/j.scienta.2018.10.044 >. Accessed: Jun. 3, 2019. doi: 10.1016/j.scienta.2018.10.044.
https://doi.org/10.1016/j.scienta.2018.1...
) and to describe the production of strawberry cultivars from different seedling origins grown on organic substrates (DIEL et al., 2018DIEL, M. I. et al. Nonlinear regression for description of strawberry (Fragaria x ananassa) production. The Journal of Horticultural Science and Biotechnology, p.1-15, 2018. Available from: <Available from: https://doi.org/10.1080/14620316.2018.1472045 >. Accessed: Dec. 18, 2018. doi: 10.1080/14620316.2018.1472045.
https://doi.org/10.1080/14620316.2018.14...
).

The parameters of the Logistic model were estimated as a function of the relations between the productive characters and the accumulated thermal sum. The parameters (a, b and c) estimated in this study can be used for simulation and prediction of growth of cvs. Gloriosa, Pira Verde, and Stella in the autumn-winter period, for research or production. However, we recommend the use of the thermal sum from the crop site, because the lettuce crop is more influenced by temperature during its vegetative phase (LOPES et al., 2004LOPES, S. J. et al. Models to estimate phytomass accumulation of hydroponic lettuce. Scientia Agricola, v.61, p.392-400, 2004. Available from: <Available from: http://dx.doi.org/10.1590/S0103-90162004000400007 >. Accessed: May, 7, 2018. doi: 10.1590/S0103-90162004000400007.
http://dx.doi.org/10.1590/S0103-90162004...
). Thus, this prediction can be used, but the values obtained will be approximated to those reported in this study and should follow its growth curve pattern. In addition, because we found no studies focusing on growing these cultivars in autumn-winter, the models developed here can become a reference for further research.

# CONCLUSION:

The growth models developed show differences between the experiments (years) and among the cultivars. The Logistic and Gompertz growth models showed a satisfactory Goodness-of-fit for the fresh and dry matter of leaves and fresh and dry matter of shoots of the lettuce cultivars Gloriosa, Pira Verde, and Stella, in autumn and winter. The Logistic model best describes the growth of the lettuce cultivars.

# ACKNOWLEDGEMENTS

To the National Council of Scientifc and Technological Development (CNPQ - Processes number 401045/2016-1 and 304652/2017-2), the Coordination of Improvement of Personnel of Superior Level (CAPES) and the Foundation of Support to Research of the state of Rio Grande do Sul (FAPERGS), by the granting of scholarship to the authors.

# REFERENCES

• CR-2019-0534.R1

# Publication Dates

• Publication in this collection
15 June 2020
• Date of issue
2020