Search results for " Programming"
showing 10 items of 1616 documents
Bayesian joint ordinal and survival modeling for breast cancer risk assessment
2016
We propose a joint model to analyze the structure and intensity of the association between longitudinal measurements of an ordinal marker and time to a relevant event. The longitudinal process is defined in terms of a proportional-odds cumulative logit model. Time-to-event is modeled through a left-truncated proportionalhazards model, which incorporates information of the longitudinal marker as well as baseline covariates. Both longitudinal and survival processes are connected by means of a common vector of random effects. General inferences are discussed under the Bayesian approach and include the posterior distribution of the probabilities associated to each longitudinal category and the …
On 1-Laplacian Elliptic Equations Modeling Magnetic Resonance Image Rician Denoising
2015
Modeling magnitude Magnetic Resonance Images (MRI) rician denoising in a Bayesian or generalized Tikhonov framework using Total Variation (TV) leads naturally to the consideration of nonlinear elliptic equations. These involve the so called $1$-Laplacian operator and special care is needed to properly formulate the problem. The rician statistics of the data are introduced through a singular equation with a reaction term defined in terms of modified first order Bessel functions. An existence theory is provided here together with other qualitative properties of the solutions. Remarkably, each positive global minimum of the associated functional is one of such solutions. Moreover, we directly …
Bayesian survival analysis with BUGS
2020
Survival analysis is one of the most important fields of statistics in medicine and biological sciences. In addition, the computational advances in the last decades have favored the use of Bayesian methods in this context, providing a flexible and powerful alternative to the traditional frequentist approach. The objective of this article is to summarize some of the most popular Bayesian survival models, such as accelerated failure time, proportional hazards, mixture cure, competing risks, multi-state, frailty, and joint models of longitudinal and survival data. Moreover, an implementation of each presented model is provided using a BUGS syntax that can be run with JAGS from the R programmin…
Reassessing Accuracy Rates of Median Decisions
2007
We show how Bruno de Finetti''s fundamental theorem of prevision has computable applications in statistical problems that involve only partial information. Specifically, we assess accuracy rates for median decision procedures used in the radiological diagnosis of asbestosis. Conditional exchangeability of individual radiologists'' diagnoses is recognized as more appropriate than independence which is commonly presumed. The FTP yields coherent bounds on probabilities of interest when available information is insufficient to determine a complete distribution. Further assertions that are natural to the problem motivate a partial ordering of conditional probabilities, extending the computation …
Assessing uncertainty of voter transitions estimated from aggregated data. Application to the 2017 French presidential election
2020
[EN] Inferring electoral individual behaviour from aggregated data is a very active research area, with ramifications in sociology and political science. A new approach based on linear programming is proposed to estimate voter transitions among parties (or candidates) between two elections. Compared to other linear and quadratic programming models previously published, our approach presents two important innovations. Firstly, it explicitly deals with new entries and exits in the election census without assuming unrealistic hypotheses, enabling a reasonable estimation of vote behaviour of young electors voting for the first time. Secondly, by exploiting the information contained in the model…
Macroscopic Dynamic Effects of Migrations in Patchy Predator-prey Systems
1997
Abstract Different mechanisms at the behaviourial or physiological levels determine many properties of predator-prey systems at the population level. In this paper, we present a method of obtaining complex predator-prey dynamic models from models at a detailed, behaviourial level of description. We consider a multi-patch predator-prey model, the dynamics of which contains two time-scales: a fast one, associated with migrations between patches, and a slow one, on which interactions, reproduction and mortality occur. We use methods of perturbation theory in order to aggregate the multi-patch system into a reduced system of two differential equations for the total prey and predator populations…
Evaluation of Insurance Products with Guarantee in Incomplete Markets
2008
Abstract Life insurance products are usually equipped with minimum guarantee and bonus provision options. The pricing of such claims is of vital importance for the insurance industry. Risk management, strategic asset allocation, and product design depend on the correct evaluation of the written options. Also regulators are interested in such issues since they have to be aware of the possible scenarios that the overall industry will face. Pricing techniques based on the Black & Scholes paradigm are often used, however, the hypotheses underneath this model are rarely met. To overcome Black & Scholes limitations, we develop a stochastic programming model to determine the fair price of the mini…
Calibration of optimal execution of financial transactions in the presence of transient market impact
2012
Trading large volumes of a financial asset in order driven markets requires the use of algorithmic execution dividing the volume in many transactions in order to minimize costs due to market impact. A proper design of an optimal execution strategy strongly depends on a careful modeling of market impact, i.e. how the price reacts to trades. In this paper we consider a recently introduced market impact model (Bouchaud et al., 2004), which has the property of describing both the volume and the temporal dependence of price change due to trading. We show how this model can be used to describe price impact also in aggregated trade time or in real time. We then solve analytically and calibrate wit…
A non-linear optimization procedure to estimate distances and instantaneous substitution rate matrices under the GTR model.
2006
Abstract Motivation: The general-time-reversible (GTR) model is one of the most popular models of nucleotide substitution because it constitutes a good trade-off between mathematical tractability and biological reality. However, when it is applied for inferring evolutionary distances and/or instantaneous rate matrices, the GTR model seems more prone to inapplicability than more restrictive time-reversible models. Although it has been previously noted that the causes for intractability are caused by the impossibility of computing the logarithm of a matrix characterised by negative eigenvalues, the issue has not been investigated further. Results: Here, we formally characterize the mathematic…
Designing and pricing guarantee options in defined contribution pension plans
2015
Abstract The shift from defined benefit (DB) to defined contribution (DC) is pervasive among pension funds, due to demographic changes and macroeconomic pressures. In DB all risks are borne by the provider, while in plain vanilla DC all risks are borne by the beneficiary. However, for DC to provide income security some kind of guarantee is required. A minimum guarantee clause can be modeled as a put option written on some underlying reference portfolio and we develop a discrete model that selects the reference portfolio to minimize the cost of a guarantee. While the relation DB–DC is typically viewed as a binary one, the model shows how to price a wide range of guarantees creating a continu…