RESUMEN

This paper discusses regression analysis of current status data with dependent censoring, a problem that often occurs in many areas such as cross-sectional studies, epidemiological investigations and tumorigenicity experiments. Copula model-based methods are commonly employed to tackle this issue. However, these methods often face challenges in terms of model and parameter identification. The primary aim of this paper is to propose a copula-based analysis for dependent current status data, where the association parameter is left unspecified. Our method is based on a general class of semiparametric linear transformation models and parametric copulas. We demonstrate that the proposed semiparametric model is identifiable under certain regularity conditions from the distribution of the observed data. For inference, we develop a sieve maximum likelihood estimation method, using Bernstein polynomials to approximate the nonparametric functions involved. The asymptotic consistency and normality of the proposed estimators are established. Finally, to demonstrate the effectiveness and practical applicability of our method, we conduct an extensive simulation study and apply the proposed method to a real data example.

RESUMEN

The immune checkpoint inhibitor anti-PD-1, commonly used in cancer immunotherapy, has not been successful as a monotherapy for the highly aggressive brain cancer glioblastoma. However, when used in conjunction with a CC-chemokine receptor-2 (CCR2) antagonist, anti-PD-1 has shown efficacy in preclinical studies. In this paper, we aim to optimize treatment regimens for this combination immunotherapy using optimal control theory. We extend a treatment-free glioblastoma-immune dynamics ODE model to include interventions with anti-PD-1 and the CCR2 antagonist. An optimized regimen increases the survival of an average mouse from 32 days post-tumor implantation without treatment to 111 days with treatment. We scale this approach to a virtual murine cohort to evaluate mortality and quality of life concerns during treatment, and predict survival, tumor recurrence, or death after treatment. A parameter identifiability analysis identifies five parameters suitable for personalizing treatment within the virtual cohort. Sampling from these five practically identifiable parameters for the virtual murine cohort reveals that personalized, optimized regimens enhance survival: 84% of the virtual mice survive to day 100, compared to 60% survival in a previously studied experimental regimen. Subjects with high tumor growth rates and low T cell kill rates are identified as more likely to die during and after treatment due to their compromised immune systems and more aggressive tumors. Notably, the MDSC death rate emerges as a long-term predictor of either disease-free survival or death.

RESUMEN

Antimicrobial heteroresistance refers to the presence of different subpopulations with heterogeneous antimicrobial responses within the same bacterial isolate, so they show reduced susceptibility compared with the main population. Though it is widely accepted that heteroresistance can play a crucial role in the outcome of antimicrobial treatments, predictive Antimicrobial Resistance (AMR) models accounting for bacterial heteroresistance are still scarce and need to be refined as the techniques to measure heteroresistance become standardised and consistent conclusions are drawn from data. In this work, we propose a multivariate Birth-Death (BD) model of bacterial heteroresistance and analyse its properties in detail. Stochasticity in the population dynamics is considered since heteroresistance is often characterised by low initial frequencies of the less susceptible subpopulations, those mediating AMR transmission and potentially leading to treatment failure. We also discuss the utility of the heteroresistance model for practical applications and calibration under realistic conditions, demonstrating that it is possible to infer the model parameters and heteroresistance distribution from time-kill data, i.e., by measuring total cell counts alone and without performing any heteroresistance test.

Asunto(s)

RESUMEN

The advent of rapid and inexpensive sequencing technologies has necessitated the development of computationally efficient methods for analyzing sequence data for many genes simultaneously in a phylogenetic framework. The coalescent process is the most commonly used model for linking the underlying genealogies of individual genes with the global species-level phylogeny, but inference under the coalescent model is computationally daunting in the typical inference frameworks (e.g., the likelihood and Bayesian frameworks) due to the dimensionality of the space of both gene trees and species trees. Here we consider estimation of the branch lengths in fixed species trees with three or four taxa, and show that these branch lengths are identifiable. We also show that for three and four taxa simple estimators for the branch lengths can be derived based on observed site pattern frequencies. Properties of these estimators, such as their asymptotic variances and large-sample distributions, are examined, and performance of the estimators is assessed using simulation. Finally, we use these estimators to develop a hypothesis test that can be used to delimit species under the coalescent model for three or four putative taxa.

RESUMEN

In competing risks data, in practice, there may be lack of information or uncertainty about the true failure type, termed as 'missing failure type', for some subjects. We consider a general pattern of missing failure type in which we observe, if not the true failure type, a set of possible failure types containing the true one. In this work, we focus on both parametric and non-parametric estimation based on current status data with two competing risks and the above-mentioned missing failure type. Here, the missing probabilities are assumed to be time-dependent, that is, dependent on both failure and monitoring time points, in addition to being dependent on the true failure type. This makes the missing mechanism non-ignorable. We carry out maximum likelihood estimation and obtain the asymptotic properties of the estimators. Simulation studies are conducted to investigate the finite sample properties of the estimators. Finally, the methods are illustrated through a data set on hearing loss.

RESUMEN

When hybridization or other forms of lateral gene transfer have occurred, evolutionary relationships of species are better represented by phylogenetic networks than by trees. While inference of such networks remains challenging, several recently proposed methods are based on quartet concordance factors-the probabilities that a tree relating a gene sampled from the species displays the possible 4-taxon relationships. Building on earlier results, we investigate what level-1 network features are identifiable from concordance factors under the network multispecies coalescent model. We obtain results on both topological features of the network, and numerical parameters, uncovering a number of failures of identifiability related to 3-cycles in the network. Addressing these identifiability issues is essential for designing statistically consistent inference methods.

Asunto(s)

RESUMEN

Metabolic fluxes are the rates of life-sustaining chemical reactions within a cell and metabolites are the components. Determining the changes in these fluxes is crucial to understanding diseases with metabolic causes and consequences. Kinetic flux profiling (KFP) is a method for estimating flux that utilizes data from isotope tracing experiments. In these experiments, the isotope-labeled nutrient is metabolized through a pathway and integrated into the downstream metabolite pools. Measurements of proportion labeled for each metabolite in the pathway are taken at multiple time points and used to fit an ordinary differential equations model with fluxes as parameters. We begin by generalizing the process of converting diagrams of metabolic pathways into mathematical models composed of differential equations and algebraic constraints. The scaled differential equations for proportions of unlabeled metabolite contain parameters related to the metabolic fluxes in the pathway. We investigate flux parameter identifiability given data collected only at the steady state of the differential equation. Next, we give criteria for valid parameter estimations in the case of a large separation of timescales with fast-slow analysis. Bayesian parameter estimation on simulated data from KFP experiments containing both irreversible and reversible reactions illustrates the accuracy and reliability of flux estimations. These analyses provide constraints that serve as guidelines for the design of KFP experiments to estimate metabolic fluxes.

RESUMEN

COVID-19 highlighted the importance of considering human behavior change when modeling disease dynamics. This led to developing various models that incorporate human behavior. Our objective is to contribute to an in-depth, mathematical examination of such models. Here, we consider a simple deterministic compartmental model with endogenous incorporation of human behavior (i.e., behavioral feedback) through transmission in a classic Susceptible-Exposed-Infectious-Recovered (SEIR) structure. Despite its simplicity, the SEIR structure with behavior (SEIRb) was shown to perform well in forecasting, especially compared to more complicated models. We contrast this model with an SEIR model that excludes endogenous incorporation of behavior. Both models assume permanent immunity to COVID-19, so we also consider a modification of the models which include waning immunity (SEIRS and SEIRSb). We perform equilibria, sensitivity, and identifiability analyses on all models and examine the fidelity of the models to replicate COVID-19 data across the United States. Endogenous incorporation of behavior significantly improves a model's ability to produce realistic outbreaks. While the two endogenous models are similar with respect to identifiability and sensitivity, the SEIRSb model, with the more accurate assumption of the waning immunity, strengthens the initial SEIRb model by allowing for the existence of an endemic equilibrium, a realistic feature of COVID-19 dynamics. When fitting the model to data, we further consider the addition of simple seasonality affecting disease transmission to highlight the explanatory power of the models.

Asunto(s)

RESUMEN

Cognitive diagnostic models (CDMs) are a popular family of discrete latent variable models that model students' mastery or deficiency of multiple fine-grained skills. CDMs have been most widely used to model categorical item response data such as binary or polytomous responses. With advances in technology and the emergence of varying test formats in modern educational assessments, new response types, including continuous responses such as response times, and count-valued responses from tests with repetitive tasks or eye-tracking sensors, have also become available. Variants of CDMs have been proposed recently for modeling such responses. However, whether these extended CDMs are identifiable and estimable is entirely unknown. We propose a very general cognitive diagnostic modeling framework for arbitrary types of multivariate responses with minimal assumptions, and establish identifiability in this general setting. Surprisingly, we prove that our general-response CDMs are identifiable under Q -matrix-based conditions similar to those for traditional categorical-response CDMs. Our conclusions set up a new paradigm of identifiable general-response CDMs. We propose an EM algorithm to efficiently estimate a broad class of exponential family-based general-response CDMs. We conduct simulation studies under various response types. The simulation results not only corroborate our identifiability theory, but also demonstrate the superior empirical performance of our estimation algorithms. We illustrate our methodology by applying it to a TIMSS 2019 response time dataset.

RESUMEN

This paper considers the practical and ethical issues related to the death of a sperm donor. It looks at whether sperm banks should check whether the donor is alive at the time the spermatozoa is made available. Knowing that the donor has died in combination with the cause of death can provide important medical information but the chance is very small. However, when the information is available, it may help to decide whether to tell previous recipients and/or to block the remaining samples for future use. A second advantage may be that the donor's offspring can be informed that contact will not be possible and that recipients who are planning to order spermatozoa from an identity-release donor can be told that the donor has died. However, these advantages presume that identifiable equals contactable while there is no strict link between these two conditions.

RESUMEN

In this paper we develop a four compartment within-host model of nutrition and HIV. We show that the model has two equilibria: an infection-free equilibrium and infection equilibrium. The infection free equilibrium is locally asymptotically stable when the basic reproduction number $ \mathcal{R}_0 < 1 $, and unstable when $ \mathcal{R}_0 > 1 $. The infection equilibrium is locally asymptotically stable if $ \mathcal{R}_0 > 1 $ and an additional condition holds. We show that the within-host model of HIV and nutrition is structured to reveal its parameters from the observations of viral load, CD4 cell count and total protein data. We then estimate the model parameters for these 3 data sets. We have also studied the practical identifiability of the model parameters by performing Monte Carlo simulations, and found that the rate of clearance of the virus by immunoglobulins is practically unidentifiable, and that the rest of the model parameters are only weakly identifiable given the experimental data. Furthermore, we have studied how the data frequency impacts the practical identifiability of model parameters.

Asunto(s)

RESUMEN

Parameter identification involves the estimation of undisclosed parameters within a system based on observed data and mathematical models. In this investigation, we employ DAISY to meticulously examine the structural identifiability of parameters of a within-host SARS-CoV-2 epidemic model, taking into account an array of observable datasets. Furthermore, Monte Carlo simulations are performed to offer a comprehensive practical analysis of model parameters. Lastly, sensitivity analysis is employed to ascertain that decreasing the replication rate of the SARS-CoV-2 virus and curbing the infectious period are the most efficacious measures in alleviating the dissemination of COVID-19 amongst hosts.

RESUMEN

Many commonly used mathematical models in the field of mathematical biology involve challenges of parameter non-identifiability. Practical non-identifiability, where the quality and quantity of data does not provide sufficiently precise parameter estimates is often encountered, even with relatively simple models. In particular, the situation where some parameters are identifiable and others are not is often encountered. In this work we apply a recent likelihood-based workflow, called Profile-Wise Analysis (PWA), to non-identifiable models for the first time. The PWA workflow addresses identifiability, parameter estimation, and prediction in a unified framework that is simple to implement and interpret. Previous implementations of the workflow have dealt with idealised identifiable problems only. In this study we illustrate how the PWA workflow can be applied to both structurally non-identifiable and practically non-identifiable models in the context of simple population growth models. Dealing with simple mathematical models allows us to present the PWA workflow in a didactic, self-contained document that can be studied together with relatively straightforward Julia code provided on GitHub . Working with simple mathematical models allows the PWA workflow prediction intervals to be compared with gold standard full likelihood prediction intervals. Together, our examples illustrate how the PWA workflow provides us with a systematic way of dealing with non-identifiability, especially compared to other approaches, such as seeking ad hoc parameter combinations, or simply setting parameter values to some arbitrary default value. Importantly, we show that the PWA workflow provides insight into the commonly-encountered situation where some parameters are identifiable and others are not, allowing us to explore how uncertainty in some parameters, and combinations of parameters, regardless of their identifiability status, influences model predictions in a way that is insightful and interpretable.

Asunto(s)

RESUMEN

This research aimed to provide experimental evidence on whether identifying an edible animal by a name and specific preferences encourages children to perceive the animal as more similar to humans, increases their willingness to befriend the animal, and makes them less willing to consume it. In two pre-registered studies involving 208 preschool children, participants were presented with pictures of pigs (Study 1) and chickens (Study 2). In the identifiability condition, one animal was depicted with individual qualities such as a name and personal preferences, while in the non-identifiability condition, animals were portrayed with characteristics representative of the entire species. The children then rated their desire to befriend and consume the animal, while in Study 2, they also rated the animal's similarity to humans. The results revealed that animal identifiability led to higher perceived similarity to humans, increased the desire to befriend it, and reduced inclination to consume the animal. These findings highlight animal identifiability's powerful and robust effect on children's attitudes toward edible animals.

Asunto(s)

RESUMEN

Quantitative models of sequence-function relationships are ubiquitous in computational biology, e.g., for modeling the DNA binding of transcription factors or the fitness landscapes of proteins. Interpreting these models, however, is complicated by the fact that the values of model parameters can often be changed without affecting model predictions. Before the values of model parameters can be meaningfully interpreted, one must remove these degrees of freedom (called "gauge freedoms" in physics) by imposing additional constraints (a process called "fixing the gauge"). However, strategies for fixing the gauge of sequence-function relationships have received little attention. Here we derive an analytically tractable family of gauges for a large class of sequence-function relationships. These gauges are derived in the context of models with all-order interactions, but an important subset of these gauges can be applied to diverse types of models, including additive models, pairwise-interaction models, and models with higher-order interactions. Many commonly used gauges are special cases of gauges within this family. We demonstrate the utility of this family of gauges by showing how different choices of gauge can be used both to explore complex activity landscapes and to reveal simplified models that are approximately correct within localized regions of sequence space. The results provide practical gauge-fixing strategies and demonstrate the utility of gauge-fixing for model exploration and interpretation.

RESUMEN

BACKGROUND: Integration of a patient's non-invasive imaging data in a digital twin (DT) of the heart can provide valuable insight into the myocardial disease substrates underlying left ventricular (LV) mechanical discoordination. However, when generating a DT, model parameters should be identifiable to obtain robust parameter estimations. In this study, we used the CircAdapt model of the human heart and circulation to find a subset of parameters which were identifiable from LV cavity volume and regional strain measurements of patients with different substrates of left bundle branch block (LBBB) and myocardial infarction (MI). To this end, we included seven patients with heart failure with reduced ejection fraction (HFrEF) and LBBB (study ID: 2018-0863, registration date: 2019-10-07), of which four were non-ischemic (LBBB-only) and three had previous MI (LBBB-MI), and six narrow QRS patients with MI (MI-only) (study ID: NL45241.041.13, registration date: 2013-11-12). Morris screening method (MSM) was applied first to find parameters which were important for LV volume, regional strain, and strain rate indices. Second, this parameter subset was iteratively reduced based on parameter identifiability and reproducibility. Parameter identifiability was based on the diaphony calculated from quasi-Monte Carlo simulations and reproducibility was based on the intraclass correlation coefficient ( ICC ) obtained from repeated parameter estimation using dynamic multi-swarm particle swarm optimization. Goodness-of-fit was defined as the mean squared error ( χ 2 ) of LV myocardial strain, strain rate, and cavity volume. RESULTS: A subset of 270 parameters remained after MSM which produced high-quality DTs of all patients ( χ 2 < 1.6), but minimum parameter reproducibility was poor ( ICC min = 0.01). Iterative reduction yielded a reproducible ( ICC min = 0.83) subset of 75 parameters, including cardiac output, global LV activation duration, regional mechanical activation delay, and regional LV myocardial constitutive properties. This reduced subset produced patient-resembling DTs ( χ 2 < 2.2), while septal-to-lateral wall workload imbalance was higher for the LBBB-only DTs than for the MI-only DTs (p < 0.05). CONCLUSIONS: By applying sensitivity and identifiability analysis, we successfully determined a parameter subset of the CircAdapt model which can be used to generate imaging-based DTs of patients with LV mechanical discoordination. Parameters were reproducibly estimated using particle swarm optimization, and derived LV myocardial work distribution was representative for the patient's underlying disease substrate. This DT technology enables patient-specific substrate characterization and can potentially be used to support clinical decision making.

Asunto(s)

RESUMEN

Practical limitations of quality and quantity of data can limit the precision of parameter identification in mathematical models. Model-based experimental design approaches have been developed to minimise parameter uncertainty, but the majority of these approaches have relied on first-order approximations of model sensitivity at a local point in parameter space. Practical identifiability approaches such as profile-likelihood have shown potential for quantifying parameter uncertainty beyond linear approximations. This research presents a genetic algorithm approach to optimise sample timing across various parameterisations of a demonstrative PK-PD model with the goal of aiding experimental design. The optimisation relies on a chosen metric of parameter uncertainty that is based on the profile-likelihood method. Additionally, the approach considers cases where multiple parameter scenarios may require simultaneous optimisation. The genetic algorithm approach was able to locate near-optimal sampling protocols for a wide range of sample number (n = 3-20), and it reduced the parameter variance metric by 33-37% on average. The profile-likelihood metric also correlated well with an existing Monte Carlo-based metric (with a worst-case r > 0.89), while reducing computational cost by an order of magnitude. The combination of the new profile-likelihood metric and the genetic algorithm demonstrate the feasibility of considering the nonlinear nature of models in optimal experimental design at a reasonable computational cost. The outputs of such a process could allow for experimenters to either improve parameter certainty given a fixed number of samples, or reduce sample quantity while retaining the same level of parameter certainty.

Asunto(s)

RESUMEN

Accurately reconstructing the reticulate histories of polyploids remains a central challenge for understanding plant evolution. Although phylogenetic networks can provide insights into relationships among polyploid lineages, inferring networks may be hindered by the complexities of homology determination in polyploid taxa. We use simulations to show that phasing alleles from allopolyploid individuals can improve phylogenetic network inference under the multispecies coalescent by obtaining the true network with fewer loci compared to haplotype consensus sequences or sequences with heterozygous bases represented as ambiguity codes. Phased allelic data can also improve divergence time estimates for networks, which is helpful for evaluating allopolyploid speciation hypotheses and proposing mechanisms of speciation. To achieve these outcomes in empirical data, we present a novel pipeline that leverages a recently developed phasing algorithm to reliably phase alleles from polyploids. This pipeline is especially appropriate for target enrichment data, where depth of coverage is typically high enough to phase entire loci. We provide an empirical example in the North American Dryopteris fern complex that demonstrates insights from phased data as well as the challenges of network inference. We establish that our pipeline (PATÉ: Phased Alleles from Target Enrichment data) is capable of recovering a high proportion of phased loci from both diploids and polyploids. These data may improve network estimates compared to using haplotype consensus assemblies by accurately inferring the direction of gene flow, but statistical non-identifiability of phylogenetic networks poses a barrier to inferring the evolutionary history of reticulate complexes.

RESUMEN

Individualizing treatment assignment can improve outcomes for diseases with patient-to-patient variability in comparative treatment effects. When a clinical trial demonstrates that some patients improve on treatment while others do not, it is tempting to assume that treatment effect heterogeneity exists. However, if outcome variability is mainly driven by factors other than variability in the treatment effect, investigating the extent to which covariate data can predict differential treatment response is a potential waste of resources. Motivated by recent meta-analyses assessing the potential of individualizing treatment for major depressive disorder using only summary statistics, we provide a method that uses summary statistics widely available in published clinical trial results to bound the benefit of optimally assigning treatment to each patient. We also offer alternate bounds for settings in which trial results are stratified by another covariate. Our upper bounds can be especially informative when they are small, as there is then little benefit to collecting additional covariate data. We demonstrate our approach using summary statistics from a depression treatment trial. Our methods are implemented in the rct2otrbounds R package.

Asunto(s)

RESUMEN

Nonparametric item response models provide a flexible framework in psychological and educational measurements. Douglas (Psychometrika 66(4):531-540, 2001) established asymptotic identifiability for a class of models with nonparametric response functions for long assessments. Nevertheless, the model class examined in Douglas (2001) excludes several popular parametric item response models. This limitation can hinder the applications in which nonparametric and parametric models are compared, such as evaluating model goodness-of-fit. To address this issue, We consider an extended nonparametric model class that encompasses most parametric models and establish asymptotic identifiability. The results bridge the parametric and nonparametric item response models and provide a solid theoretical foundation for the applications of nonparametric item response models for assessments with many items.