-
Exploiting the polyhedral geometry of stochastic linear bilevel programming
Authors:
Gonzalo Muñoz,
David Salas,
Anton Svensson
Abstract:
We study linear bilevel programming problems whose lower-level objective is given by a random cost vector with known distribution. We consider the case where this distribution is nonatomic, allowing to reformulate the problem of the leader using the Bayesian approach in the sense of Salas and Svensson (2023), with a decision-dependent distribution that concentrates on the vertices of the feasible…
▽ More
We study linear bilevel programming problems whose lower-level objective is given by a random cost vector with known distribution. We consider the case where this distribution is nonatomic, allowing to reformulate the problem of the leader using the Bayesian approach in the sense of Salas and Svensson (2023), with a decision-dependent distribution that concentrates on the vertices of the feasible set of the follower's problem. We call this a vertex-supported belief. We prove that this formulation is piecewise affine over the so-called chamber complex of the feasible set of the high-point relaxation. We propose two algorithmic approaches to solve general problems enjoying this last property. The first one is based on enumerating the vertices of the chamber complex. This approach is not scalable, but we present it as a computational baseline and for its theoretical interest. The second one is a Monte-Carlo approximation scheme based on the fact that randomly drawn points of the domain lie, with probability 1, in the interior of full-dimensional chambers, where the problem (restricted to this chamber) can be reduced to a linear program. Finally, we evaluate these methods through computational experiments showing both approaches' advantages and challenges.
△ Less
Submitted 22 May, 2024; v1 submitted 4 November, 2022;
originally announced November 2022.
-
Online Learning Models for Vehicle Usage Prediction During COVID-19
Authors:
Tobias Lindroth,
Axel Svensson,
Niklas Åkerblom,
Mitra Pourabdollah,
Morteza Haghir Chehreghani
Abstract:
Today, there is an ongoing transition to more sustainable transportation, for which an essential part is the switch from combustion engine vehicles to battery electric vehicles (BEVs). BEVs have many advantages from a sustainability perspective, but issues such as limited driving range and long recharge times slow down the transition from combustion engines. One way to mitigate these issues is by…
▽ More
Today, there is an ongoing transition to more sustainable transportation, for which an essential part is the switch from combustion engine vehicles to battery electric vehicles (BEVs). BEVs have many advantages from a sustainability perspective, but issues such as limited driving range and long recharge times slow down the transition from combustion engines. One way to mitigate these issues is by performing battery thermal preconditioning, which increases the energy efficiency of the battery. However, to optimally perform battery thermal preconditioning, the vehicle usage pattern needs to be known, i.e., how and when the vehicle will be used. This study attempts to predict the departure time and distance of the first drive each day using online machine learning models. The online machine learning models are trained and evaluated on historical driving data collected from a fleet of BEVs during the COVID-19 pandemic. Additionally, the prediction models are extended to quantify the uncertainty of their predictions, which can be used to decide whether the prediction should be used or dismissed. Based on our results, the best-performing prediction models yield an aggregated mean absolute error of 2.75 hours when predicting departure time and 13.37 km when predicting trip distance.
△ Less
Submitted 28 February, 2024; v1 submitted 28 October, 2022;
originally announced October 2022.
-
Natural Almost Hermitian Structures on Conformally Foliated 4-Dimensional Lie Groups with Minimal Leaves
Authors:
Emma Andersdotter Svensson
Abstract:
Let $(G,g)$ be a 4-dimensional Riemannian Lie group with a 2-dimensional left-invariant, conformal foliation $\mathcal{F}$ with minimal leaves. Let $J$ be an almost Hermitian structure on $G$ adapted to the foliation $\mathcal{F}$. The corresponding Lie algebra $\mathfrak{g}$ must then belong to one of 20 families $\mathfrak{g}_1,\dots,\mathfrak{g}_{20}$ according to S. Gudmundsson and M. Svensson…
▽ More
Let $(G,g)$ be a 4-dimensional Riemannian Lie group with a 2-dimensional left-invariant, conformal foliation $\mathcal{F}$ with minimal leaves. Let $J$ be an almost Hermitian structure on $G$ adapted to the foliation $\mathcal{F}$. The corresponding Lie algebra $\mathfrak{g}$ must then belong to one of 20 families $\mathfrak{g}_1,\dots,\mathfrak{g}_{20}$ according to S. Gudmundsson and M. Svensson. We classify such structures $J$ which are almost Kähler $(\mathcal{A}\mathcal{K})$, integrable $(\mathcal{I})$ or Kähler $(\mathcal{K})$. Hereby, we construct 16 multi-dimensional almost Kähler families, 18 integrable families and 11 Kähler families.
△ Less
Submitted 3 March, 2022;
originally announced March 2022.
-
Natural almost Hermitian structures on conformally foliated 4-dimensional Lie groups with minimal leaves
Authors:
Emma Andersdotter Svensson,
Sigmundur Gudmundsson
Abstract:
Let $(G,g)$ be a 4-dimensional Riemannian Lie group with a 2-dimensional left-invariant, conformal foliation $\F$ with minimal leaves. Let $J$ be an almost Hermitian structure on $G$ adapted to the foliation $\F$. We classify such structures $J$ which are almost Kähler $(\A\K)$, integrable $(\I)$ or Kähler $(\K)$. Hereby we construct several new multi-dimensional examples in each class.
Let $(G,g)$ be a 4-dimensional Riemannian Lie group with a 2-dimensional left-invariant, conformal foliation $\F$ with minimal leaves. Let $J$ be an almost Hermitian structure on $G$ adapted to the foliation $\F$. We classify such structures $J$ which are almost Kähler $(\A\K)$, integrable $(\I)$ or Kähler $(\K)$. Hereby we construct several new multi-dimensional examples in each class.
△ Less
Submitted 11 January, 2022;
originally announced January 2022.
-
Existence of solutions for deterministic bilevel games under a general Bayesian approach
Authors:
David Salas,
Anton Svensson
Abstract:
In 1996, Mallozzi and Morgan [33] proposed a new model for Stackelberg games which we refer here to as the Bayesian approach. The leader has only partial information about how followers select their reaction among possibly multiple optimal ones. This partial information is modeled as a decision-dependent distribution, the so-called belief of the leader. In this work, we formalize the setting of th…
▽ More
In 1996, Mallozzi and Morgan [33] proposed a new model for Stackelberg games which we refer here to as the Bayesian approach. The leader has only partial information about how followers select their reaction among possibly multiple optimal ones. This partial information is modeled as a decision-dependent distribution, the so-called belief of the leader. In this work, we formalize the setting of this approach for bilevel games admitting multiple leaders and we provide new results of existence of solutions. We pay particular attention to the fundamental case of linear bilevel problems, which has not been studied before, and which main difficulty is given by possible variations in the dimension of the reaction set of the follower. Our main technique to address this difficulty is based on a stronger notion of continuity for set-valued maps that we call rectangular continuity, and which is verified by the solution set of parametric linear problems. Finally, we provide some numerical experiments to address linear bilevel problems under the Bayesian approach.
△ Less
Submitted 11 May, 2023; v1 submitted 11 October, 2020;
originally announced October 2020.
-
Finite intersection property for bifunctions and existence for quasi-equilibrium problems
Authors:
John Cotrina,
Anton Svensson
Abstract:
The "finite intersection property" for bifunctions is introduced and its relationship with generalized monotonicity properties is studied. Some results concerning existence of solution for (quasi-)equilibrium problems are established and several results well-known in the literature are recovered. Furthermore, two applications are considered.
The "finite intersection property" for bifunctions is introduced and its relationship with generalized monotonicity properties is studied. Some results concerning existence of solution for (quasi-)equilibrium problems are established and several results well-known in the literature are recovered. Furthermore, two applications are considered.
△ Less
Submitted 11 February, 2020;
originally announced February 2020.
-
Coerciveness condition for quasi-equilibrium problems
Authors:
John Cotrina,
Abderrahim Hantoute,
Anton Svensson
Abstract:
A quasi-equilibrium problem is an equilibrium problem where the constraint set does depend on the reference point. It generalizes important problems such as quasi-variational inequalities and generalized Nash equilibrium problems. We study the existence of equilibria on unbounded sets under a coerciveness condition adapted from one specific for quasi-variational inequalities recently proposed by A…
▽ More
A quasi-equilibrium problem is an equilibrium problem where the constraint set does depend on the reference point. It generalizes important problems such as quasi-variational inequalities and generalized Nash equilibrium problems. We study the existence of equilibria on unbounded sets under a coerciveness condition adapted from one specific for quasi-variational inequalities recently proposed by Aussel and Sultana. We discuss the relation of our results with others that are present in the literature.
△ Less
Submitted 2 June, 2019; v1 submitted 25 January, 2019;
originally announced January 2019.
-
Data Consistency Approach to Model Validation
Authors:
Andreas Svensson,
Dave Zachariah,
Petre Stoica,
Thomas B. Schön
Abstract:
In scientific inference problems, the underlying statistical modeling assumptions have a crucial impact on the end results. There exist, however, only a few automatic means for validating these fundamental modelling assumptions. The contribution in this paper is a general criterion to evaluate the consistency of a set of statistical models with respect to observed data. This is achieved by automat…
▽ More
In scientific inference problems, the underlying statistical modeling assumptions have a crucial impact on the end results. There exist, however, only a few automatic means for validating these fundamental modelling assumptions. The contribution in this paper is a general criterion to evaluate the consistency of a set of statistical models with respect to observed data. This is achieved by automatically gauging the models' ability to generate data that is similar to the observed data. Importantly, the criterion follows from the model class itself and is therefore directly applicable to a broad range of inference problems with varying data types, ranging from independent univariate data to high-dimensional time-series. The proposed data consistency criterion is illustrated, evaluated and compared to several well-established methods using three synthetic and two real data sets.
△ Less
Submitted 20 May, 2019; v1 submitted 17 August, 2018;
originally announced August 2018.
-
On a stochastic model of epidemic spread with an application to competing infections
Authors:
Åke Svensson
Abstract:
A simple, but ``classical``, stochastic model for epidemic spread in a finite, but large, population is studied. The progress of the epidemic can be divided into three different phases that requires different tools to analyse. Initially the process is approximated by a branching process. It is discussed for how long time this approximation is valid. When a non-negligible proportion of the populati…
▽ More
A simple, but ``classical``, stochastic model for epidemic spread in a finite, but large, population is studied. The progress of the epidemic can be divided into three different phases that requires different tools to analyse. Initially the process is approximated by a branching process. It is discussed for how long time this approximation is valid. When a non-negligible proportion of the population is already infected the process can be studied using differential equations. In a final phase the spread will fade out. The results are used to investigate what happens if two strains of infectious agents, with different potential for spread, are simultaneously introduced in a totally susceptible population. It is assumed that an infection causes immunity, and that a person can only be infected by one strain. The two epidemics will initially develop approximately as independent branching processes. However, if both strains causes large epidemics they will, due to immunity, eventually interact. We will mainly be interested in the final outcome of the spread, i.e., how large proportion of the population is infected by the different strains.
△ Less
Submitted 27 May, 2018;
originally announced May 2018.
-
How consistent is my model with the data? Information-Theoretic Model Check
Authors:
Andreas Svensson,
Dave Zachariah,
Thomas B. Schön
Abstract:
The choice of model class is fundamental in statistical learning and system identification, no matter whether the class is derived from physical principles or is a generic black-box. We develop a method to evaluate the specified model class by assessing its capability of reproducing data that is similar to the observed data record. This model check is based on the information-theoretic properties…
▽ More
The choice of model class is fundamental in statistical learning and system identification, no matter whether the class is derived from physical principles or is a generic black-box. We develop a method to evaluate the specified model class by assessing its capability of reproducing data that is similar to the observed data record. This model check is based on the information-theoretic properties of models viewed as data generators and is applicable to e.g. sequential data and nonlinear dynamical models. The method can be understood as a specific two-sided posterior predictive test. We apply the information-theoretic model check to both synthetic and real data and compare it with a classical whiteness test.
△ Less
Submitted 19 December, 2017; v1 submitted 7 December, 2017;
originally announced December 2017.
-
Learning nonlinear state-space models using smooth particle-filter-based likelihood approximations
Authors:
Andreas Svensson,
Fredrik Lindsten,
Thomas B. Schön
Abstract:
When classical particle filtering algorithms are used for maximum likelihood parameter estimation in nonlinear state-space models, a key challenge is that estimates of the likelihood function and its derivatives are inherently noisy. The key idea in this paper is to run a particle filter based on a current parameter estimate, but then use the output from this particle filter to re-evaluate the lik…
▽ More
When classical particle filtering algorithms are used for maximum likelihood parameter estimation in nonlinear state-space models, a key challenge is that estimates of the likelihood function and its derivatives are inherently noisy. The key idea in this paper is to run a particle filter based on a current parameter estimate, but then use the output from this particle filter to re-evaluate the likelihood function approximation also for other parameter values. This results in a (local) deterministic approximation of the likelihood and any standard optimization routine can be applied to find the maximum of this local approximation. By iterating this procedure we eventually arrive at a final parameter estimate.
△ Less
Submitted 29 November, 2017;
originally announced November 2017.
-
A general representation of delta-normal sets to sublevels of convex functions
Authors:
Abderrahim Hantoute,
Anton Svensson
Abstract:
The (delta-) normal cone to an arbitrary intersection of sublevel sets of proper, lower semicontinuous, and convex functions is characterized, using either epsilon-subdifferentials at the nominal point or exact subdifferentials at nearby points. Our tools include (epsilon-) calculus rules for sup/max functions. The framework of this work is that of a locally convex space, however, formulas using e…
▽ More
The (delta-) normal cone to an arbitrary intersection of sublevel sets of proper, lower semicontinuous, and convex functions is characterized, using either epsilon-subdifferentials at the nominal point or exact subdifferentials at nearby points. Our tools include (epsilon-) calculus rules for sup/max functions. The framework of this work is that of a locally convex space, however, formulas using exact subdifferentials require some restriction either on the space (e.g. Banach), or on the function (e.g. epi-pointed).
△ Less
Submitted 27 October, 2017;
originally announced October 2017.
-
Probabilistic learning of nonlinear dynamical systems using sequential Monte Carlo
Authors:
Thomas B. Schön,
Andreas Svensson,
Lawrence Murray,
Fredrik Lindsten
Abstract:
Probabilistic modeling provides the capability to represent and manipulate uncertainty in data, models, predictions and decisions. We are concerned with the problem of learning probabilistic models of dynamical systems from measured data. Specifically, we consider learning of probabilistic nonlinear state-space models. There is no closed-form solution available for this problem, implying that we a…
▽ More
Probabilistic modeling provides the capability to represent and manipulate uncertainty in data, models, predictions and decisions. We are concerned with the problem of learning probabilistic models of dynamical systems from measured data. Specifically, we consider learning of probabilistic nonlinear state-space models. There is no closed-form solution available for this problem, implying that we are forced to use approximations. In this tutorial we will provide a self-contained introduction to one of the state-of-the-art methods---the particle Metropolis--Hastings algorithm---which has proven to offer a practical approximation. This is a Monte Carlo based method, where the particle filter is used to guide a Markov chain Monte Carlo method through the parameter space. One of the key merits of the particle Metropolis--Hastings algorithm is that it is guaranteed to converge to the "true solution" under mild assumptions, despite being based on a particle filter with only a finite number of particles. We will also provide a motivating numerical example illustrating the method using a modeling language tailored for sequential Monte Carlo methods. The intention of modeling languages of this kind is to open up the power of sophisticated Monte Carlo methods---including particle Metropolis--Hastings---to a large group of users without requiring them to know all the underlying mathematical details.
△ Less
Submitted 15 December, 2017; v1 submitted 7 March, 2017;
originally announced March 2017.
-
Learning of state-space models with highly informative observations: a tempered Sequential Monte Carlo solution
Authors:
Andreas Svensson,
Thomas B. Schön,
Fredrik Lindsten
Abstract:
Probabilistic (or Bayesian) modeling and learning offers interesting possibilities for systematic representation of uncertainty using probability theory. However, probabilistic learning often leads to computationally challenging problems. Some problems of this type that were previously intractable can now be solved on standard personal computers thanks to recent advances in Monte Carlo methods. In…
▽ More
Probabilistic (or Bayesian) modeling and learning offers interesting possibilities for systematic representation of uncertainty using probability theory. However, probabilistic learning often leads to computationally challenging problems. Some problems of this type that were previously intractable can now be solved on standard personal computers thanks to recent advances in Monte Carlo methods. In particular, for learning of unknown parameters in nonlinear state-space models, methods based on the particle filter (a Monte Carlo method) have proven very useful. A notoriously challenging problem, however, still occurs when the observations in the state-space model are highly informative, i.e. when there is very little or no measurement noise present, relative to the amount of process noise. The particle filter will then struggle in estimating one of the basic components for probabilistic learning, namely the likelihood $p($data$|$parameters$)$. To this end we suggest an algorithm which initially assumes that there is substantial amount of artificial measurement noise present. The variance of this noise is sequentially decreased in an adaptive fashion such that we, in the end, recover the original problem or possibly a very close approximation of it. The main component in our algorithm is a sequential Monte Carlo (SMC) sampler, which gives our proposed method a clear resemblance to the SMC^2 method. Another natural link is also made to the ideas underlying the approximate Bayesian computation (ABC). We illustrate it with numerical examples, and in particular show promising results for a challenging Wiener-Hammerstein benchmark problem.
△ Less
Submitted 13 December, 2017; v1 submitted 6 February, 2017;
originally announced February 2017.
-
A flexible state space model for learning nonlinear dynamical systems
Authors:
Andreas Svensson,
Thomas B. Schön
Abstract:
We consider a nonlinear state-space model with the state transition and observation functions expressed as basis function expansions. The coefficients in the basis function expansions are learned from data. Using a connection to Gaussian processes we also develop priors on the coefficients, for tuning the model flexibility and to prevent overfitting to data, akin to a Gaussian process state-space…
▽ More
We consider a nonlinear state-space model with the state transition and observation functions expressed as basis function expansions. The coefficients in the basis function expansions are learned from data. Using a connection to Gaussian processes we also develop priors on the coefficients, for tuning the model flexibility and to prevent overfitting to data, akin to a Gaussian process state-space model. The priors can alternatively be seen as a regularization, and helps the model in generalizing the data without sacrificing the richness offered by the basis function expansion. To learn the coefficients and other unknown parameters efficiently, we tailor an algorithm using state-of-the-art sequential Monte Carlo methods, which comes with theoretical guarantees on the learning. Our approach indicates promising results when evaluated on a classical benchmark as well as real data.
△ Less
Submitted 28 March, 2017; v1 submitted 17 March, 2016;
originally announced March 2016.
-
Nonlinear State Space Model Identification Using a Regularized Basis Function Expansion
Authors:
Andreas Svensson,
Thomas B. Schön,
Arno Solin,
Simo Särkkä
Abstract:
This paper is concerned with black-box identification of nonlinear state space models. By using a basis function expansion within the state space model, we obtain a flexible structure. The model is identified using an expectation maximization approach, where the states and the parameters are updated iteratively in such a way that a maximum likelihood estimate is obtained. We use recent particle me…
▽ More
This paper is concerned with black-box identification of nonlinear state space models. By using a basis function expansion within the state space model, we obtain a flexible structure. The model is identified using an expectation maximization approach, where the states and the parameters are updated iteratively in such a way that a maximum likelihood estimate is obtained. We use recent particle methods with sound theoretical properties to infer the states, whereas the model parameters can be updated using closed-form expressions by exploiting the fact that our model is linear in the parameters. Not to over-fit the flexible model to the data, we also propose a regularization scheme without increasing the computational burden. Importantly, this opens up for systematic use of regularization in nonlinear state space models. We conclude by evaluating our proposed approach on one simulation example and two real-data problems.
△ Less
Submitted 2 October, 2015;
originally announced October 2015.
-
Computationally Efficient Bayesian Learning of Gaussian Process State Space Models
Authors:
Andreas Svensson,
Arno Solin,
Simo Särkkä,
Thomas B. Schön
Abstract:
Gaussian processes allow for flexible specification of prior assumptions of unknown dynamics in state space models. We present a procedure for efficient Bayesian learning in Gaussian process state space models, where the representation is formed by projecting the problem onto a set of approximate eigenfunctions derived from the prior covariance structure. Learning under this family of models can b…
▽ More
Gaussian processes allow for flexible specification of prior assumptions of unknown dynamics in state space models. We present a procedure for efficient Bayesian learning in Gaussian process state space models, where the representation is formed by projecting the problem onto a set of approximate eigenfunctions derived from the prior covariance structure. Learning under this family of models can be conducted using a carefully crafted particle MCMC algorithm. This scheme is computationally efficient and yet allows for a fully Bayesian treatment of the problem. Compared to conventional system identification tools or existing learning methods, we show competitive performance and reliable quantification of uncertainties in the model.
△ Less
Submitted 15 April, 2016; v1 submitted 7 June, 2015;
originally announced June 2015.
-
Sequential Monte Carlo Methods for System Identification
Authors:
Thomas B. Schön,
Fredrik Lindsten,
Johan Dahlin,
Johan Wågberg,
Christian A. Naesseth,
Andreas Svensson,
Liang Dai
Abstract:
One of the key challenges in identifying nonlinear and possibly non-Gaussian state space models (SSMs) is the intractability of estimating the system state. Sequential Monte Carlo (SMC) methods, such as the particle filter (introduced more than two decades ago), provide numerical solutions to the nonlinear state estimation problems arising in SSMs. When combined with additional identification tech…
▽ More
One of the key challenges in identifying nonlinear and possibly non-Gaussian state space models (SSMs) is the intractability of estimating the system state. Sequential Monte Carlo (SMC) methods, such as the particle filter (introduced more than two decades ago), provide numerical solutions to the nonlinear state estimation problems arising in SSMs. When combined with additional identification techniques, these algorithms provide solid solutions to the nonlinear system identification problem. We describe two general strategies for creating such combinations and discuss why SMC is a natural tool for implementing these strategies.
△ Less
Submitted 10 March, 2016; v1 submitted 20 March, 2015;
originally announced March 2015.
-
Nonlinear state space smoothing using the conditional particle filter
Authors:
Andreas Svensson,
Thomas B. Schön,
Manon Kok
Abstract:
To estimate the smoothing distribution in a nonlinear state space model, we apply the conditional particle filter with ancestor sampling. This gives an iterative algorithm in a Markov chain Monte Carlo fashion, with asymptotic convergence results. The computational complexity is analyzed, and our proposed algorithm is successfully applied to the challenging problem of sensor fusion between ultra-w…
▽ More
To estimate the smoothing distribution in a nonlinear state space model, we apply the conditional particle filter with ancestor sampling. This gives an iterative algorithm in a Markov chain Monte Carlo fashion, with asymptotic convergence results. The computational complexity is analyzed, and our proposed algorithm is successfully applied to the challenging problem of sensor fusion between ultra-wideband and accelerometer/gyroscope measurements for indoor positioning. It appears to be a competitive alternative to existing nonlinear smoothing algorithms, in particular the forward filtering-backward simulation smoother.
△ Less
Submitted 16 September, 2015; v1 submitted 12 February, 2015;
originally announced February 2015.
-
Marginalizing Gaussian Process Hyperparameters using Sequential Monte Carlo
Authors:
Andreas Svensson,
Johan Dahlin,
Thomas B. Schön
Abstract:
Gaussian process regression is a popular method for non-parametric probabilistic modeling of functions. The Gaussian process prior is characterized by so-called hyperparameters, which often have a large influence on the posterior model and can be difficult to tune. This work provides a method for numerical marginalization of the hyperparameters, relying on the rigorous framework of sequential Mont…
▽ More
Gaussian process regression is a popular method for non-parametric probabilistic modeling of functions. The Gaussian process prior is characterized by so-called hyperparameters, which often have a large influence on the posterior model and can be difficult to tune. This work provides a method for numerical marginalization of the hyperparameters, relying on the rigorous framework of sequential Monte Carlo. Our method is well suited for online problems, and we demonstrate its ability to handle real-world problems with several dimensions and compare it to other marginalization methods. We also conclude that our proposed method is a competitive alternative to the commonly used point estimates maximizing the likelihood, both in terms of computational load and its ability to handle multimodal posteriors.
△ Less
Submitted 2 October, 2015; v1 submitted 6 February, 2015;
originally announced February 2015.
-
Identification of jump Markov linear models using particle filters
Authors:
Andreas Svensson,
Thomas B. Schön,
Fredrik Lindsten
Abstract:
Jump Markov linear models consists of a finite number of linear state space models and a discrete variable encoding the jumps (or switches) between the different linear models. Identifying jump Markov linear models makes for a challenging problem lacking an analytical solution. We derive a new expectation maximization (EM) type algorithm that produce maximum likelihood estimates of the model param…
▽ More
Jump Markov linear models consists of a finite number of linear state space models and a discrete variable encoding the jumps (or switches) between the different linear models. Identifying jump Markov linear models makes for a challenging problem lacking an analytical solution. We derive a new expectation maximization (EM) type algorithm that produce maximum likelihood estimates of the model parameters. Our development hinges upon recent progress in combining particle filters with Markov chain Monte Carlo methods in solving the nonlinear state smoothing problem inherent in the EM formulation. Key to our development is that we exploit a conditionally linear Gaussian substructure in the model, allowing for an efficient algorithm.
△ Less
Submitted 25 September, 2014;
originally announced September 2014.
-
The effect of time distribution shape on simulated epidemic models
Authors:
Martin Camitz,
Ake Svensson
Abstract:
By convention, and even more often, as an unintentional consequence of design, time distributions of latency and infectious durations in stochastic epidemic simulations are often exponential. The skewed distribtion typically leads to unrealistically short times. We examine the effects of altering the distribution latency and infectious times by comparing the key results after simulation with exp…
▽ More
By convention, and even more often, as an unintentional consequence of design, time distributions of latency and infectious durations in stochastic epidemic simulations are often exponential. The skewed distribtion typically leads to unrealistically short times. We examine the effects of altering the distribution latency and infectious times by comparing the key results after simulation with exponential and gamma distributions in a homogeneous mixing model aswell as a model with regional divisions connected by a travel intensity matrix. We show a delay in spread with more realistic latency times and offer an explanation of the effect.
△ Less
Submitted 18 December, 2007; v1 submitted 3 December, 2007;
originally announced December 2007.
-
Dynamics of Crystal Formation in the Greenland NorthGRIP Ice Core
Authors:
Joachim Mathiesen,
Jesper Ferkinghoff-Borg,
Mogens H. Jensen,
Mogens Levinsen,
Poul Olesen,
Dorthe Dahl-Jensen,
Anders Svensson
Abstract:
The North Greenland Ice Core Project (NorthGRIP) provides paleoclimatic information back to at about 120 kyr before present (Dahl-Jensen and others, 2002). Each year, precipitation on the ice sheet covers it with a new layer of snow, which gradually transforms into ice crystals as the layer sinks into the ice sheet. The size distribution of ice crystals has been measured at selected depths in th…
▽ More
The North Greenland Ice Core Project (NorthGRIP) provides paleoclimatic information back to at about 120 kyr before present (Dahl-Jensen and others, 2002). Each year, precipitation on the ice sheet covers it with a new layer of snow, which gradually transforms into ice crystals as the layer sinks into the ice sheet. The size distribution of ice crystals has been measured at selected depths in the upper 880 m of the NorthGRIP ice core (Svensson and others, 2003b), which covers a time span of 5300 years. The distributions change with time toward a universal curve, indicating a common underlying physical process in the formation of crystals. We identify this process as an interplay between fragmentation of the crystals and diffusion of their grain boundaries. The process is described by a two-parameter differential equation to which we obtain the exact solution. The solution is in excellent agreement with the experimentally observed distributions.
△ Less
Submitted 12 January, 2004; v1 submitted 17 October, 2003;
originally announced October 2003.