Parameter Estimation Research Articles

Validating a Calibrated Model of a Groundwater Pump-And-Treat System Using Robust Multiple Regression

Validation of groundwater models is relatively challenging due to the need to reserve scarce water level data for calibration targets. In addition, traditional statistical validation metrics are unintuitive for non-technical audiences and do not directly identify model behaviors that require further refinement. We developed a novel model validation method that analyzes rate change events at pump-and-treat wells and statistically compares the water level responses at nearby monitoring wells between the data and model. The method takes advantage of events that occur alongside ambient pumping, unlike parameter estimation techniques that require independent drawdown or recovery events. The ability of the model to match well connections that are evident (or not evident) in the observations is characterized statistically, leading to four decision scenarios: model matches the observed connection (1) or lack thereof (2), model exhibits a connection that is not observed (3), or model over- or understates the observed connection (4). The method is applied to an FEHM-based groundwater flow and transport model that is shown to match 84.5% of the well connections analyzed. The method provides novel perspectives on the influence of calibration targets on the flow field and suggests that although the overall effect of drawdown targets was to improve the model, the choice of target well pairs creates flow pathways that may be inconsequential during normal operational conditions. The model adequately matches the flow over short spatial scales (<800 m) and over-represents the flow over greater distances (300–1200 m), suggesting the need for “null” drawdown targets in subsequent rounds of calibration.

Population-based cross-sectional survey of cervical cancer screening prevalence and socio-demographic correlates in Bangladeshi women.

Cervical cancer, albeit preventable, is the second-most deadly gynecological cancer in developing nations. Little is known about cervical cancer among Bangladeshi women. This study aims to estimate the prevalence of cervical cancer screening and demographic correlates to identify potential variabilities in screening rates among different demographic groups and regions. This study used secondary data from the WHO STEPS 2018 Survey. We used Bayesian regression to perform the bivariate analyses between the outcome and each explanatory factor, as it generates more acceptable results and improves parameter estimates. The top-ranked socio-demographic factors were identified using a two-step cluster analysis. This method determines the relevance of predictor variables and automatically establishes the number of clusters. The prevalence of Bangladeshi women who had ever been screened for cervical cancer was 6.2%. In the adjusted model, women with the following socio-demographic factors had a higher likelihood of developing cervical cancer: being 18-29 years old (AOR = 3.3, 95% CI: 0.24, 15.27) or 45-59 years old (AOR = 2.8, 95% CI: 1.22, 6.0), currently married (AOR = 2.3, 95% CI: 1.36, 3.70), and employed (AOR = 2.4, 95% CI: 1.40, 4.06). Women in the Barisal division were found to have higher odds of being screened for cervical cancer (AOR = 21, 95% CI: 0.66, 121.97). Cluster analysis found residence status predisposes women to cervical cancer screening. There is a significant potential for substantial reductions in the burden of cervical cancer in Bangladesh by strengthening the application of cervical cancer screening. Future studies should examine how socioeconomic status, culture, and healthcare access affect cervical cancer screening trends for different divisions in Bangladesh. An independent national cancer registry is urgently needed to evaluate screening trends and outcomes.

A note on potential gains in precision of radiation risk estimates from joint analysis.

In estimating radiation-related risk of cancer and other diseases based on the RERF Life Span Study (LSS), joint analyses can be performed where multiple health outcome endpoints are combined in the same model, allowing some parameters to be estimated in common among all endpoints with possible increase in precision of radiation risk and other model parameter estimates. Using as a basis excess relative risk (ERR) and excess absolute risk (EAR) models of the type commonly used in analysis of LSS data at RERF, we use maximum likelihood theory to compute the asymptotic relative standard error of endpoint-specific radiation effect and other parameter estimates using joint analyses as compared to traditional independent analysis. We show that some gains in precision of endpoint-specific radiation risk parameter estimates can be achieved by sharing effect modifier and other model parameters, but only small or negligible gains in precision are achieved for endpoint-specific background modifying or effect modifying parameters when other model parameters are shared. The magnitude of the precision gain for radiation risk estimates depends on the number of endpoints, the baseline incidence rate of the endpoint, and the type of model being used.

A rapid multi-modal parameter estimation technique for LISA

Abstract The Laser Interferometer Space Antenna (LISA) will observe gravitational-wave signals from a wide range of sources, including massive black hole binaries. Although numerous techniques have been developed to perform Bayesian inference for LISA, they are often computationally expensive; analyses often take at least $\sim 1$ month on a single CPU, even when using accelerated techniques. Not only does this make it difficult to concurrently analyse more than one gravitational-wave signal, it also makes it challenging to rapidly produce parameter estimates for possible electromagnetic follow-up campaigns. simple-pe was recently developed to produce rapid parameter estimates for gravitational-wave signals observed with ground-based gravitational-wave detectors. In this work, we extend simple-pe to produce rapid parameter estimates for LISA sources, including the effects of higher order multipole moments. We show that simple-pe infers the source properties of massive black hole binaries in zero-noise at least $\sim 100\times$ faster than existing techniques; $\sim 12$ hours on a single CPU. We further demonstrate that simple-pe can be applied before existing Bayesian techniques to mitigate biases in multi-modal parameter estimation analyses of MBHBs.

Water quality estimates using machine learning techniques in an experimental watershed

ABSTRACT Estimation of the water quality parameters is important to enhance time and cost-effectiveness to that of the conventional approach. This study is aimed to identify the best machine learning (ML) approach to predict concentrations of biochemical oxygen demand (BOD), nitrate, and phosphate. Four ML techniques including decision tree, random forest, gradient boosting, and XGBoost were compared to estimate the water quality parameters based on biophysical (i.e., population, basin area, river slope, water level, and stream flow) and physicochemical properties (i.e., conductivity, turbidity, pH, temperature, and dissolved oxygen) as input parameters. The data were split into training and test sets, and the model performances were evaluated using coefficient of determination (R2), Nash–Sutcliffe efficiency coefficient (NSE), and root mean squared error (RMSE). The mean squared error (MSE) was used as the optimization target. The robust fivefold cross-validation, along with hyperparameter tuning, achieved R2 values of 0.76, 0.67, and 0.71 for phosphate, nitrate, and BOD, and NSE values of 0.73, 0.67, and 0.66, respectively. XGBoost yielded the lowest RMSE across all parameters, showcasing superior performance when considering all metrics performed. In conclusion, ML techniques, particularly with a robust cross-validation technique and hyperparameter optimization, showed good results in water quality parameter prediction.

Full-field Modal Analysis of a Tensegrity Column Using a Three-dimensional Scanning Laser Doppler Vibrometer with a Mirror

Abstract Tensegrity structures become important components of various engineering structures due to their high stiffness, light weight, and deployable capability. Existing studies on their dynamic analyses mainly focus on responses of their nodal points while overlook deformations of their cable and strut members. This study proposes a non-contact approach for experimental modal analysis of a tensegrity structure to identify its three-dimensional (3D) natural frequencies and full-field mode shapes, which include modes with deformations of its cable and strut members. A 3D scanning laser Doppler vibrometer is used with a mirror for extending its field of view to measure full-field vibration of a novel three-strut metal tensegrity column with free boundaries. Tensions and axial stiffnesses of its cable members are determined using natural frequencies of their transverse and longitudinal modes, respectively, to build its theoretical model for dynamic analysis and model validation purposes. Modal assurance criterion (MAC) values between experimental and theoretical mode shapes are used to identify their paired modes. Modal parameters of the first 15 elastic modes of the tensegrity column identified from the experiment, including those of the overall structure and its cable members, can be classified into five mode groups depending on their types. Modes paired between experimental and theoretical results have MAC values larger than 78%. Differences between natural frequencies of paired modes of the tensegrity column are less than 15%. The proposed non-contact 3D vibration measurement approach allows accurate estimation of 3D full-field modal parameters of the tensegrity column.

How the window of visibility varies around polar angle.

Contrast sensitivity, the amount of contrast required to discriminate an object, depends on spatial frequency (SF). The contrast sensitivity function (CSF) peaks at intermediate SFs and drops at other SFs. The CSF varies from foveal to peripheral vision, but only a couple of studies have assessed how the CSF changes with polar angle of the visual field. For many visual dimensions, sensitivity is better along the horizontal than the vertical meridian and at the lower than the upper vertical meridian, yielding polar angle asymmetries. Here, for the first time, to our knowledge, we investigate CSF attributes around polar angle at both group and individual levels and examine the relations in CSFs across locations and individual observers. To do so, we used hierarchical Bayesian modeling, which enables precise estimation of CSF parameters. At the group level, maximum contrast sensitivity and the SF at which the sensitivity peaks are higher at the horizontal than vertical meridian and at the lower than the upper vertical meridian. By analyzing the covariance across observers (n = 28), we found that, at the individual level, CSF attributes (e.g., maximum sensitivity) across locations are highly correlated. This correlation indicates that, although the CSFs differ across locations, the CSF at one location is predictive of that at another location. Within each location, the CSF attributes covary, indicating that CSFs across individuals vary in a consistent manner (e.g., as maximum sensitivity increases, so does the corresponding SF), but more so at the horizontal than the vertical meridian locations. These results show similarities and uncover some critical polar angle differences across locations and individuals, suggesting that the CSF should not be generalized across isoeccentric locations around the visual field.

Improved adaptive robust control of direct-drive pump-valve cooperative brake-by-wire unit with disturbance compensation

Abstract Aiming at the needs of high-level autonomous driving, a direct-drive pump-valve coordinated brake-by-wire unit was designed. To improve the response speed, control accuracy, and robustness of the brake-by-wire unit, an improved adaptive integral robust control method (AIRC-AD) based on a novel adaptive law and disturbance compensation was proposed. By integrating pressure tracking error information and overall parameter estimation error information, a novel parameter adaptive law with fast convergence speed was designed. A finite-time disturbance observer was designed to observe the overall system disturbance, improving the system’s anti-disturbance performance. Finally, the stability of the control method was proven using the Lyapunov stability proof method. The results show that under step target signals, sinusoidal target signals and triangular wave target signals, AIRC-AD has better response speed and control accuracy. In conditions such as increased permanent magnet temperature, increased coil resistance, or reduced flow area due to hydraulic pipeline blockage, AIRC-AD has better adaptability. The unit and the method can meet the demand for precise braking force control in high-level autonomous driving.

Real-time correction of gain nonlinearity in electrostatic actuation for whole-angle micro-shell resonator gyroscope

MEMS gyroscopes are well known for their outstanding advantages in Cost Size Weight and Power (CSWaP), which have inspired great research attention in recent years. A higher signal-to-noise ratio (SNR) for MEMS gyroscopes operating at larger vibrating amplitudes provides improved measuring resolution and ARW performance. However, the increment of amplitude causes strong nonlinear effects of MEMS gyroscopes due to their micron size, which has negative influences on the performance. This paper carries out detailed research on a general nonlinear mechanism on the sensors using parallel-plate capacitive transducers, which is called the gain nonlinearity in electrostatic actuation. The theoretical model established in this paper demonstrates the actuation gain nonlinearity causes the control-force coupling and brings extra angle-dependent bias with the 4th component for the whole-angle gyroscopes, which are verified by the experiments carried out on a micro-shell resonator gyroscope (MSRG). Furthermore, a real-time correction method is proposed to restore a linear response of the electrostatic actuation, which is realized by the gain modification with an online parameter estimation based on the harmonic-component relationship of capacitive detection. This real-time correction method could reduce the 4th component of the angle-dependent bias by over 95% from 0.003°/s to less than 0.0001°/s even under different temperatures. After the correction of actuation gain nonlinearity, the bias instability (BI) of whole-angle MSRG is improved by about 3.5 times from 0.101°/h to 0.029°/h and the scale factor nonlinearity (SFN) is reduced by almost one order of magnitude from 2.02 ppm to 0.21 ppm.

Quantum correlations and parameter estimation for two superconducting qubits interacting with a quantized field.

In the present manuscript, we introduce a quantum system of two superconducting qubits (S-Qs) interacting with a quantized field under the influence of the Kerr nonlinear medium and Ising interaction. We formulate the Hamiltonian of the quantum model and determine the density operator of whole quantum system as well as quantum subsystems. We examine the dynamics of the quantumness measures for subsequent times including the S-Qs entanglement, S-Qs-field entanglement and quantum Fisher information in relation to the system parameters. Finally, we display the connection among the measures of quantumness during the time evolution.

Jupyter Notebooks for Parameter Estimation, Uncertainty Analysis, and Optimization with PEST+.

Jupyter Notebooks for Parameter Estimation, Uncertainty Analysis, and Optimization with PEST+.

Extending a Level Three Order Polynomial Probability Model

Inverse and power extensions were derived for a three order polynomial probability distribution and studied in this paper as post-development. The two extensions serve to make up for the limitations of their baseline model. Some relevant properties—shape of the PDF, moments, survival function, hazard function, quantile function, stress-strength reliability, order statistics and parameter estimations—were studied. The hazard shapes of the distributions were found to be inverted bathtub for the former and three different shapes namely: increasing function, decreasing function and bathtub, for the latter. This implies that the distributions can altogether model many varieties of datasets emanating from different life phenomena. This statement “if S(∞)=0, then H(∞)=∞” was examined and was discovered to apply for exponential distribution but not any of the extended distributions, and the baseline distribution too. Finally, the extensions showed to be competent over both the baseline distribution and their respective counterpart distributions, with respect to the datasets used. For the latter, the superiority was accountably hinged on the extra parameters, since the dataset has an outlier.

Evaluation of spatial models in the estimation of genetic parameters for incidence of frosty pod rot and production in Theobroma cacao full-sib family trials

Evaluation of spatial models in the estimation of genetic parameters for incidence of frosty pod rot and production in Theobroma cacao full-sib family trials

Variational inference for acceleration of SN Ia photometric distance estimation with BayeSN

Abstract Type Ia supernovae (SNe Ia) are standarizable candles whose observed light curves can be used to infer their distances, which can in turn be used in cosmological analyses. As the quantity of observed SNe Ia grows with current and upcoming surveys, increasingly scalable analyses are necessary to take full advantage of these new datasets for precise estimation of cosmological parameters. Bayesian inference methods enable fitting SN Ia light curves with robust uncertainty quantification, but traditional posterior sampling using Markov Chain Monte Carlo (MCMC) is computationally expensive. We present an implementation of variational inference (VI) to accelerate the fitting of SN Ia light curves using the BayeSN hierarchical Bayesian model for time-varying SN Ia spectral energy distributions (SEDs). We demonstrate and evaluate its performance on both simulated light curves and data from the Foundation Supernova Survey with two different forms of surrogate posterior – a multivariate normal and a custom multivariate zero-lower-truncated normal distribution – and compare them with the Laplace Approximation and full MCMC analysis. To validate of our variational approximation, we calculate the pareto-smoothed importance sampling (PSIS) diagnostic, and perform variational simulation-based calibration (VSBC). The VI approximation achieves similar results to MCMC but with an order-of-magnitude speedup for the inference of the photometric distance moduli. Overall, we show that VI is a promising method for scalable parameter inference that enables analysis of larger datasets for precision cosmology.

Two New Families of Local Asymptotically Minimax Lower Bounds in Parameter Estimation

We propose two families of asymptotically local minimax lower bounds on parameter estimation performance. The first family of bounds applies to any convex, symmetric loss function that depends solely on the difference between the estimate and the true underlying parameter value (i.e., the estimation error), whereas the second is more specifically oriented to the moments of the estimation error. The proposed bounds are relatively easy to calculate numerically (in the sense that their optimization is over relatively few auxiliary parameters), yet they turn out to be tighter (sometimes significantly so) than previously reported bounds that are associated with similar calculation efforts, across many application examples. In addition to their relative simplicity, they also have the following advantages: (i) Essentially no regularity conditions are required regarding the parametric family of distributions. (ii) The bounds are local (in a sense to be specified). (iii) The bounds provide the correct order of decay as functions of the number of observations, at least in all the examples examined. (iv) At least the first family of bounds extends straightforwardly to vector parameters.

Novel family of probability generating distributions: Properties and data analysis

Abstract Current probability distributions often fall short when modeling lifespan data with non-monotonic hazard rate shapes, highlighting a gap in probability theory. To address this issue, this study introduces the novel family of probability-generating distributions (NP-G) as a solution to capture complex data patterns better. This research aims to develop and validate the NP-G distributions, which include a two-parameter continuous model based on the Weibull distribution. The study extensively explores the mathematical properties of the NP-G distributions, including series expansions of probability density functions, and evaluates various probabilistic measures such as moments, stress strength, mean deviations, entropy measures, and order statistics. Using the maximum likelihood method for parameter estimation and conducting a numerical simulation analysis, the study assesses the performance of the proposed extended exponentiated Weibull model. The findings demonstrate that this new model offers superior fit and accuracy compared to existing distributions when applied to two referenced datasets. This advancement is significant as it provides a more effective tool for analyzing lifespan data with non-monotonic hazard rates, enhancing both theoretical understanding and practical applications in fields such as reliability engineering and survival analysis.

Metal–Organic Framework-Derived CeO2/Gold Nanospheres in a Highly Sensitive Electrochemical Sensor for Uric Acid Quantification in Milk

In this work, CeBTC (a cerium(III) 1,3,5-benzene-tricarboxylate), was used as a precursor for obtaining CeO2 nanoparticles (nanoceria) with better sensor performances than CeO2 nanoparticles synthesized by the solvothermal method. Metal–organic framework-derived nanoceria (MOFdNC) were functionalized with spheric gold nanoparticles (AuNPs) to further improve non-enzymatic electrode material for highly sensitive detection of prominent biocompound uric acid (UA) at this modified carbon paste electrode (MOFdNC/AuNPs&CPE). X-ray powder diffraction (XRPD) and transmission electron microscopy (TEM) analysis were used for morphological structure characterization of the obtained nanostructures. Cyclic voltammetry and electrochemical impedance spectroscopy, both in an [Fe(CN)6]3−/4− redox system and uric acid standard solutions, were used for the characterization of material electrocatalytic performances, the selection of optimal electrode modifier, and the estimation of nature and kinetic parameters of the electrode process. Square-wave voltammetry (SWV) was chosen, and the optimal parameters of technique and experimental conditions were established for determining uric acid over MOFdNC/AuNPs&CPE. Together with the development of the sensor, the detection procedure was optimized with the following analytical parameters: linear operating ranges of 0.05 to 1 µM and 1 to 50 µM and a detection limit of 0.011 µM, with outstanding repeatability, reproducibility, and stability of the sensor surface. Anti-interference experiments yielded a stable and nearly unchanged current response with negligible or no change in peak potential. After minor sample pretreatment, the proposed electrode was successfully applied for the quantification of UA in milk.

Parameter estimation and hypothesis tests in logistic model for complex correlated data

Parameter estimation and hypothesis tests in logistic model for complex correlated data

Reconstruction of Nearshore Surface Gravity Wave Heights From Distributed Acoustic Sensing Data

AbstractDistributed Acoustic Sensing (DAS) is a photonics technology converting seafloor telecommunications and optical fiber cables into dense arrays of strain sensors, allowing to monitor various oceanic physical processes. Yet, several applications are hindered by the limited knowledge of the transfer function between geophysical variables and DAS measurements. This study investigates the quantitative relationship between surface gravity DAS‐recorded wave‐generated strain signals along the seafloor and the pressure at a colocated sensor. A remarkable linear correlation is found over various sea conditions allowing us to reliably determine significant wave heights from DAS data. Utilizing linear wave potential theory, we derive an analytical transfer function linking cable deformation and wave kinematic parameters. This transfer function provides a first quantification of the effects related to surface gravity waves and fiber responses. Our results validate DAS's potential for real‐time reconstruction of the surface gravity wave spectrum over extended coastal areas. It also enables the estimation of waves hydraulic parameters at depth without the need from offshore deployments.

Detection of QPO Soft Lag during the Outburst of Swift J1727.8-1613: Estimation of Intrinsic Parameters from Spectral Study

The recently discovered bright transient black hole candidate Swift J1727.8-1613 is studied in a broad energy range (0.5–79 keV) using combined NICER and NuSTAR data taken on 2023 August 29. A prominent type C quasiperiodic oscillation (QPO) at 0.89 ± 0.01 Hz with its harmonic was observed in NICER data of 0.5–10 keV. Interestingly, the harmonic becomes weaker in the lower energy bands (0.5–1 and 1–3 keV). We also report the first detection of a soft time lag of 0.014 ± 0.001 s at the QPO frequency between harder (3–10 keV) and softer (0.5–3 keV) band photons observed with the NICER/X-ray timing instrument. This indicates that the inclination of the accretion disk in the binary system might be high. From the detailed spectral analysis with the relxill reflection model, we found the disk inclination angle of the source to be ∼85°. We discuss how the accretion flow configuration inferred from spectral analysis can help us understand the origin of QPOs and soft lag in this source.