<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="brief-report" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2024.1426424</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Brief Research Report</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Comparison of five equations in describing the variation of leaf area distributions of <italic>Alangium chinense</italic> (Lour.) Harms</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Deng</surname>
<given-names>Linli</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>He</surname>
<given-names>Ke</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="author-notes" rid="fn003">
<sup>&#x2020;</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Niklas</surname>
<given-names>Karl J.</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/33184"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Shi</surname>
<given-names>Zhuyue</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mu</surname>
<given-names>Youying</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Shi</surname>
<given-names>Peijian</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/255240"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Department of Applied Mathematics, College of Science, Nanjing Forestry University</institution>, <addr-line>Nanjing</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Architectural Design and Research Institute, Shenzhen University</institution>, <addr-line>Shenzhen</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>School of Civil Engineering and Architecture, Xiamen University of Technology</institution>, <addr-line>Xiamen</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>School of Integrative Plant Science, Cornell University</institution>, <addr-line>Ithaca, NY</addr-line>, <country>United States</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Jie Gao, Xinjiang Normal University, China</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Hede Gong, Southwest Forestry University, China</p>
<p>Wang Haitao, University of G&#xf6;ttingen, Germany</p>
<p>Tian Dong, Central South University, China</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Peijian Shi, <email xlink:href="mailto:pjshi@njfu.edu.cn">pjshi@njfu.edu.cn</email>
</p>
</fn>
<fn fn-type="equal" id="fn003">
<p>&#x2020;These authors have contributed equally to this work</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>04</day>
<month>07</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>15</volume>
<elocation-id>1426424</elocation-id>
<history>
<date date-type="received">
<day>01</day>
<month>05</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>13</day>
<month>06</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Deng, He, Niklas, Shi, Mu and Shi</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Deng, He, Niklas, Shi, Mu and Shi</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Previous studies have validated a performance equation (PE) and its generalized version (GPE) in describing the rotated and right-shifted Lorenz curves of organ size (e.g., leaf area and fruit volume) distributions of herbaceous plants. Nevertheless, there are still two questions that have not been adequately addressed by prior work: (i) whether the PE and GPE apply to woody plant species and (ii) how do the PE and GPE perform in comparison with other Lorenz equations when fitting data. To address these deficiencies, we measured the lamina length and width of each leaf on 60 <italic>Alangium chinense</italic> saplings to compare the performance of the PE and GPE with three other Lorenz equations in quantifying the inequality of leaf area distributions across individual trees. Leaf area is shown to be the product of a proportionality coefficient (<italic>k</italic>) and leaf length and width. To determine the numerical value of <italic>k</italic>, we scanned 540 leaves to obtain the leaf area empirically. Using the estimated <italic>k</italic>, the leaf areas of 60&#xa0;A<italic>. chinense</italic> saplings were calculated. Using these data, the two performance equations and three other Lorenz equations were then compared and assessed using the root-mean-square error (RMSE) and Akaike information criterion (AIC). The PE and GPE were found to be valid in describing the rotated and right-shifted Lorenz curves of the <italic>A. chinense</italic> leaf area distributions, and GPE has the lowest RMSE and AIC values. This work validates the GPE as the best model in gauging variations in leaf area of the woody species.</p>
</abstract>
<kwd-group>
<kwd>Akaike information criterion</kwd>
<kwd>leaf area distribution</kwd>
<kwd>Lorenz curve</kwd>
<kwd>model comparison</kwd>
<kwd>nonlinear regression</kwd>
</kwd-group>
<counts>
<fig-count count="4"/>
<table-count count="0"/>
<equation-count count="8"/>
<ref-count count="31"/>
<page-count count="8"/>
<word-count count="3771"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Functional Plant Ecology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Prior work has shown that even small differences in the morphological traits of leaves can affect the photosynthetic and respiration rates (<xref ref-type="bibr" rid="B31">Zirbel et&#xa0;al., 2017</xref>) and that leaf (lamina) area is a critical functional trait for a large number of plant species (<xref ref-type="bibr" rid="B29">Westoby, 1998</xref>), providing the maximum information on plant growth and resource utilization (<xref ref-type="bibr" rid="B3">Hodgson et&#xa0;al., 1999</xref>). Therefore, accurately measuring the leaf area and the extent to which it varies across and within species is important to understand plant ecology and evolution.</p>
<p>The Montgomery equation (ME) assumes a proportional relationship between leaf area and the product of leaf length and leaf width (<xref ref-type="bibr" rid="B12">Montgomery, 1911</xref>). Previous work has shown that the ME can effectively predict the area of even intricately shaped lamina (<xref ref-type="bibr" rid="B23">Shi et&#xa0;al., 2019a</xref>, <xref ref-type="bibr" rid="B24">b</xref>; <xref ref-type="bibr" rid="B30">Yu et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B20">Schrader et&#xa0;al., 2021</xref>). In addition, <xref ref-type="bibr" rid="B30">Yu et&#xa0;al. (2020)</xref> found that the ME is superior to other mathematical models using leaf length and width based on data drawn from 15 species of vines which included some complex leaf shapes. However, prior studies delving into the utility of the ME have used computer recognition methods to identify the leaf maximum width perpendicular to the leaf length axis. This approach is difficult under field research conditions to accurately and non-destructively identify leaf maximum width. Fortunately, a recent study has shown a negligible difference between the use of computer recognition methods and manual measurements defining leaf width by simply connecting two leaf boundary landmarks (<xref ref-type="bibr" rid="B13">Mu et&#xa0;al., 2024</xref>).</p>
<p>To assess variations in leaf area, prior work has turned to using the mathematics of the Lorenz curve, which is commonly used in economics to describe the inequality of income distributions (<xref ref-type="bibr" rid="B9">Lorenz, 1905</xref>; <xref ref-type="bibr" rid="B7">Kakwani, 1977</xref>) by plotting the cumulative proportion of household income against the cumulative proportion of the number of households. When the income distributions across households are absolutely equal, the Lorenz curve overlaps with the&#xa0;straight line passing through (0, 0) and (1, 1), i.e., <italic>y = x</italic>. In turn, the Gini coefficient (GC) is used to quantify the extent to which the Lorenz curve deviates from the line of absolute equality, which equals twice the area formed by the Lorenz curve and the line of absolute equality. The closer the GC is to 0, the more equal the income distribution tends to be. More recently, the GC has been used in some botanical studies due to its ability to quantify the inequality of urban green space, plant, and plant organ size distributions (<xref ref-type="bibr" rid="B27">Taylor and Aarssen, 1989</xref>; <xref ref-type="bibr" rid="B10">Matlack, 1994</xref>; <xref ref-type="bibr" rid="B11">Metsaranta and Lieffers, 2008</xref>; <xref ref-type="bibr" rid="B2">Chen et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B15">Nero, 2017</xref>; <xref ref-type="bibr" rid="B4">Huang et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B8">Lian et&#xa0;al., 2023</xref>)&#x2014;for example, the unevenness of leaf size distributions can be quantified by the GC because the Lorenz curve can quantify the relationship between the cumulative proportion of leaf area and the cumulative proportion of the number of leaves.</p>
<p>The accuracy of this approach has been investigated using a performance equation (PE) and its generalized version (GPE) to fit the rotated and right-shifted Lorenz curves of the organ size distributions of herbaceous species (<xref ref-type="bibr" rid="B4">Huang et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B8">Lian et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B21">Shi et&#xa0;al., 2024</xref>; <xref ref-type="bibr" rid="B28">Wang et&#xa0;al., 2024</xref>)&#x2014;for example, <xref ref-type="bibr" rid="B28">Wang et&#xa0;al. (2024)</xref> compared the PE and the GPE with two other performance equations based on a temperature-dependent square root equation proposed by <xref ref-type="bibr" rid="B16">Ratkowsky et&#xa0;al. (1983)</xref> and found that the GPE had the best goodness of fit and lowest Akaike information criterion (AIC). In recent years, researchers have compared Lorenz equations in quantifying the inequality or non-uniformity of the size distributions of abiotic as well as biotic quantities (<xref ref-type="bibr" rid="B18">Sarabia, 1997</xref>; <xref ref-type="bibr" rid="B19">Sarabia et&#xa0;al., 1999</xref>; <xref ref-type="bibr" rid="B25">Sitthiyot and Holasut, 2023</xref>). However, prior work has not compared the two performance equations that describe the rotated and right-shifted Lorenz equations with other Lorenz equations, perhaps because the parameters of the Lorenz equations are estimated by minimizing the residual sum of squares (RSS) between the observed and predicted cumulative proportions of plant (or organ) size, whereas the parameters of the PE and the GPE are estimated by minimizing the residual sum of squares (RSS) between the observed and predicted cumulative proportions of plant (organ) size when rotated and right-shifted. Consequently, there is a need to develop a method that can directly compare the two performance equations and other Lorenz equations to evaluate which equation is more accurate when quantifying the inequality of plant (or organ) size distributions.</p>
<p>In the present study, we measured the leaf length (<italic>L</italic>) and width (<italic>W</italic>) data of 752 leaves from 60 <italic>Alangium chinense</italic> (Lour.) Harms saplings. The leaf area (<italic>A</italic>) was determined using the ME, which assumes a proportional relationship between <italic>A</italic> and <italic>LW</italic>. In addition, 540 leaves were scanned to determine the <italic>A</italic>, <italic>L</italic>, and <italic>W</italic> of the 540 leaves empirically. We employed five models to fit the rotated and right-shifted Lorenz curves of individual leaf area distributions (i.e., the cumulative proportion of leaf area per sapling vs. the cumulative proportion of the number of leaves per sapling). Our goals were to address two questions: (i) Can the PE and GPE be used to describe the rotated and right-shifted Lorenz curve of leaf area distributions of woody species? (ii) Are the two performance equations superior to other Lorenz equations used in nonlinear regression? If true, this approach provides a non-destructive method to quantify leaf area and inequalities in organ size distributions of both woody and herbaceous species, which could provide deeper insights into extant (and possibly extinct) plant ecology and evolution.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Leaf sampling information and data acquisition</title>
<p>Dataset 1: A total of 540 fresh, mature, and undamaged leaves of <italic>A. chinense</italic> were sampled from the middle canopy of 10 saplings with a diameter at breast height (DBH) of 2 to 7&#xa0;cm growing in the Nanjing Forestry University campus, Nanjing, China (32&#xb0;4&#x2032;48&#x2033; N, 118&#xb0;49&#x2032;12&#x2033; E) on July 31, 2019. Leaves in groups of four were placed in plastic self-sealing bags (28&#xa0;cm &#xd7; 20&#xa0;cm) to avoid tissue dehydration and taken to the laboratory within 1&#xa0;h. Laminae were scanned to bitmap images at 600-dpi resolution using a photo scanner (Aficio MP 7502; Ricoh, Tokyo, Japan). ImageJ software (<ext-link ext-link-type="uri" xlink:href="https://imagej.nih.gov/ij/index.html">https://imagej.nih.gov/ij/index.html</ext-link>) was used to measure the area (<italic>A</italic>), length (<italic>L</italic>), and width (<italic>W</italic>) of each leaf (see online <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table S1</bold>
</xref>). <italic>W</italic> was defined by the maximum distance between two lobe apices (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>).</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Freehand drawing of the above-ground part of <italic>A. chinensis</italic> and definition of leaf length (<italic>L</italic>) and width (<italic>W</italic>). <italic>L</italic> was defined as the distance between leaf tip and the connection point between the lamina and the petiole, and <italic>W</italic> was defined as the distance between the two leaf apices.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1426424-g001.tif"/>
</fig>
<p>Dataset 2: On October 17, 2023 and on May 18, 2024, we randomly selected 60&#xa0;A<italic>. chinense</italic> saplings, each with at least nine leaves, to provide a sample size of not less than 1.5 times the minimum data required to perform nonlinear regression analyses to estimate the regression parameters of <xref ref-type="disp-formula" rid="eq4">Equation 4</xref> (see below for details). The above-ground height of the saplings ranged between 30.6 and 136.4&#xa0;cm, with a mean &#xb1; standard error of 81.3 &#xb1; 28.5&#xa0;cm. We measured the <italic>L</italic> and <italic>W</italic> of each leaf of the 60 saplings without removing the leaves from the branches. A total of 752 leaves from the 60 saplings were examined (see online <xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Table S2</bold>
</xref>).</p>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Estimate of leaf area based on leaf length and width</title>
<p>
<xref ref-type="bibr" rid="B12">Montgomery (1911)</xref> proposed a formula to estimate the leaf area (<italic>A</italic>) of corn (<italic>Zea mays</italic> Linn.) using the product of a proportionality coefficient (<italic>k</italic>) and leaf length (<italic>L</italic>) and width (<italic>W</italic>):</p>
<disp-formula id="eq1">
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>k</mml:mi>
<mml:mi>L</mml:mi>
<mml:mi>W</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p>which is denoted as Montgomery equation (ME), and <italic>k</italic> is referred to as Montgomery parameter (MP) hereinafter. To normalize <italic>A</italic>, <xref ref-type="disp-formula" rid="eq1">Equation 1</xref> was log-transformed:</p>
<disp-formula id="eq2">
<label>(2)</label>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>log</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mi>A</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mo>=</mml:mo>
<mml:mi>c</mml:mi>
<mml:mo>+</mml:mo>
<mml:mtext>log</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>L</mml:mi>
<mml:mi>W</mml:mi>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>c</italic> is a constant to be estimated. It is apparent that <italic>c</italic> = log<italic>k</italic>. <xref ref-type="disp-formula" rid="eq2">Equation 2</xref> was used to fit the 540 empirical observations of <italic>A</italic> vs. <italic>LW</italic> and thereby determine the numerical value of <italic>c</italic>.</p>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Equations to describe the rotated and right-shifted Lorenz curve</title>
<p>The cumulative proportions of leaf area per sapling vs. the cumulative proportion of the number of leaves per sapling were rotated counterclockwise by 135&#xb0; around the origin and shifted to the right by a distance of <inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>. Previous studies have used two performance equations for fitting the size frequency distributions of leaf size and fruit volume data of herbaceous plants (<xref ref-type="bibr" rid="B4">Huang et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B8">Lian et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B28">Wang et&#xa0;al., 2024</xref>):</p>
<p>(i) The performance equation, denoted as PE, which was initially used to describe the jumping distance of green frogs at different body temperatures (<xref ref-type="bibr" rid="B5">Huey, 1975</xref>; <xref ref-type="bibr" rid="B6">Huey and Stevenson, 1979</xref>):</p>
<disp-formula id="eq3">
<label>(3)</label>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>c</mml:mi>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>In the context of our study, <italic>x</italic> and <italic>y</italic> represent the rotated and right-shifted cumulative proportion of the number of leaves per sapling and the cumulative proportion of leaf area per sapling, respectively; <italic>c</italic>, <italic>K</italic>
<sub>1</sub>, and <italic>K</italic>
<sub>2</sub> are constants to be estimated; <italic>x</italic>
<sub>1</sub> and <italic>x</italic>
<sub>2</sub> are equal to 0 and <inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>, respectively (<xref ref-type="bibr" rid="B8">Lian et&#xa0;al., 2023</xref>; <xref ref-type="bibr" rid="B28">Wang et&#xa0;al., 2024</xref>).</p>
<p>(ii) The generalized performance equation, denoted as GPE, which increases the flexibility of data fitting by adding two parameters, <italic>a</italic> and <italic>b</italic>, to PE (<xref ref-type="bibr" rid="B8">Lian et&#xa0;al., 2023</xref>):</p>
<disp-formula id="eq4">
<label>(4)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:mi>y</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>c</mml:mi>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>a</mml:mi>
</mml:msup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mi>e</mml:mi>
<mml:mrow>
<mml:msub>
<mml:mi>K</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>b</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>In the context of our study, <italic>x</italic> and <italic>y</italic> represent the rotated and right-shifted cumulative proportion of the number of leaves per sapling and the cumulative proportion of leaf area per sapling, respectively; <italic>x</italic>
<sub>1</sub>, and <italic>x</italic>
<sub>2</sub> are also equal to 0 and <inline-formula>
<mml:math display="inline" id="im3">
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula> (<xref ref-type="bibr" rid="B8">Lian et&#xa0;al., 2023</xref>).</p>
<p>
<xref ref-type="disp-formula" rid="eq3">Equations 3</xref> and <xref ref-type="disp-formula" rid="eq4">4</xref> were used to evaluate the performance of Lorenz curve fitting. In terms of the latter, <xref ref-type="bibr" rid="B18">Sarabia (1997)</xref> proposed a highly flexible parameter family of Lorenz curves based on the generalized Tukey &#x3bb; distribution for fitting Lorenz curves:</p>
<disp-formula id="eq5">
<label>(5)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3bb;</mml:mtext>
<mml:mo>+</mml:mo>
<mml:mtext>&#x3b7;</mml:mtext>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mtext>&#x3bb;</mml:mtext>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mrow><mml:mi>l</mml:mi></mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3b7;</mml:mtext>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:math>
</disp-formula>
<p>In the context of our study, <italic>x<sub>l</sub>
</italic> and <italic>y<sub>l</sub>
</italic> represent the cumulative proportion of the number of leaves per sapling and the cumulative proportion of leaf area per sapling, respectively, and &#x3bb;, &#x3b7;, <italic>a</italic>
<sub>1</sub>, and <italic>a</italic>
<sub>2</sub> are constants to be estimated, where <inline-formula>
<mml:math display="inline" id="im4">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>1</mml:mn>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im5">
<mml:mrow>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im6">
<mml:mrow>
<mml:mtext>&#x3b7;</mml:mtext>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>+</mml:mo>
<mml:mtext>&#x3bb;</mml:mtext>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im7">
<mml:mrow>
<mml:mtext>&#x3bb;</mml:mtext>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im8">
<mml:mrow>
<mml:mtext>&#x3b7;</mml:mtext>
<mml:msub>
<mml:mi>a</mml:mi>
<mml:mn>2</mml:mn>
</mml:msub>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. <xref ref-type="disp-formula" rid="eq5">Equation 5</xref> is denoted as SarabiaE hereinafter.</p>
<p>
<xref ref-type="bibr" rid="B19">Sarabia et&#xa0;al. (1999)</xref> proposed a general method for building parametric-functional families of Lorenz curves generated from an initial Lorenz curve, which satisfies some regularity conditions:</p>
<disp-formula id="eq6">
<label>(6)</label>
<mml:math display="block" id="M6">
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:msubsup>
<mml:mi>x</mml:mi>
<mml:mrow><mml:mi>l</mml:mi></mml:mrow>
<mml:mrow><mml:mtext>&#x3b3;</mml:mtext></mml:mrow>
</mml:msubsup>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">[</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mtext>&#x3b1;</mml:mtext>
</mml:msup>
</mml:mrow>
<mml:mo stretchy="false">]</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mtext>&#x3b2;</mml:mtext>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>In the context of our study, <italic>x<sub>l</sub>
</italic> and <italic>y<sub>l</sub>
</italic> represent the cumulative proportion of the number of leaves per sapling and the cumulative proportion of leaf area per sapling, respectively, and &#x3b1;, &#x3b2;, and &#x3b3; are constants to be estimated, where <inline-formula>
<mml:math display="inline" id="im9">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&lt;</mml:mo>
<mml:mtext>&#x3b1;</mml:mtext>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im10">
<mml:mrow>
<mml:mtext>&#x3b2;</mml:mtext>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im11">
<mml:mrow>
<mml:mtext>&#x3b3;</mml:mtext>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. <xref ref-type="disp-formula" rid="eq6">Equation 6</xref> is denoted as SCSE hereinafter.</p>
<p>Importantly, the existing parameter function forms of Lorenz curve are not suitable for extreme inequalities in size distributions. To fit sample data with typical convex segments better in a Lorenz curve, a universal function for fitting Lorenz curves was proposed (<xref ref-type="bibr" rid="B25">Sitthiyot and Holasut, 2023</xref>):</p>
<disp-formula id="eq7">
<label>(7)</label>
<mml:math display="block" id="M7">
<mml:mtable>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3c1;</mml:mtext>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mn>2</mml:mn>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>+</mml:mo>
<mml:mtext>&#x3c1;</mml:mtext>
<mml:mrow>
<mml:mo>[</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3c9;</mml:mtext>
</mml:mrow>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mi>P</mml:mi>
</mml:msup>
<mml:mo>+</mml:mo>
<mml:mtext>&#x3c9;</mml:mtext>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mrow>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mrow>
<mml:mo>(</mml:mo>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
<mml:mo>)</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mi>P</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:msup>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:mrow>
<mml:mo>]</mml:mo>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:math>
</disp-formula>
<p>when <inline-formula>
<mml:math display="inline" id="im12">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&gt;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula>
<mml:math display="inline" id="im13">
<mml:mrow>
<mml:msub>
<mml:mi>y</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mn>0</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, when <inline-formula>
<mml:math display="inline" id="im14">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>l</mml:mi>
</mml:msub>
<mml:mo>&#x2264;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
</mml:mrow>
</mml:math>
</inline-formula>. Here <italic>x<sub>l</sub>
</italic> and <italic>y<sub>l</sub>
</italic> represent the cumulative proportion of the number of leaves per sapling and the cumulative proportion of leaf area per sapling, respectively, and &#x3b4;, &#x3c1;, &#x3c9;, and <italic>P</italic> are constants to be estimated, where <inline-formula>
<mml:math display="inline" id="im15">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:mtext>&#x3b4;</mml:mtext>
<mml:mo>&lt;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im16">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:mtext>&#x3c1;</mml:mtext>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula>
<mml:math display="inline" id="im17">
<mml:mrow>
<mml:mn>0</mml:mn>
<mml:mo>&#x2264;</mml:mo>
<mml:mtext>&#x3c9;</mml:mtext>
<mml:mo>&#x2264;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>, and <inline-formula>
<mml:math display="inline" id="im18">
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#x2265;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
</mml:math>
</inline-formula>. <xref ref-type="disp-formula" rid="eq7">Equation 7</xref> is denoted as SHE hereinafter.</p>
<p>We rotated and right-shifted the three foregoing Lorenz equations [i.e., <xref ref-type="disp-formula" rid="eq5">Equations 5</xref>-<xref ref-type="disp-formula" rid="eq7">7</xref>] to evaluate the performance of the PE and GPE with these three Lorenz equations in describing the rotated and right-shifted data of the cumulative proportion of leaf area per sapling vs. the cumulative proportion of the number of leaves per sapling.</p>
</sec>
<sec id="s2_4">
<label>2.4</label>
<title>Data fitting and model evaluation</title>
<p>We used the PE, GPE, SarabiaE, SCSE, and SHE equations to fit the empirical data after the data were rotated counterclockwise by 135&#xb0; and shifted to the right by a distance of <inline-formula>
<mml:math display="inline" id="im19">
<mml:mrow>
<mml:msqrt>
<mml:mn>2</mml:mn>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>. The Nelder&#x2013;Mead optimization algorithm (<xref ref-type="bibr" rid="B14">Nelder and Mead, 1965</xref>) was then used to minimize the fitting criteria for nonlinear regression and the residual sum of squares (RSS) between empirical and predicted <italic>y</italic>-values to estimate the parameters of each model. To better evaluate the goodness of fit of nonlinear regression, the root-mean-square error (RMSE) was used to measure the prediction accuracy. The smaller the RMSE, the higher the prediction accuracy of the model. The &#x201c;fitLorenz&#x201d; function in the &#x201c;biogeom&#x201d; package (version 1.4.3; <xref ref-type="bibr" rid="B22">Shi et&#xa0;al., 2022</xref>) based on statistical software R (version 4.3.3; <xref ref-type="bibr" rid="B17">R Core Team, 2024</xref>) was used to fit the leaf area distribution:</p>
<disp-formula id="eq8">
<label>(8)</label>
<mml:math display="block" id="M8">
<mml:mrow>
<mml:mtext>RMSE</mml:mtext>
<mml:mo>=</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mtext>RSS</mml:mtext>
<mml:mo stretchy="false">/</mml:mo>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
</disp-formula>
<p>where <italic>n</italic> represents the number of leaves per sapling. The Akaike information criterion (AIC) that compared to the adjusted coefficient of determination is frequently recommended in nonlinear regression to balance the goodness of fit and the model constructive complexity (<xref ref-type="bibr" rid="B26">Spiess and Neumeyer, 2010</xref>). The model with the smallest AIC is considered as the optimal model. Paired <italic>t</italic>-test at a significance level of 0.05 was used to test the significance of differences in the RMSE or AIC values between any two of the five equations.</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<p>The number of leaves per sapling ranged between nine and 25, with a mean &#xb1; standard error of 12.5 &#xb1; 3.8. The Montgomery equation (ME) was validated in estimating <italic>A</italic> by multiplying the leaf length and width by the estimated <italic>k</italic>. The correlation coefficient (<italic>r</italic>) between <italic>A</italic> and <italic>LW</italic> was 0.9726, and the RMSE (see <xref ref-type="disp-formula" rid="eq8">Equation 8</xref>) of linear regression was 0.1557 (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). The estimated value of the proportionality coefficient <italic>k</italic> was 0.6469. Using this value of <italic>k</italic>, the area of each leaf of the 60 saplings was obtained. The total leaf area per sapling ranged between 564 and 2,875 cm<sup>2</sup>, with a mean &#xb1; standard of 1,391 &#xb1; 548 cm<sup>2</sup>. The mean leaf area of each sapling ranged between 56 and 183 cm<sup>2</sup>, with a mean &#xb1; standard error of 111 &#xb1; 28 cm<sup>2</sup>.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Results of fitting the Montgomery equation that assumes a proportional relationship between <italic>A</italic> and <italic>LW</italic>, which represent the leaf area, length, and width, respectively. Here RMSE represents the root-mean-square error; <italic>r</italic> represents the correlation coefficient; <italic>n</italic> represents the number of samples; <inline-formula>
<mml:math display="inline" id="im20">
<mml:mover accent="true">
<mml:mi>k</mml:mi>
<mml:mo>^</mml:mo>
</mml:mover>
</mml:math>
</inline-formula> represents the estimated value of the Montgomery parameter, i.e., the proportionality coefficient of the Montgomery equation; 95% CI represents 95% confidence intervals of the Montgomery parameter based on 3,000 bootstrap replicates.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1426424-g002.tif"/>
</fig>
<p>The PE, GPE, SarabiaE, SCSE, and SHE fitted the rotated and right-shifted data of the leaf area distributions of the 60 saplings well, with RMSE values &lt;0.013. However, the GPE had the lowest RMSE and AIC values among the five equations. The SarabiaE and PE worked the second best and were better than SCSE and SHE (<xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>). The RMSE and AIC values of GPE were both significantly lower than those of the other equations (<italic>p</italic> &lt; 0.05). <xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref> shows the fitted results of the five equations to the rotated and right-shifted data of leaf area distribution of one sapling example.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>
<bold>(A, B)</bold> Boxplot of the root-mean-square error between any two of the five equations (i.e., PE, GPE, SarabiaE, SCSE, and SHE) for 60 datasets. Paired <italic>t</italic>-test was used to determine the significance of differences between any two equations at the 0.05 significance level, and different letters represent a significant difference of any two equations. The vertical solid line in each box represents the median.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1426424-g003.tif"/>
</fig>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>
<bold>(A&#x2013;D)</bold> Comparison of the observed and predicted data of the rotated and right-shifted Lorenz curve for the leaf area distribution of a representative individual sapling of <italic>A. chinensis</italic>. Data points represent observations; curves represent predicted values. Specifically, in <bold>(A)</bold>, the blue curve represents the estimated PE, and the red curve represents the estimated GPE. <italic>n</italic> represents the number of leaves on this individual sapling; AIC represents the Akaike information criterion of the corresponding nonlinear equation in each panel; RMSE represents the root-mean-square error.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-15-1426424-g004.tif"/>
</fig>
</sec>
<sec id="s4" sec-type="discussion">
<label>4</label>
<title>Discussion</title>
<p>The goals of this paper were to compare equations describing the inequality of plant size distributions and to evaluate which among the contending equations perform the best. It is important, therefore, to review the criteria used to evaluate the equations investigated in this study. A secondary but equally important goal was to demonstrate the utility of the approach taken in this and other studies using the Montgomery equation. These goals are discussed in the following two sections.</p>
<sec id="s4_1">
<label>4.1</label>
<title>Criteria to evaluate competitive equations</title>
<p>When evaluating the superiority of competitive nonlinear models, researchers often use the root-mean-square error (RMSE) and Akaike information criterion (AIC) (<xref ref-type="bibr" rid="B1">Akaike, 1998</xref>). RMSE can be used not only to measure the prediction accuracy of a model on continuous data but also to measure the average degree of deviation between the predicted and observed values. A smaller RMSE indicates a more accurate prediction of the model. The AIC is based on the concept of entropy and can balance the complexity of the estimated model with the goodness of fit of the estimated model. It is apparent that increasing the number of free parameters in a model can improve the goodness of fit and enhance the flexibility of a model to fit the data. Thus, the value of the likelihood function increases as the model complexity increases, which leads to a smaller AIC. Conversely, when the increased rate of the likelihood function slows down, it can lead to an increase in AIC, i.e., AIC encourages the goodness of fit but avoids an overfitting, if the model is complex. We argue that the priority of model selection should focus on the lowest AIC value.</p>
<p>In the present study, the performance equation (PE) and its generalized version (GPE) are validated based on 60 empirical datasets (i.e., the leaves removed from 60 saplings). The RMSE and AIC values of GPE were significantly smaller than those of the four other contending equations. Specifically, GPE worked best, and PE worked second best in that it had the second lowest AIC value. Both were superior to the SarabiaE, SCSE, and SHE. Clearly, the GPE has more parameters than PE. However, PE contains fewer parameters that can reduce the probability of overfitting. Although GPE has the lowest RMSE and AIC values, PE is still recommended due to its fewer parameters and validity in fitting the rotated and right-shifted data. Only for extremely right-skewed distributions can SHE perform better than the two performance equations (<xref ref-type="bibr" rid="B21">Shi et&#xa0;al., 2024</xref>).</p>
</sec>
<sec id="s4_2">
<label>4.2</label>
<title>The general applicability of the Montgomery&#x2013;Lorenz Performance approach</title>
<p>A pivotal goal of this paper was to evaluate a very specific and targeted approach to quantifying the size frequency distributions of plants and plant organs (in this case, the leaves of saplings). The approach presented here, which can be called the Montgomery-Lorenz-Performance (MLP) approach, involves a threefold mathematical methodology. The Montgomery equation allows us to quantify leaf area (as well as many other variables of interest) non-destructively once the Montgomery parameter (i.e., the proportionality coefficient of ME) is determined (<xref ref-type="bibr" rid="B20">Schrader et&#xa0;al., 2021</xref>), which requires a comparatively small number of samples from a population of leaves. The Lorenz curve, in tandem with the Gini coefficient, as a methodology, allows us to quantify the evenness/unevenness of any size frequency distribution. Lastly, the performance equation and its generalized form allow us to evaluate the reliability of the quantification of size frequency distribution evenness/unevenness.</p>
<p>At issue is whether this overall approach has general applicability. We argue that the evidence indicates that the methodology has general applicability because the mathematics underlying each component of the overall approach has been validated here and elsewhere both in terms of its mathematical rigor and its empirical predictive ability. Because of its analytical rigor and empirical robustness, the methodology presented here and elsewhere has application to many aspects of ecological enquiry, not just to plant organs, because the metrics used to measure &#x201c;size&#x201d; (e.g., length, width, volume, and biomass) can be applied to any cellular form of life, extant or extinct.</p>
</sec>
</sec>
<sec id="s5" sec-type="conclusions">
<label>5</label>
<title>Conclusions</title>
<p>We compared the two performance equations (i.e., PE and GPE) with three other Lorenz equations (i.e., SarabiaE, SCSE, and SHE) in describing the rotated and right-shifted data of the cumulative proportion of leaf area per sapling vs. the number of leaves per sapling for a representative broad-leaved plant (i.e., <italic>Alangium chinensis</italic>). The five equations were found to be valid in describing the rotated and right-shifted Lorenz curve of the leaf area distributions of this tree species. However, the GPE was found to be superior, and PE worked the second best compared with the other three equations based on the comparison of their AIC values. This refined method provides a general protocol for quantifying the inequality of any organic size distribution, which is of great significance in analyzing the allocation of organismic resources.</p>
</sec>
<sec id="s6" sec-type="data-availability">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="supplementary-material" rid="SM1">
<bold>Supplementary Material</bold>
</xref>. Further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7" sec-type="author-contributions">
<title>Author contributions</title>
<p>LD: Formal analysis, Methodology, Writing &#x2013; original draft. KH: Investigation, Writing &#x2013; original draft. KN: Formal analysis, Methodology, Writing &#x2013; review &amp; editing. ZS: Investigation, Writing &#x2013; review &amp; editing. YM: Investigation, Writing &#x2013; review &amp; editing. PS: Formal analysis, Methodology, Writing &#x2013; review &amp; editing.</p>
</sec>
</body>
<back>
<sec id="s8" sec-type="funding-information">
<title>Funding</title>
<p>The author(s) declare that no financial support was received for the research, authorship, and/or publication of this article.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>We are grateful to the handling editor and three reviewers for their review work. We thank Ms. Mengdi Liu for her prior leaf sampling work.</p>
</ack>
<sec id="s9" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
<p>The handling editor JG declared a past co-authorship with the author(s) KN and PS.</p>
</sec>
<sec id="s10" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s11" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2024.1426424/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2024.1426424/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Table_1.xlsx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Akaike</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>1998</year>). &#x201c;<article-title>Information theory and an extension of the maximum likelihood principle</article-title>,&#x201d; in <source>Selected papers of Hirotugu Akaike</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Parzen</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Tanabe</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Kitagawa</surname> <given-names>G.</given-names>
</name>
</person-group> (<publisher-loc>New York</publisher-loc>:<publisher-name>Springer</publisher-name>), <fpage>199</fpage>&#x2013;<lpage>213</lpage>.</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>B. J. W.</given-names>
</name>
<name>
<surname>During</surname> <given-names>H. J.</given-names>
</name>
<name>
<surname>Vermeulen</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Anten</surname> <given-names>N. P. R.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>The presence of a below-ground neighbour alters within-plant seed size distribution in <italic>Phaseolus vulgaris</italic>
</article-title>. <source>Ann. Bot.</source> <volume>114</volume>, <fpage>937</fpage>&#x2013;<lpage>943</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/aob/mcu162</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hodgson</surname> <given-names>J. G.</given-names>
</name>
<name>
<surname>Wilson</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Hunt</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Grime</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Thompson</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Allocating C-S-R plan functional types: a soft approach to a hard problem</article-title>. <source>Oikos</source> <volume>85</volume>, <fpage>282</fpage>&#x2013;<lpage>294</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/3546494</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>L. C.</given-names>
</name>
<name>
<surname>Ratkowsky</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Hui</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Gielis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Lian</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>W. H.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Inequality measure of leaf area distribution for a drought-tolerant landscape plant</article-title>. <source>Plants</source> <volume>12</volume>, <elocation-id>3143</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/plants12173143</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Huey</surname> <given-names>R. B.</given-names>
</name>
</person-group> (<year>1975</year>). <source>Ecology of lizard thermoregulation</source> (<publisher-loc>Cambridge</publisher-loc>: <publisher-name>Harvard University</publisher-name>).</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huey</surname> <given-names>R. B.</given-names>
</name>
<name>
<surname>Stevenson</surname> <given-names>R. D.</given-names>
</name>
</person-group> (<year>1979</year>). <article-title>Integrating thermal physiology and ecology of ectotherms: a discussion of approaches</article-title>. <source>Am. Zool.</source> <volume>19</volume>, <fpage>357</fpage>&#x2013;<lpage>366</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/icb/19.1.357</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kakwani</surname> <given-names>N. C.</given-names>
</name>
</person-group> (<year>1977</year>). <article-title>Applications of Lorenz curves in economic analysis</article-title>. <source>Econometrica</source> <volume>45</volume>, <fpage>719</fpage>&#x2013;<lpage>728</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/1911684</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lian</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L. Y.</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>W. H.</given-names>
</name>
<name>
<surname>Gielis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Niklas</surname> <given-names>K. J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A generalized performance equation and its application in measuring the Gini index of leaf size inequality</article-title>. <source>Trees</source> <volume>37</volume>, <fpage>1555</fpage>&#x2013;<lpage>1565</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00468-023-02448-8</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lorenz</surname> <given-names>M. O.</given-names>
</name>
</person-group> (<year>1905</year>). <article-title>Methods of measuring the concentration of wealth</article-title>. <source>J. Am. Stat. Assoc.</source> <volume>9</volume>, <fpage>209</fpage>&#x2013;<lpage>219</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/15225437.1905.10503443</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Matlack</surname> <given-names>G. L.</given-names>
</name>
</person-group> (<year>1994</year>). <article-title>Vegetation dynamics of the forest edge&#x2013;trends in space and successional time</article-title>. <source>J. Ecol.</source> <volume>82</volume>, <fpage>113</fpage>&#x2013;<lpage>123</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/2261391</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Metsaranta</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Lieffers</surname> <given-names>V. J.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Inequality of size and size increment in <italic>Pinus banksiana</italic> in relation to stand dynamics and annual growth rate</article-title>. <source>Ann. Bot.</source> <volume>101</volume>, <fpage>561</fpage>&#x2013;<lpage>571</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/aob/mcm320</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Montgomery</surname> <given-names>E. G.</given-names>
</name>
</person-group> (<year>1911</year>). <article-title>Correlation studies in corn</article-title>. <source>Nebraska Agric. Exp. Stn. Annu. Rep.</source> <volume>24</volume>, <fpage>108</fpage>&#x2013;<lpage>159</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mu</surname> <given-names>Y. Y.</given-names>
</name>
<name>
<surname>He</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>L. L.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>Z. Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>comparation between computer recognition and manual measurement methods for the estimate on of leaf area</article-title>. <source>Ann. Bot</source>.  In press.  doi:&#xa0;<pub-id pub-id-type="doi">10.1093/aob/mcae090</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nelder</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Mead</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>1965</year>). <article-title>A simplex method for function minimization</article-title>. <source>Comput. J.</source> <volume>7</volume>, <fpage>308</fpage>&#x2013;<lpage>313</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/comjnl/7.4.308</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nero</surname> <given-names>B. F.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Urban green space dynamics and socio-environmental inequity: Multi-resolution and spatiotemporal data analysis of Kumasi, Ghana</article-title>. <source>Int. J. Remote. Sens.</source> <volume>38</volume>, <fpage>6993</fpage>&#x2013;<lpage>7020</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/01431161.2017.1370152</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ratkowsky</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Lowry</surname> <given-names>R. K.</given-names>
</name>
<name>
<surname>McMeekin</surname> <given-names>T. A.</given-names>
</name>
<name>
<surname>Stokes</surname> <given-names>A. N.</given-names>
</name>
<name>
<surname>Chandler</surname> <given-names>R. E.</given-names>
</name>
</person-group> (<year>1983</year>). <article-title>Model for bacterial culture growth rate throughout the entire biokinetic temperature range</article-title>. <source>J. Bacteriol.</source> <volume>154</volume>, <fpage>1222</fpage>&#x2013;<lpage>1226</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1128/jb.154.3.1222-1226.1983</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="web">
<person-group person-group-type="author">
<collab>R Core Team</collab>
</person-group> (<year>2024</year>). <article-title>R: A Language and Environment for Statistical Computing</article-title>. Available online at: <uri xlink:href="https://www.rproject.org/">https://www.rproject.org/</uri> (Accessed <access-date>1 April 2024</access-date>).</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sarabia</surname> <given-names>J. M.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>A hierarchy of Lorenz curves based on the generalized tukey's lambda distribution</article-title>. <source>Econom. Rev.</source> <volume>16</volume>, <fpage>305</fpage>&#x2013;<lpage>320</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/07474939708800389</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sarabia</surname> <given-names>J. M.</given-names>
</name>
<name>
<surname>Castillo</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Slottje</surname> <given-names>D. J.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>An ordered family of Lorenz curves</article-title>. <source>J. Econometrics</source> <volume>91</volume>, <fpage>43</fpage>&#x2013;<lpage>60</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S0304-4076(98)00048-7</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schrader</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Royer</surname> <given-names>D. L.</given-names>
</name>
<name>
<surname>Peppe</surname> <given-names>D. J.</given-names>
</name>
<name>
<surname>Gallagher</surname> <given-names>R. V.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Leaf size estimation based on leaf length, width and shape</article-title>. <source>Ann. Bot.</source> <volume>128</volume>, <fpage>395</fpage>&#x2013;<lpage>406</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/aob/mcab078</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Deng</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Niklas</surname> <given-names>K. J.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Rotated Lorenz curves of biological size distributions follow two performance equations</article-title>. <source>Symmetry</source> <volume>16</volume>, <elocation-id>565</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/sym16050565</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Gielis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Quinn</surname> <given-names>B. K.</given-names>
</name>
<name>
<surname>Niklas</surname> <given-names>K. J.</given-names>
</name>
<name>
<surname>Ratkowsky</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Schrader</surname> <given-names>J.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>&#x2018;biogeom&#x2019;: an R package for simulating and fitting natural shapes</article-title>. <source>Ann. N. Y. Acad. Sci.</source> <volume>1516</volume>, <fpage>123</fpage>&#x2013;<lpage>134</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/nyas.14862</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>M. D.</given-names>
</name>
<name>
<surname>Ratkowsky</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Gielis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Su</surname> <given-names>J. L.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>X. J.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>a). <article-title>Leaf area&#x2013;length allometry and its implications in leaf shape evolution</article-title>. <source>Trees Struct. Funct.</source> <volume>33</volume>, <fpage>1073</fpage>&#x2013;<lpage>1085</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s00468-019-01843-4</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>M. D.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>X. J.</given-names>
</name>
<name>
<surname>Gielis</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Ratkowsky</surname> <given-names>D. A.</given-names>
</name>
</person-group> (<year>2019</year>b). <article-title>Proportional relationship between leaf area and the product of leaf length width of four types of special leaf shapes</article-title>. <source>Forests</source> <volume>10</volume>, <elocation-id>178</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/f10020178</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sitthiyot</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Holasut</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>A universal model for the Lorenz curve with novel applications for datasets containing zeros and/or exhibiting extreme inequality</article-title>. <source>Sci. Rep.</source> <volume>13</volume>, <fpage>4729</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-023-31827-x</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Spiess</surname> <given-names>A.-N.</given-names>
</name>
<name>
<surname>Neumeyer</surname> <given-names>N.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>An evaluation of <italic>R</italic>
<sup>2</sup> as an inadequate measure for nonlinear models in pharmacological and biochemical research: a Monte Carlo approach</article-title>. <source>BMC Pharmacol.</source> <volume>10</volume>, <elocation-id>6</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/1471-2210-10-6</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Taylor</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Aarssen</surname> <given-names>L. W.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Neighbor effects in mast year seedlings of <italic>Acer saccharum</italic>
</article-title>. <source>Am. J. Bot.</source> <volume>76</volume>, <fpage>546</fpage>&#x2013;<lpage>554</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/j.1537-2197.1989.tb11346.x</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>He</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Hui</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Ratkowsky</surname> <given-names>D. A.</given-names>
</name>
<name>
<surname>Yao</surname> <given-names>W. H.</given-names>
</name>
<name>
<surname>Lian</surname> <given-names>M.</given-names>
</name>
<etal/>
</person-group>. (<year>2024</year>). <article-title>Comparison of four performance models in quantifying the inequality of leaf and fruit size distribution</article-title>. <source>Ecol. Evol.</source> <volume>14</volume>, <elocation-id>e11072</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/ece3.11072</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Westoby</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>A leaf-height-seed (LHS) plant ecology strategy scheme</article-title>. <source>Plant Soil</source> <volume>199</volume>, <fpage>213</fpage>&#x2013;<lpage>227</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1023/A:1004327224729</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname> <given-names>X. J.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>P. J.</given-names>
</name>
<name>
<surname>Schrader</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Niklas</surname> <given-names>K. J.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Nondestructive estimation of leaf area for 15 species of vines with different leaf shapes</article-title>. <source>Am. J. Bot.</source> <volume>107</volume>, <fpage>1481</fpage>&#x2013;<lpage>1490</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/ajb2.1560</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zirbel</surname> <given-names>C. R.</given-names>
</name>
<name>
<surname>Bassett</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Grman</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Brudvig</surname> <given-names>L. A.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Plant functional traits and environmental conditions shape community assembly and ecosystem functioning during restoration</article-title>. <source>J. Appl. Ecol.</source> <volume>54</volume>, <fpage>1070</fpage>&#x2013;<lpage>1079</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/1365-2664.12885</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>