skip to main content

Title: Bayesian Sequential Design Based on Dual Objectives for Accelerated Life Tests
Traditional accelerated life test plans are typically based on optimizing the C-optimality for minimizing the variance of an interested quantile of the lifetime distribution. These methods often rely on some specified planning values for the model parameters, which are usually unknown prior to the actual tests. The ambiguity of the specified parameters can lead to suboptimal designs for optimizing the reliability performance of interest. In this paper, we propose a sequential design strategy for life test plans based on considering dual objectives. In the early stage of the sequential experiment, we suggest allocating more design locations based on optimizing the D-optimality to quickly gain precision in the estimated model parameters. In the later stage of the experiment, we can allocate more observations based on optimizing the C-optimality to maximize the precision of the estimated quantile of the lifetime distribution. We compare the proposed sequential design strategy with existing test plans considering only a single criterion and illustrate the new method with an example on the fatigue testing of polymer composites.  more » « less
Award ID(s):
1565314 1838271
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
Statistical Quality Technologies
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Past work on optimizing fabrication plans given a carpentry design can provide Pareto-optimal plans trading off between material waste, fabrication time, precision, and other considerations. However, when developing fabrication plans, experts rarely restrict to a single design , instead considering families of design variations , sometimes adjusting designs to simplify fabrication. Jointly exploring the design and fabrication plan spaces for each design is intractable using current techniques. We present a new approach to jointly optimize design and fabrication plans for carpentered objects. To make this bi-level optimization tractable, we adapt recent work from program synthesis based on equality graphs (e-graphs), which encode sets of equivalent programs. Our insight is that subproblems within our bi-level problem share significant substructures. By representing both designs and fabrication plans in a new bag of parts (BOP) e-graph, we amortize the cost of optimizing design components shared among multiple candidates. Even using BOP e-graphs, the optimization space grows quickly in practice. Hence, we also show how a feedback-guided search strategy dubbed Iterative Contraction and Expansion on E-graphs (ICEE) can keep the size of the e-graph manageable and direct the search towards promising candidates. We illustrate the advantages of our pipeline through examples from the carpentry domain. 
    more » « less
  2. Breast cancer is the most commonly diagnosed cancer type worldwide. Given high survivorship, increased focus has been placed on long-term treatment outcomes and patient quality of life. While breast-conserving surgery (BCS) is the preferred treatment strategy for early-stage breast cancer, anticipated healing and breast deformation (cosmetic) outcomes weigh heavily on surgeon and patient selection between BCS and more aggressive mastectomy procedures. Unfortunately, surgical outcomes following BCS are difficult to predict, owing to the complexity of the tissue repair process and significant patient-to-patient variability. To overcome this challenge, we developed a predictive computational mechanobiological model that simulates breast healing and deformation following BCS. The coupled biochemical-biomechanical model incorporates multi-scale cell and tissue mechanics, including collagen deposition and remodeling, collagen-dependent cell migration and contractility, and tissue plastic deformation. Available human clinical data evaluating cavity contraction and histopathological data from an experimental porcine lumpectomy study were used for model calibration. The computational model was successfully fit to data by optimizing biochemical and mechanobiological parameters through Gaussian process surrogates. The calibrated model was then applied to define key mechanobiological parameters and relationships influencing healing and breast deformation outcomes. Variability in patient characteristics including cavity-to-breast volume percentage and breast composition were further evaluated to determine effects on cavity contraction and breast cosmetic outcomes, with simulation outcomes aligning well with previously reported human studies. The proposed model has the potential to assist surgeons and their patients in developing and discussing individualized treatment plans that lead to more satisfying post-surgical outcomes and improved quality of life. 
    more » « less
  3. Despite the success of multiscale modeling in science and engineering, embedding molecular-level information into nonlinear reactor design and control optimization problems remains challenging. In this work, we propose a computationally tractable scale-bridging approach that incorporates information from multi-product microkinetic (MK) models with thousands of rates and chemical species into nonlinear reactor design optimization problems. We demonstrate reduced-order kinetic (ROK) modeling approaches for catalytic oligomerization in shale gas processing. We assemble a library of six candidate ROK models based on literature and MK model structure. We find that three metrics—quality of fit (e.g., mean squared logarithmic error), thermodynamic consistency (e.g., low conversion of exothermic reactions at high temperatures), and model identifiability—are all necessary to train and select ROK models. The ROK models that closely mimic the structure of the MK model offer the best compromise to emulate the product distribution. Using the four best ROK models, we optimize the temperature profiles in staged reactors to maximize conversions to heavier oligomerization products. The optimal temperature starts at 630–900K and monotonically decreases to approximately 560 K in the final stage, depending on the choice of ROK model. For all models, staging increases heavier olefin production by 2.5% and there is minimal benefit to more than four stages. The choice of ROK model, i.e., model-form uncertainty, results in a 22% difference in the objective function, which is twice the impact of parametric uncertainty; we demonstrate sequential eigendecomposition of the Fisher information matrix to identify and fix sloppy model parameters, which allows for more reliable estimation of the covariance of the identifiable calibrated model parameters. First-order uncertainty propagation determines this parametric uncertainty induces less than a 10% variability in the reactor optimization objective function. This result highlights the importance of quantifying model-form uncertainty, in addition to parametric uncertainty, in multi-scale reactor and process design and optimization. Moreover, the fast dynamic optimization solution times suggest the ROK strategy is suitable for incorporating molecular information in sequential modular or equation-oriented process simulation and optimization frameworks. 
    more » « less
  4. Abstract

    Structured population models are among the most widely used tools in ecology and evolution. Integral projection models (IPMs) use continuous representations of how survival, reproduction and growth change as functions of state variables such as size, requiring fewer parameters to be estimated than projection matrix models (PPMs). Yet, almost all published IPMs make an important assumption that size‐dependent growth transitions are or can be transformed to be normally distributed. In fact, many organisms exhibit highly skewed size transitions. Small individuals can grow more than they can shrink, and large individuals may often shrink more dramatically than they can grow. Yet, the implications of such skew for inference from IPMs has not been explored, nor have general methods been developed to incorporate skewed size transitions into IPMs, or deal with other aspects of real growth rates, including bounds on possible growth or shrinkage.

    Here, we develop a flexible approach to modelling skewed growth data using a modified beta regression model. We propose that sizes first be converted to a (0,1) interval by estimating size‐dependent minimum and maximum sizes through quantile regression. Transformed data can then be modelled using beta regression with widely available statistical tools. We demonstrate the utility of this approach using demographic data for a long‐lived plant, gorgonians and an epiphytic lichen. Specifically, we compare inferences of population parameters from discrete PPMs to those from IPMs that either assume normality or incorporate skew using beta regression or, alternatively, a skewed normal model.

    The beta and skewed normal distributions accurately capture the mean, variance and skew of real growth distributions. Incorporating skewed growth into IPMs decreases population growth and estimated life span relative to IPMs that assume normally distributed growth, and more closely approximate the parameters of PPMs that do not assume a particular growth distribution. A bounded distribution, such as the beta, also avoids the eviction problem caused by predicting some growth outside the modelled size range.

    Incorporating biologically relevant skew in growth data has important consequences for inference from IPMs. The approaches we outline here are flexible and easy to implement with existing statistical tools.

    more » « less
  5. Abstract

    The study compared the life cycle environmental impacts of three coastal flood management strategies: grey infrastructure (levee), green–grey infrastructure (levee and oyster reef), and a do-nothing scenario, considering the flood damage of a single flooding event in the absence of protection infrastructure. A case study was adopted from a New Orleans, Louisiana residential area to facilitate the comparison. Hazus software, design guidelines, reports, existing projects, and literature were utilized as foreground data for modelling materials. A process-based life cycle assessment was used to assess environmental impacts. The life cycle environmental impacts included global warming, ozone depletion, acidification, eutrophication, smog formation, resource depletion, ecotoxicity, and various human health effects. The ecoinvent database was used for the selected life cycle unit processes. The mean results show green–grey infrastructure as the most promising strategy across most impact categories, reducing 47% of the greenhouse gas (GHG) emissions compared to the do-nothing strategy. Compared to grey infrastructure, green–grey infrastructure mitigates 13%–15% of the environmental impacts while providing equivalent flood protection. A flooding event with a 100-year recurrence interval in the study area is estimated at 34 million kg of CO2equivalent per kilometre of shoreline, while grey and green–grey infrastructure mitigating such flooding is estimated to be 21 and 18 million kg, respectively. This study reinforced that coastal flooding environmental impacts are primarily caused by rebuilding damaged houses, especially concrete and structural timber replacement, accounting for 90% of GHG emissions, with only 10% associated with flood debris waste treatment. The asphalt cover of the levee was identified as the primary contributor to environmental impacts in grey infrastructure, accounting for over 75% of GHG emissions during construction. We found that there is an important interplay between grey and green infrastructure and optimizing their designs can offer solutions to sustainable coastal flood protection.

    more » « less