-
Stable Coherent Systems
Authors:
Rodrigo Iglesias,
Fatemeh Mohammadi,
Patricia Pascual-Ortigosa,
Eduardo Sáenz-de-Cabezón,
Henry P. Wynn
Abstract:
We describe the notion of stability of coherent systems as a framework to deal with redundancy. We define stable coherent systems and show how this notion can help the design of reliable systems. We demonstrate that the reliability of stable systems can be efficiently computed using the algebraic versions of improved inclusion-exclusion formulas and sum of disjoint products.
We describe the notion of stability of coherent systems as a framework to deal with redundancy. We define stable coherent systems and show how this notion can help the design of reliable systems. We demonstrate that the reliability of stable systems can be efficiently computed using the algebraic versions of improved inclusion-exclusion formulas and sum of disjoint products.
△ Less
Submitted 2 June, 2024;
originally announced June 2024.
-
Active Discrimination Learning for Gaussian Process Models
Authors:
Elham Yousefi,
Luc Pronzato,
Markus Hainy,
Werner G. Müller,
Henry P. Wynn
Abstract:
The paper covers the design and analysis of experiments to discriminate between two Gaussian process models, such as those widely used in computer experiments, kriging, sensor location and machine learning. Two frameworks are considered. First, we study sequential constructions, where successive design (observation) points are selected, either as additional points to an existing design or from the…
▽ More
The paper covers the design and analysis of experiments to discriminate between two Gaussian process models, such as those widely used in computer experiments, kriging, sensor location and machine learning. Two frameworks are considered. First, we study sequential constructions, where successive design (observation) points are selected, either as additional points to an existing design or from the beginning of observation. The selection relies on the maximisation of the difference between the symmetric Kullback Leibler divergences for the two models, which depends on the observations, or on the mean squared error of both models, which does not. Then, we consider static criteria, such as the familiar log-likelihood ratios and the Fréchet distance between the covariance functions of the two models. Other distance-based criteria, simpler to compute than previous ones, are also introduced, for which, considering the framework of approximate design, a necessary condition for the optimality of a design measure is provided. The paper includes a study of the mathematical links between different criteria and numerical illustrations are provided.
△ Less
Submitted 21 November, 2022;
originally announced November 2022.
-
Propagating uncertainty in a network of energy models
Authors:
Victoria Volodina,
Nikki Sonenberg,
Jim Q. Smith,
Peter G. Challenor,
Chris J. Dent,
Henry P. Wynn
Abstract:
Computer models are widely used in decision support for energy systems operation, planning and policy. A system of models is often employed, where model inputs themselves arise from other computer models, with each model being developed by different teams of experts. Gaussian Process emulators can be used to approximate the behaviour of complex, computationally intensive models and used to generat…
▽ More
Computer models are widely used in decision support for energy systems operation, planning and policy. A system of models is often employed, where model inputs themselves arise from other computer models, with each model being developed by different teams of experts. Gaussian Process emulators can be used to approximate the behaviour of complex, computationally intensive models and used to generate predictions together with a measure of uncertainty about the predicted model output. This paper presents a computationally efficient framework for propagating uncertainty within a network of models with high-dimensional outputs used for energy planning. We present a case study from a UK county council considering low carbon technologies to transform its infrastructure to reach a net-zero carbon target. The system model considered for this case study is simple, however the framework can be applied to larger networks of more complex models.
△ Less
Submitted 19 April, 2022; v1 submitted 24 January, 2022;
originally announced January 2022.
-
Non-dimensional Newton-Puiseux Expansions
Authors:
C. J. Chapman,
H. P. Wynn,
M. A. Atherton,
R. A. Bates
Abstract:
Recent results in the theory and application of Newton-Puiseux expansions, i.e. fractional power series solutions of equations, suggest further developments within a more abstract algebraic-geometric framework, involving in particular the theory of toric varieties and ideals. Here, we present a number of such developments, especially in relation to the equations of van der Pol, Riccati, and Schröd…
▽ More
Recent results in the theory and application of Newton-Puiseux expansions, i.e. fractional power series solutions of equations, suggest further developments within a more abstract algebraic-geometric framework, involving in particular the theory of toric varieties and ideals. Here, we present a number of such developments, especially in relation to the equations of van der Pol, Riccati, and Schrödinger. Some pure mathematical concepts we are led to are Graver, Gröbner, lattice and circuit bases, combinatorial geometry and differential algebra, and algebraic-differential equations. Two techniques are coordinated: classical dimensional analysis (DA) in applied mathematics and science, and a polynomial and differential-polynomial formulation of asymptotic expansions, referred to here as Newton-Puiseux (NP) expansions. The latter leads to power series with rational exponents, which depend on the choice of the dominant part of the equations, often referred to as the method of "dominant balance". The paper shows, using a new approach to DA based on toric ideals, how dimensionally homogeneous equations may best be non-dimensionalised, and then, with examples involving differential equations, uses recent work on NP expansions to show how non-dimensional parameters affect the results. Our approach finds a natural home within computational algebraic geometry, i.e. at the interface of abstract and algorithmic mathematics.
△ Less
Submitted 9 November, 2021;
originally announced November 2021.
-
Circuits for robust designs
Authors:
Roberto Fontana,
Fabio Rapallo,
Henry P. Wynn
Abstract:
This paper continues the application of circuit theory to experimental design started by the first two authors. The theory gives a very special and detailed representation of the kernel of the design model matrix. This representation turns out to be an appropriate way to study the optimality criteria referred to as robustness: the sensitivity of the design to the removal of design points. Many exa…
▽ More
This paper continues the application of circuit theory to experimental design started by the first two authors. The theory gives a very special and detailed representation of the kernel of the design model matrix. This representation turns out to be an appropriate way to study the optimality criteria referred to as robustness: the sensitivity of the design to the removal of design points. Many examples are given, from classical combinatorial designs to two-level factorial design including interactions. The complexity of the circuit representations are useful because the large range of options they offer, but conversely require the use of dedicated software. Suggestions for speed improvement are made.
△ Less
Submitted 21 June, 2021;
originally announced June 2021.
-
Circuit bases for randomisation
Authors:
Elena Pesce,
Fabio Rapallo,
Eva Riccomagno,
Henry P. Wynn
Abstract:
After a rich history in medicine, randomisation control trials both simple and complex are in increasing use in other areas such as web-based AB testing and planning and design decisions. A main objective is to be able to measure parameters, and contrasts in particular, while guarding against biases from hidden confounders. After careful definitions of classical entities such as contrasts, an alge…
▽ More
After a rich history in medicine, randomisation control trials both simple and complex are in increasing use in other areas such as web-based AB testing and planning and design decisions. A main objective is to be able to measure parameters, and contrasts in particular, while guarding against biases from hidden confounders. After careful definitions of classical entities such as contrasts, an algebraic method based on circuits is introduced which gives a wide choice of randomisation schemes.
△ Less
Submitted 7 May, 2021;
originally announced May 2021.
-
Building A Bayesian Decision Support System for Evaluating COVID-19 Countermeasure Strategies
Authors:
Peter Strong,
Aditi Shenvi,
Xuewen Yu,
K. Nadia Papamichail,
Henry P Wynn,
Jim Q Smith
Abstract:
Decision making in the face of a disaster requires the consideration of several complex factors. In such cases, Bayesian multi-criteria decision analysis provides a framework for decision making. In this paper, we present how to construct a multi-attribute decision support system for choosing between countermeasure strategies, such as lockdowns, designed to mitigate the effects of COVID-19. Such a…
▽ More
Decision making in the face of a disaster requires the consideration of several complex factors. In such cases, Bayesian multi-criteria decision analysis provides a framework for decision making. In this paper, we present how to construct a multi-attribute decision support system for choosing between countermeasure strategies, such as lockdowns, designed to mitigate the effects of COVID-19. Such an analysis can evaluate both the short term and long term efficacy of various candidate countermeasures. The expected utility scores of a countermeasure strategy capture the expected impact of the policies on health outcomes and other measures of population well-being. The broad methodologies we use here have been established for some time. However, this application has many novel elements to it: the pervasive uncertainty of the science; the necessary dynamic shifts between regimes within each candidate suite of countermeasures; and the fast moving stochastic development of the underlying threat all present new challenges to this domain. Our methodology is illustrated by demonstrating in a simplified example how the efficacy of various strategies can be formally compared through balancing impacts of countermeasures, not only on the short term (e.g. COVID-19 deaths) but the medium to long term effects on the population (e.g increased poverty).
△ Less
Submitted 18 December, 2021; v1 submitted 12 January, 2021;
originally announced January 2021.
-
The use of multiple models within an organisation
Authors:
Chris J Dent,
Michael Goldstein,
Andrew Wright,
Henry P. Wynn
Abstract:
Organisations, whether in government, industry or commerce, are required to make decisions in a complex and uncertain environment. The way models are used is intimately connected to the way organisations make decisions and the context in which they make them. Typically, in a complex organisation, multiple related models will often be used in support of a decision. For example, engineering models m…
▽ More
Organisations, whether in government, industry or commerce, are required to make decisions in a complex and uncertain environment. The way models are used is intimately connected to the way organisations make decisions and the context in which they make them. Typically, in a complex organisation, multiple related models will often be used in support of a decision. For example, engineering models might be combined with financial models and macro-economic models in order to decide whether to invest in new production capability. Different parts of a complex organisation might operate their own related models which might then be presented to a central decision maker. Yet in practice, there is little awareness of the practical challenges of using models in a robust way to support decision making. There is significant scope to improve decision making though an enhanced understanding of the role and limitations of modelling and through the application of cutting edge methodologies and organisational best practice. This report is in the form of a 'white paper', whose purpose is to identify key issues for consideration whist postulating tentative approaches to these issues that might be worthy of further exploration, focussing on both technical and organisational aspects. It begins with a framework for consideration of how model-based decisions are made in organisations. It then looks more closely at the questions of uncertainty and multiple models. It then postulates some technical statistical and organisational approaches for managing some of these issues. Finally, it considers the way forward, and the possible focus for further work.
△ Less
Submitted 4 August, 2020;
originally announced August 2020.
-
Algebraic reliability of multi-state $k$-out-of-$n$ systems
Authors:
Patricia Pascual-Ortigosa,
Eduardo Sáenz-de-Cabezón,
Henry P. Wynn
Abstract:
In this paper we review different definitions that multi-state $k$-out-of-$n$ systems have received along the literature and study them in a unified way using the algebra of monomial ideals. We thus obtain formulas and algorithms to compute their reliability and bounds for it.
In this paper we review different definitions that multi-state $k$-out-of-$n$ systems have received along the literature and study them in a unified way using the algebra of monomial ideals. We thus obtain formulas and algorithms to compute their reliability and bounds for it.
△ Less
Submitted 11 April, 2020; v1 submitted 20 February, 2020;
originally announced February 2020.
-
Approximation of probability density functions for PDEs with random parameters using truncated series expansions
Authors:
Giacomo Capodaglio,
Max Gunzburger,
Henry P. Wynn
Abstract:
The probability density function (PDF) of a random variable associated with the solution of a partial differential equation (PDE) with random parameters is approximated using a truncated series expansion. The random PDE is solved using two stochastic finite element methods, Monte Carlo sampling and the stochastic Galerkin method with global polynomials. The random variable is a functional of the s…
▽ More
The probability density function (PDF) of a random variable associated with the solution of a partial differential equation (PDE) with random parameters is approximated using a truncated series expansion. The random PDE is solved using two stochastic finite element methods, Monte Carlo sampling and the stochastic Galerkin method with global polynomials. The random variable is a functional of the solution of the random PDE, such as the average over the physical domain. The truncated series are obtained considering a finite number of terms in the Gram-Charlier or Edgeworth series expansions. These expansions approximate the PDF of a random variable in terms of another PDF, and involve coefficients that are functions of the known cumulants of the random variable. To the best of our knowledge, their use in the framework of PDEs with random parameters has not yet been explored.
△ Less
Submitted 23 September, 2020; v1 submitted 1 October, 2018;
originally announced October 2018.
-
Polarization and depolarization of monomial ideals with application to multistate system reliability
Authors:
Fatemeh Mohammadi,
Patricia Pascual-Ortigosa,
Eduardo Sáenz-de-Cabezón,
Henry P. Wynn
Abstract:
We study the family of depolarizations of a squarefree monomial ideal $I$, i.e. all monomial ideals whose polarization is $I$. We describe a method to find all depolarizations of $I$ and study some of the properties they share and some they do not share. We then apply polarization and depolarization tools to study the reliability of multi-state coherent systems via binary systems and vice versa.
We study the family of depolarizations of a squarefree monomial ideal $I$, i.e. all monomial ideals whose polarization is $I$. We describe a method to find all depolarizations of $I$ and study some of the properties they share and some they do not share. We then apply polarization and depolarization tools to study the reliability of multi-state coherent systems via binary systems and vice versa.
△ Less
Submitted 2 December, 2018; v1 submitted 16 July, 2018;
originally announced July 2018.
-
Efficient multicut enumeration of k-out-of-n:F and consecutive k-out-of-n:F systems
Authors:
Fatemeh Mohammadi,
Eduardo Saenz-de-Cabezon,
Henry P. Wynn
Abstract:
We study multiple simultaneous cut events for k-out-of-n:F and linear consecutive k-out-of-n:F systems in which each component has a constant failure probability. We list the multicuts of these systems and describe the structural differences between them. Our approach, based on combinatorial commutative algebra, allows complete enumeration of the sets of multicuts for both kinds of systems. We als…
▽ More
We study multiple simultaneous cut events for k-out-of-n:F and linear consecutive k-out-of-n:F systems in which each component has a constant failure probability. We list the multicuts of these systems and describe the structural differences between them. Our approach, based on combinatorial commutative algebra, allows complete enumeration of the sets of multicuts for both kinds of systems. We also analyze the computational issues of multicut enumeration and reliability computations.
△ Less
Submitted 21 December, 2017;
originally announced December 2017.
-
Experimental Design Issues in Big Data. The Question of Bias
Authors:
Elena Pesce,
Eva Riccomagno,
Henry P. Wynn
Abstract:
Data can be collected in scientific studies via a controlled experiment or passive observation. Big data is often collected in a passive way, e.g. from social media. In studies of causation great efforts are made to guard against bias and hidden confounders or feedback which can destroy the identification of causation by corrupting or omitting counterfactuals (controls). Various solutions of these…
▽ More
Data can be collected in scientific studies via a controlled experiment or passive observation. Big data is often collected in a passive way, e.g. from social media. In studies of causation great efforts are made to guard against bias and hidden confounders or feedback which can destroy the identification of causation by corrupting or omitting counterfactuals (controls). Various solutions of these problems are discussed, including randomization.
△ Less
Submitted 20 November, 2018; v1 submitted 19 December, 2017;
originally announced December 2017.
-
Optimal experimental design that minimizes the width of simultaneous confidence bands
Authors:
Satoshi Kuriki,
Henry P. Wynn
Abstract:
We propose an optimal experimental design for a curvilinear regression model that minimizes the band-width of simultaneous confidence bands. Simultaneous confidence bands for curvilinear regression are constructed by evaluating the volume of a tube about a curve that is defined as a trajectory of a regression basis vector (Naiman, 1986). The proposed criterion is constructed based on the volume of…
▽ More
We propose an optimal experimental design for a curvilinear regression model that minimizes the band-width of simultaneous confidence bands. Simultaneous confidence bands for curvilinear regression are constructed by evaluating the volume of a tube about a curve that is defined as a trajectory of a regression basis vector (Naiman, 1986). The proposed criterion is constructed based on the volume of a tube, and the corresponding optimal design that minimizes the volume of tube is referred to as the tube-volume optimal (TV-optimal) design. For Fourier and weighted polynomial regressions, the problem is formalized as one of minimization over the cone of Hankel positive definite matrices, and the criterion to minimize is expressed as an elliptic integral. We show that the Möbius group keeps our problem invariant, and hence, minimization can be conducted over cross-sections of orbits. We demonstrate that for the weighted polynomial regression and the Fourier regression with three bases, the tube-volume optimal design forms an orbit of the Möbius group containing D-optimal designs as representative elements.
△ Less
Submitted 30 March, 2019; v1 submitted 13 April, 2017;
originally announced April 2017.
-
Adaptive Experimental Design for Path-following Performance Assessment of Unmanned Vehicles
Authors:
Eleonora Saggini,
Eva Riccomagno,
Massimo Caccia,
Henry P. Wynn
Abstract:
The definition of Good Experimental Methodologies (GEMs) in robotics is a topic of widespread interest due also to the increasing employment of robots in everyday civilian life. The present work contributes to the ongoing discussion on GEMs for Unmanned Surface Vehicles (USVs). It focuses on the definition of GEMs and provides specific guidelines for path-following experiments. Statistically desig…
▽ More
The definition of Good Experimental Methodologies (GEMs) in robotics is a topic of widespread interest due also to the increasing employment of robots in everyday civilian life. The present work contributes to the ongoing discussion on GEMs for Unmanned Surface Vehicles (USVs). It focuses on the definition of GEMs and provides specific guidelines for path-following experiments. Statistically designed experiments (DoE) offer a valid basis for developing an empirical model of the system being investigated. A two-step adaptive experimental procedure for evaluating path-following performance and based on DoE, is tested on the simulator of the Charlie USV. The paper argues the necessity of performing extensive simulations prior to the execution of field trials.
△ Less
Submitted 14 November, 2016;
originally announced November 2016.
-
Types of signature analysis in reliability based on Hilbert series
Authors:
Fatemeh Mohammadi,
Eduardo Sáenz-de-Cabezón,
Henry P. Wynn
Abstract:
The present paper studies multiple failure and signature analysis of coherent systems using the theory of monomial ideals. While system reliability has been studied using Hilbert series of monomial ideals, this is not enough to understand in a deeper sense the ideal structure features that reflect the behavior of the system under multiple simultaneous failures and signature. Therefore, we introduc…
▽ More
The present paper studies multiple failure and signature analysis of coherent systems using the theory of monomial ideals. While system reliability has been studied using Hilbert series of monomial ideals, this is not enough to understand in a deeper sense the ideal structure features that reflect the behavior of the system under multiple simultaneous failures and signature. Therefore, we introduce the lcm-filtration of a monomial ideal, and we study the Hilbert series and resolutions of the corresponding ideals. Given a monomial ideal, we explicitly compute the resolutions for all ideals in the associated lcm-filtration, and we apply this to study coherent systems. Some computational results are shown in examples to demonstrate the usefulness of this approach and the computational issues that arise. We also study the failure distribution from a statistical point of view by means of the algebraic tools described.
△ Less
Submitted 15 October, 2015;
originally announced October 2015.
-
The algebraic method in tree percolation
Authors:
Fatemeh Mohammadi,
Eduardo Sáenz-de-Cabezón,
Henry P. Wynn
Abstract:
We apply the methods of algebraic reliability to the study of percolation on trees. To a complete $k$-ary tree $T_{k,n}$ of depth $n$ we assign a monomial ideal $I_{k,n}$ on $\sum_{i=1}^n k^i$ variables and $k^n$ minimal monomial generators. We give explicit recursive formulae for the Betti numbers of $I_{k,n}$ and their Hilbert series, which allow us to study explicitly percolation on $T_{k,n}$.…
▽ More
We apply the methods of algebraic reliability to the study of percolation on trees. To a complete $k$-ary tree $T_{k,n}$ of depth $n$ we assign a monomial ideal $I_{k,n}$ on $\sum_{i=1}^n k^i$ variables and $k^n$ minimal monomial generators. We give explicit recursive formulae for the Betti numbers of $I_{k,n}$ and their Hilbert series, which allow us to study explicitly percolation on $T_{k,n}$. We study bounds on this percolation and study its asymptotical behavior with the mentioned commutative algebra techniques.
△ Less
Submitted 31 March, 2016; v1 submitted 14 October, 2015;
originally announced October 2015.
-
"Building" exact confidence nets
Authors:
Andrew R. Francis,
Milan Stehlik,
Henry P. Wynn
Abstract:
Confidence nets, that is, collections of confidence intervals that fill out the parameter space and whose exact parameter coverage can be computed, are familiar in nonparametric statistics. Here, the distributional assumptions are based on invariance under the action of a finite reflection group. Exact confidence nets are exhibited for a single parameter, based on the root system of the group. The…
▽ More
Confidence nets, that is, collections of confidence intervals that fill out the parameter space and whose exact parameter coverage can be computed, are familiar in nonparametric statistics. Here, the distributional assumptions are based on invariance under the action of a finite reflection group. Exact confidence nets are exhibited for a single parameter, based on the root system of the group. The main result is a formula for the generating function of the coverage interval probabilities. The proof makes use of the theory of "buildings" and the Chevalley factorization theorem for the length distribution on Cayley graphs of finite reflection groups.
△ Less
Submitted 9 March, 2016; v1 submitted 31 July, 2014;
originally announced July 2014.
-
Empirical geodesic graphs and CAT(k) metrics for data analysis
Authors:
Kei Kobayashi,
Henry P. Wynn
Abstract:
A methodology is developed for data analysis based on empirically constructed geodesic metric spaces. For a probability distribution, the length along a path between two points can be defined as the amount of probability mass accumulated along the path. The geodesic, then, is the shortest such path and defines a geodesic metric. Such metrics are transformed in a number of ways to produce parametri…
▽ More
A methodology is developed for data analysis based on empirically constructed geodesic metric spaces. For a probability distribution, the length along a path between two points can be defined as the amount of probability mass accumulated along the path. The geodesic, then, is the shortest such path and defines a geodesic metric. Such metrics are transformed in a number of ways to produce parametrised families of geodesic metric spaces, empirical versions of which allow computation of intrinsic means and associated measures of dispersion. These reveal properties of the data, based on geometry, such as those that are difficult to see from the raw Euclidean distances. Examples of application include clustering and classification. For certain parameter ranges, the spaces become CAT(0) spaces and the intrinsic means are unique. In one case, a minimal spanning tree of a graph based on the data becomes CAT(0). In another, a so-called "metric cone" construction allows extension to CAT($k$) spaces. It is shown how to empirically tune the parameters of the metrics, making it possible to apply them to a number of real cases.
△ Less
Submitted 14 March, 2019; v1 submitted 13 January, 2014;
originally announced January 2014.
-
Computational algebraic methods in efficient estimation
Authors:
Kei Kobayashi,
Henry P. Wynn
Abstract:
A strong link between information geometry and algebraic statistics is made by investigating statistical manifolds which are algebraic varieties. In particular it it shown how first and second order efficient estimators can be constructed, such as bias corrected Maximum Likelihood and more general estimators, and for which the estimating equations are purely algebraic. In addition it is shown how…
▽ More
A strong link between information geometry and algebraic statistics is made by investigating statistical manifolds which are algebraic varieties. In particular it it shown how first and second order efficient estimators can be constructed, such as bias corrected Maximum Likelihood and more general estimators, and for which the estimating equations are purely algebraic. In addition it is shown how Gröbner basis technology, which is at the heart of algebraic statistics, can be used to reduce the degrees of the terms in the estimating equations. This points the way to the feasible use, to find the estimators, of special methods for solving polynomial equations, such as homotopy continuation methods. Simple examples are given showing both equations and computations. *** The proof of Theorem 2 was corrected by the latest version. Some minor errors were also corrected.
△ Less
Submitted 10 January, 2014; v1 submitted 24 October, 2013;
originally announced October 2013.
-
Dimensional analysis using toric ideals
Authors:
M. A. Atherton,
R. A. Bates,
H. P. Wynn
Abstract:
Classical dimensional analysis is one of the cornerstones of qualitative physics and is also used in the analysis of engineering systems, for example in engineering design. The basic power product relationship in dimensional analysis is identical to one way of defining toric ideals in algebraic geometry, a large and growing field. This paper exploits the toric representation to provide a method fo…
▽ More
Classical dimensional analysis is one of the cornerstones of qualitative physics and is also used in the analysis of engineering systems, for example in engineering design. The basic power product relationship in dimensional analysis is identical to one way of defining toric ideals in algebraic geometry, a large and growing field. This paper exploits the toric representation to provide a method for automatic dimensional analysis for engineering systems. In particular all "primitive", invariants for a particular problem, in a well defined sense, can be found using such methods.
△ Less
Submitted 12 April, 2013;
originally announced April 2013.
-
Subgroup Majorization
Authors:
Andrew R. Francis,
Henry P. Wynn
Abstract:
The extension of majorization (also called the rearrangement ordering), to more general groups than the symmetric (permutation) group, is referred to as $G$-majorization. There are strong results in the case that $G$ is a reflection group and this paper builds on this theory in the direction of subgroups, normal subgroups, quotient groups and extensions. The implications for fundamental cones and…
▽ More
The extension of majorization (also called the rearrangement ordering), to more general groups than the symmetric (permutation) group, is referred to as $G$-majorization. There are strong results in the case that $G$ is a reflection group and this paper builds on this theory in the direction of subgroups, normal subgroups, quotient groups and extensions. The implications for fundamental cones and order-preserving functions are studied. The main example considered is the hyperoctahedral group, which, acting on a vector in $\mathbb R^n$, permutes and changes the signs of components.
△ Less
Submitted 27 November, 2013; v1 submitted 11 March, 2013;
originally announced March 2013.
-
The algebraic method in experimental design
Authors:
Hugo Maruri-Aguilar,
Henry P. Wynn
Abstract:
The algebraic method provides useful techniques to identify models in designs and to understand aliasing of polynomial models. The present note surveys the topic of Gröbner bases in experimental design and then describes the notion of confounding and the algebraic fan of a design. The ideas are illustrated with a variety of design examples ranging from Latin squares to screening designs.
The algebraic method provides useful techniques to identify models in designs and to understand aliasing of polynomial models. The present note surveys the topic of Gröbner bases in experimental design and then describes the notion of confounding and the algebraic fan of a design. The ideas are illustrated with a variety of design examples ranging from Latin squares to screening designs.
△ Less
Submitted 12 July, 2012;
originally announced July 2012.
-
(U,V)-Ordering and a Duality Theorem for Risk Aversion and Lorenz-type Orderings
Authors:
Alessandra Giovagnoli,
Henry P. Wynn
Abstract:
There is a duality theory connecting certain stochastic orderings between cumulative distribution functions F_1,F_2 and stochastic orderings between their inverses F_1^(-1),F_2^(-1). This underlies some theories of utility in the case of the cdf and deprivation indices in the case of the inverse. Under certain conditions there is an equivalence between the two theories. An example is the equivalen…
▽ More
There is a duality theory connecting certain stochastic orderings between cumulative distribution functions F_1,F_2 and stochastic orderings between their inverses F_1^(-1),F_2^(-1). This underlies some theories of utility in the case of the cdf and deprivation indices in the case of the inverse. Under certain conditions there is an equivalence between the two theories. An example is the equivalence between second order stochastic dominance and the Lorenz ordering. This duality is generalised to include the case where there is "distortion" of the cdf of the form v(F) and also of the inverse. A comprehensive duality theorem is presented in a form which includes the distortions and links the duality to the parallel theories of risk and deprivation indices. It is shown that some well-known examples are special cases of the results, including some from the Yaari social welfare theory and the theory of majorization.
△ Less
Submitted 4 August, 2011;
originally announced August 2011.
-
Differential cumulants, hierachical models and monomial ideals
Authors:
Daniel Bruynooghe,
Henry P. Wynn
Abstract:
For a joint probability density function f(x) of a random vector X the mixed partial derivatives of log f(x) can be interpreted as limiting cumulants in an infinitesimally small open neighborhood around x. Moreover, setting them to zero everywhere gives independence and conditional independence conditions. The latter conditions can be mapped, using an algebraic differential duality, into monomial…
▽ More
For a joint probability density function f(x) of a random vector X the mixed partial derivatives of log f(x) can be interpreted as limiting cumulants in an infinitesimally small open neighborhood around x. Moreover, setting them to zero everywhere gives independence and conditional independence conditions. The latter conditions can be mapped, using an algebraic differential duality, into monomial ideal conditions. This provides an isomorphism between hierarchical models and monomial ideals. It is thus shown that certain monomial ideals are associated with particular classes of hierarchical models.
△ Less
Submitted 10 February, 2011;
originally announced February 2011.
-
Smooth supersaturated models
Authors:
Ron A. Bates,
Hugo Maruri-Aguilar,
Henry P. Wynn
Abstract:
In areas such as kernel smoothing and non-parametric regression there is emphasis on smooth interpolation and smooth statistical models. Splines are known to have optimal smoothness properties in one and higher dimensions. It is shown, with special attention to polynomial models, that smooth interpolators can be constructed by first extending the monomial basis and then minimising a measure of s…
▽ More
In areas such as kernel smoothing and non-parametric regression there is emphasis on smooth interpolation and smooth statistical models. Splines are known to have optimal smoothness properties in one and higher dimensions. It is shown, with special attention to polynomial models, that smooth interpolators can be constructed by first extending the monomial basis and then minimising a measure of smoothness with respect to the free parameters in the extended basis. Algebraic methods are a help in choosing the extended basis which can also be found as a saturated basis for an extended experimental design with dummy design points. One can get arbitrarily close to optimal smoothing for any dimension and over any region, giving a simple alternative models of spline type. The relationship to splines is shown in one and two dimensions. A case study is given which includes benchmarking against kriging methods.
△ Less
Submitted 26 September, 2008;
originally announced September 2008.
-
Asymptotic behaviour of a family of gradient algorithms in R^d and Hilbert spaces
Authors:
Luc Pronzato,
Henry P. Wynn,
Anatoly A. Zhigljavsky
Abstract:
The asymptotic behaviour of a family of gradient algorithms (including the methods of steepest descent and minimum residues) for the optimisation of bounded quadratic operators in R^d and Hilbert spaces is analyzed. The results obtained generalize those of Akaike (1959) in several directions. First, all algorithms in the family are shown to have the same asymptotic behaviour (convergence to a tw…
▽ More
The asymptotic behaviour of a family of gradient algorithms (including the methods of steepest descent and minimum residues) for the optimisation of bounded quadratic operators in R^d and Hilbert spaces is analyzed. The results obtained generalize those of Akaike (1959) in several directions. First, all algorithms in the family are shown to have the same asymptotic behaviour (convergence to a two-point attractor), which implies in particular that they have similar asymptotic convergence rates. Second, the analysis also covers the Hilbert space case. A detailed analysis of the stability property of the attractor is provided.
△ Less
Submitted 29 February, 2008;
originally announced February 2008.
-
Monomial ideals and the Scarf complex for coherent systems in reliability theory
Authors:
Beatrice Giglio,
Henry P. Wynn
Abstract:
A certain type of integer grid, called here an echelon grid, is an object found both in coherent systems whose components have a finite or countable number of levels and in algebraic geometry. If α=(α_1,...,α_d) is an integer vector representing the state of a system, then the corresponding algebraic object is a monomial x_1^{α_1}... x_d^{α_d} in the indeterminates x_1,..., x_d. The idea is to r…
▽ More
A certain type of integer grid, called here an echelon grid, is an object found both in coherent systems whose components have a finite or countable number of levels and in algebraic geometry. If α=(α_1,...,α_d) is an integer vector representing the state of a system, then the corresponding algebraic object is a monomial x_1^{α_1}... x_d^{α_d} in the indeterminates x_1,..., x_d. The idea is to relate a coherent system to monomial ideals, so that the so-called Scarf complex of the monomial ideal yields an inclusion-exclusion identity for the probability of failure, which uses many fewer terms than the classical identity. Moreover in the ``general position'' case we obtain via the Scarf complex the tube bounds given by Naiman and Wynn [J. Inequal. Pure Appl. Math. (2001) 2 1-16].
Examples are given for the binary case but the full utility is for general multistate coherent systems and a comprehensive example is given.
△ Less
Submitted 25 June, 2004;
originally announced June 2004.