Parameter Estimation Research Articles

Removal of acidic dyes; acid yellow 25 and acid red 4 from wastewater by degassed activated carbon

Abstract Activated carbon was prepared at 300 °C and 600 °C, characterized by SEM, EDX and XRD, and was then used as an adsorbent for the removal of acidic dyes; acid yellow 25 and acid red 4. The activated carbon prepared at high temperature (600 °C) due to its high carbon contents and surface area was subsequently used as adsorbent for the selected dyes adsorption using batch adsorption approaches to estimate different adsorption parameters. For the estimation of kinetics and equilibrium parameters a number of kinetics and isotherm models were employed. Dyes were adsorbed on activated carbon surface at a high rate for the first 15 min, after which it began to diffuse into the micro pores and thus the process became steady. The rate constant was estimated for first and second order kinetics models. The maximum adsorption capacities recorded were 526.32 mg g−1 for acid red 4 and 555.55 mg g−1 for acid yellow 25. The enthalpy change values recorded were; 19.44 kJ mol−1 for acid yellow 25 adsorption and 16 kJ mol−1 for acid red 4 adsorption, meant that the process is endothermic. The negative values of Gibbs free energy change (−393.28, −1,515.48, −2,634.68 J mol−1) of acid red 4 and acid yellow 25 (−251.72, −1,058.06, −2,367.84 J mol−1) at tested temperatures, confirmed the feasibility and spontaneity of the adsorption processes. The adsorption of dyes on the carbon surface was diffusion-controlled process, as demonstrated by the linear graph of intraparticle diffusion model.

A novel statistical approach to COVID-19 variability using the Weibull-Inverse Nadarajah Haghighi distribution

Researchers have devoted decades to striving to create a plethora of distinctive distributions in order to meet specific objectives. The argument is that traditional distributions have typically been found to lack fit in real-world situations, which include pharmaceutical studies, the field of engineering, hydrology, environmental science, and a number of others. The Weibull-inverse Nadarajah Haghighi (WINH) distribution is developed by combining the Weibull and inverse Nadarajah Haghighi distributions. The proposed distribution's fundamental characteristics have been established and analyzed. Several plots of the distributional properties, notably probability density function (PDF) with corresponding cumulative distribution function (CDF) are displayed. The estimation of model parameter is performed via the MLE procedure. Simulation-based research is conducted to demonstrate the performance of proposed estimator’s using some measure, like the average bias, variance, and associated mean square error (MSE). Two real datasets represent the morality due to COVID 19 in France and Canada are illustrated to see the practicality of the recommended model.

The potential of mixing model of wind speed distribution in Algerian High Plateaus

The evaluation of wind energy relies primarily on the probability density function PDF, which corresponds well with the wind speed data. Single PDFs are widely used in the assessment of wind. In contrast, homogeneous or heterogeneous mixed models are rarely used, especially in Algeria, where the bimodal wind speed distribution is expected. This research aims to investigate the potential of heterogeneous PDFs Generalized Extreme Value-Weibul and Normal-Extreme Value PDFs in assessing wind energy at three meteorological stations in the high plateau against the single widespread PDFs Weibull and GEV by analyzing five years of archived wind speed data. The estimation of mixed model parameters is obtained by applying the Expectation-Maximization algorithm, and the identification of the appropriate PDF is obtained by four goodness-of-fit (GOF) criteria and compared with the widespread single distributions. The results show that the mixed model surpasses widely the single model for all the GOF criteria used at the three selected sites. The proposed mixed model fits all the wind speed distributions related to unimodal and bimodal regimes.

Landslide Analysis with Incomplete Data: A Framework for Critical Parameter Estimation

Landslides are one of the most common geohazards, posing significant risks to infrastructure, recreation, and human life. Slope stability analyses rely on detailed data, accurate materials testing, and careful model parameter selection. These factors are not always readily available, and estimations must be made, introducing uncertainty and error to the final slope stability analysis results. The most critical slope stability parameters that are often missing or incompletely constrained include slope topography, depth to water table, depth to failure plane, and material property parameters. Though estimation of these values is common practice, there is limited guidance or best practice instruction for this important step in the analysis. Guidance is provided for the estimation of: original and/or post-failure slope topography via traditional methods as well as the use of open-source digital elevation models, water table depth across variable hydrologic settings, and the iterative estimation of depth to failure plane and slope material properties. Workflows are proposed for the systematic estimation of critical parameters based primarily on slide type and scale. The efficacy of the proposed estimation techniques, uncertainty quantification, and final parameter estimation protocol for data-sparse landslide analysis is demonstrated via application at a landslide in Colorado, USA.

Bayesian Federated Learning with Hamiltonian Monte Carlo: Algorithm and Theory

This work introduces a novel and efficient Bayesian federated learning algorithm, namely, the Federated Averaging stochastic Hamiltonian Monte Carlo (FA-HMC), for parameter estimation and uncertainty quantification. We establish rigorous convergence guarantees of FA-HMC on non-iid distributed datasets, under the strong convexity and Hessian smoothness assumptions. Our analysis investigates the effects of parameter space dimension, noise on gradients and momentum, and the frequency of communication (between the central node and local nodes) on the convergence and communication costs of FA-HMC. Beyond that, we establish the tightness of our analysis by showing that the convergence rate cannot be improved even for continuous FA-HMC process. Moreover, extensive empirical studies demonstrate that FA-HMC outperforms the existing Federated Averaging-Langevin Monte Carlo (FA-LD) algorithm. Supplementary materials for this article are available online.

Extension of the characterization of non-Gaussianity in gravitational wave detectors with a statistical hypothesis test

Abstract In gravitational wave (GW) astronomy, non-Gaussian noise, such as scattered light noise disturbs stable interferometer operation, limiting the interferometer’s sensitivity, and reducing the reliability of the analyses. In scattered light noise, the non-Gaussian noise dominates the sensitivity in a low frequency range of less than a few hundred Hz, which is sensitive to GWs from compact binary coalescence. This non-Gaussian noise prevents reliable parameter estimation, since several analysis methods are optimized only for Gaussian noise. Therefore, identifying data contaminated by non-Gaussian noise is important. In this work, we extended the conventional method to evaluate non-Gaussian noise, the Rayleigh statistic, by using a statistical hypothesis test to determine a threshold for non-Gaussian noise. First, we estimated the distribution of the Rayleigh statistic against Gaussian noise, called the background distribution, and validated that our extension serves as the hypothetical test. The threshold on the Rayleigh statistic is estimated at 0.73 and 1.28 when the significance level is 0.05, and the sample size is 39. Moreover, we investigated the detection efficiency by assuming two non-Gaussian noise models. For example, for the model with strong scattered light noise, the true positive rate (TPR) was always above 0.7 when the significance level was 0.05. For the demonstration, we applied our extension with the estimated thresholds to the real data. We confirmed our extension provides the binarized outputs of the statistical tests and contributes to finding the non-Gaussian noise in the real data. The results showed that our extension can contribute to an initial detection of non-Gaussian noise and lead to further investigation of the origin of the non-Gaussian noise.

Comparison of estimation limits for quantum two-parameter estimation

Measurement estimation bounds for local quantum multiparameter estimation, which provide lower bounds on possible measurement uncertainties, have so far been formulated in two ways: by extending the classical Cramér-Rao bound (e.g., the quantum Cramér-Rao bound and the Nagaoka Cramér-Rao bound) and by incorporating the parameter estimation framework with the uncertainty principle, as in the Lu-Wang uncertainty relation. In this work, we present a general framework that allows a direct comparison between these different types of estimation limits. Specifically, we compare the attainability of the Nagaoka Cramér-Rao bound and the Lu-Wang uncertainty relation, using analytical and numerical techniques. We show that these two limits can provide different information about the physically attainable precision. We present an example where both limits provide the same attainable precision and an example where the Lu-Wang uncertainty relation is not attainable even for pure states. We further demonstrate that the unattainability in the latter case arises because the figure of merit underpinning the Lu-Wang uncertainty relation (the difference between the quantum and classical Fisher information matrices) does not necessarily agree with the conventionally used figure of merit (mean-squared error). The results offer insights into the general attainability and applicability of the Lu-Wang uncertainty relation. Furthermore, our proposed framework for comparing bounds of different types may prove useful in other settings. Published by the American Physical Society 2024

Flexible Modeling of Nonstationary Extremal Dependence using Spatially Fused LASSO and Ridge Penalties

Statistical modeling of a nonstationary spatial extremal dependence structure is challenging. Parametric max-stable processes (MSPs) are common choices for modeling spatially indexed block maxima, where an assumption of stationarity is usual to make inference feasible. However, this assumption is unrealistic for data observed over a large or complex domain. We develop a computationally efficient method for estimating extremal dependence using a globally nonstationary but locally stationary MSP construction, with the spatial domain divided into a fine grid of subregions, each with its own dependence parameters. We use LASSO ( L 1 ) or ridge ( L 2 ) penalties to obtain spatially smooth parameter estimates. We then develop a novel data-driven algorithm to merge homogeneous neighboring subregions. The algorithm facilitates model parsimony and interpretability. To make our model suitable for high-dimensional data, we exploit a pairwise likelihood to perform inference and discuss its computational and statistical efficiency. We apply our proposed method to model monthly maximum temperature data at over 1400 sites in Nepal and the surrounding Himalayan and sub-Himalayan regions; we show significant improvements in model fit compared to a stationary model. Furthermore, we demonstrate that the estimated merged partition is interpretable from a geographic perspective and leads to better model diagnostics by adequately reducing the number of parameters.

A New Multimodal Modification of the Skew Family of Distributions: Properties and Applications to Medical and Environmental Data

The skew distribution has the characteristic of appropriately modeling asymmetric unimodal data. However, in practice, there are several cases in which the data present more than one mode. In the literature, it is possible to find a large number of authors who have studied extensions based on the skew distribution to model this type of data. In this article, a new family is introduced, consisting of a multimodal modification to the family of skew distributions. Using the methodology of the weighted version of a function, we perform the product of the density function of a family of skew distributions with a polynomial of degree 4, thus obtaining a more flexible model that allows modeling data sets, whose distribution contains at most three modes. The density function, some properties, moments, skewness coefficients, and kurtosis of this new family are presented. This study focuses on the particular cases of skew-normal and Laplace distributions, although it can be applied to any other distribution. A simulation study was carried out, to study the behavior of the model parameter estimates. Illustrations with real data, referring to medicine and environmental data, show the practical performance of the proposed model in the two particular cases presented.

Copula Approximate Bayesian Computation Using Distribution Random Forests

Ongoing modern computational advancements continue to make it easier to collect increasingly large and complex datasets, which can often only be realistically analyzed using models defined by intractable likelihood functions. This Stats invited feature article introduces and provides an extensive simulation study of a new approximate Bayesian computation (ABC) framework for estimating the posterior distribution and the maximum likelihood estimate (MLE) of the parameters of models defined by intractable likelihoods, that unifies and extends previous ABC methods proposed separately. This framework, copulaABCdrf, aims to accurately estimate and describe the possibly skewed and high-dimensional posterior distribution by a novel multivariate copula-based meta-t distribution based on univariate marginal posterior distributions that can be accurately estimated by distribution random forests (drf), while performing automatic summary statistics (covariates) selection, based on robustly estimated copula dependence parameters. The copulaABCdrf framework also provides a novel multivariate mode estimator to perform MLE and posterior mode estimation and an optional step to perform model selection from a given set of models using posterior probabilities estimated by drf. The posterior distribution estimation accuracy of the ABC framework is illustrated and compared with previous standard ABC methods through several simulation studies involving low- and high-dimensional models with computable posterior distributions, which are either unimodal, skewed, or multimodal; and exponential random graph and mechanistic network models, each defined by an intractable likelihood from which it is costly to simulate large network datasets. This paper also proposes and studies a new solution to the simulation cost problem in ABC involving the posterior estimation of parameters from datasets simulated from the given model that are smaller compared to the potentially large size of the dataset being analyzed. This proposal is motivated by the fact that, for many models defined by intractable likelihoods, such as the network models when they are applied to analyze massive networks, the repeated simulation of large datasets (networks) for posterior-based parameter estimation can be too computationally costly and vastly slow down or prohibit the use of standard ABC methods. The copulaABCdrf framework and standard ABC methods are further illustrated through analyses of large real-life networks of sizes ranging between 28,000 and 65.6 million nodes (between 3 million and 1.8 billion edges), including a large multilayer network with weighted directed edges. The results of the simulation studies show that, in settings where the true posterior distribution is not highly multimodal, copulaABCdrf usually produced similar point estimates from the posterior distribution for low-dimensional parametric models as previous ABC methods, but the copula-based method can produce more accurate estimates from the posterior distribution for high-dimensional models, and, in both dimensionality cases, usually produced more accurate estimates of univariate marginal posterior distributions of parameters. Also, posterior estimation accuracy was usually improved when pre-selecting the important summary statistics using drf compared to ABC employing no pre-selection of the subset of important summaries. For all ABC methods studied, accurate estimation of a highly multimodal posterior distribution was challenging. In light of the results of all the simulation studies, this article concludes by discussing how the copulaABCdrf framework can be improved for future research.

Driving low-carbon mechanisms through smart investments in renewable resources and green finance initiatives among G20 nations

In light of the escalating concerns regarding climate change and environmental decline, major nations are actively exploring strategies to mitigate environmental harm and achieve future sustainability. The surge in economic expansion in developed economies is linked to an increase in CO2 emissions. Consequently, in their pursuit of carbon-neutral policies, these countries are increasingly turning towards renewable energy as a means to enhance resource conservation and efficiency. This research investigates the varied impacts of renewable energy investment, green finance, geopolitical risk, GDP growth, foreign direct investment, and gross fixed capital formation on the carbon emissions of G20 countries. The study uses the CUP-FM (Continuously Updated Fully Modified) and CUP-BC (Continuously Updated Bias-Corrected) estimators, which are sophisticated econometric approaches designed to handle non-stationary panel data and cross-sectional dependency, to produce robust long-term parameter estimates. The CUP-FM estimator adjusts for potential endogeneity and serial correlation, improving the accuracy of long-run relationships in panel data. The CUP-BC estimator provides bias-corrected estimates to further enhance the precision of these long-run connections.The long-term parameter estimates reveal a negative correlation between renewable energy investment, green finance, and carbon emissions. In contrast, foreign direct investment, gross fixed capital formation, GDP growth, and geopolitical risk are positively associated with CO2 emissions. This suggests that financial stability often leads to investments in carbon-heavy economic ventures, thereby implicating economic growth as a contributing factor to environmental degradation in G20 countries.

Beef-on-dairy: current and potential meat quality of dairy-based beef production with culled Holstein cows and Danish Blue × Holstein crossbred calves

ABSTRACT This study characterizes rosé veal from 412 Danish Blue × Holstein heifer and bull calves and beef from 183 highly marbled culled Holstein dairy cows. Furthermore, it investigates the potential for breeding higher meat quality in dairy crossbreeds via genetic parameter estimation. Carcasses were split between the 5th and 6th thoracic vertebra and the meat quality of M. longissimus thoracis was evaluated for pH, color, cooking loss, Warner-Bratzler shear force (WBSF), and chemical intramuscular fat concentration (IMF%). Genetic parameters were estimated using an animal model. The crossbred bull calves outperformed heifers in carcass characteristics, while heifers were superior in redness, cooking loss, WBSF and IMF%. Heritability was significant for WBSF, cooking loss, and IMF%, indicating that meat quality can be improved through selective breeding. The meat quality of highly marbled Holstein cows did not vary across parities, highlighting them as an under-utilized potential source of high-quality beef from a beef-on-dairy system.

Modular Many-Body Quantum Sensors.

Quantum many-body systems undergoing phase transitions have been proposed as probes enabling beyond-classical enhancement of sensing precision. However, this enhancement is usually limited to a very narrow region around the critical point. Here, we systematically develop a modular approach for introducing multiple phase transitions in a many-body system. This naturally allows us to enlarge the region of quantum-enhanced precision by encompassing the newly created phase boundaries. Our approach is general and can be applied to both symmetry-breaking and topological quantum sensors. In symmetry-breaking sensors, we show that the newly created critical points inherit the original universality class and a simple total magnetization measurement already suffices to locate them. In topological sensors, our modular construction creates multiple bands which leads to a rich phase diagram. In both cases, Heisenberg scaling for Hamiltonian parameter estimation is achieved at all the phase boundaries. This can be exploited to create a global sensor which significantly outperforms a uniform probe.

Likelihood and pivotal-based reliability inference for inverse exponentiated Rayleigh distribution under block progressive censoring

ABSTRACT In this paper, focus is on the estimation of survival characteristics, specifically the reliability function, hazard rate function, median time to failure, and differences in different test facilities, using block progressive censored data. Estimations are carried out through both maximum likelihood and pivotal methods, assuming the lifetime distribution of test units follows an inverse exponentiated Rayleigh distribution. The paper derives maximum likelihood estimators for unknown parameters, exploring their existence and uniqueness properties. Approximate confidence intervals for survival characteristics are constructed using the delta method and likelihood theory. Moreover, point and generalized confidence interval estimators are developed through a pivotal quantity-based method. A simulation study is conducted to compare the performance of the proposed approaches, revealing that the pivotal quantity-based approach yields superior estimation results. Finally, the proposed estimates are applied to analyze two real datasets, illustrating their practical applicability.

Bootstrap inference for linear time-varying coefficient models in locally stationary time series

Time-varying coefficient models can capture evolving relationships. However, constructing asymptotic confidence bands for coefficient curves in these models is challenging due to slow convergence rates and the presence of various nuisance parameters. A residual-based sieve bootstrap method has recently been proposed to address these issues. While it successfully produces confidence bands with accurate empirical coverage, its applicability is restricted to strictly stationary processes. We introduce a new bootstrap scheme, the local blockwise wild bootstrap (LBWB), that allows for locally stationary processes. The LBWB can replicate the distribution of the parameter estimates while automatically accounting for nuisance parameters. An extensive simulation study reveals the superior performance of the LBWB compared to various benchmark approaches. It also shows the potential applicability of the LBWB in broader scenarios, including time-varying cointegrating models. We then examine herding effects in the Chinese renewable energy market using the proposed methods. Our findings strongly support the presence of herding behaviors before 2016, aligning with earlier studies. However, contrary to previous research, we find no significant evidence of herding between around 2018 and 2021. Online supplementary materials are available for this article.

Individual participant data meta-analysis to examine linear or non-linear treatment-covariate interactions at multiple time-points for a continuous outcome.

Individual participant data (IPD) meta-analysis projects obtain, harmonise, and synthesise original data from multiple studies. Many IPD meta-analyses of randomised trials are initiated to identify treatment effect modifiers at the individual level, thus requiring statistical modelling of interactions between treatment effect and participant-level covariates. Using a two-stage approach, the interaction is estimated in each trial separately and combined in a meta-analysis. In practice, two complications often arise with continuous outcomes: examining non-linear relationships for continuous covariates and dealing with multiple time-points. We propose a two-stage multivariate IPD meta-analysis approach that summarises non-linear treatment-covariate interaction functions at multiple time-points for continuous outcomes. A set-up phase is required to identify a small set of time-points; relevant knot positions for a spline function, at identical locations in each trial; and a common reference group for each covariate. Crucially, the multivariate approach can include participants or trials with missing outcomes at some time-points. In the first stage, restricted cubic spline functions are fitted and their interaction with each discrete time-point is estimated in each trial separately. In the second stage, the parameter estimates defining these multiple interaction functions are jointly synthesised in a multivariate random-effects meta-analysis model accounting for within-trial and across-trial correlation. These meta-analysis estimates define the summary non-linear interactions at each time-point, which can be displayed graphically alongside confidence intervals. The approach is illustrated using an IPD meta-analysis examining effect modifiers for exercise interventions in osteoarthritis, which shows evidence of non-linear relationships and small gains in precision by analysing all time-points jointly.

Rapid sloshing-free transport of liquids in arbitrarily shaped containers

AbstractWe present a model-based feedforward control strategy suitable for designing swift rest-to-rest maneuvers for liquids in arbitrarily shaped containers. We employ the commonly used equivalent pendulum model to represent the sloshing dynamics and suggest a novel parameter identification scheme suitable for arbitrary container shapes and any number of sloshing modes. By computing natural modes and fluid reaction forces and torques for imposed harmonic container motions via a finite element model, we obtain data for the identification scheme. A fitting procedure then yields highly accurate parameters for a physical pendulum model, where each pendulum represents one sloshing mode. We also provide a thorough analysis of parameter identifiability and guidelines for obtaining robust parameter estimates. The proposed feedforward control method uses a virtual tray pendulum on which we place the container (in the form of its equivalent pendulum model). Designing the virtual tray such that the fluid’s dominant sloshing mode cannot be excited by horizontally moving the tray pendulum pivot effectively zeros out any sloshing motion in this mode. We then exploit the flatness property of the resulting system to design rest-to-rest maneuvers where any residual sloshing motion (in higher modes) can be exactly stopped at the end of the maneuver. The effectiveness of the proposed method is demonstrated through extensive simulations and experimental results using a Martini cocktail glass, whose shape is challenging in terms of sloshing. The experimental results show the successful, accurate suppression of sloshing, validating the efficacy of the proposed concept.

Mathematical Modeling of Electrical Power Systems for Fault Detection

Control networks are critical offices that have to be exceptionally solid and productive. Finding issues in these frameworks is exceptionally imperative for keeping them steady, keeping gear from breaking, and maintaining a strategic distance from long power blackouts. It has gotten to be clear that scientific modeling may be a awesome way to think about and discover issues in control frameworks. Taking after this strategy requires making a scientific show of the control framework, which incorporates motors, transformers, transmission lines, and loads. A few physical laws, like Ohm's Law and Kirchhoff's laws, are frequently utilized within the models to accurately show how electrical circuits work when they are working legitimately and when they aren't. On the off chance that you need to discover issues, you'll utilize numerical models to form flaws like brief circuits, open circuits, and ground deficiencies appear genuine. Through examining how the framework responds to these fake issues, engineers can make programs that can discover genuine flaws and tell where they are when they happen. These calculations frequently utilize strategies like evaluating the state, evaluating the parameters, and recognizing designs. Utilizing real-time information to figure the current state of the framework is what state estimation is all approximately. On the other hand, parameter estimation is all almost finding changes in framework parameters that might cruel there is a problem. Design acknowledgment employments past information and machine learning to discover peculiarities and sort diverse sorts of flaws. Differential conditions, framework polynomial math, and numerical strategies are fair some of the progressed science instruments that are frequently utilized to unravel the complicated conditions that depict how the power framework works. Putting these models beside unused computer innovations like high-performance computing and real-time simulations makes issue location strategies more exact and quicker. As green vitality sources and shrewd network advances are included to control frameworks, they will continue to alter. Mathematical modeling will become more and more important to make sure they work safely and efficiently. It talks about how to mathematically model electrical power systems to find faults, what they can be used for, and what the future holds for this field.

The Comparison of Inverse Gaussian and Gamma Regression: Application on Stunting Data in Jepara

Many research data have distributions other than the normal distribution, called exponential family distributions. The exponential family of distributions includes the inverse Gaussian and Gamma distributions. There are parallels between these two distributions in terms of the kind of random variable and how well they work. Finding the optimal model using inverse Gaussian and Gamma regression on stunting data in Jepara is the goal of this study. Maximum Likelihood Estimation is used for parameter estimation, Maximum Likelihood Ratio Test is used for simultaneous parameter testing, and Wald testing is used for partial parameter testing. For this case, the best model is inverse Gaussian regression. Exclusive breastfeeding, low birth weight babies, clean drinking water facilities, and the number of Integrated Service Post (Posyandu) influence the percentage of stunting in Jepara..

Economic Evaluation of Conservation through Use of an Araucaria angustifolia Provenance and Progeny Test.

Araucaria angustifolia is a species known for its valuable wood and nuts, but it is threatened with extinction. The plantation of forests for genetic resource conservation is a complementary strategy designed to reduce the species' genetic variability loss. This study aimed to evaluate the technical and economic viability of A. angustifolia for genetic conservation through use. The analyzed provenance and progeny trial was established in 1982 in Itapeva, Brazil. It was structured using a compact family blocks design with 110 open-pollinated progenies from five natural populations, three replicates, ten plants per subplot, and 3.0 m × 2.0 m spacing. After 33 years, the trial was evaluated for total height, diameter at breast height, wood volume, and survival. The variance components and genetic parameter estimates were performed using Restricted Maximum Likelihood/Best Linear Unbiased Prediction methods (REML/BLUP) methods with the Selegen software (version 2014). The production and management scenarios were obtained using the SisAraucaria software (version 2003). Sensitivity analysis and economic parameter estimates were obtained through various economic evaluation methods using the Planin software (version 1995). In general, the genetic parameters indicated that the population has enough variability for both conservation and breeding purposes, suggesting technical viability for the establishment of a seed orchard. The economic parameters indicated that the commercialization of wood and araucaria nuts proved to be more profitable than wood production by itself. In conclusion, araucaria genetic conservation through use is a technically and economically viable ex situ conservation strategy.