BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Raphael Hauser

Raphael Hauser appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
math.OCmath.PRq-fin.PMq-fin.PRComputational Engineering, Finance, and ScienceComputer Visionmath.STnlin.AOphysics.med-phq-fin.CPq-fin.RMStatistics Theory

Catalog footprint

What is connected

17works
12topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2016arXiv

3D Image Reconstruction from X-Ray Measurements with Overlap

3D image reconstruction from a set of X-ray projections is an important image reconstruction problem, with applications in medical imaging, industrial inspection and airport security. The innovation of X-ray emitter arrays allows for a novel type of X-ray scanners with multiple simultaneously emitting sources. However, two or more sources emitting at the same time can yield measurements from overlapping rays, imposing a new type of image reconstruction problem based on nonlinear constraints. Using traditional linear reconstruction methods, respective scanner geometries have to be implemented such that no rays overlap, which severely restricts the scanner design. We derive a new type of 3D image reconstruction model with nonlinear constraints, based on measurements with overlapping X-rays. Further, we show that the arising optimization problem is partially convex, and present an algorithm to solve it. Experiments show highly improved image reconstruction results from both simulated and real-world measurements.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A New Approach to Model Free Option Pricing

In this paper we introduce a new approach to model-free path-dependent option pricing. We first introduce a general duality result for linear optimisation problems over signed measures introduced in [3] and show how the the problem of model-free option pricing can be formulated in the new framework. We then introduce a model to solve the problem numerically when the only information provided is the market data of vanilla call or put option prices. Compared to the common approaches in the literature, e.g. [4], the model does not require the marginal distributions of the stock price for different maturities. Though the experiments are carried out for simple path-dependent options on a single stock, the model is easy to generalise for multi-asset framework.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Strong Duality of Linear Optimisation Problems over Measure Spaces

In this work we present two particular cases of the general duality result for linear optimisation problems over signed measures with infinitely many constraints in the form of integrals of functions with respect to the decision variables (the measure in question) for which strong duality holds. In the first case the optimisation problems are over measures with $L^p$ density functions with $1 < p < \infty$. In the second case we consider a semi-infinite optimisation problem where finitely many constraints are given in form of bounds on integrals. The latter case has a particular importance in practice where the model can be applied in robust risk management and model-free option pricing.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A General Duality Relation with Applications in Quantitative Risk Management

A fundamental problem in risk management is the robust aggregation of different sources of risk in a situation where little or no data are available to infer information about their dependencies. A popular approach to solving this problem is to formulate an optimization problem under which one maximizes a risk measure over all multivariate distributions that are consistent with the available data. In several special cases of such models, there exist dual problems that are easier to solve or approximate, yielding robust bounds on the aggregated risk. In this chapter we formulate a general optimization problem, which can be seen as a doubly infinite linear programming problem, and we show that the associated dual generalizes several well known special cases and extends to new risk management models we propose.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

An Upper Bound on the Convergence Rate of a Second Functional in Optimal Sequence Alignment

Consider finite sequences $X_{[1,n]}=X_1\dots X_n$ and $Y_{[1,n]}=Y_1\dots Y_n$ of length $n$, consisting of i.i.d.\ samples of random letters from a finite alphabet, and let $S$ and $T$ be chosen i.i.d.\ randomly from the unit ball in the space of symmetric scoring functions over this alphabet augmented by a gap symbol. We prove a probabilistic upper bound of linear order in $n^{0.75}$ for the deviation of the score relative to $T$ of optimal alignments with gaps of $X_{[1,n]}$ and $Y_{[1,n]}$ relative to $S$. It remains an open problem to prove a lower bound. Our result contributes to the understanding of the microstructure of optimal alignments relative to one given scoring function, extending a theory begun by the first two authors.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Calculation of a power price equilibrium

In this paper we propose a tractable quadratic programming formulation for calculating the equilibrium term structure of electricity prices. We rely on a theoretical model described in [21], but extend it so that it reflects actually traded electricity contracts, transaction costs and liquidity considerations. Our numerical simulations examine the properties of the term structure and its dependence on various parameters of the model. The proposed quadratic programming formulation is applied to calculate the equilibrium term structure of electricity prices in the UK power grid consisting of a few hundred power plants. The impact of ramp up and ramp down constraints are also studied.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The existence and uniqueness of a power price equilibrium

We propose a term structure power price model that, in contrast to widely accepted no-arbitrage based approaches, accounts for the non-storable nature of power. It belongs to a class of equilibrium game theoretic models with players divided into producers and consumers. The consumers' goal is to maximize a mean-variance utility function subject to satisfying an inelastic demand of their own clients (e.g households, businesses etc.) to whom they sell the power. The producers, who own a portfolio of power plants each defined by a running fuel (e.g. gas, coal, oil...) and physical characteristics (e.g. efficiency, capacity, ramp up/down times...), similarly, seek to maximize a mean-variance utility function consisting of power, fuel, and emission prices subject to production constraints. Our goal is to determine the term structure of the power price at which production matches consumption. In this paper we show that in such a setting the equilibrium price exists and also discuss the conditions for its uniqueness.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The impact of startup costs and the grid operator on the power price equilibrium

In this paper we propose a quadratic programming model that can be used for calculating the term structure of electricity prices while explicitly modeling startup costs of power plants. In contrast to other approaches presented in the literature, we incorporate the startup costs in a mathematically rigorous manner without relying on ad hoc heuristics. Moreover, we propose a tractable approach for estimating the startup costs of power plants based on their historical production. Through numerical simulations applied to the entire UK power grid, we demonstrate that the inclusion of startup costs is necessary for the modeling of electricity prices in realistic power systems. Numerical results show that startup costs make electricity prices very spiky. In the second part of the paper, we extend the initial model by including the grid operator who is responsible for managing the grid. Numerical simulations demonstrate that robust decision making of the grid operator can significantly decrease the number and severity of spikes in the electricity price and improve the reliability of the power grid.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Letter Change Bias and Local Uniqueness in Optimal Sequence Alignments

Considering two optimally aligned random sequences, we investigate the effect on the alignment score caused by changing a random letter in one of the two sequences. Using this idea in conjunction with large deviations theory, we show that in alignments with a low proportion of gaps the optimal alignment is locally unique in most places with high probability. This has implications in the design of recently pioneered alignment methods that use the local uniqueness as a homology indicator.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Regression techniques for Portfolio Optimisation using MOSEK

Regression is widely used by practioners across many disciplines. We reformulate the underlying optimisation problem as a second-order conic program providing the flexibility often needed in applications. Using examples from portfolio management and quantitative trading we solve regression problems with and without constraints. Several Python code fragments are given. The code and data are available online at http://www.github.com/tschm/MosekRegression.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Seven Sins in Portfolio Optimization

Although modern portfolio theory has been in existence for over 60 years, fund managers often struggle to get its models to produce reliable portfolio allocations without strongly constraining the decision vector by tight bands of strategic allocation targets. The two main root causes to this problem are inadequate parameter estimation and numerical artifacts. When both obstacles are overcome, portfolio models yield excellent allocations. In this paper, which is primarily aimed at practitioners, we discuss the most common mistakes in setting up portfolio models and in solving them algorithmically.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The S-Procedure via Dual Cone Calculus

Given a quadratic function $h$ that satisfies a Slater condition, Yakubovich's S-Procedure (or S-Lemma) gives a characterization of all other quadratic functions that are copositive with $h$ in a form that is amenable to numerical computations. In this paper we present a deep-rooted connection between the S-Procedure and the dual cone calculus formula $(K_1\cap K_2)^*= K_1^*+K_2^*$, which holds for closed convex cones in $\R^2$. To establish the link with the S-Procedure, we generalize the dual cone calculus formula to a situation where $K_1$ is nonclosed, nonconvex and nonconic but exhibits sufficient mathematical resemblance to a closed convex cone. As a result, we obtain a new proof of the S-Lemma and an extension to Hilbert space kernels.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

A Monte Carlo Approach to the Fluctuation Problem in Optimal Alignments of Random Strings

The problem of determining the correct order of fluctuation of the optimal alignment score of two random strings of length $n$ has been open for several decades. It is known that the biased expected effect of a random letter-change on the optimal score implies an order of fluctuation linear in $\sqrt{n}$. However, in many situations where such a biased effect is observed empirically, it has been impossible to prove analytically. The main result of this paper shows that when the rescaled-limit of the optimal alignment score increases in a certain direction, then the biased effect exists. On the basis of this result one can quantify a confidence level for the existence of such a biased effect and hence of an order $\sqrt{n}$ fluctuation based on simulation of optimal alignments scores.This is an important step forward, as the correct order of fluctuation was previously known only for certain special distributions. To illustrate the usefulness of our new methodology, we apply it to optimal alignments of strings written in the DNA-alphabet. As scoring function, we use the BLASTZ default-substitution matrix together with a realistic gap penalty. BLASTZ is one of the most widely used sequence alignment methodologies in bioinformatics. For this DNA-setting, we show that with a high level of confidence, the fluctuation of the optimal alignment score is of order $Θ(\sqrt{n})$. An important special case of optimal alignment score is the Longest Common Subsequence (LCS) of random strings. For binary sequences with equiprobable symbols, the question of the fluctuation of the LCS remains open. The symmetry in that case does not allow for our method. On the other hand, in real-life DNA sequences, it is not the case that all letters occur with the same frequency. Thus, for many real life situations, our method allows to determine the order of the fluctuation up to a high confidence level.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Distribution of Aligned Letter Pairs in Optimal Alignments of Random Sequences

Considering the optimal alignment of two i.i.d. random sequences of length $n$, we show that when the scoring function is chosen randomly, almost surely the empirical distribution of aligned letter pairs in all optimal alignments converges to a unique limiting distribution as $n$ tends to infinity. This result is interesting because it helps understanding the microscopic path structure of a special type of last passage percolation problem with correlated weights, an area of long-standing open problems. Characterizing the microscopic path structure yields furthermore a robust alternative to optimal alignment scores for testing the relatedness of genetic sequences.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Inferring the Composition of a Trader Population in a Financial Market

We discuss a method for predicting financial movements and finding pockets of predictability in the price-series, which is built around inferring the heterogeneity of trading strategies in a multi-agent trader population. This work explores extensions to our previous framework (arXiv:physics/0506134). Here we allow for more intelligent agents possessing a richer strategy set, and we no longer constrain the estimate for the heterogeneity of the agents to a probability space. We also introduce a scheme which allows the incorporation of models with a wide variety of agent types, and discuss a mechanism for the removal of bias from relevant parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2001arXiv

Self-scaled barrier functions on symmetric cones and their classification

Self-scaled barrier functions on self-scaled cones were introduced through a set of axioms in 1994 by Y.E. Nesterov and M.J. Todd as a tool for the construction of long-step interior point algorithms. This paper provides firm foundation for these objects by exhibiting their symmetry properties, their intimate ties with the symmetry groups of their domains of definition, and subsequently their decomposition into irreducible parts and algebraic classification theory. In a first part we recall the characterisation of the family of self-scaled cones as the set of symmetric cones and develop a primal-dual symmetric viewpoint on self-scaled barriers, results that were first discovered by the second author. We then show in a short, simple proof that any pointed, convex cone decomposes into a direct sum of irreducible components in a unique way, a result which can also be of independent interest. We then show that any self-scaled barrier function decomposes in an essentially unique way into a direct sum of self-scaled barriers defined on the irreducible components of the underlying symmetric cone. Finally, we present a complete algebraic classification of self-scaled barrier functions using the correspondence between symmetric cones and Euclidean Jordan algebras.

0 citations0 reviewsNo code linkedOpen access