Refine
Year of publication
Language
- English (1103) (remove)
Keywords
- optimal control (27)
- stability (14)
- integer programming (11)
- Stochastic programming (9)
- finite elements (9)
- mixed integer programming (9)
- Hamiltonian matrix (8)
- finite element method (8)
- model reduction (8)
- state constraints (8)
Riemann surfaces naturally appear in the analysis of complex functions that are branched over the complex plane. However, they usually possess a complicated topology and are thus hard to understand. We present an algorithm for constructing Riemann surfaces as meshes in R3 from explicitly given branch points with corresponding branch indices. The constructed surfaces cover the complex plane by the canonical
projection onto R2 and can therefore be considered as multivalued graphs
over the plane – hence they provide a comprehensible visualization of the
topological structure. Complex functions are elegantly visualized using domain coloring on
a subset of C. By applying domain coloring to the automatically constructed Riemann surface models, we generalize this approach to deal with functions which cannot be entirely visualized in the complex plane.
Diffusion weighted imaging is a magnetic resonance based method to investigate
tissue micro-structure especially in the human brain via water diffusion.
Since the standard diffusion tensor model for the acquired data failes in
large portion of the brain voxel more sophisticated models have bee developed.
Here, we report on the package dti and how some of these models
can be used with the package.
The package fmri is provided for analysis of single run functional
Magnetic Resonance Imaging data. It implements structural adaptive smoothing
methods with signal detection for adaptive noise reduction which avoids blurring
of edges of activation areas. fmri provides fmri analysis from time series
modeling to signal detection and publication-ready images.
Modeling the orientation distribution function by mixtures of angular central Gaussian distributions
(2010)
In this paper we develop a tensor mixture model for diffusion weighted imaging
data using an automatic model selection criterion for the order of tensor
components in a voxel. We show that the weighted orientation distribution
function for this model can be expanded into a mixture of angular central
Gaussian distributions. We show properties of this model in extensive
simulations and in a high angular resolution experimental data set. The results
suggest that the model may improve imaging of cerebral fiber tracts. We
demonstrate how inference on canonical model parameters may give rise to new
clinical applications.
In this work, we study the spectra and eigenmodes of the Hessian of various discrete surface energies and discuss applications to shape analysis. In particular, we consider a physical model that describes the vibration modes and frequencies of a surface through the eigenfunctions and eigenvalues of the Hessian of a deformation energy, and we
derive a closed form representation for the Hessian (at the rest state of the energy) for a general class of deformation energies. Furthermore, we design a quadratic energy, such that the eigenmodes of the Hessian of
this energy are sensitive to the extrinsic curvature of the surface. Based on these spectra and eigenmodes, we derive two shape signatures. One that measures the similarity of points on a surface, and another that
can be used to identify features of the surface. In addition, we discuss a
spectral quadrangulation scheme for surfaces.
We discuss the perturbation analysis for
eigenvalues and eigenvectors of structured homogeneous matrix polynomials with
Hermitian, skew-Hermitian, H-even and H-odd structure.
We construct minimal structured perturbations (structured backward errors) such that an
approximate eigenpair is an exact eigenpair of an appropriate perturbed structured matrix
polynomial. We present various comparisons with unstructured backward
errors and previous error bounds derived for the non-homogeneous case
and show that our bounds present a significant improvement.
Many applications give rise to matrix polynomials whose coefficients have
a kind of reversal symmetry, a structure we call palindromic.
Several properties of scalar palindromic polynomials are derived,
and together with properties of compound matrices, used to
establish the Smith form of regular and singular T-palindromic matrix polynomials,
over arbitrary fields.
The invariant polynomials are shown to
inherit palindromicity,
and their structure is described in detail.
Jordan structures of palindromic matrix polynomials are characterized,
and necessary conditions for the
existence of structured linearizations established.
In the odd degree case, a constructive procedure for building
palindromic linearizations shows that the necessary conditions are sufficient as well.
The Smith form for *-palindromic polynomials is also analyzed. Finally, results for palindromic matrix polynomials over fields of
characteristic two are presented.
We define a risk averse nonanticipative feasible policy for multistage stochastic programs and propose a methodology to implement it. The approach is based on dynamic programming equations written for a risk averse formulation of the problem.
This formulation relies on a new class of multiperiod risk functionals called extended polyhedral risk measures. Dual representations of such risk functionals are given and used to derive conditions of coherence. In the one-period case, conditions for convexity and consistency with second order stochastic dominance are also provided. The risk averse dynamic programming equations are specialized considering convex combinations of one-period extended polyhedral risk measures such as spectral risk measures.
To implement the proposed policy, the approximation of the risk averse recourse functions for stochastic linear programs is discussed. In this context, we detail a stochastic dual dynamic programming algorithm which converges to the optimal value of the risk averse problem.
We derive a Crooks-Jarzynski-type identity for computing free energy differences between metastable states that is based on nonequilibrium diffusion processes. Furthermore we outline a brief derivation of an infinite-dimensional stochastic partial differential equation that can be used to efficiently generate the ensemble of trajectories connecting the metastable states.
To address the plurality of interpretations of the subjective notion of risk, we describe it by means of a risk order and concentrate on the context invariant features of diversification and monotonicity. Our main results are uniquely characterized robust representations of lower semicontinuous risk orders on vector spaces and convex sets. This representation covers most instruments related to risk and allow for a differentiated interpretation depending on the underlying context which is illustrated in different settings: For random variables, risk perception can be interpreted as model risk, and we compute among others the robust representation of the economic index of riskiness. For lotteries, risk perception can be viewed as distributional risk and we study the "Value at Risk". For consumption patterns, which excerpt an intertemporality dimension in risk perception, we provide an interpretation in terms of discounting risk and discuss some examples.
We study the risk assessment of uncertain cash flows in terms of dynamic convex risk measures for processes as introduced in Cheridito, Delbaen, and Kupper (2006). These risk measures take into account not only the amounts but also the timing of a cash flow. We discuss their robust representation in terms of suitably penalized probability measures on the optional $\sigma$-field. This yields an explicit analysis both of model and discounting ambiguity. We focus on supermartingale criteria for time consistency. In particular we show how ``bubbles'' may appear in the dynamic penalization, and how they cause a breakdown of asymptotic safety of the risk assessment procedure.
Dynamic risk measures
(2010)
This paper gives an overview of the theory of dynamic convex risk measures for random variables in discrete time setting. We summarize robust representation results of conditional convex risk measures, and we characterize various time consistency properties of dynamic risk measures in terms of acceptance sets, penalty functions, and by supermartingale properties of risk processes and penalty functions.
The classical valuation of an uncertain cash flow in discrete time consists in taking the expectation of the sum of the discounted future payoffs under a fixed probability measure, which is assumed to be known. Here we discuss the valuation problem in the context of Knightian uncertainty. Using results from the theory of convex risk measures, but without assuming the existence of a global reference measure, we derive a robust representation of concave valuations with an infinite time horizon, which specifies the interplay between model uncertainty and uncertainty about the time value of money.
Recently, there is a growing trend to offer guarantee products where the investor is allowed to shift her account/investment value between multiple funds. The switching right is granted a finite number per year, i.e. it is American style with multiple exercise possibilities. In consequence, the pricing and the risk management is based on the switching strategy which maximizes the value of the guarantee put option. We analyze the optimal stopping problem in the case of one switching right within different model classes and compare the exact price with the lower price bound implied by the optimal deterministic switching time. We show that, within the class of log-price processes with independent increments, the stopping problem is solved by a deterministic stopping time if (and only if) the price process is in addition continuous. Thus, in a sense, the Black & Scholes model is the only (meaningful) pricing model where the lower price bound gives the exact price. It turns out that even moderate deviations from the Black & Scholes model assumptions give a lower price bound which is really below the exact price. This is illustrated by means of a stylized stochastic volatility model setup.
In this paper we consider the optimal stopping problem for general dynamic monetary utility functionals. Sufficient conditions for the Bellman principle and the existence of optimal stopping times are provided. Particular attention is payed to representations which allow for a numerical treatment in real situations. To this aim, generalizations of standard evaluation methods like policy iteration, dual and consumption based approaches are developed in the context of general dynamic monetary utility functionals. As a result, it turns out that the possibility of a particular generalization depends on specific properties of the utility functional under consideration.
Boolean modeling frameworks have long since proved their worth for capturing and analyzing essential characteristics of complex systems.
Hybrid approaches aim at exploiting the advantages of Boolean formalisms while refining expressiveness. In this paper, we present a formalism that augments Boolean models with stochastic aspects. More specifically, biological reactions effecting a system in a given state are associated
with probabilities, resulting in dynamical behavior represented as a Markov chain. Using this approach, we model and analyze the cytokinin
response network of Arabidopsis thaliana with a focus on clarifying the character of an important feedback mechanism.
Bovine fertility is the subject of extensive research in animal sciences,
especially because fertility of dairy cows has declined during the last
decades. The regulation of estrus is controlled by the complex interplay
of various organs and hormones. Mathematical modeling of the bovine
estrous cycle could help in understanding the dynamics of this complex
biological system. In this paper we present a mechanistic mathematical
model of the bovine estrous cycle that includes the processes of follicle
and corpus luteum development and the key hormones that interact to
control these processes. The model generates successive estrous cycles of
21 days, with three waves of follicle growth per cycle. The model contains
12 differential equations and 54 parameters. Focus in this paper is on
development of the model, but also some simulation results are presented,
showing that a set of equations and parameters is obtained that describes
the system consistent with empirical knowledge. Even though the majority
of the mechanisms that are included in the model are based on relations
that in literature have only been described qualitatively (i.e. stimulation
and inhibition), the output of the model is surprisingly well in line with
empirical data. This model of the bovine estrous cycle could be used
as a basis for more elaborate models with the ability to study effects of
external manipulations and genetic differences.
We consider simple models of financial markets with less and better
informed investors described by a smaller and a larger filtration on a
general stochastic basis that describes the market dynamics, including
continuous and jump components. We study the relation between different forms of non existance of arbitrage and the characteristics of the stochastic basis under the different filtrations. This is achieved through the analysis of the properties of the numéraire portfolio. Furthermore, we focus on the problem of calculating the additional logarithmic utility of the better informed investor in terms of the Shannon antropy of is additional information. The information drift, i.e. the drift to eliminate in order to preserved the martingale property in the larger filtration terms out to be the crucial quantity needed to tackle these problems. We show that the expected
ed logarithmic utility increment due to better information equals its Shannon
entropy also in case of a pure jump basis with jumps that are quadratically
hedgeable, and so extend a similar result known for bases consisting of
continuous semimartingales. An example illustrates that the equality may
not persist if both continuous and jump components are present in the
underlying.