On Complexity of Multistage Stochastic Programs
|
|
- Douglas Osborne
- 5 years ago
- Views:
Transcription
1 On Complexity of Multistage Stochastic Programs Alexander Shapiro School of Industrial and Systems Engineering, Georgia Institute of Technology, Atlanta, Georgia , USA Abstract In this paper we derive estimates of the sample sizes required to solve a multistage stochastic programming problem with a given accuracy by the (conditional sampling) sample average approximation method. The presented analysis is self contained and is based on a, relatively elementary, one dimensional Cramér s Large Deviations Theorem. Key words: stochastic programming, Monte Carlo sampling, sample average method, large deviations exponential bounds, complexity.
2 Introduction Consider the following stochastic programming problem } f(x) :=E[F (x, ξ)], (.) Min x X where ξ is a random vector supported on a set Ξ R d, the expectation in (.) is taken with respect to a (known) probability distribution of ξ, X is a nonempty subset of R n and F : X Ξ R. In the case of two-stage stochastic programming, the function F (x, ξ) is given as the optimal value of a corresponding second stage problem. In that case the assumption that F (x, ξ) is real valued for all x X and ξ Ξ can only hold if the corresponding recourse is relatively complete. Only in very specific situations the expected value function f(x) can be written in a closed form. Therefore, it should be calculated by a numerical integration. Already for the number of random variables d 5 it is typically impossible to evaluate the corresponding multidimensional integral (expectation) with a high accuracy. This makes stochastic programming problems of the form (.) really difficult. A way of estimating the expected value function is suggested by the Monte Carlo method. That is, a random sample ξ,..., ξ N of N realizations of ξ is generated, and the expected value function f(x) is approximated by the sample average function ˆf N (x) :=N N i= F (x, ξi ). This is the basic idea of the so-called sample average approximation (SAA) method. It is possible to show, under mild regularity conditions, that for >0 and α (0, ) the sample size N O()σ2 2 [ n log ( ) ( )] DL O() + log α (.2) guarantees that any /2-optimal solution of the SAA problem is an -optimal solution of the true problem with probability at least α (see [3, 7, 8]). Here O() is a generic constant, D is the diameter of the set X (assumed to be finite), L is Lipschitz constant of f(x) and σ 2 is a certain constant measuring variability of the objective function F (x, ξ). Recall that for >0 it is said that x is an -optimal solution of problem (.) if x X and f( x) inf x X f(x)+. In a sense the estimate (.2) of the sample size gives a bound on complexity of solving, with a specified probability, the (true) problem (.) by using the sample average approximation. Note that the estimated sample size grows linearly in the dimension n of the first stage problem and is proportional to the squared ratio of the variability coefficient σ to the desired accuracy. (The following Example shows that this estimate cannot be significantly improved.) This indicates that one may expect to solve the true problem (.) with a manageable sample size to a reasonable accuracy by using the SSA method. And, indeed, this was verified in various numerical experiments (cf., [4, 5, 9]). Example Consider problem (.) with F (x, ξ) := x 2k 2k ξ,x, where k is a positive integer, X := x R n : x }, x, y denotes the standard scalar product of two
3 vectors x, y R n and x = x, x. Suppose, further, that random vector ξ has normal distribution N(0,σ 2 I n ), where σ 2 is a positive constant, i.e., components ξ i of ξ are independent and ξ i N(0,σ 2 ), i =,..., n. It follows that f(x) = x 2k, and hence for [0, ] the set of -optimal solutions of the true problem (.) is x : x 2k }. Now let ξ,..., ξ N be an iid random sample of ξ and ξ N := (ξ ξ N )/N. The corresponding sample average function is ˆf N (x) = x 2k 2k ξ N,x, and the optimal solution ˆx N of the SAA problem is ˆx N = ξ N γ ξn, where γ := 2k 2 if ξ 2k N, and γ =if ξ N >. It follows that, for (0, ), the optimal solution of the corresponding SAA problem is an -optimal solution of the true problem iff ξ N ν, where ν := 2k. 2k We have that ξ N N(0,σ 2 N I n ), and hence N ξ N 2 /σ 2 has the chi-square distribution with n degrees of freedom. Consequently, the probability that ξ N ν >is equal to the probability P ( χ 2 n >N 2/ν /σ 2). Moreover, E[χ 2 n]=n and P(χ 2 n >n) increases and tends to /2 asn increases, e.g., P(χ 2 > )=0.373, P(χ 2 2 > 2)=0.3679, P(χ 2 3 > 3)=0.396 and etc. Consequently, for α (0, 0.3) and (0, ), for example, the sample size N should satisfy N> nσ2 (.3) 2/ν in order to have the property: with probability α an (exact) optimal solution of the SAA problem is an -optimal solution of the true problem. Compared with (.2), the lower bound (.3) also grows linearly in n and is proportional to σ 2 / 2/ν. It remains to note that the constant ν decreases to one as k increases. The aim of this paper is to extend this analysis of the SAA method to the multistage stochastic programming (MSP) setting. A discussion of complexity of the MSP can be found in [8]. It was already argued there that complexity of the SAA method, when applied to the MSP, grows fast with increase of the number of stages and seemingly simple MSP problems can be computationally unmanageable. We estimate sample sizes, required to solve the true problem with a given accuracy, by using tools of the Large Deviations (LD) theory (see, e.g., [2] for a thorough discussion of the LD theory). In that respect our analysis is self contained and rather elementary since we only employ the upper bound of the (one dimensional) Cramér s LD Theorem. That is, if X,..., X N is a sequence of iid realizations of random variable X and X N := N N i= X i is the corresponding average, then Here P(A) denotes probability of event A, } I(z) :=sup tz log M(t) t R P( X N a) e NI(a). (.4) is the so-called rate function, and M(t) :=E[e tx ] is the moment generating function of random variable X. 2
4 Let us make the following simple observation which will be used in our derivations. Let X and Y be two random variables and R. We have that if X and Y 2, where + 2 =, then X + Y, and hence ω : X(ω)+Y (ω) >} ω : X(ω) > } ω : Y (ω) > 2 }. This implies the following inequality for the corresponding probabilities P(X + Y>) P(X > )+P(Y > 2 ). (.5) 2 Sample average approximations of multistage stochastic programs Consider the following T stage stochastic programming problem [ [ Min F (x )+E inf F 2 (x 2,ξ 2 )+E + E [ inf F T (x T,ξ T ) ]]] (2.) x 2 X 2 (x,ξ 2 ) x T X T (x T,ξ T ) driven by the random data process ξ 2,..., ξ T. Here x t R nt, t =,..., T, are decision variables, F t : R nt R dt R are continuous functions and X t : R n t R dt R nt, t = 2,..., T, are measurable multifunctions, the function F : R n R and the set X R n are deterministic. We assume that the set X is nonempty. For example, in linear case F t (x t,ξ t ):= c t,x t, X := x : A x = b,x 0}, X t (x t,ξ t ):=x t : B t x t + A t x t = b t,x t 0},t=2,..., T, ξ := (c,a,b ) is known at the first stage (and hence is nonrandom), and ξ t := (c t,b t,a t,b t ) R dt, t =2,..., T, are data vectors some (all) elements of which can be random. In the sequel we use ξ t to denote random data vector and its particular realization. Which one of these two meanings will be used in a particular situation will be clear from the context. If we denote by Q 2 (x,ξ 2 ) the optimal value of the (T ) stage problem: Min x 2 X 2 (x,ξ 2 ) F 2(x 2,ξ 2 )+E [ + E [ min x T X T (x T,ξ T ) F T (x T,ξ T ) ]], (2.2) then we can write the T stage problem (2.) in the following form of two-stage programming problem Min F (x )+E [ Q 2 (x,ξ 2 ) ]. (2.3) Note, however, that if T > 2, then problem (2.2) in itself is a stochastic programming problem. Consequently, if the number of scenarios involved in (2.2) is very large, or 3
5 infinite, then the optimal value Q 2 (x,ξ 2 ) can be calculated only approximately, say by sampling. For the sake of simplicity we make the following derivations for the 3-stage problem, i.e., we assume that T = 3 (it will be clear how the obtained results can be extended to an analysis of T>3). In that case Q 2 (x,ξ 2 ) is given by the optimal value of the problem Min F 2 (x 2,ξ 2 )+E[Q 3 (x 2,ξ 3 ) ξ 2 ], (2.4) x 2 X 2 (x,ξ 2 ) where the expectation is taken with respect to the conditional distribution of ξ 3 given ξ 2 and Q 3 (x 2,ξ 3 ) := We make the following assumption: inf F 3 (x 3,ξ 3 ). x 3 X 3 (x 2,ξ 3 ) For every the expectation E [ Q 2 (x,ξ 2 ) ] is well defined and finite valued. Of course, finite valuedness of E [ Q 2 (x,ξ 2 ) ] can only holds if Q 2 (x,ξ 2 ) is finite valued for a.e. ξ 2, which in turn implies that X 2 (x,ξ 2 ) is nonemty for a.e. ξ 2 etc. That is, the above assumption implies that the recourse is relatively complete. Now let ξ2, i i =,..., N, be a random sample of independent realizations of the random vector ξ 2. We can approximate problem (2.3) by the following SAA problem } Min ˆf N (x ):=F (x )+ N Q 2 (x,ξ i N 2). (2.5) Since Q 2 (x,ξ2) i are not given explicitly we need to estimate these values by conditional sampling (note that in order for the SAA method to produce consistent estimators, conditional sampling is required, see [6]). That is, we generate random sample ξ ij 3, j =,..., N 2,ofN 2 independent realizations according to conditional distribution of ξ 3 given ξ2, i i =,..., N. Consequently, we approximate Q 2 (x,ξ2)by i } ˆQ 2,N2 (x,ξ2) i := inf F 2 (x 2,ξ2)+ i N 2 Q 3 (x 2,ξ ij x 2 X 2 (x,ξ2 i ) 3 ). (2.6) N 2 Finally, we approximate the true (expected value) problem (2.3) by the following socalled Sample Average Approximating (SAA) problem } Min f N,N 2 (x ):=F (x )+ N ˆQ 2,N2 (x,ξ i N 2). (2.7) The above SAA problem is obtained by approximating the objective function f(x ):=F (x )+E [ Q 2 (x,ξ 2 ) ] of problem (2.3) with f N,N 2 (x ). 4 i= i= j=
6 3 Sample size estimates In order to proceed with our analysis we need to estimate the probability P sup f(x ) f N,N 2 (x ) } > (3.) for an arbitrary constant >0. To this end we use the following result about Large Deviations (LD) bounds for the uniform convergence of sample average approximations. Consider a function h : X Ξ R and the corresponding expected value function φ(x) := E[h(x, ξ)], where the expectation is taken with respect to the probability distribution P of the random vector ξ = ξ(ω), X is a nonempty closed subset of R n and Ξ R d is the support of the probability distribution P. Assume that for every x X the expectation φ(x) is well defined, i.e., h(x, ) is measurable and P -integrable. Let ξ,..., ξ N be an iid sample of the random vector ξ(ω), and ˆφ N (x) := N N i= h(x, ξi ) be the corresponding sample average function. Theorem Suppose that the set X has finite diameter D, and the following conditions hold: (i) there exists a constant σ>0 such that M x (t) exp σ 2 t 2 /2 }, t R, x X, (3.2) where M x (t) is the moment generating function of the random variable h(x, ξ) φ(x), (ii) there exists a constant L>0 such that Then for any >0, P h(x,ξ) h(x, ξ) L x x, ξ Ξ, x,x X. (3.3) sup x X ˆφN (x) φ(x) } ( ) n DL O() exp N2 6σ 2 }. (3.4) To make the paper self contained we give a proof of this theorem in the appendix. We can apply the LD bound (3.4) to obtain estimates of the probability (3.). We have that sup f(x ) f N,N 2 (x ) sup f(x ) ˆf N (x ) + sup ˆfN (x ) f N,N 2 (x ), and hence, by (.5), P sup f(x ) f N,N 2 (x ) } > P sup f(x ) ˆf N (x ) } >/2 + P sup ˆfN (x ) f N,N 2 (x ) } >/2. (3.5) 5
7 Note that f(x ) ˆf N (x )=E[Q 2 (x,ξ 2 )] N Q 2 (x,ξ N 2), i i= and ˆf N (x ) f N,N 2 (x )= N [ Q 2 (x,ξ N 2) i ˆQ ] 2,N2 (x,ξ2) i. i= Let us assume, for the sake of simplicity, the between stages independence of the random process. That is, suppose that the following condition holds. (A) Random vectors ξ 2 and ξ 3 are independent. Of course, under this condition the conditional expectation in formula (2.4) does not depend on ξ 2. Also in that case the conditional sample ξ ij 3 has the (marginal) distribution of ξ 3 and is independent of ξ2, i and can be generated in two ways. Namely, we can either generate the same random sample ξ ij 3 for each i =,..., N, or these samples can be generated independently of each other. Let us make further the following assumptions. (A2) The set X has finite diameter D. (A3) There is a constant L > 0 such that Q 2 (x,ξ 2 ) Q 2 (x,ξ 2 ) L x x for all x, and a.e. ξ 2. (A4) There exists constant σ > 0 such that for any it holds that M,x (t) exp σ 2 t 2 /2 }, t R, where M,x (t) is the moment generating function of Q 2 (x,ξ 2 ) E[Q 2 (x,ξ 2 )]. (A5) There is a positive constant D 2 such that for every and a.e. ξ 2 the set X 2 (x,ξ 2 ) has a finite diameter less than or equal to D 2. (A6) There is a constant L 2 > 0 such that F2 (x 2,ξ 2 ) F 2 (x 2,ξ 2 )+Q 3 (x 2,ξ 3 ) Q 3 (x 2,ξ 3 ) L2 x 2 x 2 for all x 2,x 2 X 2 (x,ξ 2 ), and a.e. ξ 2 and ξ 3. 6
8 (A7) There exists constant σ 2 > 0 such that for any x 2 X 2 (x,ξ 2 ) and all and a.e. ξ 2 it holds that M 2,x2 (t) exp σ 2 2t 2 /2 }, t R, where M 2,x2 (t) is the moment generating function of Q 3 (x 2,ξ 3 ) E[Q 3 (x 2,ξ 3 )]. By Theorem, under assumptions (A2) (A4), we have that P sup f(x ) ˆf N (x ) } ( ) n D L >/2 O() exp O()N } 2. (3.6) σ 2 For ξ 2 and random sample ξ3,..., ξ N 2 3 of N 2 independent replications of ξ 3, consider function ˆψ N2 (x 2,ξ 2 ):=F 2 (x 2,ξ 2 )+ N 2 Q 3 (x 2,ξ N 3) j 2 and its expected value j= ψ(x 2,ξ 2 )=F 2 (x 2,ξ 2 )+E[Q 3 (x 2,ξ 3 )]. By Theorem, under assumptions (A5) (A7), we have that for any, P sup ˆψN2 (x 2,ξ 2 ) ψ(x 2,ξ 2 ) } >/2 C (N 2 ), (3.7) x 2 X 2 (x,ξ 2 ) where ( ) n2 D2 L 2 C (N 2 )=O() exp O()N } 2 2. σ2 2 It follows that } P inf ˆψ N2 (x 2,ξ 2 ) inf ψ(x 2,ξ 2 ) x 2 X 2 (x,ξ 2 ) x 2 X 2 (x,ξ 2 ) >/2 C (N 2 ). (3.8) Note that inf x 2 X 2 (x,ξ 2 ) ψ(x 2,ξ 2 )=Q 2 (x,ξ 2 ), and for ξ 2 = ξ2, i inf ˆψ N2 (x 2,ξ 2 )= ˆQ 2,N2 (x,ξ2). i x 2 X 2 (x,ξ 2 ) It follows from (3.8) that (for both strategies of using the same or independent samples for each ξ2) i the following inequality holds P ˆfN (x ) f N,N 2 (x ) } >/2 C (N 2 ). (3.9) Suppose further that: 7
9 There is L 3 > 0 such that ˆf N ( ) f N,N 2 ( ) is Lipschitz continuous on X with constant L 3. Then by constructing a ν-net in X and using (3.9) it can be shown (in a way similar to the proof of Theorem in the Appendix) that P sup ˆfN (x ) f N,N 2 (x ) } >/2 O() ( D L 3 ) n ( D2 L 2 ) } n2 exp O()N 2 2. σ 2 x X 2 (3.0) Combining (3.5) with estimates (3.6) and (3.0) gives an upper bound for the probability (3.). Let us also observe that if ˆx is an /2-optimal solution of the SAA problem (2.7) and sup x X f(x ) f N,N 2 (x ) /2, then ˆx is an -optimal solution of the true problem (2.3). Therefore, we obtain the following result. Theorem 2 Under the specified assumptions and for >0and α (0, ), and the sample sizes N and N 2 satisfying [ ( O() D ) } L n exp O()N 2 + ( ) D L n ( 3 D2 ) }] L n2 2 exp O()N 2 2 α, (3.) σ 2 we have that any /2-optimal solution of the SAA problem (2.7) is an -optimal solution of the true problem (2.3) with probability at least α. In particular, suppose that N = N 2. Then for L := maxl,l 2,L 3 }, D := maxd,d 2 } and σ := maxσ,σ 2 } we can use the following estimate of the required sample size N = N 2 : ( ) n +n DL 2 O() exp O()N } 2 α, (3.2) σ 2 which is equivalent to 4 Discussion N O()σ2 2 [ (n + n 2 ) log ( ) ( DL O() + log α σ 2 2 )]. (3.3) The estimate (3.3), for 3-stage programs, looks similar to the estimate (.2) for two-stage programs. Note, however, that if we use the SAA method with conditional sampling and respective sample sizes N and N 2, then the total number of scenarios is N = N N 2. Therefore, our analysis seems to indicate that for 3-stage problems we need random samples with the total number of scenarios of order of the square of the corresponding sample size for two-stage problems. This analysis can be extended to T -stage problems with the conclusion that the total number of scenarios needed to solve the true problem with a 8
10 reasonable accuracy grows exponentially with increase of the number of stages T. Some numerical experiments seem to confirm this conclusion (cf., []). Of course, it should be mentioned that the above analysis does not prove in a rigorous mathematical sense that complexity of multistage programming grows exponentially with increase of the number of stages. It only indicates that the SAA method, which showed a considerable promise for solving two stage problems, could be practically inapplicable for solving multistage problems with a large (say greater than 5) number of stages. Our analysis was performed under several simplifying assumptions. In particular, we consider a 3-stage setting and assumed the between stages independence condition. An extension of the analysis from 3 to a higher number of stages is straightforward. Removing the between stages independence assumption may create technical difficulties and requires a further investigation. 5 Appendix Proof of Theorem. By the LD bound (.4) we have that for any x X and >0it holds that } P ˆφN (x) φ(x) exp NI x ()}, (5.) where I x (z) :=sup zt log Mx (t) } (5.2) t R is the rate function of random variable h(x, ξ) φ(x). Similarly } P ˆφN (x) φ(x) exp NI x ( )}, and hence P ˆφN (x) φ(x) } exp NI x ()} + exp NI x ( )}. (5.3) For a constant ν>0, let x,..., x M X be such that for every x X there exists x l, l,..., M}, such that x x l ν. Such set x,..., x M } is called a ν-net in X.We can choose this net in such a way that M O()(D/ν) n, where D := sup x,x X x x is the diameter of X and O() is a generic constant. By (3.3) we have that φ(x ) φ(x) L x x (5.4) and ˆφN (x ) ˆφ N (x) L x x, (5.5) 9
11 for any x, x X. It follows by (5.3) that ( P max ˆφN ( x l ) φ( x l ) ) ( M P l= ˆφN ( x l ) φ( x l ) }) l M M ( P ˆφN ( x l ) φ( x l ) ) 2 M exp N[I xl () I xl ( )]}. l= l= (5.6) For an x X consider l(x) arg min l M x x l. By construction of the ν-net we have that x x l(x) ν for every x X. Then ˆφN (x) φ(x) ˆφN (x) ˆφ N ( x l(x) ) + ˆφN ( x l(x) ) φ( x l(x) ) + φ( x l(x) ) φ(x) Lν + ˆφN ( x l(x) ) φ( x l(x) ) + Lν. Let us take now a ν-net with such ν that Lν = /4, i.e., ν := [/(4L)]. Then P ˆφN (x) φ(x) } P ˆφN ( x l ) φ( x l ) } /2, sup x X max l M which together with (5.6) implies that P ˆφN (x) φ(x) } 2 M exp N [ I xl (/2) I xl ( /2) ]}. (5.7) sup x X l= Moreover, because of the condition (i) we have that log M x (t) σ 2 t 2 /2, and hence I x (z) z2, z R. (5.8) 2σ2 It follows from (5.7) and (5.8) that P sup x X ˆφN (x) φ(x) } } N 2M exp 2. (5.9) 6σ 2 Finally, since M O()(D/ν) n = O()(DL/) n, we obtain that (5.9) implies (3.4), and hence the proof is complete. References [] J. Blomvall and A. Shapiro, Solving multistage asset investment problems by Monte Carlo based optimization, E-print available at: [2] A. Dembo and O. Zeitouni, Large Deviations Techniques and Applications, Springer- Verlag, New York, NY,
12 [3] A. J. Kleywegt, A. Shapiro, and T. Homem-De-Mello, The sample average approximation method for stochastic discrete optimization, SIAM Journal of Optimization, 2 (200), [4] J. Linderoth, A. Shapiro, and S. Wright, The empirical behavior of sampling methods for stochastic programming, Annals of Operations Research, to appear. [5] W. K. Mak, D. P. Morton, and R. K. Wood, Monte Carlo bounding techniques for determining solution quality in stochastic programs, Operations Research Letters, 24 (999), [6] A. Shapiro, Inference of statistical bounds for multistage stochastic programming problems, Mathematical Methods of Operations Research, 58 (2003), [7] A. Shapiro, Monte Carlo sampling methods. In: A. Rusczyński and A. Shapiro (editors), Stochastic Programming, volume 0 of Handbooks in Operations Research and Management Science, North-Holland, [8] A. Shapiro and A. Nemirovski, On complexity of stochastic programming problems, E-print available at: [9] B. Verweij, S. Ahmed, A.J. Kleywegt, G. Nemhauser, and A. Shapiro, The sample average approximation method applied to stochastic routing problems: a computational study, Computational Optimization and Applications, 24 (2003),
High Dimensional Edgeworth Expansion. Applications to Bootstrap and Its Variants
With Applications to Bootstrap and Its Variants Department of Statistics, UC Berkeley Stanford-Berkeley Colloquium, 2016 Francis Ysidro Edgeworth (1845-1926) Peter Gavin Hall (1951-2016) Table of Contents
More informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
More informationScenario Generation and Sampling Methods
Scenario Generation and Sampling Methods Güzin Bayraksan Tito Homem-de-Mello SVAN 2016 IMPA May 9th, 2016 Bayraksan (OSU) & Homem-de-Mello (UAI) Scenario Generation and Sampling SVAN IMPA May 9 1 / 30
More informationarxiv: v1 [math.st] 18 Sep 2018
Gram Charlier and Edgeworth expansion for sample variance arxiv:809.06668v [math.st] 8 Sep 08 Eric Benhamou,* A.I. SQUARE CONNECT, 35 Boulevard d Inkermann 900 Neuilly sur Seine, France and LAMSADE, Universit
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationApproximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications
Approximation of Continuous-State Scenario Processes in Multi-Stage Stochastic Optimization and its Applications Anna Timonina University of Vienna, Abraham Wald PhD Program in Statistics and Operations
More informationThe ruin probabilities of a multidimensional perturbed risk model
MATHEMATICAL COMMUNICATIONS 231 Math. Commun. 18(2013, 231 239 The ruin probabilities of a multidimensional perturbed risk model Tatjana Slijepčević-Manger 1, 1 Faculty of Civil Engineering, University
More informationAsymptotic results discrete time martingales and stochastic algorithms
Asymptotic results discrete time martingales and stochastic algorithms Bernard Bercu Bordeaux University, France IFCAM Summer School Bangalore, India, July 2015 Bernard Bercu Asymptotic results for discrete
More informationAssessing Policy Quality in Multi-stage Stochastic Programming
Assessing Policy Quality in Multi-stage Stochastic Programming Anukal Chiralaksanakul and David P. Morton Graduate Program in Operations Research The University of Texas at Austin Austin, TX 78712 January
More informationApproximations of Stochastic Programs. Scenario Tree Reduction and Construction
Approximations of Stochastic Programs. Scenario Tree Reduction and Construction W. Römisch Humboldt-University Berlin Institute of Mathematics 10099 Berlin, Germany www.mathematik.hu-berlin.de/~romisch
More informationarxiv: v1 [q-fin.pm] 13 Mar 2014
MERTON PORTFOLIO PROBLEM WITH ONE INDIVISIBLE ASSET JAKUB TRYBU LA arxiv:143.3223v1 [q-fin.pm] 13 Mar 214 Abstract. In this paper we consider a modification of the classical Merton portfolio optimization
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationNon replication of options
Non replication of options Christos Kountzakis, Ioannis A Polyrakis and Foivos Xanthos June 30, 2008 Abstract In this paper we study the scarcity of replication of options in the two period model of financial
More informationMultistage risk-averse asset allocation with transaction costs
Multistage risk-averse asset allocation with transaction costs 1 Introduction Václav Kozmík 1 Abstract. This paper deals with asset allocation problems formulated as multistage stochastic programming models.
More informationWorst-case-expectation approach to optimization under uncertainty
Worst-case-expectation approach to optimization under uncertainty Wajdi Tekaya Joint research with Alexander Shapiro, Murilo Pereira Soares and Joari Paulo da Costa : Cambridge Systems Associates; : Georgia
More informationMartingales. by D. Cox December 2, 2009
Martingales by D. Cox December 2, 2009 1 Stochastic Processes. Definition 1.1 Let T be an arbitrary index set. A stochastic process indexed by T is a family of random variables (X t : t T) defined on a
More informationSolving real-life portfolio problem using stochastic programming and Monte-Carlo techniques
Solving real-life portfolio problem using stochastic programming and Monte-Carlo techniques 1 Introduction Martin Branda 1 Abstract. We deal with real-life portfolio problem with Value at Risk, transaction
More informationOptimal Security Liquidation Algorithms
Optimal Security Liquidation Algorithms Sergiy Butenko Department of Industrial Engineering, Texas A&M University, College Station, TX 77843-3131, USA Alexander Golodnikov Glushkov Institute of Cybernetics,
More informationDASC: A DECOMPOSITION ALGORITHM FOR MULTISTAGE STOCHASTIC PROGRAMS WITH STRONGLY CONVEX COST FUNCTIONS
DASC: A DECOMPOSITION ALGORITHM FOR MULTISTAGE STOCHASTIC PROGRAMS WITH STRONGLY CONVEX COST FUNCTIONS Vincent Guigues School of Applied Mathematics, FGV Praia de Botafogo, Rio de Janeiro, Brazil vguigues@fgv.br
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationOn Existence of Equilibria. Bayesian Allocation-Mechanisms
On Existence of Equilibria in Bayesian Allocation Mechanisms Northwestern University April 23, 2014 Bayesian Allocation Mechanisms In allocation mechanisms, agents choose messages. The messages determine
More informationForecast Horizons for Production Planning with Stochastic Demand
Forecast Horizons for Production Planning with Stochastic Demand Alfredo Garcia and Robert L. Smith Department of Industrial and Operations Engineering Universityof Michigan, Ann Arbor MI 48109 December
More informationIntroduction to Probability Theory and Stochastic Processes for Finance Lecture Notes
Introduction to Probability Theory and Stochastic Processes for Finance Lecture Notes Fabio Trojani Department of Economics, University of St. Gallen, Switzerland Correspondence address: Fabio Trojani,
More informationThe rth moment of a real-valued random variable X with density f(x) is. x r f(x) dx
1 Cumulants 1.1 Definition The rth moment of a real-valued random variable X with density f(x) is µ r = E(X r ) = x r f(x) dx for integer r = 0, 1,.... The value is assumed to be finite. Provided that
More informationOptimal Allocation of Policy Limits and Deductibles
Optimal Allocation of Policy Limits and Deductibles Ka Chun Cheung Email: kccheung@math.ucalgary.ca Tel: +1-403-2108697 Fax: +1-403-2825150 Department of Mathematics and Statistics, University of Calgary,
More information4 Martingales in Discrete-Time
4 Martingales in Discrete-Time Suppose that (Ω, F, P is a probability space. Definition 4.1. A sequence F = {F n, n = 0, 1,...} is called a filtration if each F n is a sub-σ-algebra of F, and F n F n+1
More informationConvergence Analysis of Monte Carlo Calibration of Financial Market Models
Analysis of Monte Carlo Calibration of Financial Market Models Christoph Käbe Universität Trier Workshop on PDE Constrained Optimization of Certain and Uncertain Processes June 03, 2009 Monte Carlo Calibration
More informationSTAT 830 Convergence in Distribution
STAT 830 Convergence in Distribution Richard Lockhart Simon Fraser University STAT 830 Fall 2013 Richard Lockhart (Simon Fraser University) STAT 830 Convergence in Distribution STAT 830 Fall 2013 1 / 31
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
More informationEquivalence between Semimartingales and Itô Processes
International Journal of Mathematical Analysis Vol. 9, 215, no. 16, 787-791 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/1.12988/ijma.215.411358 Equivalence between Semimartingales and Itô Processes
More information4: SINGLE-PERIOD MARKET MODELS
4: SINGLE-PERIOD MARKET MODELS Marek Rutkowski School of Mathematics and Statistics University of Sydney Semester 2, 2016 M. Rutkowski (USydney) Slides 4: Single-Period Market Models 1 / 87 General Single-Period
More informationLesson 3: Basic theory of stochastic processes
Lesson 3: Basic theory of stochastic processes Dipartimento di Ingegneria e Scienze dell Informazione e Matematica Università dell Aquila, umberto.triacca@univaq.it Probability space We start with some
More informationStochastic Dual Dynamic Programming Algorithm for Multistage Stochastic Programming
Stochastic Dual Dynamic Programg Algorithm for Multistage Stochastic Programg Final presentation ISyE 8813 Fall 2011 Guido Lagos Wajdi Tekaya Georgia Institute of Technology November 30, 2011 Multistage
More informationMonte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50)
Monte Carlo and Empirical Methods for Stochastic Inference (MASM11/FMSN50) Magnus Wiktorsson Centre for Mathematical Sciences Lund University, Sweden Lecture 5 Sequential Monte Carlo methods I January
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Other Miscellaneous Topics and Applications of Monte-Carlo Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationOutline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0.
Outline Coordinate Minimization Daniel P. Robinson Department of Applied Mathematics and Statistics Johns Hopkins University November 27, 208 Introduction 2 Algorithms Cyclic order with exact minimization
More informationFast Convergence of Regress-later Series Estimators
Fast Convergence of Regress-later Series Estimators New Thinking in Finance, London Eric Beutner, Antoon Pelsser, Janina Schweizer Maastricht University & Kleynen Consultants 12 February 2014 Beutner Pelsser
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationFinancial Risk Forecasting Chapter 9 Extreme Value Theory
Financial Risk Forecasting Chapter 9 Extreme Value Theory Jon Danielsson 2017 London School of Economics To accompany Financial Risk Forecasting www.financialriskforecasting.com Published by Wiley 2011
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationNumerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More informationInformation, Interest Rates and Geometry
Information, Interest Rates and Geometry Dorje C. Brody Department of Mathematics, Imperial College London, London SW7 2AZ www.imperial.ac.uk/people/d.brody (Based on work in collaboration with Lane Hughston
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationAdvanced Probability and Applications (Part II)
Advanced Probability and Applications (Part II) Olivier Lévêque, IC LTHI, EPFL (with special thanks to Simon Guilloud for the figures) July 31, 018 Contents 1 Conditional expectation Week 9 1.1 Conditioning
More informationOPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF FINITE
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 WeA11.6 OPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF
More informationHomework Assignments
Homework Assignments Week 1 (p. 57) #4.1, 4., 4.3 Week (pp 58 6) #4.5, 4.6, 4.8(a), 4.13, 4.0, 4.6(b), 4.8, 4.31, 4.34 Week 3 (pp 15 19) #1.9, 1.1, 1.13, 1.15, 1.18 (pp 9 31) #.,.6,.9 Week 4 (pp 36 37)
More informationROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices
ROM SIMULATION Exact Moment Simulation using Random Orthogonal Matrices Bachelier Finance Society Meeting Toronto 2010 Henley Business School at Reading Contact Author : d.ledermann@icmacentre.ac.uk Alexander
More information1 Rare event simulation and importance sampling
Copyright c 2007 by Karl Sigman 1 Rare event simulation and importance sampling Suppose we wish to use Monte Carlo simulation to estimate a probability p = P (A) when the event A is rare (e.g., when p
More informationAmerican Option Pricing Formula for Uncertain Financial Market
American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn
More informationMartingale Pricing Theory in Discrete-Time and Discrete-Space Models
IEOR E4707: Foundations of Financial Engineering c 206 by Martin Haugh Martingale Pricing Theory in Discrete-Time and Discrete-Space Models These notes develop the theory of martingale pricing in a discrete-time,
More information3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors
3.4 Copula approach for modeling default dependency Two aspects of modeling the default times of several obligors 1. Default dynamics of a single obligor. 2. Model the dependence structure of defaults
More information3.2 No-arbitrage theory and risk neutral probability measure
Mathematical Models in Economics and Finance Topic 3 Fundamental theorem of asset pricing 3.1 Law of one price and Arrow securities 3.2 No-arbitrage theory and risk neutral probability measure 3.3 Valuation
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Further Variance Reduction Methods Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com Outline
More informationPricing Volatility Derivatives with General Risk Functions. Alejandro Balbás University Carlos III of Madrid
Pricing Volatility Derivatives with General Risk Functions Alejandro Balbás University Carlos III of Madrid alejandro.balbas@uc3m.es Content Introduction. Describing volatility derivatives. Pricing and
More informationCONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES
CONVERGENCE OF OPTION REWARDS FOR MARKOV TYPE PRICE PROCESSES MODULATED BY STOCHASTIC INDICES D. S. SILVESTROV, H. JÖNSSON, AND F. STENBERG Abstract. A general price process represented by a two-component
More informationThe value of foresight
Philip Ernst Department of Statistics, Rice University Support from NSF-DMS-1811936 (co-pi F. Viens) and ONR-N00014-18-1-2192 gratefully acknowledged. IMA Financial and Economic Applications June 11, 2018
More informationMATH 5510 Mathematical Models of Financial Derivatives. Topic 1 Risk neutral pricing principles under single-period securities models
MATH 5510 Mathematical Models of Financial Derivatives Topic 1 Risk neutral pricing principles under single-period securities models 1.1 Law of one price and Arrow securities 1.2 No-arbitrage theory and
More informationDRAFT. 1 exercise in state (S, t), π(s, t) = 0 do not exercise in state (S, t) Review of the Risk Neutral Stock Dynamics
Chapter 12 American Put Option Recall that the American option has strike K and maturity T and gives the holder the right to exercise at any time in [0, T ]. The American option is not straightforward
More informationStratified Sampling in Monte Carlo Simulation: Motivation, Design, and Sampling Error
South Texas Project Risk- Informed GSI- 191 Evaluation Stratified Sampling in Monte Carlo Simulation: Motivation, Design, and Sampling Error Document: STP- RIGSI191- ARAI.03 Revision: 1 Date: September
More informationMTH6154 Financial Mathematics I Stochastic Interest Rates
MTH6154 Financial Mathematics I Stochastic Interest Rates Contents 4 Stochastic Interest Rates 45 4.1 Fixed Interest Rate Model............................ 45 4.2 Varying Interest Rate Model...........................
More informationc 2014 CHUAN XU ALL RIGHTS RESERVED
c 2014 CHUAN XU ALL RIGHTS RESERVED SIMULATION APPROACH TO TWO-STAGE BOND PORTFOLIO OPTIMIZATION PROBLEM BY CHUAN XU A thesis submitted to the Graduate School New Brunswick Rutgers, The State University
More informationEstimation of Value at Risk and ruin probability for diffusion processes with jumps
Estimation of Value at Risk and ruin probability for diffusion processes with jumps Begoña Fernández Universidad Nacional Autónoma de México joint work with Laurent Denis and Ana Meda PASI, May 21 Begoña
More informationSPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin
SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market
More informationAsymptotic methods in risk management. Advances in Financial Mathematics
Asymptotic methods in risk management Peter Tankov Based on joint work with A. Gulisashvili Advances in Financial Mathematics Paris, January 7 10, 2014 Peter Tankov (Université Paris Diderot) Asymptotic
More informationRisk Minimization Control for Beating the Market Strategies
Risk Minimization Control for Beating the Market Strategies Jan Večeř, Columbia University, Department of Statistics, Mingxin Xu, Carnegie Mellon University, Department of Mathematical Sciences, Olympia
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More informationViability, Arbitrage and Preferences
Viability, Arbitrage and Preferences H. Mete Soner ETH Zürich and Swiss Finance Institute Joint with Matteo Burzoni, ETH Zürich Frank Riedel, University of Bielefeld Thera Stochastics in Honor of Ioannis
More informationIEOR E4602: Quantitative Risk Management
IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationHIGHER ORDER BINARY OPTIONS AND MULTIPLE-EXPIRY EXOTICS
Electronic Journal of Mathematical Analysis and Applications Vol. (2) July 203, pp. 247-259. ISSN: 2090-792X (online) http://ejmaa.6te.net/ HIGHER ORDER BINARY OPTIONS AND MULTIPLE-EXPIRY EXOTICS HYONG-CHOL
More information1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:
1 Mathematics in a Pill The purpose of this chapter is to give a brief outline of the probability theory underlying the mathematics inside the book, and to introduce necessary notation and conventions
More informationEVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS. Rick Katz
1 EVA Tutorial #1 BLOCK MAXIMA APPROACH IN HYDROLOGIC/CLIMATE APPLICATIONS Rick Katz Institute for Mathematics Applied to Geosciences National Center for Atmospheric Research Boulder, CO USA email: rwk@ucar.edu
More informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationVariable-Number Sample-Path Optimization
Noname manuscript No. (will be inserted by the editor Geng Deng Michael C. Ferris Variable-Number Sample-Path Optimization the date of receipt and acceptance should be inserted later Abstract The sample-path
More informationYao s Minimax Principle
Complexity of algorithms The complexity of an algorithm is usually measured with respect to the size of the input, where size may for example refer to the length of a binary word describing the input,
More informationOn Using Shadow Prices in Portfolio optimization with Transaction Costs
On Using Shadow Prices in Portfolio optimization with Transaction Costs Johannes Muhle-Karbe Universität Wien Joint work with Jan Kallsen Universidad de Murcia 12.03.2010 Outline The Merton problem The
More informationMonte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationGAME THEORY. Department of Economics, MIT, Follow Muhamet s slides. We need the following result for future reference.
14.126 GAME THEORY MIHAI MANEA Department of Economics, MIT, 1. Existence and Continuity of Nash Equilibria Follow Muhamet s slides. We need the following result for future reference. Theorem 1. Suppose
More informationOn the Distribution of Multivariate Sample Skewness for Assessing Multivariate Normality
On the Distribution of Multivariate Sample Skewness for Assessing Multivariate Normality Naoya Okamoto and Takashi Seo Department of Mathematical Information Science, Faculty of Science, Tokyo University
More informationS t d with probability (1 p), where
Stochastic Calculus Week 3 Topics: Towards Black-Scholes Stochastic Processes Brownian Motion Conditional Expectations Continuous-time Martingales Towards Black Scholes Suppose again that S t+δt equals
More informationA Continuity Correction under Jump-Diffusion Models with Applications in Finance
A Continuity Correction under Jump-Diffusion Models with Applications in Finance Cheng-Der Fuh 1, Sheng-Feng Luo 2 and Ju-Fang Yen 3 1 Institute of Statistical Science, Academia Sinica, and Graduate Institute
More informationPORTFOLIO OPTIMIZATION AND EXPECTED SHORTFALL MINIMIZATION FROM HISTORICAL DATA
PORTFOLIO OPTIMIZATION AND EXPECTED SHORTFALL MINIMIZATION FROM HISTORICAL DATA We begin by describing the problem at hand which motivates our results. Suppose that we have n financial instruments at hand,
More informationELEMENTS OF MONTE CARLO SIMULATION
APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the
More informationinduced by the Solvency II project
Asset Les normes allocation IFRS : new en constraints assurance induced by the Solvency II project 36 th International ASTIN Colloquium Zürich September 005 Frédéric PLANCHET Pierre THÉROND ISFA Université
More informationLecture Notes 6. Assume F belongs to a family of distributions, (e.g. F is Normal), indexed by some parameter θ.
Sufficient Statistics Lecture Notes 6 Sufficiency Data reduction in terms of a particular statistic can be thought of as a partition of the sample space X. Definition T is sufficient for θ if the conditional
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationSupplementary Material for Combinatorial Partial Monitoring Game with Linear Feedback and Its Application. A. Full proof for Theorems 4.1 and 4.
Supplementary Material for Combinatorial Partial Monitoring Game with Linear Feedback and Its Application. A. Full proof for Theorems 4.1 and 4. If the reader will recall, we have the following problem-specific
More informationOptimal Dam Management
Optimal Dam Management Michel De Lara et Vincent Leclère July 3, 2012 Contents 1 Problem statement 1 1.1 Dam dynamics.................................. 2 1.2 Intertemporal payoff criterion..........................
More informationConvergence of trust-region methods based on probabilistic models
Convergence of trust-region methods based on probabilistic models A. S. Bandeira K. Scheinberg L. N. Vicente October 24, 2013 Abstract In this paper we consider the use of probabilistic or random models
More informationChapter 6. Importance sampling. 6.1 The basics
Chapter 6 Importance sampling 6.1 The basics To movtivate our discussion consider the following situation. We want to use Monte Carlo to compute µ E[X]. There is an event E such that P(E) is small but
More informationNo-arbitrage theorem for multi-factor uncertain stock model with floating interest rate
Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer
More informationModel-independent bounds for Asian options
Model-independent bounds for Asian options A dynamic programming approach Alexander M. G. Cox 1 Sigrid Källblad 2 1 University of Bath 2 CMAP, École Polytechnique University of Michigan, 2nd December,
More information12 The Bootstrap and why it works
12 he Bootstrap and why it works For a review of many applications of bootstrap see Efron and ibshirani (1994). For the theory behind the bootstrap see the books by Hall (1992), van der Waart (2000), Lahiri
More informationWeek 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals
Week 2 Quantitative Analysis of Financial Markets Hypothesis Testing and Confidence Intervals Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg :
More informationPosterior Inference. , where should we start? Consider the following computational procedure: 1. draw samples. 2. convert. 3. compute properties
Posterior Inference Example. Consider a binomial model where we have a posterior distribution for the probability term, θ. Suppose we want to make inferences about the log-odds γ = log ( θ 1 θ), where
More informationStochastic Optimization with cvxpy EE364b Project Final Report
Stochastic Optimization with cvpy EE364b Project Final Report Alnur Ali alnurali@cmu.edu June 5, 2015 1 Introduction A stochastic program is a conve optimization problem that includes random variables,
More informationBlackwell Optimality in Markov Decision Processes with Partial Observation
Blackwell Optimality in Markov Decision Processes with Partial Observation Dinah Rosenberg and Eilon Solan and Nicolas Vieille April 6, 2000 Abstract We prove the existence of Blackwell ε-optimal strategies
More informationPerformance of Stochastic Programming Solutions
Performance of Stochastic Programming Solutions Operations Research Anthony Papavasiliou 1 / 30 Performance of Stochastic Programming Solutions 1 The Expected Value of Perfect Information 2 The Value of
More informationRegression estimation in continuous time with a view towards pricing Bermudan options
with a view towards pricing Bermudan options Tagung des SFB 649 Ökonomisches Risiko in Motzen 04.-06.06.2009 Financial engineering in times of financial crisis Derivate... süßes Gift für die Spekulanten
More information