A Study on the Series Expansion of Gaussian Quadratic Forms
|
|
- Jade Parker
- 6 years ago
- Views:
Transcription
1 c WISRL, KAIST, APRIL, 202 A Study on the Seres Expanson of Gaussan Quadratc Forms Techncal Report WISRL-202-APR- KAIST Juho Park, Youngchul Sung, Donggun Km, and H. Vncent Poor All rghts reserved c WISRL, KAIST Jan. 4, 200
2 c WISRL, KAIST, APRIL, Ths report s a supplementary document to the paper Outage Probablty and Outage-Based Robust Beamformng for MIMO Interference Channels wth Imperfect Channel State Informaton, by J. Park, Y. Sung, D. Km and H. V. Poor [Parket2TWC], submtted to IEEE Transactons on Wreless Communcatons. I. Dstrbuton of a Non-Central Gaussan Quadratc Form A. Prevous work and lterature survey There exst extensve lterature about the probablty dstrbuton and statstcal propertes of a quadratc form of non-central complex Gaussan random varables n the communcatons area and the probablty and statstcs communty. Through a lterature survey, we found that the man technque to compute the dstrbuton of a central or a non-central Gaussan quadratc form s based on seres fttng, whch was concretely unfed and developed by S. Kotz [Kotz-67a, Kotz-67b], and most of other works are ts varants, e.g., [Nabar-05]. Frst, we brefly explan ths seres fttng method here. Consder a Gaussan quadratc form x H Qx, where x CNµ,Σ wth sze n and Q = QH. The frst step of the seres fttng method s to convert the non-central Gaussan quadratc form nto a lnear combnaton of ch-square random varables: n x H Qx = z + δ 2 = = n [Rez + δ 2 + Imz + δ 2 ], = ndependent where z CN0, 2 for =,, n, and {δ, } are constants determned by Q, µ and Σ. Note that Rez N0, and Rez N0,. Thus, the non-central Gaussan quadratc form s equvalent to a weghted sum of non-central Ch-square random varables of whch moment generatng functon MGF s known. The MGF of a weghted sum of n ndependent non-central χ 2 random varables wth degrees of freedom 2m and non-centralty parameter µ 2 s gven by Φs = exp{ 2 n µ = n = µ 2 n 2 s } 2 s m. 2 = Note here that Φ s s nothng but the Laplace transform of the probablty densty functon PDF of x H Qx or equvalently n = z + δ 2. Now, the seres fttng method expresses the PDF as an nfnte seres composed of a set of known bass functons and tres to fnd the lnear combnaton coeffcents so that the Laplace transform of ths seres s the same as the known Jan. 4, 200
3 c WISRL, KAIST, APRIL, Fg.. Computaton of the dstrbuton of a Gaussan quadratc form Φ s. Specfcally, let the PDF be g n Q, µ,σ; y = c k h k y, 3 where {h k y, k = 0,, } s the set of known bass functons and {c k, k = 0,, } s the set of lnear combnaton coeffcents to be determned. Here, to make the problem tractable, n most cases, the followng condtons are mposed. Frst, the sequence {h k y} of bass functons s chosen among measurable complex-valued functons on [0, ] such that c k h k y Ae by, y [0, ] almost everywhere, 4 where A and b are real constants. Second, the Laplace transform ĥks of h k y has a specal form: ĥ k s = ξsη k s, 5 where ξs s a non-vanshng, analytc functon for Res > b, and ηs s analytc for Res > b and has an nverse functon. The frst condton s for the exstence of Laplace transform and the second condton s to make the problem tractable. Fnally, wth the pre-determned {h k y} Jan. 4, 200
4 c WISRL, KAIST, APRIL, wth the condtons, the coeffcents {c k } are computed so that Lg n Q, µ,σ; y = c k ĥ k s = Φ s, 6 where L denote the Laplace transform of a functon. Wdely used {h k y} for the seres expanson of the PDF of a quadratc form of non-central Gaussan random varables s as follows. [Kotz-67a, Kotz-67b]. Power seres: h k y = k y/2 n/2+k 2Γn/2+k. 2. Laguerre polynomals: Γn/2 h k y = gn; y/β[k! βγn/2 + k ]Ln/2 k y/2β, 7 where gn; y s the central χ 2 densty wth n degrees of freedom and L n/2 k x s the generalzed Laguerre polynomal defned by Rodrges formula for a > and a postve control parameter β. L n/2 k x = k! ex n/2 dk x dx k e x x k+ For the detal computaton of {c k }, please refer to [Kotz-67a, Kotz-67b, Matha-92]. The whole procedure s summarzed n Fg.. Reference group [Kotz-67a] S. Kotz, N. L. Johnson, and D. W. Boyd, Seres representaton of dstrbutons of quadratc forms n normal varables. I. Central Case, Ann. Math. Statst., vol, 38, pp , Jun [Kotz-67b] S. Kotz, N. L. Johnson, and D. W. Boyd, Seres representaton of dstrbutons of quadratc forms n normal varables. II. Non-central Case, Ann. Math. Statst., vol. 38, pp , Jun [Matha-92] A. M. Matha and S. B. Provost, Quadratc forms n random varables: Theory and applcatons, New York:M. Dekker, 992. [Nabar-05] R. Nabar, H. Bolcske, and A. Paulraj, Dversty and Outage Performance of Space- Tme Block Coded Rcean MIMO Channels, IEEE Trans. on Wreless Commun., vol. 4, no. 5, Sept Reference group 2 Jan. 4, 200
5 c WISRL, KAIST, APRIL, [Pachares-55] J. Pachares, Note on the dstrbuton of a defnte quadratc form, Ann. Math. Statst., vol. 26, pp. 28 3, Mar Power seres representaton of quadratc form of central Gaussan random varables. [Shah-6] B. K. Shah and C. G. Khatr, Dstrbuton of a defnte quadratc form for noncentral normal varates, Ann. Math. Statst., vol. 32, pp , Sep. 96. Power seres representaton of quadratc form of non-central Gaussan random varables. [Shah-63] B. K. Shah, Dstrbuton of defnte and of ndefnte quadratc forms from a noncentral normal dstrbuton, Ann. Math. Stats., vol. 34, pp , Mar Extends [Gurland-55] to derve a representaton of quadratc form of non-central Gaussan random vector wth Laguerre polynomal. Double seres of Laguerre polynomals s requred. [Gurland-55] J. Gurland, Dstrbuton of defnte and ndefnte quadratc forms, Ann. Math. Statst., vol. 26, pp , Jan Provdes a smple representaton of quadratc form of central Gaussan random vector n Laguerre polynomal. [Gurland-56] J. Gurland, Quadratc forms n normally dstrbuted random varables, Sankhya: The Indan Journal of Statstcs vol. 7, pp , Jan CDF for the ndefnte quadratc form of central random varable. [Ruben-63] H. Ruben, A new result on the dstrbuton of quadratc forms, Ann. Math. Statst., vol. 34, pp , Dec Represents the CDF of quadratc form of central and non-central Gaussan random vector wth central/non-central χ 2 dstrbuton functon. [Tku-65] M. L. Tku, Laguerre seres forms of non central χ 2 and F dstrbutons, Bometrka, vol. 52, pp , Dec Another seres representaon wth Laguerre polynomals. [Davs-77] A. W. Davs, A dfferental equaton approach to lnear combnatons of ndependent ch-squares, J. of the Ame. Statst. Assoc. vol. 72, pp , Mar Provdes another seres representaton wth power seres. [Imhof-6] J. P. Imhof, Computng the dstrbuton of quadratc forms n normal varables, Bometrka vol. 48, pp , Dec. 96. Provdes a numercal method of computng the dstrbuton [Rce-80] S. O. Rce, Dstrbuton of quadratc forms n normal varables - Evaluaton by numercal ntegraton, SIAM J. Scent. Statst. Comput., vol., no. 4, pp , 980. Another numercal method of computng dstrbuton. [Byar-93] K. H. Byar and W. C. Lndsey, Statstcal dstrbuton of Hermtan quadratc forms Jan. 4, 200
6 c WISRL, KAIST, APRIL, n complex Gaussan varables, IEEE Trans. Inform. Theory, vol. 39, pp , Mar Seres expanson of mult-varate complex Gaussan random varables. Ths paper deals wth the case that the Hermtan matrx n the quadratc form s a specal block-dagonal matrx. Reference group 3 [Raphael-96] D. Raphael, Dstrbuton of noncentral ndefnte quadratc forms n complex normal varables, IEEE Trans. Inf. Theory, vol. 42, pp , May 996. [Al-Naffour-09] T. Al-Naffour and B. Hassb, On the dstrbuton of ndefnte quadratc forms n Gaussan random varables, n Proc. of IEEE Int. Symp. Inform. Theory, Seoul, Korea, Jun. Jul B. The dfference of our work from the prevous works Frst, let us remnd our outage event n MIMO nterference channels. From equatons 5, 6 and 7 n [Parket2TWC], we have K d Pr{outage} = Pr where X mj k = j= X mjh k X mj k umh k Ĥ kk v m k 2 2 Rm k σ 2 =: τ, 8 s a non zero-mean Gaussan random varable. Note that the outage probablty s an upper tal probablty of the dstrbuton of the Gaussan quadratc form d j= XmjH k X mj k. However, as seen n Fg. 2, the most wdely-used seres fttng method explaned n the prevous subsecton yelds a good approxmaton of the dstrbuton at the lower tal not at the upper tal. The dscrepancy between the seres and the true PDF s large at the upper tal for a truncated seres. On the other hand, our approach yelds a good approxmaton to the true dstrbuton at the upper tal. Thus, the proposed seres s more relevant to our problem than the seres fttng method. Our approach to the upper tal approxmaton s based on the recent works by Raphael [Raphael-96] and by Al-Naffour and Hassb [Al-Naffour-09]. Frst, let us explan Raphael s method. The procedure n Fg. up to obtanng the MGF of the Gaussan quadratc form s common to both the sequence fttng method and Raphael s method. However, Raphael s method obtans the PDF by drect nverse Laplace transform of the MGF Φs. Typcally, the In the case of the problem consdered n [Nabar-05], the outage defned n [Nabar-05] s assocated wth the lower tal of the dstrbuton and thus the seres fttng method s well suted to that case. However, our system setup and consdered problem are dfferent from those n [Nabar-05]. Jan. 4, 200
7 c WISRL, KAIST, APRIL, [Proposed] 5 [Proposed] 0 [Laguerre] 5 [Laguerre] [Proposed] 5 [Proposed] 30 [Laguerre] 25 [Laguerre] 20 [Laguerre] Pr{Y y} [Proposed] Pr{Y y} [Proposed] Emprcal dstrbuton 20 Proposed expresson [Proposed, Laguerre] Laguerre polynomal y Emprcal dstrbuton 20 Proposed expresson [Proposed] Laguerre polynomal y a b Fg. 2. Seres fttng method versus drect nverse Laplace transform method: number of varables = 4, µ = 0.5, Q = [,0.5,0,0;0.5,,0,0;0,0,,0;0,0,0,], and Σ = 0.3I. a β = and b β = 2. β s the control parameter for the Laguerre polynomals n 7. Note that the convergent speed of the seres fttng method based on the Laguerre polynomals depends much on β. In the case of β = 2, the seres fttng method based on the Laguerre polynomals yelds large errors at the upper tal. It s not smple how to choose β and an effcent method s not known. One cannot run smulatons for emprcal dstrbutons for all cases. The seres fttng method based on the power seres shows bad performance, and t cannot be used n practce. nverse Laplace transform of the MGF s represented as a complex contour ntegral and then the complex contour ntegral s computed as an nfnte seres by the resdue theorem. However, to obtan the cumulatve dstrbuton functon CDF, whch s actually necessary to compute the tal probablty, Raphael s method requres one more step, the ntegraton of the PDF, to obtan the CDF snce the MGF Φs s the Laplace transform of the PDF. In [Parket2TWC], to obtan the CDF of a general Gaussan quadratc form, we dd not use the MGF Φs, whch s a bt complcated and requres an addtonal step, lke Raphael, but nstead we drectly used a smple contour ntegral for the CDF, eq. 2 n [Parket2TWC], obtaned by Al-Naffour and Hassb [Al-Naffour-09]. Then, the contour ntegral was computed as an In [Al-Naffour-09], Al-Naffour and Hassb obtaned the contour ntegral, eq. 2 n [Parket2TWC] for the CDF of a Gaussan quadratc form. However, they dd not obtan closed-form seres expressons for the contour ntegral n general cases except a few smple cases. The man goal of [Al-Naffour-09] was to derve a nce and smple contour ntegral form for the CDF. Jan. 4, 200
8 c WISRL, KAIST, APRIL, nfnte seres by the resdue theorem. Usng the resdue theorem s borrowed from Raphael s work. Thus, our result s smpler than Raphael s approach and does not requre the ntegraton of a PDF for the CDF. As mentoned already, the seres expanson n [Parket2TWC] has a partcular advantage over the seres fttng method consdered n [Nabar-05] for the outage event defned n [Parket2TWC]; The seres n [Parket2TWC] fts the upper tal of the dstrbuton well wth a few number of terms. We shall provde a detaled proof for ths n a specal case n the next subsecton. Thus, our seres expressons for outage probablty n MIMO nterference channels are meanngful and relevant. II. Computatonal Issues and Convergence of the Obtaned Seres A. Computng hgher order dervatves Recall the general outage expresson n Theorem n [Parket2TWC]: Pr{outage} = Pr{log 2 + SINR m k R m κ = = e τ + κ j= χj 2 k } n=κ n! gn 0 n κ +! κ j= χj 2 n κ + 9, where g s = eτs exp s / To compute 9, we need to compute s / p κp p +s / p q= χq p 2 κp. 0 + s / p p { } the egenvalues of the Kd Kd covarance matrx Σ = ΨΛΨ H, {χ j } the elements of Kd vector χ = Λ /2 Ψ H µ, where µ s the mean vector of the Gaussan dstrbuton, and the hgher order dervatves of g s. The computaton of { } and {χ j } s smple snce the szes of the mean vector and the covarance matrx are Kd and Kd Kd, respectvely. Furthermore, the hgher order dervatves of g s can also be computed effcently based on recurson [Matha-92],[Raphael-96]. Note that g s = Jan. 4, 200
9 c WISRL, KAIST, APRIL, e log g s. Thus, the dervatve of g s can be wrtten as g s = g s[log g s], g 2 s = g s[log g s] + g s[log g s] 2,. g n s = n n g l s[log g s] n l, n l l=0 where g l s and [log g s] l denote the l-th dervatves of g s and log g s, respectvely. Here, [log g s] n can be computed from 0 as [log g s] n n! n = τδ n s / n p n! n n p + p s / n+ κ p q= χ q p 2 n! n κ p n p + p s / n p where δ n s Kronecker delta functon. Thus, for gven g s and [log g s] l, we can compute g l s effcently n a recursve way, as shown n. B. Convergence analyss In ths subsecton, we provde some convergence analyss on the derved seres expanson n [Parket2TWC]. Consder the general result n Theorem of [Parket2TWC] for the CDF of a Gaussan quadratc form: Pr{Y y} = + where κ = e y + κ j= χj 2 g s, y = esy s n=κ exp n! gn 0, y n κ +! s / p κp p +s / p q= χq p 2 κp. + s / p p κ j= χj 2 n κ + Here, we explctly use the varable y as an nput parameter of the functon g s for later explanaton. g n s, y denotes the n-th partal dervatve of g s, y wth respect to s. Here, κ s the number of dstnct egenvalues of the Kd Kd covarance matrx Σ and κ s the geometrc order of egenvalue. κ = κ = Kd. The resdual error caused by truncatng the nfnte seres after the frst N terms s gven by 2 R N y = κ = e y + κ j= χj 2 n! gn 0, y n κ +! κ j= χj 2 n κ +, 3 Jan. 4, 200
10 c WISRL, KAIST, APRIL, and we have Pr{Y y; nfnte sum} = Pr{Y y; truncaton at N} + R N y. The truncaton error R N y can be expressed as κ R N y = RNy, 4 where R Ny = e y + κ j= χj 2 = n! gn 0, y n κ +! κ j= χj 2 n κ + for each κ. Then, the magntude of each term RN y n the truncaton error s bounded as RNy { y exp + κ j= χ j 2 } g n 0, y n! 5 κ j= χj 2 n κ +. n κ +! As seen n Fg. 2, our seres expanson fts the upper tal dstrbuton frst. Now, to assess the overall convergence speed of our seres, for the same step as n Fg. 2, we ran some smulatons to obtan an emprcal dstrbuton, and computed the overall mean square error MSE between the truncated seres and the emprcal dstrbuton over 0 y 0 as CDF MSE = Pr{Y y ; N, type of seres} Pr{Y y ; emprcal} 2, = where {y } are the unform samples of [0, 0]. Fg. 3 shows the CDF MSE of the three methods n Fg. 2: the proposed seres, the seres fttng method wth β = and the seres fttng method wth β = 2. It s seen n Fg. 3 that the overall convergence of the proposed seres can be worse than the seres fttng method at the small values for the number of summaton terms for the settng n Fg. 2. The bad overall convergence s due to worse fttng at the lower tal of the dstrbuton, but the bad lower tal approxmaton s not mportant to our outage computaton. Please see Fg. 2. Fg. 4 shows another case. In ths case, the proposed seres outperforms the seres fttng method both n the overall convergence and n the upper tal convergence. It s seen numercally that the proposed seres fts the upper tal dstrbuton frst. Now, we shall prove ths property of the proposed seres. However, t s a dffcult problem to prove ths property n general cases. Thus, n the next subsecton, we provde a proof of ths property when the number of dstnct egenvalues of the covarance matrx Σ s one, e.g., n the..d. case. 6 Jan. 4, 200
11 c WISRL, KAIST, APRIL, CDF MSE Proposed expresson Laguerre polynomal β= Laguerre polynomal β= Order Fg. 3. CDF MSE of the CDFs n Fg [Laguerre] 0 0 Proposed expresson Laguerre polynomal β= [Laguerre] Pr{Y y} [Proposed] CDF MSE [Proposed] Emprcal dstrbuton 0.4 Proposed expresson Laguerre polynomal y Order a b Fg. 4. number of varables = 4, µ = 0.5, Q = I, and Σ = [ ; ; ; ]. In ths case egenvalues are.0000, 0.638, 0.258, and wth β =. a CDF, b CDF MSE. Unform sample of y s taken over [0,5.9]. B. The dentty covarance matrx case Suppose that there s only one egenvalue, > 0, wth multplcty κ for the covarance matrx Σ. Ths case corresponds to Corollary 4 n [Parket2TWC], and the outage probablty Jan. 4, 200
12 c WISRL, KAIST, APRIL, s gven by where Pr{Y y} = + exp η2 exp y g n 0, y η2 / n κ+ n!n κ +!, 7 n=κ gs, y = eys s 8 and η 2 = κ j= χj 2. The resdual error caused by truncatng the nfnte seres after the frst N terms s gven by R N y = exp η2 exp y g n 0, y η2 / n κ+ n!n κ +!. 9 Before we proceed, we frst obtan the n-th dervatve of gs, y at s = 0, whch s gven n the followng lemma. Lemma : For n 0, g n 0, y = n n! n k! k y n k. 20 Proof: Proof s gven by nducton. The valdty of the clam for n = 0, and 2 s shown by drecton computaton: g 0 0, y = = yeys 0 0! s / = = s=0 0 k! k y 0 k, g 0, y = = yeys s / e ys! s / 2 = y + = s=0 k! k y k, g 2 0, y = yeys ys y/ + e ys y s 2 2e ys ys y/ s s / 4 2 = y 2 + 2y ! = 2 k! k y 2 k. Now, suppose that 20 holds up to the n -th dervatve of gs, y. From the recursve formula n, g n 0, y s obtaned as g n 0,y n n = k n = 0 g k 0,ylog g0,y n k g 0 0,ylog g0,y n + n g 0,ylog g0,y n + + s=0 n g n 0,ylog g0,y. n 2 Jan. 4, 200
13 c WISRL, KAIST, APRIL, Snce [log gs] = ys logs /, we can easly see that [log g0] = y + and [log g0] n = n! n for n 2. Therefore, 2 can be rewrtten as n g n 0, y =n!g0, y n + n g 0, yn 2! n + g 2 0, yn 3! n n g n 2 0, y 2 + g n 0, yy + =n!g0, y n + n!g 0, y n + + n g n 2 0, y 2 + g n 0, y + yg n 0, y [ n l n n! l! l! l k! k y l k n l +y l=0 m=0 [ n l n a = = l=0 n! l! l! l k! k y l k n l + n! g 2 0, y n m=0 n! n m! m y n m n! n m! m y n m ] ] 22 where a holds snce 20 holds for all g 0 0, y,, g n 0, y by the nducton assumpton. Here, consder the coeffcent of each y n 22 for = 0,, n. y n s obtaned only when m = 0. The coeffcent of y n from 22 s therefore gven by. It corresponds to the coeffcent of y n n 20. For 0 < p n, the coeffcent of y n p s obtaned by consderng all l, k that satsfes l k = n p due to the frst term n the rght-hand sde RHS of 22, and m = p due to the second term of the RHS of 22. In the frst case, we obtan y n p wth the followng pars l, k = n, p, n 2, p 2,, n p, 0. For these l, k pars, we have n l=n p n! l! l! n p! l n+p y n p n l = In the second case of m = p, we have n l=n p n! n p! p y n p n! = p n p! p y n p. 23 n! n p! p y n p. 24 Fnally, the coeffcent of y n q s gven by addng 23 and 24: n!! + pn p y n p n p! n p! n! = + p p y n p n p! n p n! = n p! p y n p, Jan. 4, 200
14 c WISRL, KAIST, APRIL, whch s equvalent to the coeffcent for y n p n 20 0 < p n. Thus, 20 holds for g n 0, y. Note that g n 0, y < 0 for all n 0 from 20. Therefore, R N y 0 for all N and y and g n 0, y = g n 0, y. Now, consder the resdual error term R N y n 9. The magntude of the resdual error can be upper bounded as follows: R N y = exp η2 = exp η2 = exp µ2 = exp η2 exp y g n 0, y η2 / n κ+ n!n κ +! exp y g n 0, y η2 / n κ+ n!n κ +! exp y exp y = 2 κ exp η2 exp g n 0, y η2 / n κ+ n!n κ +! n2η 2 n κ+ 2 κ n! gn 0, y 2 n κ +! y a 2 κ exp η2 exp y = 2κ expη2 exp y b 2κ expη2 exp y c = 2κ expη2 exp y g d = 2κ expη2 exp y = 2 κ expη 2 exp y 2 n 2η 2 n κ+ n! gn 0, y 2 n κ +! n n! gn 0, y exp2η 2 2 n! gn 0, y 2 n! gn 0, y 2 n=0 2, y expy/2 /2 n n 25 where a s from γk k! expγ = p=0 γp /p! for any γ > 0, b s from the fact that summand s negatve, c s by usng the Taylor seres expanson, and d s from 8. Snce η s a fxed constant, from 25, for any N 0 lm R Ny = y Thus, t s clear that the proposed seres converges from the upper tal dstrbuton! Jan. 4, 200
15 c WISRL, KAIST, APRIL, Now, let us consder the resdual error magntude as a functon of y for gven N. From 20, we have g n 0, y y Dfferentatng R N y wth respect to y yelds R N y y = exp η2 + exp η2 = exp η2 = ng n 0, y. 27 exp y g n 0, y η2 / n κ+ n!n κ +! exp y g n 0, y y exp y Furthermore, from 20 we have η 2 / n κ+ n!n κ +! η 2 / n κ+ n!n κ +! gn 0, y + ng n 0, y. 28 gn 0, y + ng n 0, y = y n. 29 By substtutng 29 nto 28, we have R N y y = exp η2 exp y η 2 / n κ+ y n n!n κ +!, 30 whch s postve. Snce R N y 0, lm y R N y = 0 and R Ny y > 0, the resdual error magntude monotoncally decreases as y ncreases and the maxmum error occurs at y = 0 for any gven N. Now, let us compute the worst truncaton error R N 0, whch s gven by R N 0 = exp η2 From 20, we have g n 0, 0 = n! n+. Therefore, R N 0 = exp η2 = exp η2 = exp η2 = exp η 2 g n 0, 0 η2 / n κ+ n!n κ +!. 3 n! n+ η2 / n κ+ n!n κ +! n+η2 / n κ+ n κ +! η 2 n κ+ n κ +! κ η 2 n κ+ n κ +!. 32 Jan. 4, 200
16 c WISRL, KAIST, APRIL, From 7, N κ 2. For general N κ 2, let m = n κ +. Then, R N 0 = exp η 2 m=n κ+2 η 2 m. m! Note that η 2 m m=n κ+2 m! s the resdual error of the Taylor seres expanson of expx after the frst N κ + terms. By the Taylor theorem, m=n κ+2 η 2 m m! = η2 N κ+2 N κ + 2! expαη2 33 where some α [0, ]. Therefore, the worst truncaton error s gven by R N 0 = exp α η 2 η2 N κ+2 N κ + 2! η2 N κ+2 N κ + 2!, 34 where the nequalty holds snce expα η 2 for 0 α. Furthermore, the resdual error magntude s a strctly decreasng functon of N for any y, Ths can be shown easly as follows. R N y = exp η2 = exp η2 =R N+ y + exp η2 R N y > R N+ y. 35 exp y g n 0, y η2 / n κ+ n!n κ +! exp y { g n 0, y η2 / n κ+ n!n κ +! + gn+ 0, y n=n+2 exp y g N+ η 2 / N κ+2 0, y N +!N κ + 2!. } η 2 / N κ+2 N +!N κ + 2! Snce R N y < 0 and g N+ y < 0 for all y 0 and N, we have 35. Now, based on 34 and 35, wth gven χ k and σh 2, we can compute the requred number N of terms n the seres to acheve the desred level of accuracy snce η 2 s known. Fnally, consder the worst case of N = κ 2 and y = 0: R κ 2 0 = exp η 2 η 2 n κ+ n κ +! = exp η2 n=κ m=0 η 2 m m! =, where the second equalty s by replacng m = n κ +. It s easy to see that the worst case error s - n the dentty covarance matrx case. Fg. 5 shows the performance of the proposed seres expanson n the case of the dentty covarance matrx. The numercal results well match our theoretcal analyss n ths subsecton. From the fgure, t seems reasonable to choose N for accurate outage probablty computaton. Jan. 4, 200
17 c WISRL, KAIST, APRIL, N=5 0 N=20 Pr{Y y} N=0 N=5 0.2 Emprcal dstrbuton N=20 Proposed expresson y RNy 0.2 N= N=0 0.8 N=5 0 2 y 3 4 a CDF b Resdual error Fg. 5. number of varables = 4, µ = 0.5, Q = I, and Σ = 0.I. Jan. 4, 200
Appendix for Solving Asset Pricing Models when the Price-Dividend Function is Analytic
Appendx for Solvng Asset Prcng Models when the Prce-Dvdend Functon s Analytc Ovdu L. Caln Yu Chen Thomas F. Cosmano and Alex A. Hmonas January 3, 5 Ths appendx provdes proofs of some results stated n our
More informationAppendix - Normally Distributed Admissible Choices are Optimal
Appendx - Normally Dstrbuted Admssble Choces are Optmal James N. Bodurtha, Jr. McDonough School of Busness Georgetown Unversty and Q Shen Stafford Partners Aprl 994 latest revson September 00 Abstract
More informationMultifactor Term Structure Models
1 Multfactor Term Structure Models A. Lmtatons of One-Factor Models 1. Returns on bonds of all maturtes are perfectly correlated. 2. Term structure (and prces of every other dervatves) are unquely determned
More informationRandom Variables. b 2.
Random Varables Generally the object of an nvestgators nterest s not necessarly the acton n the sample space but rather some functon of t. Techncally a real valued functon or mappng whose doman s the sample
More information3/3/2014. CDS M Phil Econometrics. Vijayamohanan Pillai N. Truncated standard normal distribution for a = 0.5, 0, and 0.5. CDS Mphil Econometrics
Lmted Dependent Varable Models: Tobt an Plla N 1 CDS Mphl Econometrcs Introducton Lmted Dependent Varable Models: Truncaton and Censorng Maddala, G. 1983. Lmted Dependent and Qualtatve Varables n Econometrcs.
More informationMgtOp 215 Chapter 13 Dr. Ahn
MgtOp 5 Chapter 3 Dr Ahn Consder two random varables X and Y wth,,, In order to study the relatonshp between the two random varables, we need a numercal measure that descrbes the relatonshp The covarance
More information15-451/651: Design & Analysis of Algorithms January 22, 2019 Lecture #3: Amortized Analysis last changed: January 18, 2019
5-45/65: Desgn & Analyss of Algorthms January, 09 Lecture #3: Amortzed Analyss last changed: January 8, 09 Introducton In ths lecture we dscuss a useful form of analyss, called amortzed analyss, for problems
More information3: Central Limit Theorem, Systematic Errors
3: Central Lmt Theorem, Systematc Errors 1 Errors 1.1 Central Lmt Theorem Ths theorem s of prme mportance when measurng physcal quanttes because usually the mperfectons n the measurements are due to several
More informationIntroduction to PGMs: Discrete Variables. Sargur Srihari
Introducton to : Dscrete Varables Sargur srhar@cedar.buffalo.edu Topcs. What are graphcal models (or ) 2. Use of Engneerng and AI 3. Drectonalty n graphs 4. Bayesan Networks 5. Generatve Models and Samplng
More informationParallel Prefix addition
Marcelo Kryger Sudent ID 015629850 Parallel Prefx addton The parallel prefx adder presented next, performs the addton of two bnary numbers n tme of complexty O(log n) and lnear cost O(n). Lets notce the
More informationCOS 511: Theoretical Machine Learning. Lecturer: Rob Schapire Lecture #21 Scribe: Lawrence Diao April 23, 2013
COS 511: Theoretcal Machne Learnng Lecturer: Rob Schapre Lecture #21 Scrbe: Lawrence Dao Aprl 23, 2013 1 On-Lne Log Loss To recap the end of the last lecture, we have the followng on-lne problem wth N
More informationTCOM501 Networking: Theory & Fundamentals Final Examination Professor Yannis A. Korilis April 26, 2002
TO5 Networng: Theory & undamentals nal xamnaton Professor Yanns. orls prl, Problem [ ponts]: onsder a rng networ wth nodes,,,. In ths networ, a customer that completes servce at node exts the networ wth
More informationECE 586GT: Problem Set 2: Problems and Solutions Uniqueness of Nash equilibria, zero sum games, evolutionary dynamics
Unversty of Illnos Fall 08 ECE 586GT: Problem Set : Problems and Solutons Unqueness of Nash equlbra, zero sum games, evolutonary dynamcs Due: Tuesday, Sept. 5, at begnnng of class Readng: Course notes,
More informationFoundations of Machine Learning II TP1: Entropy
Foundatons of Machne Learnng II TP1: Entropy Gullaume Charpat (Teacher) & Gaétan Marceau Caron (Scrbe) Problem 1 (Gbbs nequalty). Let p and q two probablty measures over a fnte alphabet X. Prove that KL(p
More information4. Greek Letters, Value-at-Risk
4 Greek Letters, Value-at-Rsk 4 Value-at-Rsk (Hull s, Chapter 8) Math443 W08, HM Zhu Outlne (Hull, Chap 8) What s Value at Rsk (VaR)? Hstorcal smulatons Monte Carlo smulatons Model based approach Varance-covarance
More informationSupplementary material for Non-conjugate Variational Message Passing for Multinomial and Binary Regression
Supplementary materal for Non-conjugate Varatonal Message Passng for Multnomal and Bnary Regresson October 9, 011 1 Alternatve dervaton We wll focus on a partcular factor f a and varable x, wth the am
More informationEDC Introduction
.0 Introducton EDC3 In the last set of notes (EDC), we saw how to use penalty factors n solvng the EDC problem wth losses. In ths set of notes, we want to address two closely related ssues. What are, exactly,
More information/ Computational Genomics. Normalization
0-80 /02-70 Computatonal Genomcs Normalzaton Gene Expresson Analyss Model Computatonal nformaton fuson Bologcal regulatory networks Pattern Recognton Data Analyss clusterng, classfcaton normalzaton, mss.
More informationLinear Combinations of Random Variables and Sampling (100 points)
Economcs 30330: Statstcs for Economcs Problem Set 6 Unversty of Notre Dame Instructor: Julo Garín Sprng 2012 Lnear Combnatons of Random Varables and Samplng 100 ponts 1. Four-part problem. Go get some
More informationII. Random Variables. Variable Types. Variables Map Outcomes to Numbers
II. Random Varables Random varables operate n much the same way as the outcomes or events n some arbtrary sample space the dstncton s that random varables are smply outcomes that are represented numercally.
More information4.4 Doob s inequalities
34 CHAPTER 4. MARTINGALES 4.4 Doob s nequaltes The frst nterestng consequences of the optonal stoppng theorems are Doob s nequaltes. If M n s a martngale, denote M n =max applen M. Theorem 4.8 If M n s
More informationCorrelations and Copulas
Correlatons and Copulas Chapter 9 Rsk Management and Fnancal Insttutons, Chapter 6, Copyrght John C. Hull 2006 6. Coeffcent of Correlaton The coeffcent of correlaton between two varables V and V 2 s defned
More informationTests for Two Correlations
PASS Sample Sze Software Chapter 805 Tests for Two Correlatons Introducton The correlaton coeffcent (or correlaton), ρ, s a popular parameter for descrbng the strength of the assocaton between two varables.
More information1 A Primer on Linear Models. 2 Chapter 1 corrections. 3 Chapter 2 corrections. 4 Chapter 3 corrections. 1.1 Corrections 23 May 2015
A Prmer on Lnear Models. Correctons 23 May 205 2 Chapter correctons Fx page: 9 lne -7 denomnator ( ρ 2 ) s mssng, should read Cov(e t, e s ) = V ar(a t )/( ρ 2 ) Fx page: lne -5 sgn needs changng: cos(a
More informationOn the Moments of the Traces of Unitary and Orthogonal Random Matrices
Proceedngs of Insttute of Mathematcs of NAS of Ukrane 2004 Vol. 50 Part 3 1207 1213 On the Moments of the Traces of Untary and Orthogonal Random Matrces Vladmr VASILCHU B. Verkn Insttute for Low Temperature
More informationarxiv: v1 [math.nt] 29 Oct 2015
A DIGITAL BINOMIAL THEOREM FOR SHEFFER SEQUENCES TOUFIK MANSOUR AND HIEU D. NGUYEN arxv:1510.08529v1 [math.nt] 29 Oct 2015 Abstract. We extend the dgtal bnomal theorem to Sheffer polynomal sequences by
More informationFast Laplacian Solvers by Sparsification
Spectral Graph Theory Lecture 19 Fast Laplacan Solvers by Sparsfcaton Danel A. Spelman November 9, 2015 Dsclamer These notes are not necessarly an accurate representaton of what happened n class. The notes
More informationCHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS
CHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS QUESTIONS 9.1. (a) In a log-log model the dependent and all explanatory varables are n the logarthmc form. (b) In the log-ln model the dependent varable
More informationProblem Set 6 Finance 1,
Carnege Mellon Unversty Graduate School of Industral Admnstraton Chrs Telmer Wnter 2006 Problem Set 6 Fnance, 47-720. (representatve agent constructon) Consder the followng two-perod, two-agent economy.
More informationASSESSING GOODNESS OF FIT OF GENERALIZED LINEAR MODELS TO SPARSE DATA USING HIGHER ORDER MOMENT CORRECTIONS
ASSESSING GOODNESS OF FIT OF GENERALIZED LINEAR MODELS TO SPARSE DATA USING HIGHER ORDER MOMENT CORRECTIONS S. R. PAUL Department of Mathematcs & Statstcs, Unversty of Wndsor, Wndsor, ON N9B 3P4, Canada
More informationInformation Flow and Recovering the. Estimating the Moments of. Normality of Asset Returns
Estmatng the Moments of Informaton Flow and Recoverng the Normalty of Asset Returns Ané and Geman (Journal of Fnance, 2000) Revsted Anthony Murphy, Nuffeld College, Oxford Marwan Izzeldn, Unversty of Lecester
More informationChapter 5 Student Lecture Notes 5-1
Chapter 5 Student Lecture Notes 5-1 Basc Busness Statstcs (9 th Edton) Chapter 5 Some Important Dscrete Probablty Dstrbutons 004 Prentce-Hall, Inc. Chap 5-1 Chapter Topcs The Probablty Dstrbuton of a Dscrete
More informationComparison of Singular Spectrum Analysis and ARIMA
Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 0, Dubln (Sesson CPS009) p.99 Comparson of Sngular Spectrum Analss and ARIMA Models Zokae, Mohammad Shahd Behesht Unverst, Department of Statstcs
More informationTests for Two Ordered Categorical Variables
Chapter 253 Tests for Two Ordered Categorcal Varables Introducton Ths module computes power and sample sze for tests of ordered categorcal data such as Lkert scale data. Assumng proportonal odds, such
More informationGraphical Methods for Survival Distribution Fitting
Graphcal Methods for Survval Dstrbuton Fttng In ths Chapter we dscuss the followng two graphcal methods for survval dstrbuton fttng: 1. Probablty Plot, 2. Cox-Snell Resdual Method. Probablty Plot: The
More informationLikelihood Fits. Craig Blocker Brandeis August 23, 2004
Lkelhood Fts Crag Blocker Brandes August 23, 2004 Outlne I. What s the queston? II. Lkelhood Bascs III. Mathematcal Propertes IV. Uncertantes on Parameters V. Mscellaneous VI. Goodness of Ft VII. Comparson
More informationUnderstanding Annuities. Some Algebraic Terminology.
Understandng Annutes Ma 162 Sprng 2010 Ma 162 Sprng 2010 March 22, 2010 Some Algebrac Termnology We recall some terms and calculatons from elementary algebra A fnte sequence of numbers s a functon of natural
More informationPrice and Quantity Competition Revisited. Abstract
rce and uantty Competton Revsted X. Henry Wang Unversty of Mssour - Columba Abstract By enlargng the parameter space orgnally consdered by Sngh and Vves (984 to allow for a wder range of cost asymmetry,
More informationS yi a bx i cx yi a bx i cx 2 i =0. yi a bx i cx 2 i xi =0. yi a bx i cx 2 i x
LEAST-SQUARES FIT (Chapter 8) Ft the best straght lne (parabola, etc.) to a gven set of ponts. Ths wll be done by mnmzng the sum of squares of the vertcal dstances (called resduals) from the ponts to the
More informationMaximum Likelihood Estimation of Isotonic Normal Means with Unknown Variances*
Journal of Multvarate Analyss 64, 183195 (1998) Artcle No. MV971717 Maxmum Lelhood Estmaton of Isotonc Normal Means wth Unnown Varances* Nng-Zhong Sh and Hua Jang Northeast Normal Unversty, Changchun,Chna
More informationDependent jump processes with coupled Lévy measures
Dependent jump processes wth coupled Lévy measures Naoufel El-Bachr ICMA Centre, Unversty of Readng May 6, 2008 ICMA Centre Dscusson Papers n Fnance DP2008-3 Copyrght 2008 El-Bachr. All rghts reserved.
More information2.1 Rademacher Calculus... 3
COS 598E: Unsupervsed Learnng Week 2 Lecturer: Elad Hazan Scrbe: Kran Vodrahall Contents 1 Introducton 1 2 Non-generatve pproach 1 2.1 Rademacher Calculus............................... 3 3 Spectral utoencoders
More informationECONOMETRICS - FINAL EXAM, 3rd YEAR (GECO & GADE)
ECONOMETRICS - FINAL EXAM, 3rd YEAR (GECO & GADE) May 17, 2016 15:30 Frst famly name: Name: DNI/ID: Moble: Second famly Name: GECO/GADE: Instructor: E-mal: Queston 1 A B C Blank Queston 2 A B C Blank Queston
More informationA Bootstrap Confidence Limit for Process Capability Indices
A ootstrap Confdence Lmt for Process Capablty Indces YANG Janfeng School of usness, Zhengzhou Unversty, P.R.Chna, 450001 Abstract The process capablty ndces are wdely used by qualty professonals as an
More informationMidterm Exam. Use the end of month price data for the S&P 500 index in the table below to answer the following questions.
Unversty of Washngton Summer 2001 Department of Economcs Erc Zvot Economcs 483 Mdterm Exam Ths s a closed book and closed note exam. However, you are allowed one page of handwrtten notes. Answer all questons
More informationMode is the value which occurs most frequency. The mode may not exist, and even if it does, it may not be unique.
1.7.4 Mode Mode s the value whch occurs most frequency. The mode may not exst, and even f t does, t may not be unque. For ungrouped data, we smply count the largest frequency of the gven value. If all
More informationEvaluating Performance
5 Chapter Evaluatng Performance In Ths Chapter Dollar-Weghted Rate of Return Tme-Weghted Rate of Return Income Rate of Return Prncpal Rate of Return Daly Returns MPT Statstcs 5- Measurng Rates of Return
More informationThe convolution computation for Perfectly Matched Boundary Layer algorithm in finite differences
The convoluton computaton for Perfectly Matched Boundary Layer algorthm n fnte dfferences Herman Jaramllo May 10, 2016 1 Introducton Ths s an exercse to help on the understandng on some mportant ssues
More informationComparative analysis of CDO pricing models
Comparatve analyss of CDO prcng models ICBI Rsk Management 2005 Geneva 8 December 2005 Jean-Paul Laurent ISFA, Unversty of Lyon, Scentfc Consultant BNP Parbas laurent.jeanpaul@free.fr, http://laurent.jeanpaul.free.fr
More informationNotes on Debye-Hückel Theory
GEOS 0300/30300 Autumn 004 Thermodynamcs and Phase Change MS Ghorso Notes on Debye-Hückel Theory We seek: µ µ o + RT ln X + µ ex where µ ex RT ln λ or the devaton from deal behavor due to long range nteracton
More informationAlternatives to Shewhart Charts
Alternatves to Shewhart Charts CUSUM & EWMA S Wongsa Overvew Revstng Shewhart Control Charts Cumulatve Sum (CUSUM) Control Chart Eponentally Weghted Movng Average (EWMA) Control Chart 2 Revstng Shewhart
More informationExplaining Movements of the Labor Share in the Korean Economy: Factor Substitution, Markups and Bargaining Power
Explanng Movements of the abor Share n the Korean Economy: Factor Substtuton, Markups and Barganng ower Bae-Geun, Km January 2, 26 Appendx A. Dervaton of the dervatve of et us start from eq. (). For notatonal
More informationAnalysis of Variance and Design of Experiments-II
Analyss of Varance and Desgn of Experments-II MODULE VI LECTURE - 4 SPLIT-PLOT AND STRIP-PLOT DESIGNS Dr. Shalabh Department of Mathematcs & Statstcs Indan Insttute of Technology Kanpur An example to motvate
More informationLecture 7. We now use Brouwer s fixed point theorem to prove Nash s theorem.
Topcs on the Border of Economcs and Computaton December 11, 2005 Lecturer: Noam Nsan Lecture 7 Scrbe: Yoram Bachrach 1 Nash s Theorem We begn by provng Nash s Theorem about the exstance of a mxed strategy
More informationEconomics 1410 Fall Section 7 Notes 1. Define the tax in a flexible way using T (z), where z is the income reported by the agent.
Economcs 1410 Fall 2017 Harvard Unversty Yaan Al-Karableh Secton 7 Notes 1 I. The ncome taxaton problem Defne the tax n a flexble way usng T (), where s the ncome reported by the agent. Retenton functon:
More informationDOUBLE IMPACT. Credit Risk Assessment for Secured Loans. Jean-Paul Laurent ISFA Actuarial School University of Lyon & BNP Paribas
DOUBLE IMPACT Credt Rsk Assessment for Secured Loans Al Chabaane BNP Parbas Jean-Paul Laurent ISFA Actuaral School Unversty of Lyon & BNP Parbas Julen Salomon BNP Parbas julen.salomon@bnpparbas.com Abstract
More informationRaising Food Prices and Welfare Change: A Simple Calibration. Xiaohua Yu
Rasng Food Prces and Welfare Change: A Smple Calbraton Xaohua Yu Professor of Agrcultural Economcs Courant Research Centre Poverty, Equty and Growth Unversty of Göttngen CRC-PEG, Wlhelm-weber-Str. 2 3773
More informationEconomic Design of Short-Run CSP-1 Plan Under Linear Inspection Cost
Tamkang Journal of Scence and Engneerng, Vol. 9, No 1, pp. 19 23 (2006) 19 Economc Desgn of Short-Run CSP-1 Plan Under Lnear Inspecton Cost Chung-Ho Chen 1 * and Chao-Yu Chou 2 1 Department of Industral
More informationOption pricing and numéraires
Opton prcng and numérares Daro Trevsan Unverstà degl Stud d Psa San Mnato - 15 September 2016 Overvew 1 What s a numerare? 2 Arrow-Debreu model Change of numerare change of measure 3 Contnuous tme Self-fnancng
More informationApplications of Myerson s Lemma
Applcatons of Myerson s Lemma Professor Greenwald 28-2-7 We apply Myerson s lemma to solve the sngle-good aucton, and the generalzaton n whch there are k dentcal copes of the good. Our objectve s welfare
More informationoccurrence of a larger storm than our culvert or bridge is barely capable of handling? (what is The main question is: What is the possibility of
Module 8: Probablty and Statstcal Methods n Water Resources Engneerng Bob Ptt Unversty of Alabama Tuscaloosa, AL Flow data are avalable from numerous USGS operated flow recordng statons. Data s usually
More informationPrinciples of Finance
Prncples of Fnance Grzegorz Trojanowsk Lecture 6: Captal Asset Prcng Model Prncples of Fnance - Lecture 6 1 Lecture 6 materal Requred readng: Elton et al., Chapters 13, 14, and 15 Supplementary readng:
More informationQuiz on Deterministic part of course October 22, 2002
Engneerng ystems Analyss for Desgn Quz on Determnstc part of course October 22, 2002 Ths s a closed book exercse. You may use calculators Grade Tables There are 90 ponts possble for the regular test, or
More informationarxiv: v2 [math.co] 6 Apr 2016
On the number of equvalence classes of nvertble Boolean functons under acton of permutaton of varables on doman and range arxv:1603.04386v2 [math.co] 6 Apr 2016 Marko Carć and Modrag Žvkovć Abstract. Let
More informationNotes on experimental uncertainties and their propagation
Ed Eyler 003 otes on epermental uncertantes and ther propagaton These notes are not ntended as a complete set of lecture notes, but nstead as an enumeraton of some of the key statstcal deas needed to obtan
More informationBasket options and implied correlations: a closed form approach
Basket optons and mpled correlatons: a closed form approach Svetlana Borovkova Free Unversty of Amsterdam CFC conference, London, January 7-8, 007 Basket opton: opton whose underlyng s a basket (.e. a
More informationMeasures of Spread IQR and Deviation. For exam X, calculate the mean, median and mode. For exam Y, calculate the mean, median and mode.
Part 4 Measures of Spread IQR and Devaton In Part we learned how the three measures of center offer dfferent ways of provdng us wth a sngle representatve value for a data set. However, consder the followng
More information>1 indicates country i has a comparative advantage in production of j; the greater the index, the stronger the advantage. RCA 1 ij
69 APPENDIX 1 RCA Indces In the followng we present some maor RCA ndces reported n the lterature. For addtonal varants and other RCA ndces, Memedovc (1994) and Vollrath (1991) provde more thorough revews.
More informationISyE 512 Chapter 9. CUSUM and EWMA Control Charts. Instructor: Prof. Kaibo Liu. Department of Industrial and Systems Engineering UW-Madison
ISyE 512 hapter 9 USUM and EWMA ontrol harts Instructor: Prof. Kabo Lu Department of Industral and Systems Engneerng UW-Madson Emal: klu8@wsc.edu Offce: Room 317 (Mechancal Engneerng Buldng) ISyE 512 Instructor:
More informationCentre for International Capital Markets
Centre for Internatonal Captal Markets Dscusson Papers ISSN 1749-3412 Valung Amercan Style Dervatves by Least Squares Methods Maro Cerrato No 2007-13 Valung Amercan Style Dervatves by Least Squares Methods
More informationProbability distribution of multi-hop-distance in one-dimensional sensor networks q
Computer etworks (7) 77 79 www.elsever.com/locate/comnet Probablty dstrbuton of mult-hop-dstance n one-dmensonal sensor networks q Serdar Vural *, Eylem Ekc Department of Electrcal and Computer Engneerng,
More informationReal Exchange Rate Fluctuations, Wage Stickiness and Markup Adjustments
Real Exchange Rate Fluctuatons, Wage Stckness and Markup Adjustments Yothn Jnjarak and Kanda Nakno Nanyang Technologcal Unversty and Purdue Unversty January 2009 Abstract Motvated by emprcal evdence on
More informationChapter 3 Descriptive Statistics: Numerical Measures Part B
Sldes Prepared by JOHN S. LOUCKS St. Edward s Unversty Slde 1 Chapter 3 Descrptve Statstcs: Numercal Measures Part B Measures of Dstrbuton Shape, Relatve Locaton, and Detectng Outlers Eploratory Data Analyss
More informationProduction and Supply Chain Management Logistics. Paolo Detti Department of Information Engeneering and Mathematical Sciences University of Siena
Producton and Supply Chan Management Logstcs Paolo Dett Department of Informaton Engeneerng and Mathematcal Scences Unversty of Sena Convergence and complexty of the algorthm Convergence of the algorthm
More information2) In the medium-run/long-run, a decrease in the budget deficit will produce:
4.02 Quz 2 Solutons Fall 2004 Multple-Choce Questons ) Consder the wage-settng and prce-settng equatons we studed n class. Suppose the markup, µ, equals 0.25, and F(u,z) = -u. What s the natural rate of
More informationA Comparison of Statistical Methods in Interrupted Time Series Analysis to Estimate an Intervention Effect
Transport and Road Safety (TARS) Research Joanna Wang A Comparson of Statstcal Methods n Interrupted Tme Seres Analyss to Estmate an Interventon Effect Research Fellow at Transport & Road Safety (TARS)
More informationreferences Chapters on game theory in Mas-Colell, Whinston and Green
Syllabus. Prelmnares. Role of game theory n economcs. Normal and extensve form of a game. Game-tree. Informaton partton. Perfect recall. Perfect and mperfect nformaton. Strategy.. Statc games of complete
More informationElton, Gruber, Brown and Goetzmann. Modern Portfolio Theory and Investment Analysis, 7th Edition. Solutions to Text Problems: Chapter 4
Elton, Gruber, Brown and Goetzmann Modern ortfolo Theory and Investment Analyss, 7th Edton Solutons to Text roblems: Chapter 4 Chapter 4: roblem 1 A. Expected return s the sum of each outcome tmes ts assocated
More informationSurvey of Math: Chapter 22: Consumer Finance Borrowing Page 1
Survey of Math: Chapter 22: Consumer Fnance Borrowng Page 1 APR and EAR Borrowng s savng looked at from a dfferent perspectve. The dea of smple nterest and compound nterest stll apply. A new term s the
More informationNote on Cubic Spline Valuation Methodology
Note on Cubc Splne Valuaton Methodology Regd. Offce: The Internatonal, 2 nd Floor THE CUBIC SPLINE METHODOLOGY A model for yeld curve takes traded yelds for avalable tenors as nput and generates the curve
More informationFinal Exam. 7. (10 points) Please state whether each of the following statements is true or false. No explanation needed.
Fnal Exam Fall 4 Econ 8-67 Closed Book. Formula Sheet Provded. Calculators OK. Tme Allowed: hours Please wrte your answers on the page below each queston. (5 ponts) Assume that the rsk-free nterest rate
More informationCOMPARISON OF THE ANALYTICAL AND NUMERICAL SOLUTION OF A ONE-DIMENSIONAL NON-STATIONARY COOLING PROBLEM. László Könözsy 1, Mátyás Benke 2
COMPARISON OF THE ANALYTICAL AND NUMERICAL SOLUTION OF A ONE-DIMENSIONAL NON-STATIONARY COOLING PROBLEM László Könözsy 1, Mátyás Benke Ph.D. Student 1, Unversty Student Unversty of Mskolc, Department of
More informationEfficient calculation of expected shortfall contributions in large credit portfolios
Effcent calculaton of expected shortfall contrbutons n large credt portfolos Mchael Kalkbrener 1, Anna Kennedy 1, Monka Popp 2 November 26, 2007 Abstract In the framework of a standard structural credt
More informationData Mining Linear and Logistic Regression
07/02/207 Data Mnng Lnear and Logstc Regresson Mchael L of 26 Regresson In statstcal modellng, regresson analyss s a statstcal process for estmatng the relatonshps among varables. Regresson models are
More informationR Functions to Symbolically Compute the Central Moments of the Multivariate Normal Distribution
R Functons to Symbolcally Compute the Central Moments of the Multvarate Normal Dstrbuton Kem Phllps Sourland Bostatstcs Abstract The central moments of the multvarate normal dstrbuton are functons of ts
More informationDiscounted Cash Flow (DCF) Analysis: What s Wrong With It And How To Fix It
Dscounted Cash Flow (DCF Analyss: What s Wrong Wth It And How To Fx It Arturo Cfuentes (* CREM Facultad de Economa y Negocos Unversdad de Chle June 2014 (* Jont effort wth Francsco Hawas; Depto. de Ingenera
More informationOPERATIONS RESEARCH. Game Theory
OPERATIONS RESEARCH Chapter 2 Game Theory Prof. Bbhas C. Gr Department of Mathematcs Jadavpur Unversty Kolkata, Inda Emal: bcgr.umath@gmal.com 1.0 Introducton Game theory was developed for decson makng
More informationElements of Economic Analysis II Lecture VI: Industry Supply
Elements of Economc Analyss II Lecture VI: Industry Supply Ka Hao Yang 10/12/2017 In the prevous lecture, we analyzed the frm s supply decson usng a set of smple graphcal analyses. In fact, the dscusson
More informationQuadratic Games. First version: February 24, 2017 This version: December 12, Abstract
Quadratc Games Ncolas S. Lambert Gorgo Martn Mchael Ostrovsky Frst verson: February 24, 2017 Ths verson: December 12, 2017 Abstract We study general quadratc games wth mult-dmensonal actons, stochastc
More informationCorrected Maximum Likelihood Estimators in Linear Heteroskedastic Regression Models *
Corrected Maxmum Lkelhood Estmators n Lnear Heteroskedastc Regresson Models * Gauss M. Cordero ** Abstract The lnear heteroskedastc regresson model, for whch the varance of the response s gven by a sutable
More informationCreating a zero coupon curve by bootstrapping with cubic splines.
MMA 708 Analytcal Fnance II Creatng a zero coupon curve by bootstrappng wth cubc splnes. erg Gryshkevych Professor: Jan R. M. Röman 0.2.200 Dvson of Appled Mathematcs chool of Educaton, Culture and Communcaton
More informationCollective Motion from Consensus with Cartesian Coordinate Coupling - Part II: Double-integrator Dynamics
Proceedngs of the 47th IEEE Conference on Decson Control Cancun Mexco Dec. 9-8 TuB. Collectve Moton from Consensus wth Cartesan Coordnate Couplng - Part II: Double-ntegrator Dynamcs We Ren Abstract Ths
More informationAn Approximate E-Bayesian Estimation of Step-stress Accelerated Life Testing with Exponential Distribution
Send Orders for Reprnts to reprnts@benthamscenceae The Open Cybernetcs & Systemcs Journal, 25, 9, 729-733 729 Open Access An Approxmate E-Bayesan Estmaton of Step-stress Accelerated Lfe Testng wth Exponental
More informationNumerical Analysis ECIV 3306 Chapter 6
The Islamc Unversty o Gaza Faculty o Engneerng Cvl Engneerng Department Numercal Analyss ECIV 3306 Chapter 6 Open Methods & System o Non-lnear Eqs Assocate Pro. Mazen Abualtaye Cvl Engneerng Department,
More informationTo Rebalance or Not to Rebalance? Edward Qian, PhD, CFA PanAgora Asset Management
To Rebalance or Not to Rebalance? Edward Qan, PhD, CFA PanAgora Asset anagement To Rebalance or Not to Rebalance It s not THE QUESTION but a very mportant one»to rebalance fxed-weght (FW); Not to Buy and
More informationarxiv: v2 [q-fin.pr] 12 Oct 2013
Lower Bound Approxmaton to Basket Opton Values for Local Volatlty Jump-Dffuson Models Guopng Xu and Harry Zheng arxv:1212.3147v2 [q-fn.pr 12 Oct 213 Abstract. In ths paper we derve an easly computed approxmaton
More informationPricing American Options using Monte Carlo Method
Prcng Amercan Optons usng Monte Carlo Method Zhemn Wu St Catherne s College Unversty of Oxford A thess submtted for the degree of Master of Scence n Mathematcal and Computatonal Fnance June 21, 2012 Acknowledgements
More informationQuadratic Games. First version: February 24, 2017 This version: August 3, Abstract
Quadratc Games Ncolas S. Lambert Gorgo Martn Mchael Ostrovsky Frst verson: February 24, 2017 Ths verson: August 3, 2018 Abstract We study general quadratc games wth multdmensonal actons, stochastc payoff
More informationME 310 Numerical Methods. Differentiation
M 0 Numercal Metods fferentaton Tese presentatons are prepared by r. Cuneyt Sert Mecancal ngneerng epartment Mddle ast Tecncal Unversty Ankara, Turkey csert@metu.edu.tr Tey can not be used wtout te permsson
More informationProject Management Project Phases the S curve
Project lfe cycle and resource usage Phases Project Management Project Phases the S curve Eng. Gorgo Locatell RATE OF RESOURCE ES Conceptual Defnton Realzaton Release TIME Cumulated resource usage and
More information