-
Increased accuracy of planning tools for optimization of dynamic multileaf collimator delivery of radiotherapy through reformulated objective functions
Authors:
Lovisa Engberg,
Kjell Eriksson,
Anders Forsgren
Abstract:
The purpose of this study is to examine in a clinical setting a novel formulation of objective functions for intensity-modulated radiotherapy treatment plan multicriteria optimization (MCO) that we suggested in a recent study. The proposed objective functions are extended with dynamic multileaf collimator (DMLC) delivery constraints from the literature, and a tailored interior point method is desc…
▽ More
The purpose of this study is to examine in a clinical setting a novel formulation of objective functions for intensity-modulated radiotherapy treatment plan multicriteria optimization (MCO) that we suggested in a recent study. The proposed objective functions are extended with dynamic multileaf collimator (DMLC) delivery constraints from the literature, and a tailored interior point method is described to efficiently solve the resulting optimization formulation. In a numerical planning study involving three patient cases, DMLC plans Pareto optimal to the MCO formulation with the proposed objective functions are generated. Evaluated based on pre-defined plan quality indices, these DMLC plans are compared to conventionally generated DMLC plans. Comparable or superior plan quality is observed. Supported by these results, the proposed objective functions are argued to have a potential to streamline the planning process, since they are designed to overcome the methodological shortcomings associated with the conventional penalty-based objective functions assumed to cause the current need for time-consuming trial-and-error parameter tuning. In particular, the increased accuracy of the planning tools imposed by the proposed objective functions has the potential to make the planning process less complicated. These conclusions position the proposed formulation as an alternative to existing methods for automated planning.
△ Less
Submitted 23 October, 2018; v1 submitted 2 February, 2018;
originally announced February 2018.
-
An exponential limit shape of random $q$-proportion Bulgarian solitaire
Authors:
Kimmo Eriksson,
Markus Jonsson abd Jonas Sjöstrand
Abstract:
We introduce \emph{$p_n$-random $q_n$-proportion Bulgarian solitaire} ($0<p_n,q_n\le 1$), played on $n$ cards distributed in piles. In each pile, a number of cards equal to the proportion $q_n$ of the pile size rounded upward to the nearest integer are candidates to be picked. Each candidate card is picked with probability $p_n$, independently of other candidate cards. This generalizes Popov's ran…
▽ More
We introduce \emph{$p_n$-random $q_n$-proportion Bulgarian solitaire} ($0<p_n,q_n\le 1$), played on $n$ cards distributed in piles. In each pile, a number of cards equal to the proportion $q_n$ of the pile size rounded upward to the nearest integer are candidates to be picked. Each candidate card is picked with probability $p_n$, independently of other candidate cards. This generalizes Popov's random Bulgarian solitaire, in which there is a single candidate card in each pile. Popov showed that a triangular limit shape is obtained for a fixed $p$ as $n$ tends to infinity. Here we let both $p_n$ and $q_n$ vary with $n$. We show that under the conditions $q_n^2 p_n n/{\log n}\rightarrow \infty$ and $p_n q_n \rightarrow 0$ as $n\to\infty$, the $p_n$-random $q_n$-proportion Bulgarian solitaire has an exponential limit shape.
△ Less
Submitted 21 March, 2017;
originally announced March 2017.
-
Limit shapes of stable configurations of a generalized Bulgarian solitaire
Authors:
Kimmo Eriksson,
Markus Jonsson,
Jonas Sjöstrand
Abstract:
Bulgarian solitaire is played on $n$ cards divided into several piles; a move consists of picking one card from each pile to form a new pile. In a recent generalization, $σ$-Bulgarian solitaire, the number of cards you pick from a pile is some function $σ$ of the pile size, such that you pick $σ(h)\le h$ cards from a pile of size $h$. Here we consider a special class of such functions. Let us call…
▽ More
Bulgarian solitaire is played on $n$ cards divided into several piles; a move consists of picking one card from each pile to form a new pile. In a recent generalization, $σ$-Bulgarian solitaire, the number of cards you pick from a pile is some function $σ$ of the pile size, such that you pick $σ(h)\le h$ cards from a pile of size $h$. Here we consider a special class of such functions. Let us call $σ$ well-behaved if $σ(1)=1$ and if both $σ(h)$ and $h-σ(h)$ are non-decreasing functions of $h$. Well-behaved $σ$-Bulgarian solitaire has a geometric interpretation in terms of layers at certain levels being picked in each move. It also satisfies that if a stable configuration of $n$ cards exists it is unique. Moreover, if piles are sorted in order of decreasing size ($λ_1 \ge λ_2\ge \dots$) then a configuration is convex if and only if it is a stable configuration of some well-behaved $σ$-Bulgarian solitaire. If sorted configurations are represented by Young diagrams and scaled down to have unit height and unit area, the stable configurations corresponding to an infinite sequence of well-behaved functions ($σ_1, σ_2, \dots$) may tend to a limit shape $φ$. We show that every convex $φ$ with certain properties can arise as the limit shape of some sequence of well-behaved $σ_n$. For the special case when $σ_n(h)=\lceil q_n h \rceil$ for $0 < q_n \le 1$, these limit shapes are triangular (in case $q_n^2 n\rightarrow 0$), or exponential (in case $q_n^2 n\rightarrow \infty$), or interpolating between these shapes (in case $q_n^2 n\rightarrow C>0$).
△ Less
Submitted 21 March, 2017;
originally announced March 2017.
-
Towards Robust Adaptive Radiation Therapy Strategies
Authors:
Michelle Böck,
Anders Forsgren,
Kjell Eriksson,
Björn Hårdemark
Abstract:
The authors propose robust adaptive strategies based on stochastic minimax optimization for a series of simulated treatments on a one-dimensional patient phantom. The plan applied during the first fractions should be able to handle anticipated systematic and random errors. At scheduled fractions, the impact of the measured errors on the delivered dose distribution is evaluated. For a patient recei…
▽ More
The authors propose robust adaptive strategies based on stochastic minimax optimization for a series of simulated treatments on a one-dimensional patient phantom. The plan applied during the first fractions should be able to handle anticipated systematic and random errors. At scheduled fractions, the impact of the measured errors on the delivered dose distribution is evaluated. For a patient receiving a dose that does not satisfy clinical goals, the plan is reoptimized based on these individually measured errors. The adapted plan is then applied during subsequent fractions until a new scheduled adaptation becomes necessary. In the first adaptive strategy, the measured systematic and random error scenarios and their assigned probabilities are updated to guide the robust reoptimization. In the second strategy, the grade of conservativeness is adapted in response to the measured dose delivery errors. In the third strategy, the uncertainty margins around the target are recalculated based on the measured errors. The simulated systematic and random errors are either similar to the anticipated errors or unpredictably larger, in order to critically evaluate the performance of the adaptive strategies. The simulations show that robustly optimized treatment plans provide sufficient treatment quality for those treatment error scenarios similar to the anticipated error scenarios. Combining robust planning with adaptation leads to improved OAR protection. In case of unpredictably larger treatment errors, the first strategy performs best at improving treatment quality in terms of target coverage and OAR protection. According to the simulations, these robust adaptive treatment strategies are able to identify candidates for an adaptive treatment, thus giving the opportunity to provide individualized plans, and improve their treatment quality through adaptation.
△ Less
Submitted 4 November, 2016;
originally announced November 2016.
-
Explicit optimization of plan quality measures in intensity-modulated radiation therapy treatment planning
Authors:
Lovisa Engberg,
Kjell Eriksson,
Anders Forsgren,
Björn Hårdemark
Abstract:
Conventional planning objectives in optimization of intensity-modulated radiotherapy treatment (IMRT) plans are designed to minimize the violation of dose-volume histogram (DVH) thresholds using penalty functions. Although successful in guiding the DVH curve towards these thresholds, conventional planning objectives offer limited control of the individual points on the DVH curve (doses-at-volume)…
▽ More
Conventional planning objectives in optimization of intensity-modulated radiotherapy treatment (IMRT) plans are designed to minimize the violation of dose-volume histogram (DVH) thresholds using penalty functions. Although successful in guiding the DVH curve towards these thresholds, conventional planning objectives offer limited control of the individual points on the DVH curve (doses-at-volume) used to evaluate plan quality. In this study, we abandon the usual penalty-function framework and propose planning objectives that more explicitly relate to DVH statistics. The proposed planning objectives are based on mean-tail-dose, resulting in convex optimization. We also demonstrate how to adapt a standard optimization method to the proposed formulation in order to obtain a substantial reduction in computational cost.
We investigate the potential of the proposed planning objectives as tools for optimizing DVH statistics through juxtaposition with the conventional planning objectives on two patient cases. Sets of treatment plans with differently balanced planning objectives are generated using either the proposed or the conventional approach. Dominance in the sense of better distributed doses-at-volume is observed in plans optimized within the proposed framework, indicating that the DVH statistics are better optimized and more efficiently balanced using the proposed planning objectives.
△ Less
Submitted 25 January, 2017; v1 submitted 5 July, 2016;
originally announced July 2016.
-
Markov chains on graded posets: Compatibility of up-directed and down-directed transition probabilities
Authors:
Kimmo Eriksson,
Markus Jonsson,
Jonas Sjöstrand
Abstract:
We consider two types of discrete-time Markov chains where the state space is a graded poset and the transitions are taken along the covering relations in the poset. The first type of Markov chain goes only in one direction, either up or down in the poset (an \emph{up chain} or \emph{down chain}). The second type toggles between two adjacent rank levels (an \emph{up-and-down chain}).
We introduc…
▽ More
We consider two types of discrete-time Markov chains where the state space is a graded poset and the transitions are taken along the covering relations in the poset. The first type of Markov chain goes only in one direction, either up or down in the poset (an \emph{up chain} or \emph{down chain}). The second type toggles between two adjacent rank levels (an \emph{up-and-down chain}).
We introduce two compatibility concepts between the up-directed transition probabilities (an \emph{up rule}) and the down-directed (a \emph{down rule}), and we relate these to compatibility between up-and-down chains. This framework is used to prove a conjecture about a limit shape for a process on Young's lattice.
Finally, we settle the questions whether the reverse of an up chain is a down chain for some down rule and whether there exists an up or down chain at all if the rank function is not bounded.
△ Less
Submitted 11 December, 2015;
originally announced December 2015.
-
The Limit Shape of a Stochastic Bulgarian Solitaire
Authors:
Kimmo Eriksson,
Markus Jonsson,
Jonas Sjöstrand
Abstract:
We consider a stochastic version of Bulgarian solitaire: A number of cards are distributed in piles; in every round a new pile is formed by cards from the old piles, and each card is picked independently with a fixed probability. This game corresponds to a multi-square birth-and-death process on Young diagrams of integer partitions. We prove that this process converges in a strong sense to an expo…
▽ More
We consider a stochastic version of Bulgarian solitaire: A number of cards are distributed in piles; in every round a new pile is formed by cards from the old piles, and each card is picked independently with a fixed probability. This game corresponds to a multi-square birth-and-death process on Young diagrams of integer partitions. We prove that this process converges in a strong sense to an exponential limit shape as the number of cards tends to infinity. Furthermore, we bound the probability of deviation from the limit shape and relate this to the number of rounds played in the solitaire.
△ Less
Submitted 11 December, 2015; v1 submitted 10 September, 2013;
originally announced September 2013.
-
Conjugacy of Coxeter elements
Authors:
Henrik Eriksson,
Kimmo Eriksson
Abstract:
For a Coxeter group (W,S), a permutation of the set S is called a Coxeter word and the group element represented by the product is called a Coxeter element. Moving the first letter to the end of the word is called a rotation and two Coxeter elements are rotation equivalent if their words can be transformed into each other through a sequence of rotations and legal commutations.
We prove that Coxe…
▽ More
For a Coxeter group (W,S), a permutation of the set S is called a Coxeter word and the group element represented by the product is called a Coxeter element. Moving the first letter to the end of the word is called a rotation and two Coxeter elements are rotation equivalent if their words can be transformed into each other through a sequence of rotations and legal commutations.
We prove that Coxeter elements are conjugate if and only if they are rotation equivalent. This was known for some special cases but not for Coxeter groups in general.
△ Less
Submitted 12 February, 2013;
originally announced February 2013.
-
Explicit Time-Stepping for Stiff ODEs
Authors:
Kenneth Eriksson,
Claes Johnson,
Anders Logg
Abstract:
We present a new strategy for solving stiff ODEs with explicit methods. By adaptively taking a small number of stabilizing small explicit time steps when necessary, a stiff ODE system can be stabilized enough to allow for time steps much larger than what is indicated by classical stability analysis. For many stiff problems the cost of the stabilizing small time steps is small, so the improvement i…
▽ More
We present a new strategy for solving stiff ODEs with explicit methods. By adaptively taking a small number of stabilizing small explicit time steps when necessary, a stiff ODE system can be stabilized enough to allow for time steps much larger than what is indicated by classical stability analysis. For many stiff problems the cost of the stabilizing small time steps is small, so the improvement is large. We illustrate the technique on a number of well-known stiff test problems.
△ Less
Submitted 12 May, 2012;
originally announced May 2012.
-
Words with intervening neighbours in infinite Coxeter groups are reduced
Authors:
Henrik Eriksson,
Kimmo Eriksson
Abstract:
Consider a graph with vertex set S. A word in the alphabet S has the intervening neighbours property if any two occurrences of the same letter are separated by all its graph neighbours. For a Coxeter graph, words represent group elements. Speyer recently proved that words with the intervening neighbours property are irreducible if the group is infinite and irreducible. We present a new and short…
▽ More
Consider a graph with vertex set S. A word in the alphabet S has the intervening neighbours property if any two occurrences of the same letter are separated by all its graph neighbours. For a Coxeter graph, words represent group elements. Speyer recently proved that words with the intervening neighbours property are irreducible if the group is infinite and irreducible. We present a new and shorter proof using the root automaton for recognition of irreducible words.
△ Less
Submitted 26 November, 2008;
originally announced November 2008.
-
The numbers game and Dynkin diagram classification results
Authors:
Robert G. Donnelly,
Kimmo Eriksson
Abstract:
The numbers game is a one-player game played on a finite simple graph with certain "amplitudes" assigned to its edges and with an initial assignment of real numbers to its nodes. The moves of the game successively transform the numbers at the nodes using the amplitudes in a certain way. Combinatorial reasoning is used to show that those connected graphs with negative integer amplitudes for which…
▽ More
The numbers game is a one-player game played on a finite simple graph with certain "amplitudes" assigned to its edges and with an initial assignment of real numbers to its nodes. The moves of the game successively transform the numbers at the nodes using the amplitudes in a certain way. Combinatorial reasoning is used to show that those connected graphs with negative integer amplitudes for which the numbers game meets a certain finiteness requirement are precisely the Dynkin diagrams associated with the finite-dimensional complex simple Lie algebras. This strengthens a result originally due to the second author. A more general result is obtained when certain real number amplitudes are allowed. The resulting graphs are in families, each family corresponding to a finite irreducible Coxeter group. These results are used to demonstrate that the only generalized Cartan matrices for which there exist finite edge-colored ranked posets enjoying a certain structure property are the Cartan matrices for the finite-dimensional complex semisimple Lie algebras. In this setting, classifications of the finite-dimensional Kac--Moody algebras and of the finite Coxeter and Weyl groups are re-derived.
△ Less
Submitted 29 October, 2008;
originally announced October 2008.
-
Optimal stopping in a two-sided secretary problem
Authors:
Kimmo Eriksson,
Jonas Sjostrand,
Pontus Strimling
Abstract:
In the "secretary problem", well-known in the theory of optimal stopping, an employer is about to interview a maximum of N secretaries about which she has no prior information. Chow et al. proved that with an optimal strategy the expected rank of the chosen secretary tends to approximately 3.87.
We study a two-sided game-theoretic version of this optimal stopping problem, where men search for…
▽ More
In the "secretary problem", well-known in the theory of optimal stopping, an employer is about to interview a maximum of N secretaries about which she has no prior information. Chow et al. proved that with an optimal strategy the expected rank of the chosen secretary tends to approximately 3.87.
We study a two-sided game-theoretic version of this optimal stopping problem, where men search for a woman to marry at the same time as women search for a man to marry. We find that in the unique subgame perfect equilibrium, the expected rank grows as the square root of N and that, surprisingly, the leading coefficient is exactly 1. We also discuss some possible variations.
△ Less
Submitted 9 November, 2004;
originally announced November 2004.
-
Note on the lamp lighting problem
Authors:
Henrik Eriksson,
Kimmo Eriksson,
Jonas Sjostrand
Abstract:
We answer some questions concerning the so called sigma-game of Sutner. It is played on a graph where each vertex has a lamp, the light of which is toggled by pressing any vertex with an edge directed to the lamp.
For example, we show that every configuration of lamps can be lit if and only if the number of complete matchings in the graph is odd. In the special case of an orthogonal grid one g…
▽ More
We answer some questions concerning the so called sigma-game of Sutner. It is played on a graph where each vertex has a lamp, the light of which is toggled by pressing any vertex with an edge directed to the lamp.
For example, we show that every configuration of lamps can be lit if and only if the number of complete matchings in the graph is odd. In the special case of an orthogonal grid one gets a criterion for whether the number of monomer-dimer tilings of an m times n grid is odd or even.
△ Less
Submitted 9 November, 2004;
originally announced November 2004.
-
Exact expectations for random graphs and assignments
Authors:
Henrik Eriksson,
Kimmo Eriksson,
Jonas Sjostrand
Abstract:
For a random graph on n vertices where the edges appear with individual rates, we give exact formulas for the expected time at which the number of components has gone down to k and the expected length of the corresponding minimal spanning forest.
For a random bipartite graph we give a formula for the expected time at which a k-assignment appears. This result has bearing upon the random assignm…
▽ More
For a random graph on n vertices where the edges appear with individual rates, we give exact formulas for the expected time at which the number of components has gone down to k and the expected length of the corresponding minimal spanning forest.
For a random bipartite graph we give a formula for the expected time at which a k-assignment appears. This result has bearing upon the random assignment problem.
△ Less
Submitted 9 November, 2004;
originally announced November 2004.
-
Expected number of inversions after a sequence of random adjacent transpositions
Authors:
Henrik Eriksson,
Kimmo Eriksson,
Jonas Sjostrand
Abstract:
In the evolution of a genome, the gene sequence is sometimes rearranged, for example by transposition of two adjacent gene blocks. In biocombinatorics, one tries to reconstruct these rearrangement incidents from the resulting permutation. It seems that the algorithms used are too effective and find a shorter path than the real one. For the simplified case of adjacent transpositions, we give expr…
▽ More
In the evolution of a genome, the gene sequence is sometimes rearranged, for example by transposition of two adjacent gene blocks. In biocombinatorics, one tries to reconstruct these rearrangement incidents from the resulting permutation. It seems that the algorithms used are too effective and find a shorter path than the real one. For the simplified case of adjacent transpositions, we give expressions for the expected number of inversions after t random moves. This average can be much smaller than t, a fact that has largely been neglected so far.
△ Less
Submitted 9 November, 2004;
originally announced November 2004.
-
Conjectures on three-dimensional stable matching
Authors:
Kimmo Eriksson,
Jonas Sjostrand,
Pontus Strimling
Abstract:
We consider stable three-dimensional matchings of three categories of agents, such as women, men and dogs. This was suggested long ago by Knuth (1976), but very little seems to have been published on this problem. Based on computer experiments, we present a couple of conjectures as well as a few counter-examples to other natural but discarded conjectures. In particular, a circular 3D matching is…
▽ More
We consider stable three-dimensional matchings of three categories of agents, such as women, men and dogs. This was suggested long ago by Knuth (1976), but very little seems to have been published on this problem. Based on computer experiments, we present a couple of conjectures as well as a few counter-examples to other natural but discarded conjectures. In particular, a circular 3D matching is one where women only care about the man, men only care about the dog, and dogs only care about the woman they are matched with. We conjecture that a stable outcome always exists for any circular 3D matching market, and we prove it for markets with at most four agents of each category.
△ Less
Submitted 1 November, 2004;
originally announced November 2004.
-
Stable matching in a common generalization of the marriage and assignment models
Authors:
Kimmo Eriksson,
Johan Karlander
Abstract:
In the theory of two-sided matching markets there are two well-known models: the marriage model (where no money is involved) and the assignment model (where payments are involved). Roth and Sotomayor (1990) asked for an explanation for the similarities in behavior between those two models. We address this question by introducing a common generalization that preserves the two important features:…
▽ More
In the theory of two-sided matching markets there are two well-known models: the marriage model (where no money is involved) and the assignment model (where payments are involved). Roth and Sotomayor (1990) asked for an explanation for the similarities in behavior between those two models. We address this question by introducing a common generalization that preserves the two important features: the existence of a stable outcome and the lattice property of the set of stable outcomes.
△ Less
Submitted 21 January, 1998;
originally announced January 1998.