Estimates Of Model Parameters Research Articles

Modeling Glucose, Insulin, C-Peptide, and Lactate Interplay in Adolescents During an Oral Glucose Tolerance Test.

Lactate is not considered just a "waste product" of anaerobic glycolysis anymore. It has been proved to play a key role in several metabolic diseases, such as in the metabolic dysfunction-associated steatotic liver disease, obesity, and diabetes. The capability of simulating glucose-insulin-lactate interaction would be useful to design and test drugs targeting lactate metabolism in such pathological conditions. Minimal models are available, which describe and quantify glucose-lactate interaction but models to simulate postprandial glucose-insulin-C-peptide-lactate time courses are missing. The aim of this study is to fill this gap. Starting from the Padova Type 2 Diabetes Simulator (T2DS), we first added a description of glucose-lactate kinetics and then created a population of 100 in silico subjects to match glucose-insulin-C-peptide-lactate data of 44 adolescents with/without obesity who underwent a standard oral glucose tolerance test (OGTT) of 75 g. The developed model accurately predicts all molecules time courses, guaranteeing precise model parameter estimates (percent coefficient of variation [CV%] median [25th-75th percentile] = 19 [9-29]%). The generated in silico population shows good agreement with the clinical data in terms of area under the curve (AUC) (P = .6, .6, .9, .6 for glucose, insulin, C-peptide, and lactate, respectively) and parameter distributions (P > .1). We have developed a simulator to describe glucose, insulin, C-peptide, and lactate kinetics during an OGTT, which captures the behavior of a real population of adolescents with/without obesity both in terms of average and intersubject variability. Such simulator can be used to investigate the pharmacodynamics of drugs targeting lactate metabolic pathway in various pathological conditions.

Achieving Occam's razor: Deep learning for optimal model reduction.

All fields of science depend on mathematical models. Occam's razor refers to the principle that good models should exclude parameters beyond those minimally required to describe the systems they represent. This is because redundancy can lead to incorrect estimates of model parameters from data, and thus inaccurate or ambiguous conclusions. Here, we show how deep learning can be powerfully leveraged to apply Occam's razor to model parameters. Our method, FixFit, uses a feedforward deep neural network with a bottleneck layer to characterize and predict the behavior of a given model from its input parameters. FixFit has three major benefits. First, it provides a metric to quantify the original model's degree of complexity. Second, it allows for the unique fitting of data. Third, it provides an unbiased way to discriminate between experimental hypotheses that add value versus those that do not. In three use cases, we demonstrate the broad applicability of this method across scientific domains. To validate the method using a known system, we apply FixFit to recover known composite parameters for the Kepler orbit model and a dynamic model of blood glucose regulation. In the latter, we demonstrate the ability to fit the latent parameters to real data. To illustrate how the method can be applied to less well-established fields, we use it to identify parameters for a multi-scale brain model and reduce the search space for viable candidate mechanisms.

Different Statistical Inference Algorithms for the New Pareto Distribution Based on Type-II Progressively Censored Competing Risk Data with Applications

In this research, the statistical inference of unknown lifetime parameters is proposed in the presence of independent competing risks using a progressive Type-II censored dataset. The lifetime distribution associated with a failure mode is assumed to follow the new Pareto distribution, with consideration given to two distinct competing failure reasons. Maximum likelihood estimators (MLEs) for the unknown model parameters, as well as reliability and hazard functions, are derived, noting that they are not expressible in closed form. The Newton–Raphson, expectation maximization (EM), and stochastic expectation maximization (SEM) methods are employed to generate maximum likelihood (ML) estimations. Approximate confidence intervals for the unknown parameters, reliability, and hazard rate functions are constructed using the normal approximation of the MLEs and the normal approximation of the log-transformed MLEs. Additionally, the missing information principle is utilized to derive the closed form of the Fisher information matrix, which, in turn, is used with the delta approach to calculate confidence intervals for reliability and hazards. Bayes estimators are derived under both symmetric and asymmetric loss functions, with informative and non-informative priors considered, including independent gamma distributions for informative priors. The Monte Carlo Markov Chain sampling approach is employed to obtain the highest posterior density credible intervals and Bayesian point estimates for unknown parameters and reliability characteristics. A Monte Carlo simulation is conducted to assess the effectiveness of the proposed techniques, with the performances of the Bayes and maximum likelihood estimations examined using average values and mean squared errors as benchmarks. Interval estimations are compared in terms of average lengths and coverage probabilities. Real datasets are considered and examined for each topic to provide illustrative examples.

Modeling and Optimal Control of Infectious Diseases

We propose a stochastic model of infectious disease transmission that is more realistic than those found in the literature. The model is based on jump-diffusion processes. However, it is defined in such a way that the number of people susceptible to be infected decreases over time, which is the case for a population of fixed size. Next, we consider the problem of finding the optimal control of the proposed model. The dynamic programming equation satisfied by the value function is derived. Estimators of the various model parameters are obtained.

Simultaneous Inference of Multiple Binary Endpoints in Biomedical Research: Small Sample Properties of Multiple Marginal Models and a Resampling Approach.

In biomedical research, the simultaneous inference of multiple binary endpoints may be of interest. In such cases, an appropriate multiplicity adjustment is required that controls the family-wise error rate, which represents the probability of making incorrect test decisions. In this paper, we investigate two approaches that perform single-step -value adjustments that also take into account the possible correlation between endpoints. A rather novel and flexible approach known as multiple marginal models is considered, which is based on stacking of the parameter estimates of the marginal models and deriving their joint asymptotic distribution. We also investigate a nonparametric vector-based resampling approach, and we compare both approaches with the Bonferroni method by examining the family-wise error rate and power for different parameter settings, including low proportions and small sample sizes. The results show that the resampling-based approach consistently outperforms the other methods in terms of power, while still controlling the family-wise error rate. The multiple marginal models approach, on the other hand, shows a more conservative behavior. However, it offers more versatility in application, allowing for more complex models or straightforward computation of simultaneous confidence intervals. The practical application of the methods is demonstrated using a toxicological dataset from the National ToxicologyProgram.

Comparação das distribuições Normal, Logística, Laplace e t de Student para o erro experimental na descrição bayesiana do acúmulo de matéria seca de Allium sativum

ABSTRACT: This study assessed distributions associated with Bayesian nonlinear modeling error in the description of total plant dry matteraccumulation (TDMA) of Allium sativumas a function of days after planting (DAP). According to the DIC criterion, Logistic and Gompertzmodels that use student’s t distribution error exhibited the highest DIC with logistic error distribution. In general, the difference of DIC in all the scenarios was not more than 5.The Bayes factor (BF) criterion showed no difference in the Logistic and Gompertzmodel fit when four distributions are used for the errors, where BF values do not exceed 2. Posterior distributions and the usual estimators of Logistic and Gompertz model parameters were similar even forvaried error distribution. In summary, there was no difference in the use of 4 distributions associated with the modeling error of garlic plant growth by the Bayes factor, whereby the results showed that alternating between error distributions significantly changes the number of Markov Chain Monte Carlo (MCMC) iterations.

Evaluating the robustness of parameter estimates in cognitive models: A meta-analytic review of multinomial processing tree models across the multiverse of estimation methods.

Researchers have become increasingly aware that data-analysis decisions affect results. Here, we examine this issue systematically for multinomial processing tree (MPT) models, a popular class of cognitive models for categorical data. Specifically, we examine the robustness of MPT model parameter estimates that arise from two important decisions: the level of data aggregation (complete-pooling, no-pooling, or partial-pooling) and the statistical framework (frequentist or Bayesian). These decisions span a multiverse of estimation methods. We synthesized the data from 13,956 participants (164 published data sets) with a meta-analytic strategy and analyzed the magnitude of divergence between estimation methods for the parameters of nine popular MPT models in psychology (e.g., process-dissociation, source monitoring). We further examined moderators as potential sources of divergence. We found that the absolute divergence between estimation methods was small on average (<.04; with MPT parameters ranging between 0 and 1); in some cases, however, divergence amounted to nearly the maximum possible range (.97). Divergence was partly explained by few moderators (e.g., the specific MPT model parameter, uncertainty in parameter estimation), but not by other plausible candidate moderators (e.g., parameter trade-offs, parameter correlations) or their interactions. Partial-pooling methods showed the smallest divergence within and across levels of pooling and thus seem to be an appropriate default method. Using MPT models as an example, we show how transparency and robustness can be increased in the field of cognitive modeling. (PsycInfo Database Record (c) 2024 APA, all rights reserved).

Cerebrospinal Fluid Pharmacokinetics of Nicardipine Following Intrathecal Administration in Subarachnoid Hemorrhage Patients.

Subarachnoid hemorrhage (SAH) is a devastating type of stroke, leading to high mortality and morbidity rates. Cerebral vasospasm and delayed cerebral ischemia (DCI) are common complications following SAH that contribute significantly to the poor outcomes observed in these patients. Intrathecal (IT) nicardipine delivered via an existing external ventricular drain is an off-label intervention that has been shown to be correlated with reduced DCI and improved patient outcomes. The current study aims to characterize the population pharmacokinetic (popPK) properties of intermittent IT nicardipine. Following informed consent, serial cerebrospinal fluid (CSF) samples were obtained from 16 SAH patients (50.4 ± 9.3 years old; 13 females) treated with IT nicardipine every 6 h (q6h, n = 8) or every 8 h (q8h, n = 8) for an average of 72 ± 21 doses. High-performance liquid chromatography was used to quantify CSF concentration from each sample. Our popPK analysis showed that the CSF pharmacokinetics of IT nicardipine in the cohort was adequately described by a two-compartment model with a lag time. Model parameter estimates were reliable (relative standard error <50%). Intracranial pressure influenced both the total clearance and the central volume of nicardipine (i.e., negative correlation, P <-.001). Calculated PK parameters were similar between q6h and q8h dosing regimens. Despite a small cohort of SAH patients, we successfully developed a popPK model to describe the nicardipine disposition kinetics in the CSF following IT administration. These findings may help inform future clinical trials designed to examine the optimal dosing of IT nicardipine.

Bayesian estimation and prediction under progressive-stress accelerated life test for a log-logistic model

Accelerated life tests play a very critical role in reliability analysis because highly reliable products are being produced with recent advanced technologies to sustain the market demand and competition. A progressive-stress accelerated life test is one of the kinds of accelerated life tests that allows applied stress to change continuously. Considering the importance of progressive-stress accelerated life tests, this paper deals with progressive-stress accelerated life tests when the lifetimes of units follow the log-logistic distribution and are progressively type-II censored where the associated scale parameter conforms to the inverse power law. Different estimators of model parameters are derived using maximum likelihood and Bayesian methods. Interval estimation is also considered. In the sequel, approximate confidence, bootstrap, and Bayes credible intervals are constructed. Bayes estimators are obtained under squared error, LINEX, and entropy loss functions using proper and improper prior distributions. A simulation study is conducted based on various censoring schemes. The coverage percentages and interval widths are computed via Monte Carlo simulations. Bayes predictive estimates and intervals are also obtained. Finally, two different accelerated life test data are analyzed for illustration purposes.

Lomax tangent generalized family of distributions: Characteristics, simulations, and applications on hydrological-strength data

This article proposes and discusses a novel approach for generating trigonometric G-families using hybrid generalizers of distributions. The proposed generalizer is constructed by utilizing the tangent trigonometric function and distribution function of base model G(x). The newly proposed family of uni-variate continuous distributions is named the “Lomax Tangent Generalized Family of Distributions (LT-G)” and structural-mathematical-statistical properties are derived. Some special and sub-models of the proposed family are also presented.A Weibull-based model, ‘The Lomax Tangent Weibull (LT-W) Distribution,” is discussed and the plots of density (pdf) and hazard (hrf) functions are also explained. Model parameter estimates are estimated by employing the maximum likelihood estimation (MLE) procedure. The accuracy of the MLEs is evaluated through Monte Carlo simulation. Last but not least, to demonstrate the flexibility and potential of the proposed distribution, two actual hydrological and strength data sets are analyzed. The obtained results are compared with well-known, competitive, and related existing distributions.

Reliability inference for a family of inverted exponentiated distributions under block progressive Type II censoring

In this article, estimation of survival characteristics is discussed when life tests are performed against different test facilities using block progressive Type-II censoring. Under the assumption that the lifetimes of test units follow a family of inverted exponentiated distributions, estimates are derived under classical and hierarchical Bayesian frameworks when the differences in different test facilities are taken into consideration. Maximum likelihood estimators are obtained and their existence and uniqueness properties are also established. In addition, Bayes estimators are provided under a hierarchical framework, and a hybrid Metropolis-Hastings sampling procedure is proposed for complex posterior computation. The performance of all the proposed approaches is compared via extensive simulation studies and two real-life data are analyzed for illustration. From numerical results we observe that the hierarchical approach provides relatively better estimates for model parameters.

Estimation for generalized linear cointegration regression models through composite quantile regression approach

This paper introduces a meaningful approach employing composite quantile regression (CQR) to estimate generalized linear cointegration regression models. We elucidate the fundamental structure of the proposed model by presenting its underlying expressions and derive the asymptotic distribution of the estimates of model parameters. Through extensive simulations, our findings demonstrate the superior robustness and precision of the CQR method compared to ordinary least squares (OLS) and quantile regression (QR) approaches. The application of the model to economic and financial variables highlights its significant academic and practical value.

Statistical Modeling of Right-Censored Spatial Data Using Gaussian Random Fields

Consider a fixed number of clustered areas identified by their geographical coordinates that are monitored for the occurrences of an event such as a pandemic, epidemic, or migration. Data collected on units at all areas include covariates and environmental factors. We apply a probit transformation to the time to event and embed an isotropic spatial correlation function into our models for better modeling as compared to existing methodologies that use frailty or copula. Composite likelihood technique is employed for the construction of a multivariate Gaussian random field that preserves the spatial correlation function. The data are analyzed using counting process and geostatistical formulation that led to a class of weighted pairwise semiparametric estimating functions. The estimators of model parameters are shown to be consistent and asymptotically normally distributed under infill-type asymptotic spatial statistics. Detailed small sample numerical studies that are in agreement with theoretical results are provided. The foregoing procedures are applied to the leukemia survival data in Northeast England. A comparison to existing methodologies provides improvement.

A constrained maximum likelihood approach to developing well-calibrated models for predicting binary outcomes

The added value of candidate predictors for risk modeling is routinely evaluated by comparing the performance of models with or without including candidate predictors. Such comparison is most meaningful when the estimated risk by the two models are both unbiased in the target population. Very often data for candidate predictors are sourced from nonrepresentative convenience samples. Updating the base model using the study data without acknowledging the discrepancy between the underlying distribution of the study data and that in the target population can lead to biased risk estimates and therefore an unfair evaluation of candidate predictors. To address this issue assuming access to a well-calibrated base model, we propose a semiparametric method for model fitting that enforces good calibration. The central idea is to calibrate the fitted model against the base model by enforcing suitable constraints in maximizing the likelihood function. This approach enables unbiased assessment of model improvement offered by candidate predictors without requiring a representative sample from the target population, thus overcoming a significant practical challenge. We study theoretical properties for model parameter estimates, and demonstrate improvement in model calibration via extensive simulation studies. Finally, we apply the proposed method to data extracted from Penn Medicine Biobank to inform the added value of breast density for breast cancer risk assessment in the Caucasian woman population.

Abstract PO4-02-12: MDA iLobular Prognostic Tool: A Novel Approach for Risk Stratification in Invasive Lobular Carcinoma

Abstract Background: Invasive lobular carcinoma (ILC) is a distinct subtype of breast cancer that differs from the more common invasive ductal carcinoma (IDC) by its unique characteristics that influence prognosis. Recent studies have shown that ILC is associated with worse long-term outcomes compared to IDC. Risk stratification is essential in breast cancer for tailoring treatment plans to the individual patient’s needs, maximizing the chances of successful outcomes and minimizing unnecessary interventions. Although there are multiple tools available for breast cancer, there is currently not a risk stratification tool that is specific for ILC. The MDA iLobular prognostic tool is the first risk stratification tool that has been developed specifically for ILC patients. Methods: We retrospectively searched for patients treated at MD Anderson Cancer Center with a diagnosis of stage I-III ILC in our prospectively collected electronic database. The study focused on two primary endpoints: Overall Survival (OS) and Distant Recurrence Free Survival (DRFS). We utilized univariate and multivariate Cox Proportional Hazard (PH) regression models to assess the statistical significance of all variables. The univariate Cox analysis identified prognostic factors, which were further analyzed using backward and stepwise multivariate Cox proportional hazards regression analysis. This process helped identify statistically significant prognostic factors that were included in the final multivariate models. We estimated hazard ratios and 95% confidence intervals for each factor, considering a P-value of &amp;lt; 0.05 as statistically significant. To evaluate the performance of the fitted multivariate Cox PH regression models for OS and DRFS, we randomly divided two-thirds of the data points into a training dataset, while the remaining one-third constituted the test dataset. We assessed the discrimination capacity of each model using Harrell's C-index. Results: The study included a total of 4,216 female patients, with a median age of 56 years. The median pathological tumor size was 2 mm, and the median number of lymph nodes was one. Among these patients, 1,376 were pre-menopausal, while 2,837 were post-menopausal. The training cohort was a subset of 2,950 patients and the test cohort was a subset of 1,266 patients. After evaluating various prognostic models, we identified the model with the highest prognostic accuracy for OS and DRFS. This selected model demonstrated a Harrell's C-index of 0.704 for OS and 0.718 for DRFS on the training dataset and a Harrell's C-index of 0.702 for OS and 0.671 for DRFS on the test dataset. The model incorporated several covariates, including age at the time of diagnosis, number of lymph nodes, pathological tumor size (mm), ER status, tumor grade, ILC histology, and the presence or absence of concomitant LCIS (Table 1). Conclusion: In conclusion, the MDA iLobular prognostic tool represents a significant advancement as the first dedicated tool for risk stratifying ILC, providing valuable guidance for tailoring therapy to individual ILC patients based on their specific risk profiles. Table 1. Multivariate Cox proportional hazard model parameter estimates for OS and DRFS Citation Format: Jason Mouabbi, Sarah Pasyar, Roland Bassett, Akshara Singareeka Raghavendra, Amy Hassan, Rachel Layman, Debu Tripathy. MDA iLobular Prognostic Tool: A Novel Approach for Risk Stratification in Invasive Lobular Carcinoma [abstract]. In: Proceedings of the 2023 San Antonio Breast Cancer Symposium; 2023 Dec 5-9; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2024;84(9 Suppl):Abstract nr PO4-02-12.

Pedestrian injury severities resulting from vehicle/pedestrian intersection crashes: An assessment of COVID-contributing temporal shifts

Pedestrian mobility has become an increasingly important concern in transportation system analysis because of its positive impacts on the environment and healthy lifestyles. However, pedestrian safety in a vehicle-dominated transportation network remains a concern and potential barrier to pedestrian mobility, with pedestrian intersection safety being of particular concern. In addition, it is important to understand how pedestrian safety has been affected by the COVID-19 pandemic, perhaps permanently shifting pedestrian injury risks. This research seeks to provide insight into how pedestrian injury risks at intersections have changed as a result of the pandemic by estimating a series pedestrian injury severity models. To do so, unconstrained and partially constrained random parameters multinomial logit models with heterogeneity in the means of random parameters were estimated. Using Florida data, two one-year periods (one year before and one year after the COVID-19 pandemic) were defined based on vehicle miles traveled. The sample includes 3,780 single pedestrian-single vehicle crashes (2,348 from daytime and 1,432 from nighttime). A broad range of variables was considered to assess how the parameters may have shifted between the before and after periods. A series of likelihood ratio tests were conducted to examine the stability of model parameter estimates across the predefined time periods as well as to determine the differences between the daytime and nighttime crash injury severity outcomes. The results show that the nighttime crashes experienced more temporal shifts relative to daytime crashes. The findings also showed that both pedestrian and driver behavior played key temporally-shifting roles before and after the COVID-19 pandemic period. Finally, the out-of-sample simulations suggest that pedestrian injuries have become more severe after the pandemic.

Bayesian modeling of quantiles of body mass index among under-five children in Ethiopia

BackgroundBody Mass Index (BMI) is a measurement of nutritional status, which is a vital pre-condition for good health. The prevalence of childhood malnutrition and the potential long-term health risks associated with obesity in Ethiopia have recently increased globally. The main objective of this study was to investigate the factors associated with the quantiles of under-five children’s BMI in Ethiopia.MethodsData on 5,323 children, aged between 0-59 months from March 21, 2019, to June 28, 2019, were obtained from the Ethiopian Mini Demographic Health Survey (EMDHS, 2019), based on the standards set by the World Health Organization. The study used a Bayesian quantile regression model to investigate the association of factors with the quantiles of under-five children’s body mass index. Markov Chain Monte Carlo (MCMC) with Gibbs sampling was used to estimate the country-specific marginal posterior distribution estimates of model parameters, using the Brq R package.ResultsOut of a total of 5323 children included in this study, 5.09% were underweight (less than 12.92 BMI), 10.05% were overweight (BMI: 17.06 – 18.27), and 5.02% were obese (greater than or equal to 18.27 BMI) children’s. The result of the Bayesian quantile regression model, including marginal posterior credible intervals (CIs), showed that for the prediction of the 0.05 quantile of BMI, the current age of children [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document}= -0.007, 95% CI :(-0.01, -0.004)], the region Afar [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = - 0.32, 95% CI: (-0.57, -0.08)] and Somalia[β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = -0.72, 95% CI: (-0.96, -0.49)] were negatively associated with body mass index while maternal age [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = 0.01, 95% CI: (0.005, 0.02)], mothers primary education [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document}= 0.19, 95% CI: (0.08, 0.29)], secondary and above [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = 0.44, 95% CI: (0.29, 0.58)], and family follows protestant [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = 0.22, 95% CI: (0.07, 0.37)] were positively associated with body mass index. In the prediction of the 0.95 (or 0.85?) quantile of BMI, in the upper quantile, still breastfeeding [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = -0.25, 95% CI: (-0.41, -0.10)], being female [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = -0.13, 95% CI: (-0.23, -0.03)] were negatively related while wealth index [β\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\upbeta$$\\end{document} = 0.436, 95% CI: (0.25, 0.62)] was positively associated with under-five children’s BMI.ConclusionsIn conclusion, the research findings indicate that the percentage of lower and higher BMI for under-five children in Ethiopia is high. Factors such as the current age of children, sex of children, maternal age, religion of the family, region and wealth index were found to have a significant impact on the BMI of under-five children both at lower and upper quantile levels. Thus, these findings highlight the need for administrators and policymakers to devise and implement strategies aimed at enhancing the normal or healthy weight status among under-five children in Ethiopia.

How robust are estimates of key parameters in standard viral dynamic models?

Mathematical models of viral infection have been developed, fitted to data, and provide insight into disease pathogenesis for multiple agents that cause chronic infection, including HIV, hepatitis C, and B virus. However, for agents that cause acute infections or during the acute stage of agents that cause chronic infections, viral load data are often collected after symptoms develop, usually around or after the peak viral load. Consequently, we frequently lack data in the initial phase of viral growth, i.e., when pre-symptomatic transmission events occur. Missing data may make estimating the time of infection, the infectious period, and parameters in viral dynamic models, such as the cell infection rate, difficult. However, having extra information, such as the average time to peak viral load, may improve the robustness of the estimation. Here, we evaluated the robustness of estimates of key model parameters when viral load data prior to the viral load peak is missing, when we know the values of some parameters and/or the time from infection to peak viral load. Although estimates of the time of infection are sensitive to the quality and amount of available data, particularly pre-peak, other parameters important in understanding disease pathogenesis, such as the loss rate of infected cells, are less sensitive. Viral infectivity and the viral production rate are key parameters affecting the robustness of data fits. Fixing their values to literature values can help estimate the remaining model parameters when pre-peak data is missing or limited. We find a lack of data in the pre-peak growth phase underestimates the time to peak viral load by several days, leading to a shorter predicted growth phase. On the other hand, knowing the time of infection (e.g., from epidemiological data) and fixing it results in good estimates of dynamical parameters even in the absence of early data. While we provide ways to approximate model parameters in the absence of early viral load data, our results also suggest that these data, when available, are needed to estimate model parameters more precisely.

CovXtreme : MATLAB software for non-stationary penalised piecewise constant marginal and conditional extreme value models

The covXtreme software provides functionality for estimation of marginal and conditional extreme value models, non-stationary with respect to covariates, and environmental design contours. Generalised Pareto (GP) marginal models of peaks over threshold are estimated, using a piecewise-constant representation for the variation of GP threshold and scale parameters on the (potentially multidimensional) covariate domain of interest. The conditional variation of one or more associated variates, given a large value of a single conditioning variate, is described using the conditional extremes model of Heffernan and Tawn (2004), the slope term of which is also assumed to vary in a piecewise constant manner with covariates. Optimal smoothness of marginal and conditional extreme value model parameters with respect to covariates is estimated using cross-validated roughness-penalised maximum likelihood estimation. Uncertainties in model parameter estimates due to marginal and conditional extreme value threshold choice, and sample size, are quantified using a bootstrap resampling scheme. Estimates of environmental contours using various schemes, including the direct sampling approach of Huseby et al. 2013, are calculated by simulation or numerical integration under fitted models. The software was developed in MATLAB for metocean applications, but is applicable generally to multivariate samples of peaks over threshold data. The software and case study data can be downloaded from GitHub, with an accompanying user guide.

Statistical inference of comparative generalized inverted exponential populations under joint adaptive progressive type-II censored samples

Design of efficient test plans to assessing product reliability rapidly, is a critical step to obtain product reliability accurately. In this paper, we are adopted the problem of statistical inference of two generalized inverted exponential populations in a competing duration. To save the balance between the optimal test time and the number of failures needing for statistical inference, joint adaptive type-II hybrid progressive censoring scheme is applied. The model parameters are estimated with classical methods (maximum likelihood and bootstrap) for point estimate and the corresponding confidence intervals. Also, Bayes estimators of model parameters are computed under importance sample technique with the corresponding Bayes credible intervals. A real data sets are analyzed for illustrating purposes. The results are compared and assessed through Monte Carlo simulation study.