# Examination of the Effects of Curvature in Geometrical Space on Accuracy of Scaling Derived Projections of Plant Biomass Units: Applications to the Assessment of Average Leaf Biomass in Eelgrass Shoots

Authors: Héctor Echavarría-Heras; Cecilia Leal-Ramírez; Enrique Villa-Diharce; Abelardo Montesinos-López
Journal: BioMed Research International
 (2019)

---

## Abstract

Conservation of eelgrass relies on transplants and evaluation of success depends on nondestructive measurements of average leaf biomass in shoots among other variables. Allometric proxies offer a convenient way to assessments. Identifying surrogates via log transformation and linear regression can set biased results. Views conceive this approach to be meaningful, asserting that curvature in geometrical space explains bias. Inappropriateness of correction factor of retransformation bias could also explain inconsistencies. Accounting for nonlinearity of the log transformed response relied on a generalized allometric model. Scaling parameters depend continuously on the descriptor. Joining correction factor is conceived as the partial sum of series expansion of mean retransformed residuals leading to highest reproducibility strength. Fits of particular characterizations of the generalized curvature model conveyed outstanding reproducibility of average eelgrass leaf biomass in shoots. Although nonlinear heteroscedastic regression resulted also to be suitable, only log transformation approaches can unmask a size related differentiation in growth form of the leaf. Generally, whenever structure of regression error is undetermined, choosing a suitable form of retransformation correction factor becomes elusive. Compared to customary nonparametric characterizations of this correction factor, present form proved more efficient. We expect that offered generalized allometric model along with proposed correction factor form provides a suitable analytical arrangement for the general settings of allometric examination.


---

## Body

## 1. Introduction

The model of relative growth of Huxley [1] is formally stated by means of a scaling relationship of the form(1)w=βaα,where w and a are measurable traits and the parameter α is designated as the allometric exponent, while β is identified as the normalization constant. This model, also termed equation of simple allometry, has been extensively used in research problems in biology [1–5], physics [6], economics [7], earth sciences [8], resource management, and conservation [9, 10], among other fields.Eelgrass provides nursery for waterfowl and fish species. By trapping sediment and stumping wave energy, this seagrass promotes shoreline stabilization. Eelgrass services also include nutrient recycling, water filtration, and carbon dioxide removal. Current anthropogenic influences threaten eelgrass permanence. Conservation efforts rely on plot transplanting in a fundamental way. Monitoring effectiveness depends on measurements of standing stock and productivity through time. This makes the assessment of average leaf biomass in shoots a necessary input. But traditional estimation of eelgrass leaf biomass units relies on destructive methods. This could alter shoot density in a developing transplant. Thus evaluation renders indirect assessment methods necessary [10, 11]. Results show that an allometric scaling of the form of (1) for eelgrass leaf biomass w and associated area a is consistent [11]. Derived projections of individual leaf biomass convey useful surrogates for mean leaf biomass in shoots. Moreover, estimates of the parameters α and β are invariant within a given geographical region [10–12]. Hence, estimates fitted at site can endow suitable projections of leaf biomass values currently observed in other places of the region. This bears the referred allometric projections of a convenient nondestructive feature ([11, 13]).Simplicity of (1) makes allometric projection of eelgrass leaf biomass units convenient. But there are caveats on dependability. For instance, even though α and β are invariant, environmental influences can induce a relative extent of variability on local estimates ([12, 14]). Besides, the response in the power function-like scaling of (1) is very sensitive to variation of parameter estimates. Then, accuracy of derived proxies is subject of error propagation of estimates. In addition, there are factors of biological scaling which can influence precision of estimates (e.g., [10, 11, 13, 15]). Packard [16] questioned results of Mascaro et al. [17] on allometric examination, and Mascaro et al. [18] responded to criticism. Going over this exchange highlights the relevance of procedural factors in determining precision of parameter estimates of allometric scaling. It also offers a convenient framework for the aims of the present research.An important factor influencing precision of estimates of allometric parameters is analysis method. A widespread approach is the traditional analysis method of allometry (TAMA hereafter). It relies on a log transformation of data in arithmetical scale in order to contemplate a linear regression model in geometrical scale. Then, the fitted line is back-transformed to yield a two-parameter power function in the original scale. But embracing this procedure fuels a vivid unresolved debate. Views assert that this protocol can lead to biased results (e.g., [16, 19–29]). And other practitioners consistently wed to the idea that logarithmic transformations are necessary (e.g., [18, 30–39]). An alternative to the TAMA approach is using nonlinear regression methods in the direct scale of the data [16]. Echavarria-Heras et al. [11] concluded that producing allometric projections of average leaf biomass in eelgrass shoots must rely on this protocol. Yet a direct nonlinear regression approach in allometry is also not unfailing. For instance, inadequate identification of the inherent error structure can lead to significant bias [18]. Besides, Lai et al. [30] found that estimates of allometric parameters fitted by nonlinear regression can exhibit a high sensitivity to the largest values of the covariate. Therefore, evaluation of analysis method suitability in acquiring consistent eelgrass leaf biomass proxies needs revision.The adoption of methods of curvature in geometrical space could offer a way to overcome inadequacy of the TAMA procedure ([27, 40–42]). In particular, it is pertinent to examine if taking curvature into account leads to improved accuracy of eelgrass leaf biomass proxies. But, according to Mascaro et al [18], curvature could manifest because of methodological factors of data gathering. Thus, an examination of the effect of curvature in eelgrass leaf biomass allometry must also take into account a possible participation of data quality effects. Mascaro et al. [18] reminds on three ways of handling curvilinearity in geometrical space. One is by separating data to contemplate different local linear models to account for heterogeneity of effects of the covariate [43–45]. A second one is by fitting a polynomial model [46–49]. A third approach endorses direct nonlinear regression assuming a heteroscedastic error structure as contemplated by Mascaro et al. [17]. Either approach above bears complexity beyond the linear model in geometrical space that associates with the customary bivariate power function of allometry. This suggested putting forward a generalized allometric model intended to deal with curvature in geometrical space. This paradigm incorporates parameters that change as continuous functions of the log transformed covariate ([27, 39, 46]). As we explain further on, the curvature arrangements recommended by Mascaro et al. [18] can be all derived from the offered formalization. Moreover, a nonzero intercept power function that Packard [50] recommends to handle curvilinearity in geometrical space also derives from the presented generalized scaling model.But any scheme addressing curvature in geometrical space depends on a factor for correction of bias of retransformation of the regression error. In the general settings, ifϵ stands for the regression error, then the said correction factor, through denoted using the symbol δϵ, is given by the mean of the exponentiated error random variable; that is, δϵ=Ee∈ [51–53]. Furthermore, the TAMA approach relies on the essential assumption that ϵ is additive, normally distributed, and homoscedastic [33]. When this happens, δϵ takes its lognormal-mean form [51–53]. But if ϵ fails to be normally distributed, there are two possibilities. If the distribution of ϵ is known, we could derive a closed form for δϵ. In turn, if the error distribution is not identified a priori, a widespread approach is taking δϵ in the nonparametric form given by the smearing estimator of bias of Duan [54]. Still, there are provisions on this. A smearing estimate form can fail to compensate the downward retransformation bias of logged data ([53, 55, 56]). Thus, in a circumstance where ϵ is unspecified, characterizing δϵ seems elusive. Here, we put forward an arrangement for δϵ aimed to get around this circumstance. Zeng and Tang [57] proposed a nonparametric alternate to the smearing form. It matches the first three terms’ partial sum of a power series expression of Ee∈, assuming Eϵ=0. Form suggested here corresponds to a generalization of this construct. It does not abide the restriction Eϵ=0 and matches an n-terms partial sum approximation of the exponential series representation of Ee∈. The partial sum that maximizes reproducibility strength of retransformed mean response sets criterion to choose n.Present results show that a consideration of curvature in geometrical space, as well as a suitable characterization of the correction factor of retransformation bias, offers consistent allometric proxies of observed mean leaf biomass in eelgrass shoots. Hence, contrary to views asserting direct nonlinear regression as mandatory in allometric examination, our findings validate a parallel reliability of log transformation based methods. This is well in line with claims of Mascaro et al. [18] and many others about blaming the use of logarithms of incongruent results in allometric analysis. Moreover, keeping the analysis in geometrical space unraveled heterogeneity in the inherent leaf biomass scaling pattern. This could not be achieved by clinging to direct nonlinear regression in arithmetic space as the only valid approach of allometric examination. Offered analytical arrangement is expected to be applicable in the general settings of allometry.


## 2. Materials and Methods

### 2.1. Data

For the aims of the present research, we relied on an extensive eelgrass data set collected in San Quintin Bay, a coastal lagoon on the Pacific side of the Baja California Peninsula, México (30°30’ N – 116°00′W), and through a 13 months’ long sampling period covering a whole-year cycle. Data composes measurements of length (mm), width (mm), and dry weightw (g) of a total of 10412 individual eelgrass leaves taken from 20 randomly thrown 400cm2 quadrats every monthly visit to the site. A sampling visit will be further referred to as “sampling time” in the text. The length times width proxy [11] provided estimations of leaf area a (mm2). In order to test for methodological influences of data gathering, we processed raw data set according to a mean plus or minus two standard deviations outlier’s removal procedure [58]. Appendix A presents results of an exploratory analysis of data.


### 2.2. Models

As above specified symbols,w and a stand for the biomass of an individual eelgrass leaf and its respective area one to one. Echavarría-Heras et al. [11] assert that these variables can be related through the bivariate allometric model of (1). One procedure to acquire estimates for the parameters α and β is fitting directly in arithmetical scale a nonlinear homoscedastic regression model. Besides, we can use a TAMA approach, that is, fitting the linear regression model(2)v=ln⁡β+αu+ϵ,where v=ln⁡w, u=ln⁡a, and ϵ is a random error term assumed to be normally distributed with zero mean and variance σ2, that is, ϵ~N0,σ2.We conceive curvature in geometrical space as a circumstance where fitting results of regression model of (2) are inconsistent. Dealing with this situation amounts for considering complexity beyond incorporated by (1). One possible approach to address curvature is assuming that scaling parameters α and β in (2) depend continuously on the covariate ([27, 39, 46]). This is consistent with the generalized allometric model,(3)w=βaaαa,with βa and α(a) intended to be continuous and differentiable functions defined on R+ and with βa being positive. Certainly a log transformation v=ln⁡w, u=ln⁡a of (3) establishes the regression model(4)v=vCa,u+ϵ,where(5)vCa,u=ln⁡βau+αauu,where ϵ, a residual error term, is conceived as a random variable that in the general settings is ψ-distributed with mean μ and variance set by a function σ2u of the covariate u, that is, ϵ~ψμ,σ2u.Settingβa(u)=β and αau=α with α and β constants reduces (4) to the regression model of (2). In Appendix B, we explain that (3) accommodates all curvature paradigms suggested by Mascaro et al. [18]. These include a biphasic and a polynomial model in geometrical space, as well as the nonlinear heteroscedastic model referred to by Mascaro et al. [17] in direct arithmetical space. Moreover, as shown in Appendix B, the three-parameter power function chosen as an alternate standard for curvature [59] also derives from (3).


#### 2.2.1. Biphasic Model in Geometrical Space

In order to characterize the model of (3) in a biphasic mode, we let βa=βBa and αa=αBa, such that(6)αBa=∑12ϑiuaαiβBa=exp⁡∑12ϑiualn⁡βiincluding parameters βi and αi and the function ϑi(ua) given by(7)ϑiua=2-i1-Hu-ubfori=1i-1Hu-ubfori=2,for i=1,2. Hu-ub is a Heaviside function Hz [60], evaluated at z=u-ub and correspondingly ub=ln⁡ab, with amin≤ab≤amax being a point separating growth phases u≤ub and u>ub. Then, denoting by means of vBa,u the resulting form of vCa,u from (4), we get the biphasic regression model(8)v=vBa,u+ϵ,where ϵ is a random error term as defined in (4) and(9)vBa,u=∑12ϑiuafiuawith(10)fiua=ln⁡βi+αiua,where βi and αi for i=1,2 parameters are to be estimated from data.


#### 2.2.2. Polynomial Model in Geometrical Space

Similarly, assume thatαa=αP(a) and βa=βP(a), with(11)αPa=∑0nαkuakβPa=e∑0nβkuak.αk and βk for 1≤k≤m are coefficients; one can acquire a polynomial representation vPa,u for the generalized mean response function in geometrical space vCa,u. This way the polynomial form of regression (4) becomes(12)v=vPa,u+ϵ,where ϵ is a random error term as defined in (4), with(13)vPa,u=∑0mpkuak,and pk, for k=0,2,…,m parameters.


#### 2.2.3. Nonlinear Heteroscedastic Model in Arithmetical Space

As we explain ahead, direct algebraic manipulation of (3) leads to the consideration of the nonlinear heteroscedastic regression model addressed by Mascaro et al. [17]; namely,(14)w=βθaαθ+aθϵ,with αθ, βθ, and θ being parameters and ϵ being a zero mean normally distributed error term with a covariate dependent variance σ2a=σ2a2θ, that is, ϵ~N0,σ2a2θ.A nonlinear homoscedastic form derives from (14) by setting θ=0; that is,(15)w=βoaαo+ϵ,And, again ϵ is an additive error term assumed to be normally distributed with zero mean and variance σ2, that is, ϵ~N0,σ2.AppendixA deals with exploratory analysis of data Appendix B presents notation convention and also explains how all addressed paradigms derive from the generalized model of (3). Appendix C explains the addressed forms of correction factor for bias of retransformation of the regression error. Fitting results of the geometrical space based models appear in Appendix D. Those corresponding to the nonlinear heteroscedastic and homoscedastic models pertain to Appendix E. Agreement between observed and projected values is commonly evaluated by analyzing values of Lin’s Concordance Correlation Coefficient (CCC) [61]. This correspondence index is commonly denoted by means of the symbol ρ. Agreement will be defined as poor whenever ρ<0.90, moderate for 0.90≤ρ<0.95, good for 0.95≤ρ<0.99, or excellent for ρ≥0.99 [62]. Besides CCC values, we assessed reproducibility by comparing goodness-of-fit statistics, such as the coefficient of determination (CD), standard error of estimate (SEE), mean prediction error (MPE), total relative error (TRE), average systematic error (ASE), and mean percent standard error (MPSE) ([63–65]). For statistical tasks, we relied on the R package release 3.5.


## 2.1. Data

For the aims of the present research, we relied on an extensive eelgrass data set collected in San Quintin Bay, a coastal lagoon on the Pacific side of the Baja California Peninsula, México (30°30’ N – 116°00′W), and through a 13 months’ long sampling period covering a whole-year cycle. Data composes measurements of length (mm), width (mm), and dry weightw (g) of a total of 10412 individual eelgrass leaves taken from 20 randomly thrown 400cm2 quadrats every monthly visit to the site. A sampling visit will be further referred to as “sampling time” in the text. The length times width proxy [11] provided estimations of leaf area a (mm2). In order to test for methodological influences of data gathering, we processed raw data set according to a mean plus or minus two standard deviations outlier’s removal procedure [58]. Appendix A presents results of an exploratory analysis of data.


## 2.2. Models

As above specified symbols,w and a stand for the biomass of an individual eelgrass leaf and its respective area one to one. Echavarría-Heras et al. [11] assert that these variables can be related through the bivariate allometric model of (1). One procedure to acquire estimates for the parameters α and β is fitting directly in arithmetical scale a nonlinear homoscedastic regression model. Besides, we can use a TAMA approach, that is, fitting the linear regression model(2)v=ln⁡β+αu+ϵ,where v=ln⁡w, u=ln⁡a, and ϵ is a random error term assumed to be normally distributed with zero mean and variance σ2, that is, ϵ~N0,σ2.We conceive curvature in geometrical space as a circumstance where fitting results of regression model of (2) are inconsistent. Dealing with this situation amounts for considering complexity beyond incorporated by (1). One possible approach to address curvature is assuming that scaling parameters α and β in (2) depend continuously on the covariate ([27, 39, 46]). This is consistent with the generalized allometric model,(3)w=βaaαa,with βa and α(a) intended to be continuous and differentiable functions defined on R+ and with βa being positive. Certainly a log transformation v=ln⁡w, u=ln⁡a of (3) establishes the regression model(4)v=vCa,u+ϵ,where(5)vCa,u=ln⁡βau+αauu,where ϵ, a residual error term, is conceived as a random variable that in the general settings is ψ-distributed with mean μ and variance set by a function σ2u of the covariate u, that is, ϵ~ψμ,σ2u.Settingβa(u)=β and αau=α with α and β constants reduces (4) to the regression model of (2). In Appendix B, we explain that (3) accommodates all curvature paradigms suggested by Mascaro et al. [18]. These include a biphasic and a polynomial model in geometrical space, as well as the nonlinear heteroscedastic model referred to by Mascaro et al. [17] in direct arithmetical space. Moreover, as shown in Appendix B, the three-parameter power function chosen as an alternate standard for curvature [59] also derives from (3).


### 2.2.1. Biphasic Model in Geometrical Space

In order to characterize the model of (3) in a biphasic mode, we let βa=βBa and αa=αBa, such that(6)αBa=∑12ϑiuaαiβBa=exp⁡∑12ϑiualn⁡βiincluding parameters βi and αi and the function ϑi(ua) given by(7)ϑiua=2-i1-Hu-ubfori=1i-1Hu-ubfori=2,for i=1,2. Hu-ub is a Heaviside function Hz [60], evaluated at z=u-ub and correspondingly ub=ln⁡ab, with amin≤ab≤amax being a point separating growth phases u≤ub and u>ub. Then, denoting by means of vBa,u the resulting form of vCa,u from (4), we get the biphasic regression model(8)v=vBa,u+ϵ,where ϵ is a random error term as defined in (4) and(9)vBa,u=∑12ϑiuafiuawith(10)fiua=ln⁡βi+αiua,where βi and αi for i=1,2 parameters are to be estimated from data.


### 2.2.2. Polynomial Model in Geometrical Space

Similarly, assume thatαa=αP(a) and βa=βP(a), with(11)αPa=∑0nαkuakβPa=e∑0nβkuak.αk and βk for 1≤k≤m are coefficients; one can acquire a polynomial representation vPa,u for the generalized mean response function in geometrical space vCa,u. This way the polynomial form of regression (4) becomes(12)v=vPa,u+ϵ,where ϵ is a random error term as defined in (4), with(13)vPa,u=∑0mpkuak,and pk, for k=0,2,…,m parameters.


### 2.2.3. Nonlinear Heteroscedastic Model in Arithmetical Space

As we explain ahead, direct algebraic manipulation of (3) leads to the consideration of the nonlinear heteroscedastic regression model addressed by Mascaro et al. [17]; namely,(14)w=βθaαθ+aθϵ,with αθ, βθ, and θ being parameters and ϵ being a zero mean normally distributed error term with a covariate dependent variance σ2a=σ2a2θ, that is, ϵ~N0,σ2a2θ.A nonlinear homoscedastic form derives from (14) by setting θ=0; that is,(15)w=βoaαo+ϵ,And, again ϵ is an additive error term assumed to be normally distributed with zero mean and variance σ2, that is, ϵ~N0,σ2.AppendixA deals with exploratory analysis of data Appendix B presents notation convention and also explains how all addressed paradigms derive from the generalized model of (3). Appendix C explains the addressed forms of correction factor for bias of retransformation of the regression error. Fitting results of the geometrical space based models appear in Appendix D. Those corresponding to the nonlinear heteroscedastic and homoscedastic models pertain to Appendix E. Agreement between observed and projected values is commonly evaluated by analyzing values of Lin’s Concordance Correlation Coefficient (CCC) [61]. This correspondence index is commonly denoted by means of the symbol ρ. Agreement will be defined as poor whenever ρ<0.90, moderate for 0.90≤ρ<0.95, good for 0.95≤ρ<0.99, or excellent for ρ≥0.99 [62]. Besides CCC values, we assessed reproducibility by comparing goodness-of-fit statistics, such as the coefficient of determination (CD), standard error of estimate (SEE), mean prediction error (MPE), total relative error (TRE), average systematic error (ASE), and mean percent standard error (MPSE) ([63–65]). For statistical tasks, we relied on the R package release 3.5.


## 2.2.1. Biphasic Model in Geometrical Space

In order to characterize the model of (3) in a biphasic mode, we let βa=βBa and αa=αBa, such that(6)αBa=∑12ϑiuaαiβBa=exp⁡∑12ϑiualn⁡βiincluding parameters βi and αi and the function ϑi(ua) given by(7)ϑiua=2-i1-Hu-ubfori=1i-1Hu-ubfori=2,for i=1,2. Hu-ub is a Heaviside function Hz [60], evaluated at z=u-ub and correspondingly ub=ln⁡ab, with amin≤ab≤amax being a point separating growth phases u≤ub and u>ub. Then, denoting by means of vBa,u the resulting form of vCa,u from (4), we get the biphasic regression model(8)v=vBa,u+ϵ,where ϵ is a random error term as defined in (4) and(9)vBa,u=∑12ϑiuafiuawith(10)fiua=ln⁡βi+αiua,where βi and αi for i=1,2 parameters are to be estimated from data.


## 2.2.2. Polynomial Model in Geometrical Space

Similarly, assume thatαa=αP(a) and βa=βP(a), with(11)αPa=∑0nαkuakβPa=e∑0nβkuak.αk and βk for 1≤k≤m are coefficients; one can acquire a polynomial representation vPa,u for the generalized mean response function in geometrical space vCa,u. This way the polynomial form of regression (4) becomes(12)v=vPa,u+ϵ,where ϵ is a random error term as defined in (4), with(13)vPa,u=∑0mpkuak,and pk, for k=0,2,…,m parameters.


## 2.2.3. Nonlinear Heteroscedastic Model in Arithmetical Space

As we explain ahead, direct algebraic manipulation of (3) leads to the consideration of the nonlinear heteroscedastic regression model addressed by Mascaro et al. [17]; namely,(14)w=βθaαθ+aθϵ,with αθ, βθ, and θ being parameters and ϵ being a zero mean normally distributed error term with a covariate dependent variance σ2a=σ2a2θ, that is, ϵ~N0,σ2a2θ.A nonlinear homoscedastic form derives from (14) by setting θ=0; that is,(15)w=βoaαo+ϵ,And, again ϵ is an additive error term assumed to be normally distributed with zero mean and variance σ2, that is, ϵ~N0,σ2.AppendixA deals with exploratory analysis of data Appendix B presents notation convention and also explains how all addressed paradigms derive from the generalized model of (3). Appendix C explains the addressed forms of correction factor for bias of retransformation of the regression error. Fitting results of the geometrical space based models appear in Appendix D. Those corresponding to the nonlinear heteroscedastic and homoscedastic models pertain to Appendix E. Agreement between observed and projected values is commonly evaluated by analyzing values of Lin’s Concordance Correlation Coefficient (CCC) [61]. This correspondence index is commonly denoted by means of the symbol ρ. Agreement will be defined as poor whenever ρ<0.90, moderate for 0.90≤ρ<0.95, good for 0.95≤ρ<0.99, or excellent for ρ≥0.99 [62]. Besides CCC values, we assessed reproducibility by comparing goodness-of-fit statistics, such as the coefficient of determination (CD), standard error of estimate (SEE), mean prediction error (MPE), total relative error (TRE), average systematic error (ASE), and mean percent standard error (MPSE) ([63–65]). For statistical tasks, we relied on the R package release 3.5.


## 3. Results

Exploratory analysis in AppendixA identifies maximum, minimum, and sample mean values for observed leaf area values a and associated dry weights w. We also explain distribution of variables in terms of quantiles of probability 0.1, 0.25, 0.50, 0.75, and 0.90, for both crude and processed data. Statistical exploration extends to log transformed values of these variables. We present Q-Q plots (quantile-quantile) for comparison of distribution patterns, as well as boxplots for the 13 months’ long sampling scheme for both crude and processed data. We can learn that, from month 2 to month 6, a reduction in the values of weight and area occurred; this perhaps is explained by an increase in temperature during the period. We can be also aware that a similar variation pattern over time is shown in both raw and processed data sets.


### 3.1. Fitting Results of Geometrical Space Models

In order to validate curvature in geometrical space, we compared the linear model derived from (1) as well as biphasic and polynomial alternates derived from the generalized model of (3). Appendix B explains formal matters. Tables 1 and 2 summarize notation convention. Equations numbered beyond (15) belong to the appendices. Appendix D explains corresponding regression protocols.Table 1
Summary of notation convention of the bivariate scaling model of (1) and its generalization to account for curvature as given by (3). GS stands for geometrical space, AS stands for arithmetical space, and CF means correction factor of retransformation bias.
Typical bivariate Eq. Generalized for curvature Eq. Basic form w = β a α (1) w = β a a α ( a ) (3)  Regression equation in GS v = ln ⁡ β + α u + ϵ (2), (D.1) v = v C ( a , u ) + ϵ (4)  Mean response GS E T v ∣ u = ln ⁡ β + ln ⁡ α (B.4) E C v ∣ u = ln ⁡ β a + α a u a (B.1)  Back transformation AS w =  e E T v ∣ u e ϵ (B.5) w = e E C v ∣ u e ϵ (B.2)  Mean response AS E T δ w ∣ a = β a α δ ϵ δϵ=Eeϵ (B.6) E C δ w ∣ a = β a a α ( a ) δ ( ϵ ) δϵ=Eeϵ (B.3)  CF Baskerville E T B w ∣ a = β a α δ B ( ϵ ) δB(ϵ)=eσ2/2 (B.6), (C.1) E C B w ∣ a = β a a α ( a ) δ B ( ϵ ) (B.3), (C.1)  CF Duan E T D w ∣ a =  β a α δ D ϵ δDϵ=∑1meϵjm (B.6), (C.2) E C w ∣ a = β a a α ( a ) δ D ϵ (B.3), (C.2)  CF Zeng and Tang E T Z T w ∣ a =  β a α δ Z T ϵ δZTϵ=1+σ22 (B.6), (C.3) E C Z T w ∣ a = β a a α ( a ) δ Z T ϵ (B.3), (C.3)  CF n-partial sum E T n w ∣ a = β a α δ n ϵ δnϵ=∑0nE(ϵk)k! (B.6), (C.4) E C n w ∣ a = β a a α ( a ) δ n ϵ (B.3), (C.4)Table 2
Notation convention for curvature models in geometrical space. We include the biphasic and polynomial characterizations of the generalized allometric model of (3).GS stands for geometrical space, AS stands for arithmetical space, and CF means correction factor of retransformation bias of the regression error.
Biphasic model Eq. Polynomial model Eq. Basic form w = β B a a α B a (B.11) w = β P a a α P a (B.27)  Regression equation v = v B a , u + ϵ (8), (D.4) v = v P a , u + ϵ (B.32), (D.8)  Mean response GS E B v ∣ u =  ∑ 1 2 ϑ i u a f i u a (B.18) E P v ∣ u =  ∑ 0 m p k u ( a ) k (B.33)  Back transformation AS w = β B a a α B a e ϵ (B.21) w=βPaaαPaeϵ (B.34)  Mean response AS E B δ w ∣ a = β B a a α B a δ ( ϵ ) δϵ=Eeϵ (B.22) E P w ∣ a = β P a a α P a δ ( ϵ ) δϵ=Eeϵ (B.35)   CF Baskerville E B B w ∣ a = β B a a α B a δ B ( ϵ ) (B.22), (C.1) EPBw∣a=βPaaαPaδB(ϵ) (B.35), (C.1)  CF Duan E B D w ∣ a = β B a a α B a δ D ϵ (B.22), (C.2) EPDw∣a=βPaaαPaδDϵ (B.35), (C.2)  CF Zeng and Tang E B Z T w ∣ a = β B a a α B a δ Z T ϵ (B.22), (C.3) EPZTw∣a=βPaaαPaδZTϵ (B.35), (C.3)  CF n-partial sum E B n w ∣ a = β B a a α B a δ n ϵ (B.22), (C.4) E P n w ∣ a = β P a a α P a δ n ϵ (B.35), (C.4)Fitting results of the TAMA arrangement of (2) appear in Appendix D. Figure 1 shows the spread about TAMA’s linear mean response function ETv∣u. We can visually ascertain that deviations from the linear mean response function ETv∣u suggest curvature (red dots). Thus, data processing removed inconsistent replicates but shown spread still deviates from a linear mean response. This suggests that curvature in geometrical space could not be explained by methodological factors related to data gathering.Figure 1
Spread about the TAMA’s linear mean response functionET(v∣u). Deviations about ET(v∣u) shown by red dots suggest curvature.Fitting results of the biphasic protocol of (8) are summarized in Appendix D. Figure 2 displays the spread about mean response function EBv∣u in geometrical space. Compared with Figure 1, we can ascertain that the biphasic fit provides a consistent account of different variation patterns among smaller and larger leaves. We can visually ascertain that fit produced consistent results. This confirms a judgement that identified curvature might be due to intrinsic factors of leaf growth rather than methodological influences related to data gathering.Figure 2
Spread about the biphasic mean response functionEBv∣u. Compared with the plot in Figure 1, we observe that the biphasic choice offers a better account of variability of the log transformed response than the TAMA alternate.AppendixD presents fitting results of the polynomial model (m=6). Figure 3 displays dispersion about the polynomial mean response function in geometrical space EPv∣u. A polynomial representation also exhibits higher consistency than the TAMA arrangement. Recalling the biphasic scheme, the polynomial suggests a smooth transition between two growing phases.Figure 3
Spread about the polynomial mean response functionEPv∣u. Compared with the plot in Figure 1, we can be aware that, as opposite to the TAMA protocol, the polynomial scheme offers a consistent account of curvature.


### 3.2. Model Selection in Geometrical Space

Assessment of models fitted on geometrical space relied on goodness-of-fit statistics, that is, the coefficient of determination, standard error of estimate, mean prediction error, total relative error, average systematic error, and mean percent standard error ([63–65]). Besides, we took into account concordance correlation coefficient [61] and Akaike’s information index [66]. Table 3 presents results. Goodness-of-fit statistics and ρ and AIC values disfavored the TAMA protocol. On the contrary, comparison indices favored the biphasic model. Moreover, differences among indices but TRE and ASE for this scheme and the polynomial (m=6) are slight. Particularly, the highest AIC is associated with the TAMA protocol (AIC=15069.9,∆AIC=2491.9). Therefore, this model bears the less support. The biphasic choice delivered the smallest AIC’s value (AIC=12578.0,∆AIC=0). Nevertheless, difference in AIC is just barely relative to the m=6 polynomial model, since this choice conveyed (AIC=12615.0 and ∆AIC=37). Thus, model confrontation shows that the TAMA protocol is unsuited, thus backing the assertion that whatever model aims to be consistent with the present data, it ought to be nonlinear in geometrical space.Table 3
Assessment of geometrical space fitted models. Comparison took into account goodness-of-fit statistics, that is, the coefficient of determination (R2) standard error of estimate (SEE), mean prediction error (MPE), total relative error (TRE), average systematic error (ASE), and mean percent standard error (MPSE) ([63–65]). Besides, we considered concordance correlation coefficient (ρ) [61] and Akaike’s information index (AIC) [66].
Method AIC ρ R2 SEE TRE ASE MPE MPSE TAMA 15069.9 0.9505 0.9045 0.5135 -0.0286 -0.2462 -0.1879 5.7877 Biphasic 12578.0 0.9614 0.9256 0.4530 0.0037 0.0286 -0.1658 5.1446 Polynomial 12615.0 0.9614 0.9241 0.4576 0.8497 1.2008 -0.1675 5.3187


### 3.3. Retransformation Results

The TAMA protocol was not supported by the model selection criteria. Anyway, for comparison, corresponding retransformation results are included in AppendixD. Related with the TAMA protocol, fitting results of the biphasic model display a relatively improved distribution of residuals about the zero line. Nevertheless, normal Q-Q plot still shows heavier tails than those expected for a normal distribution. And, again both test statistics and p values of an Anderson-Darling test [67] provide evidence against normality of residuals. This justifies choosing the nonparametric forms δDϵ,δZTϵ, or δnϵ for compensation of downward bias induced by retransformation of the regression error ([11, 52, 53]). Table 4 displays comparison statistics for the reproducibility strength of the biphasic mean response EBδ(w∣a) as shaped by the different forms of δϵ.Table 4
Comparison of reproducibility strength statistic for the biphasic mean responseEBδ(w∣a) as calculated by different forms of the correction factor for bias of retransformation of the regression error ϵ.The ρδϵ symbol stands for the concordance correlation value associated with δϵ.
CF ρ ( δ ϵ ) R2 SEE TRE ASE MPE MPSE δ D ( ϵ ) 0.9664 0.9256 0.0049 -9.2196 4.48e-13 0.8133 31.2839 δ Z T ( ϵ ) 0.9687 0.9319 0.0047 -7.5530 1.8359 0.7780 31.4388 δ n ( ϵ ) 0.9727 0.9464 0.0042 1.9516 12.3058 0.6903 33.9525We can learn that agreement between the biphasic mean response and leaf biomass data is best forδnϵ. Figure 4 shows spread of processed leaf biomass values about the biphasic mean response function EBδw∣a as shaped by the considered forms of δ(ϵ). We can observe that both δDϵ and δZTϵ overcompensate the bias correction by δnϵ. Moreover results show that as opposed to the TAMA a biphasic protocol along with the δnϵ form offers consistent proxies of individual leaf biomass. But it is worth mentioning that in spite of the fact that model selection favored the biphasic scheme, examination of the polynomial model output reveals similar predictive strength to the biphasic alternate.Figure 4
Spread about the biphasic mean response functionEBδw∣a in arithmetical space. Black lines relate to δDϵ, green to δZTϵ, and red ones to δnϵ.


### 3.4. Assessing Curvature by Direct Nonlinear Regression

As suggested by Mascaro et al. [18], effects of curvature in geometrical space can be analyzed by means of the direct nonlinear heteroscedastic regression model of (14). In Appendix B, we explain that such a protocol also derives from the generalized bivariate allometric model of (3). Table 5 presents pertinent notation convention. For comparison, we also present results for the associated homoscedastic case.Table 5
Notation convention for the nonlinear heteroscedastic and homoscedastic models of (14) and (15), respectively.
Nonlinear model Heteroscedastic Eq. Homoscedastic Eq. Regression equation w = β θ  a α θ +  a θ ϵ ϵ=N0,a2θσ2 (14), (E.1) w =  β o  a α o +  ϵ ϵ~N0,σ2 (15), (E.5)  Mean response function E θ w ∣ a = β θ a α θ (B.45) E o w ∣ a =  β o  a α o (B.46)Fitting results of the heteroscedastic and homoscedastic models appear in AppendixE. We can learn that estimates for the normalization constant and scaling exponent parameters are very similar. Certainly, corresponding 95% confidence intervals display some overlap. As a result, we can expect similar reproducibility features for both models. Table 6 presents comparison statistics.Table 6
Assessment of models fitted in arithmetical space. This includes the nonlinear homoscedastic and heteroscedastic protocols.
Method AIC ρ R2 SEE TRE ASE MPE MPSE Heteroscedastic -92761.13 0.972 0.9467 0.0041 0.0311 26.1073 0.6882 50.2781 Homoscedastic -81386.51 0.972 0.9467 0.0041 0.2954 28.4527 0.6881 51.7400We can be aware that model assessment backs the heteroscedastic model. But selection here is mainly on qualitative grounds. It actually concerns the ability of the heteroscedastic model to identify an expected dependence of variance in the covariate. Certainly, the reproducibility strengths of both paradigms are equivalent. Indeed, Figure5 shows that mean response curves Eθw∣a and Eow∣a differ just barely.Figure 5
Spread about the mean response, curvesEθw∣a and Eow∣a associated with the nonlinear heteroscedastic and homoscedastic regression models of (14) and (15) one to one. The mean response Eθw∣a is shown in black lines and those corresponding to Eow∣a in red.Results show that as it occurred for models fitted in geometrical space, data cleaning failed to correct a heavy tails problem for the nonlinear fits. This can be ascertained from the normal Q-Q plot of residuals. This strengthens our point on the consideration of a different error structure from the one assumed here. Exploring the effects of error structure in the fitting of models for curvature addressed here will be a matter of further research. Interestingly, both the homoscedastic and heteroscedastic models seem to induce the same reproducibility strengths.


### 3.5. Model Assessment in Arithmetical Space

The model selection assay in geometrical space summarized in Table3 favored the biphasic protocol. Correspondingly, statistics in Table 6 support the nonlinear heteroscedastic model. Results of Table 4 endure δnϵ as required for largest reproducibility of retransformation output. Table 7 allows assessment of these models. We can learn that half the number of comparison indices coincide (ρ, R2, SEE, and MPE). In addition, the biphasic model is favored by AIC, ASE, and MPSE. This sets criterion for selection of curvature in geometric space as a consistent paradigm for the present data. Accordingly, the biphasic model bears adequate.Table 7
Assessment of models in arithmetical space. This includes the nonlinear heteroscedastic and biphasic protocols.
Method AIC ρ R2 SEE TRE ASE MPE MPSE Heteroscedastic -92761.13 0.972 0.946 0.004 0.031 26.107 0.690 50.278 Biphasic -96879.30 0.972 0.946 0.004 1.951 12.305 0.690 33.952


### 3.6. Implications for Allometric Proxies of Mean Leaf Biomass in Eelgrass Shoots

We in turn consider allometric proxies for average leaf biomass in eelgrass shoots. In getting these surrogates, we aggregate allometric projections of individual leaf biomass conforming a shoot. For comparison, we consider individual leaf biomass surrogates produced by the different projection methods. Table8 compares resulting reproducibility strengths.Table 8
Comparison of reproducibility strengths of proxies of average leaf biomass in shoots resulting from the biphasic polynomial or nonlinear heteroscedastic protocols.
Method ρ R2 SEE TRE ASE MPE MPSE Biphasic 0.997 0.994 8.729e-04 2.408 2.229 3.486 5.080 Polynomial (m=6) 0.996 0.992 0.001 -2.456 -1.033 4.807 5.366 Heteroscedastic 0.997 0.994 7.722e-04 0.983 -0.726 3.084 5.645Results in [11] stablished that proxies derived from the TAMA protocol are inconsistent with observed values. This endorsed nonlinear regression in the direct scale as a requirement for reliability of allometric projections of mean leaf biomass in eelgrass shoots. But Table 8 shows that a curvature model fitted in geometrical space can offer proxies entailing similar predictive power to a nonlinear regression protocol. Plots in Figure 6 allow getting a glimpse of this assertion.Figure 6
Average leaf biomass in shoots calculated from observed-processed data compared with their allometric projections. Projection lines resulting from the biphasic and polynomial models relied on theδnϵ form.
(a) (b) (c)


## 3.1. Fitting Results of Geometrical Space Models

In order to validate curvature in geometrical space, we compared the linear model derived from (1) as well as biphasic and polynomial alternates derived from the generalized model of (3). Appendix B explains formal matters. Tables 1 and 2 summarize notation convention. Equations numbered beyond (15) belong to the appendices. Appendix D explains corresponding regression protocols.Table 1
Summary of notation convention of the bivariate scaling model of (1) and its generalization to account for curvature as given by (3). GS stands for geometrical space, AS stands for arithmetical space, and CF means correction factor of retransformation bias.
Typical bivariate Eq. Generalized for curvature Eq. Basic form w = β a α (1) w = β a a α ( a ) (3)  Regression equation in GS v = ln ⁡ β + α u + ϵ (2), (D.1) v = v C ( a , u ) + ϵ (4)  Mean response GS E T v ∣ u = ln ⁡ β + ln ⁡ α (B.4) E C v ∣ u = ln ⁡ β a + α a u a (B.1)  Back transformation AS w =  e E T v ∣ u e ϵ (B.5) w = e E C v ∣ u e ϵ (B.2)  Mean response AS E T δ w ∣ a = β a α δ ϵ δϵ=Eeϵ (B.6) E C δ w ∣ a = β a a α ( a ) δ ( ϵ ) δϵ=Eeϵ (B.3)  CF Baskerville E T B w ∣ a = β a α δ B ( ϵ ) δB(ϵ)=eσ2/2 (B.6), (C.1) E C B w ∣ a = β a a α ( a ) δ B ( ϵ ) (B.3), (C.1)  CF Duan E T D w ∣ a =  β a α δ D ϵ δDϵ=∑1meϵjm (B.6), (C.2) E C w ∣ a = β a a α ( a ) δ D ϵ (B.3), (C.2)  CF Zeng and Tang E T Z T w ∣ a =  β a α δ Z T ϵ δZTϵ=1+σ22 (B.6), (C.3) E C Z T w ∣ a = β a a α ( a ) δ Z T ϵ (B.3), (C.3)  CF n-partial sum E T n w ∣ a = β a α δ n ϵ δnϵ=∑0nE(ϵk)k! (B.6), (C.4) E C n w ∣ a = β a a α ( a ) δ n ϵ (B.3), (C.4)Table 2
Notation convention for curvature models in geometrical space. We include the biphasic and polynomial characterizations of the generalized allometric model of (3).GS stands for geometrical space, AS stands for arithmetical space, and CF means correction factor of retransformation bias of the regression error.
Biphasic model Eq. Polynomial model Eq. Basic form w = β B a a α B a (B.11) w = β P a a α P a (B.27)  Regression equation v = v B a , u + ϵ (8), (D.4) v = v P a , u + ϵ (B.32), (D.8)  Mean response GS E B v ∣ u =  ∑ 1 2 ϑ i u a f i u a (B.18) E P v ∣ u =  ∑ 0 m p k u ( a ) k (B.33)  Back transformation AS w = β B a a α B a e ϵ (B.21) w=βPaaαPaeϵ (B.34)  Mean response AS E B δ w ∣ a = β B a a α B a δ ( ϵ ) δϵ=Eeϵ (B.22) E P w ∣ a = β P a a α P a δ ( ϵ ) δϵ=Eeϵ (B.35)   CF Baskerville E B B w ∣ a = β B a a α B a δ B ( ϵ ) (B.22), (C.1) EPBw∣a=βPaaαPaδB(ϵ) (B.35), (C.1)  CF Duan E B D w ∣ a = β B a a α B a δ D ϵ (B.22), (C.2) EPDw∣a=βPaaαPaδDϵ (B.35), (C.2)  CF Zeng and Tang E B Z T w ∣ a = β B a a α B a δ Z T ϵ (B.22), (C.3) EPZTw∣a=βPaaαPaδZTϵ (B.35), (C.3)  CF n-partial sum E B n w ∣ a = β B a a α B a δ n ϵ (B.22), (C.4) E P n w ∣ a = β P a a α P a δ n ϵ (B.35), (C.4)Fitting results of the TAMA arrangement of (2) appear in Appendix D. Figure 1 shows the spread about TAMA’s linear mean response function ETv∣u. We can visually ascertain that deviations from the linear mean response function ETv∣u suggest curvature (red dots). Thus, data processing removed inconsistent replicates but shown spread still deviates from a linear mean response. This suggests that curvature in geometrical space could not be explained by methodological factors related to data gathering.Figure 1
Spread about the TAMA’s linear mean response functionET(v∣u). Deviations about ET(v∣u) shown by red dots suggest curvature.Fitting results of the biphasic protocol of (8) are summarized in Appendix D. Figure 2 displays the spread about mean response function EBv∣u in geometrical space. Compared with Figure 1, we can ascertain that the biphasic fit provides a consistent account of different variation patterns among smaller and larger leaves. We can visually ascertain that fit produced consistent results. This confirms a judgement that identified curvature might be due to intrinsic factors of leaf growth rather than methodological influences related to data gathering.Figure 2
Spread about the biphasic mean response functionEBv∣u. Compared with the plot in Figure 1, we observe that the biphasic choice offers a better account of variability of the log transformed response than the TAMA alternate.AppendixD presents fitting results of the polynomial model (m=6). Figure 3 displays dispersion about the polynomial mean response function in geometrical space EPv∣u. A polynomial representation also exhibits higher consistency than the TAMA arrangement. Recalling the biphasic scheme, the polynomial suggests a smooth transition between two growing phases.Figure 3
Spread about the polynomial mean response functionEPv∣u. Compared with the plot in Figure 1, we can be aware that, as opposite to the TAMA protocol, the polynomial scheme offers a consistent account of curvature.


## 3.2. Model Selection in Geometrical Space

Assessment of models fitted on geometrical space relied on goodness-of-fit statistics, that is, the coefficient of determination, standard error of estimate, mean prediction error, total relative error, average systematic error, and mean percent standard error ([63–65]). Besides, we took into account concordance correlation coefficient [61] and Akaike’s information index [66]. Table 3 presents results. Goodness-of-fit statistics and ρ and AIC values disfavored the TAMA protocol. On the contrary, comparison indices favored the biphasic model. Moreover, differences among indices but TRE and ASE for this scheme and the polynomial (m=6) are slight. Particularly, the highest AIC is associated with the TAMA protocol (AIC=15069.9,∆AIC=2491.9). Therefore, this model bears the less support. The biphasic choice delivered the smallest AIC’s value (AIC=12578.0,∆AIC=0). Nevertheless, difference in AIC is just barely relative to the m=6 polynomial model, since this choice conveyed (AIC=12615.0 and ∆AIC=37). Thus, model confrontation shows that the TAMA protocol is unsuited, thus backing the assertion that whatever model aims to be consistent with the present data, it ought to be nonlinear in geometrical space.Table 3
Assessment of geometrical space fitted models. Comparison took into account goodness-of-fit statistics, that is, the coefficient of determination (R2) standard error of estimate (SEE), mean prediction error (MPE), total relative error (TRE), average systematic error (ASE), and mean percent standard error (MPSE) ([63–65]). Besides, we considered concordance correlation coefficient (ρ) [61] and Akaike’s information index (AIC) [66].
Method AIC ρ R2 SEE TRE ASE MPE MPSE TAMA 15069.9 0.9505 0.9045 0.5135 -0.0286 -0.2462 -0.1879 5.7877 Biphasic 12578.0 0.9614 0.9256 0.4530 0.0037 0.0286 -0.1658 5.1446 Polynomial 12615.0 0.9614 0.9241 0.4576 0.8497 1.2008 -0.1675 5.3187


## 3.3. Retransformation Results

The TAMA protocol was not supported by the model selection criteria. Anyway, for comparison, corresponding retransformation results are included in AppendixD. Related with the TAMA protocol, fitting results of the biphasic model display a relatively improved distribution of residuals about the zero line. Nevertheless, normal Q-Q plot still shows heavier tails than those expected for a normal distribution. And, again both test statistics and p values of an Anderson-Darling test [67] provide evidence against normality of residuals. This justifies choosing the nonparametric forms δDϵ,δZTϵ, or δnϵ for compensation of downward bias induced by retransformation of the regression error ([11, 52, 53]). Table 4 displays comparison statistics for the reproducibility strength of the biphasic mean response EBδ(w∣a) as shaped by the different forms of δϵ.Table 4
Comparison of reproducibility strength statistic for the biphasic mean responseEBδ(w∣a) as calculated by different forms of the correction factor for bias of retransformation of the regression error ϵ.The ρδϵ symbol stands for the concordance correlation value associated with δϵ.
CF ρ ( δ ϵ ) R2 SEE TRE ASE MPE MPSE δ D ( ϵ ) 0.9664 0.9256 0.0049 -9.2196 4.48e-13 0.8133 31.2839 δ Z T ( ϵ ) 0.9687 0.9319 0.0047 -7.5530 1.8359 0.7780 31.4388 δ n ( ϵ ) 0.9727 0.9464 0.0042 1.9516 12.3058 0.6903 33.9525We can learn that agreement between the biphasic mean response and leaf biomass data is best forδnϵ. Figure 4 shows spread of processed leaf biomass values about the biphasic mean response function EBδw∣a as shaped by the considered forms of δ(ϵ). We can observe that both δDϵ and δZTϵ overcompensate the bias correction by δnϵ. Moreover results show that as opposed to the TAMA a biphasic protocol along with the δnϵ form offers consistent proxies of individual leaf biomass. But it is worth mentioning that in spite of the fact that model selection favored the biphasic scheme, examination of the polynomial model output reveals similar predictive strength to the biphasic alternate.Figure 4
Spread about the biphasic mean response functionEBδw∣a in arithmetical space. Black lines relate to δDϵ, green to δZTϵ, and red ones to δnϵ.


## 3.4. Assessing Curvature by Direct Nonlinear Regression

As suggested by Mascaro et al. [18], effects of curvature in geometrical space can be analyzed by means of the direct nonlinear heteroscedastic regression model of (14). In Appendix B, we explain that such a protocol also derives from the generalized bivariate allometric model of (3). Table 5 presents pertinent notation convention. For comparison, we also present results for the associated homoscedastic case.Table 5
Notation convention for the nonlinear heteroscedastic and homoscedastic models of (14) and (15), respectively.
Nonlinear model Heteroscedastic Eq. Homoscedastic Eq. Regression equation w = β θ  a α θ +  a θ ϵ ϵ=N0,a2θσ2 (14), (E.1) w =  β o  a α o +  ϵ ϵ~N0,σ2 (15), (E.5)  Mean response function E θ w ∣ a = β θ a α θ (B.45) E o w ∣ a =  β o  a α o (B.46)Fitting results of the heteroscedastic and homoscedastic models appear in AppendixE. We can learn that estimates for the normalization constant and scaling exponent parameters are very similar. Certainly, corresponding 95% confidence intervals display some overlap. As a result, we can expect similar reproducibility features for both models. Table 6 presents comparison statistics.Table 6
Assessment of models fitted in arithmetical space. This includes the nonlinear homoscedastic and heteroscedastic protocols.
Method AIC ρ R2 SEE TRE ASE MPE MPSE Heteroscedastic -92761.13 0.972 0.9467 0.0041 0.0311 26.1073 0.6882 50.2781 Homoscedastic -81386.51 0.972 0.9467 0.0041 0.2954 28.4527 0.6881 51.7400We can be aware that model assessment backs the heteroscedastic model. But selection here is mainly on qualitative grounds. It actually concerns the ability of the heteroscedastic model to identify an expected dependence of variance in the covariate. Certainly, the reproducibility strengths of both paradigms are equivalent. Indeed, Figure5 shows that mean response curves Eθw∣a and Eow∣a differ just barely.Figure 5
Spread about the mean response, curvesEθw∣a and Eow∣a associated with the nonlinear heteroscedastic and homoscedastic regression models of (14) and (15) one to one. The mean response Eθw∣a is shown in black lines and those corresponding to Eow∣a in red.Results show that as it occurred for models fitted in geometrical space, data cleaning failed to correct a heavy tails problem for the nonlinear fits. This can be ascertained from the normal Q-Q plot of residuals. This strengthens our point on the consideration of a different error structure from the one assumed here. Exploring the effects of error structure in the fitting of models for curvature addressed here will be a matter of further research. Interestingly, both the homoscedastic and heteroscedastic models seem to induce the same reproducibility strengths.


## 3.5. Model Assessment in Arithmetical Space

The model selection assay in geometrical space summarized in Table3 favored the biphasic protocol. Correspondingly, statistics in Table 6 support the nonlinear heteroscedastic model. Results of Table 4 endure δnϵ as required for largest reproducibility of retransformation output. Table 7 allows assessment of these models. We can learn that half the number of comparison indices coincide (ρ, R2, SEE, and MPE). In addition, the biphasic model is favored by AIC, ASE, and MPSE. This sets criterion for selection of curvature in geometric space as a consistent paradigm for the present data. Accordingly, the biphasic model bears adequate.Table 7
Assessment of models in arithmetical space. This includes the nonlinear heteroscedastic and biphasic protocols.
Method AIC ρ R2 SEE TRE ASE MPE MPSE Heteroscedastic -92761.13 0.972 0.946 0.004 0.031 26.107 0.690 50.278 Biphasic -96879.30 0.972 0.946 0.004 1.951 12.305 0.690 33.952


## 3.6. Implications for Allometric Proxies of Mean Leaf Biomass in Eelgrass Shoots

We in turn consider allometric proxies for average leaf biomass in eelgrass shoots. In getting these surrogates, we aggregate allometric projections of individual leaf biomass conforming a shoot. For comparison, we consider individual leaf biomass surrogates produced by the different projection methods. Table8 compares resulting reproducibility strengths.Table 8
Comparison of reproducibility strengths of proxies of average leaf biomass in shoots resulting from the biphasic polynomial or nonlinear heteroscedastic protocols.
Method ρ R2 SEE TRE ASE MPE MPSE Biphasic 0.997 0.994 8.729e-04 2.408 2.229 3.486 5.080 Polynomial (m=6) 0.996 0.992 0.001 -2.456 -1.033 4.807 5.366 Heteroscedastic 0.997 0.994 7.722e-04 0.983 -0.726 3.084 5.645Results in [11] stablished that proxies derived from the TAMA protocol are inconsistent with observed values. This endorsed nonlinear regression in the direct scale as a requirement for reliability of allometric projections of mean leaf biomass in eelgrass shoots. But Table 8 shows that a curvature model fitted in geometrical space can offer proxies entailing similar predictive power to a nonlinear regression protocol. Plots in Figure 6 allow getting a glimpse of this assertion.Figure 6
Average leaf biomass in shoots calculated from observed-processed data compared with their allometric projections. Projection lines resulting from the biphasic and polynomial models relied on theδnϵ form.
(a) (b) (c)


## 4. Discussion

The customary bivariate allometric model of (1) offers nondestructive surrogates for average leaf biomass in eelgrass shoots [11]. But there are methodological factors that could influence dependability. Views assert that parameter identification based on logarithmic transformations leads to biased projections [20–29]. But other practitioners clung to this approach as meaningful and necessary in allometric examination [30–39]. This going over suggests that surpassing this controversy amounts to considering curvature in geometrical space. For that aim, we proposed the generalized model of (3). Approaches such as direct nonlinear heteroscedastic regression, as well as biphasic and polynomial protocols in geometrical space [18], became logical resultants from this construct. For present data model selection validated maintenance of the analysis in geometrical space. Nevertheless, at an empirical level, addressed protocols produced allometric projections of individual leaf biomass of correspondent precision. This was also verified for concomitant projections of average leaf biomass in shoots. But, from a qualitative standpoint, the nonlinear regression protocol mainly contributed by identifying expected dependence of the variance on covariate. Moreover. Figure 7(a) depicts manifest differences in mean response trends between the polynomial fit and the nonlinear heteroscedastic model. Nonetheless, those in Figure 7(b) corresponding to this and the biphasic models differ but only barely. Then, a nonlinear regression scheme at best shaped a reasonable approximation of the mean response function resultant from curvature methods.Figure 7
Trends of mean response functions calculated from the retransformed outputs of polynomial (a) and biphasic models (b) involving theδn(ϵ) form of δϵ. We can observe a differentiation in trends relative to that of a power function fitted in directly by means of the nonlinear heteroscedastic protocol of (14). Although relative deviations manifest for the biphasic model, differences in mean response patterns are more clearly depicted for the polynomial choice.
(a) (b)Differences in patterns of the biphasic and polynomial mean response functions relative to the nonlinear protocol exhibit that clinging to this last paradigm could impair detection of the true allometric relationship. Moreover, relying on direct nonlinear regression impairs identification of heterogeneity in the log transformed response as covariate changes. This further stresses on limitations of this device as a tool for allometric examination ([18, 30]). Oppositely, output of the selected biphasic model shown in Figure 2 suggests differentiation of growth patterns among smaller and larger leaves. Besides, the polynomial mean response in Figure 3 suggests a gradual transition between different growth phases. Thus, as opposed to direct nonlinear regression, a consideration of curvature in geometrical space could elucidate an inherent leaf growth pattern. This strengthens a judgement that the log transformation step, essential to traditional allometric examination, cannot be thrown away without losing relevant information ([33, 37]).Mascaro et al. [18] conceived curvature in geometrical space as related to methodological factors of data gathering. But present examination corroborated consistency of curvature models for processed data. This suggests that manifestation of curvature is rather explained by intrinsic factors in leaf growth. Additionally, data processing failed to amend the heavy tails problem detected on Q-Q plots. This indicates departure of residuals from an assumed error structure. As a result, numerical values of the addressed correction factor forms turned out to be different, thus conveying ambiguity in selection of suitable mean response of models fitted in geometrical space. This could entail the only advantage of nonlinear regression method over log-transformation-curvature paradigms. But, also for this analysis method, an inadequate postulation of inherent error structure can lead to significant bias [18]. It seems then reasonable considering that a suitable characterization of error structure could lead to robustness of built allometric proxies, even when they derive from crude data. Steering to an error structure different from what is assumed here is a worthwhile subject of further research.When dealing with similar data we suggest taking into account recommendations that come up from this examination. First, it is highly advisable to perform a preliminary examination of the spread around the straight line in geometrical space resulting from the model of (1). If further statistical exploration confirms that linearity and assumed error structure are consistent with data, Huxley’s bivariate allometric model could suit. Otherwise, the arrangement of curvature, error structure, and correction factor form such as proposed here could be called into account for the analysis. The use of data cleaning procedures in order to achieve a better fit is controversial [11]. Instead of performing data processing a posteriori, it is highly advisable to rely on standardized data gathering procedures. This will prevent proliferation of inconsistent replicates that could exacerbate a heavy tails problem on Q-Q plots.


## 5. Conclusion

Failure to perform both a preliminary exploration of spread of log transformed allometric data and a sound evaluation of model adequacy could impair detecting a possible manifestation of curvature. As a consequence, the output of a traditional analysis method could set biased predictions of observed values. This circumstance could result in dismissal of a log transformation step in the analysis, giving way to contemplation of direct nonlinear regression as the only protocol to acquire reliable parameter estimates [11]. Results of this examination suggest that consideration of curvature in geometrical space as set by the model of (3) could offer dependable allometric proxies of average leaf biomass in eelgrass shoots.From a general perspective, complexity as encompassed by the model of (3) can stand for curvature as conceived in allometric examination. Particularly, biphasic or polynomial protocols in geometrical space, as well as a direct nonlinear heteroscedastic regression model, derive as particular characterizations of this paradigm. Moreover all statistical models for accurate estimates of relative growth contemplated by Bervian et, al., [68] can be also accomodated by the present generalization of the model of simple allometry of Huxley. But empirical convenience on its own does not validate adoption of this paradigm as a general tool. Certainly, the Weierstrass approximation theorem [69] backs a polynomial regression model as a reasonable identification device for the generalized allometric model expressed in geometrical space. But suitability of retransformation results will sensibly depend on correction factor form. And a mean function resulting from a polynomial fitted in geometrical space will not enable characterization of functions α(a) and β(a) one to one. Furthermore, complexity of (3) could pose significant difficulties while attempting its identification through direct nonlinear regression methods. Needless to say, biological interpretation of the scaling functions α(a) and βa is also pending. A quest for efficient tools of nondestructive assessment of plant biomass units justifies addressing these examinations in a further research.



---

*Source: 3613679-2019-04-23.xml*