Article from: Risk Management August 2008 – Issue No. 13 measure. The purpose of this articleTail is to Expectation provide a high(CTE, level summary of theExpected paper “Variance of theorCTE measure. However, the Conditional also called Shortfall Estimator” B. article John Manistre and Geoffrey H. desirable Hancock that appeared in theease North American Tail-VaR) is The becoming its and purpose increasingly ofbythis isprevalent to provide adue hightolevel summary ofproperties the paper “Variance ofofthe CTE August Management The Value-at-Risk (VaR) estimator (i.e., percentile or quantile value) is still often2008 usedw Risk as a risk Actuarialby Journal inManistre 2005. We also expandH. onHancock some of the results. Estimator” John and Geoffrey in the CTE North is American interpretation. A tool thatB.can quantify the statistical precisionthat of appeared an estimated therefore measure. However, theresults. Conditional Tail Expectation (CTE, also called Expected Shortfall or Actuarial Journal in 2005. We also expand on some of the important. That the sampling error in for theestimating estimatereserves can be and placed in using perspective with other Anis, actuary who is responsible capital stochastic methods Tail-VaR) is becoming increasingly prevalent due tomust its desirable properties and ease of modeling issues, including model and assumption risk. withwho a wide range of issues. Not only must theand actuary worry about the underlying Andeal actuary isparameter, responsible forinterpretation. estimating reserves capital using stochastic methods must precision of an estimated CTE is therefore A tool that can quantify the statistical stochastic model, it’sofparameterization, assumptions and calculation formulae but he/she deal with a wide range issues. important. Not onlydata, must the actuary abouterror the underlying That is,quantify the worry sampling the estimate can be placed in perspective with other now hasmodel, (orin should have) the new issue of assumptions trying to thewhat precision ofinthe estimated If all we werestochastic interested was a regular mean then we would know to do. We draw it’s parameterization, data, and calculation formulae but he/shearisk modeling issues, including parameter, model andrisk assumption risk. measure. now has (or should have) the new issue of trying to quantify the precision of the estimated 1 . sample x1 , xmeasure. ¦i asxai risk 2 ,..., x n of size n from our model and calculate the sample average P̂ n then The Value-at-Risk (VaR) estimator percentile or quantile value) still often used If all(i.e., we were interested in was a is regular mean we would know what to do. We draw a John Manistre and GeoffreyThe Hancock measure. theestimator Conditional Expectation (CTE,value) also called Shortfall Value-at-Risk (VaR) (i.e.,Tail percentile or quantile is stillExpected often used as a riskor Statistical theory tells usHowever, three things 1 Tail-VaR) is becoming increasingly prevalent ease and ofor calculate the sample average P̂ x1 , x2due ofdesirable ,..., to x(CTE, sample sizecalled n properties from ourand model measure. However, the Conditional Tail Expectation also Expected Shortfall ¦ xi . n its n i interpretation. A toolincreasingly that can quantify the statistical precision properties of an estimated CTE is therefore Tail-VaR) is becoming prevalent due to its desirable and ease of ˆ E [ P ] P the true mean. We expect to 1. The sample average is an unbiased estimator i.e. important. That is, the errorthe instatistical the estimate canusbethree placed in perspective other Statistical theory tells things interpretation. A tool thatsampling can quantify precision of an estimated CTE iswith therefore get the right answer. modeling issues, including parameter, model and assumption risk.in perspective with other important. That is, the sampling error in the estimate can be placed he purpose of this articleissues, is toincluding provide 2.model The variance of the average is modeling parameter, assumption risk.sample 1. Theand sample average is an unbiased estimator i.e. E[ Pˆ ] P the true mean. We expect to If all we were interested in was a regular mean then2we would know2 what to do. We draw a ˆ a high level summary paper is the truevariance. variance. VAR [ P ] V / n V 2. The variance ofofthethe sample average is where is the true get the answer. 1 a If all we were interested in was a regular mean thenright we would know what to do. We draw Summary of “Variance of the CTE Estimator” T sample x1 , x2 ,..., xn of size n from our model and calculate the sample average P̂ “Variance of thesample CTE byn from our model and calculate the sample average P̂ 1 n x¦i . xi . xEstimator” , x2 ,..., xlarge ¦i i [Pˆthe 1size n of size 3. If the sample is enough, and a few variance other technical conditions are satisfied, n VAR StatisticalH. theory tells us three things ] V 2 / n where V 2 is the true variance. theis sample is B. John Manistre and Geoffrey Hancock 3.2.IfThe the sample of size largeaverage enough, and a few P̂ has an approximate normal distribution. estimator Statistical theory tells us three things E[conditions Pˆ ] P the trueare mean. We expectthe to 1. The sample average is an unbiased estimator that appeared in the North American Actuarial other technical satisfied, 1 i.e.size 2 enough, and a few other technical conditions are satisfied, the 3.from Ifestimator the is large Vˆ 2 sample Pˆ )true and the actuary would The distribution’s can be estimated [ Pˆ ]( xiPthe mean. We expect to 1. variance The average is an unbiased i.e. E¦ getsample the right answer. Journal in 2005. We also expand on some of estimator has anapproximate approximate normal disnormal distribution. estimatorn P̂ 1has i an get the right answer. 2 2 Vˆ average 1 ˆ VAR [ P ] V / n V 2. The variance of the sample is where is the true variance. 2 the results. then report the result tribution. giving any potential user sense of howfrom largeVˆthe of the work as Pˆ rThe distribution’s ( x Pˆ ) 2 and the actuary would variance canabe 2 2 estimated ¦ i 2. The variance of the sample average is VAR[ Pˆ ] V / n where V is the true variance. n n 1 i 3. If the sample size is large enough, and a few other technical conditions are satisfied, the sampling error (i.e., statistical uncertainty) might be. Users can then judge whether large ˆ this is the Vsatisfied, 3. If the enough, and a fewdistribution. other of technical conditions arecan P̂size hasisanlarge approximate normal estimator any potential user a sense of how large the then report result work as Pˆ r An actuary who is responsible forsample estimating The distribution’s variance be estior small relative to estimator other model asthe a change inthe lapse assumptions forgiving example, P̂ hassensitivities, an approximatesuch normal distribution. n 1 Vˆ 2 1 the ¦ ( xi Pˆ ) 2 of and reactusing appropriately. In methods particular, users can judge whether precision the is and theestimate actuary would The distribution’s variance can bemated estimated from reserves and capital stochastic from and the actuary sampling statistical mightwould be. Users can then judge whether this is large n 1 ( xi uncertainty) ˆ2 Pˆ ) 2 and the actuary can be estimated error from V(i.e., ¦imodel i high enough The for distribution’s the intendedvariance application. n 1 or small relative to other sensitivities, such as change in lapse assumptions for example, must deal with a wide range of issues. Not only would Vˆ then report the result of the worka as rˆ giving any potential user a sense of how large the then report the result of the work as Pˆ V and react appropriately. In particular, users can judge whether the precision of the estimate is n r then report result of the work as Pˆ estimating giving anyany potential user sense large of the how must the actuary about thethe underlying giving potential userof ahow sense How worry does this process change when we start Conditional Taila Expectations? n might enough for the application. sampling error (i.e., statistical high uncertainty) be. intended Users can then judge whether this is large stochastic model, it’s parameterization, data, largemight thesuch sampling (i.e., statistical uncersampling (i.e.,tostatistical uncertainty) be. as Users canerror then judge whether thisfor is large or smallerror relative other model sensitivities, a change in lapse assumptions example, First, the badornews. There is other no general set of formulas that are guaranteed to work inestimate all small relative to model sensitivities, such a change in lapse assumptions for example, and react appropriately. In particular, users canas judge whether thewhen precision of the is Conditional Tail Expectations? How does this process change we start estimating assumptions and calculation formulae butofhe/ tainty) might can then judge whether circumstances. The distribution CTE estimator depends onUsers a wide range of the variables such as and react appropriately. Inaparticular, users can judgebe. whether the precision of estimate is high enough for the intended application. high enough for the intended application. sample size,have) the actual sampling from and the relative method oftoestimation itself. she now has (or should the distribution new issue you of are this is large or small other model senthestart badestimating news. There is no general set of formulas that are guaranteed to work in all when we Conditional Tail Expectations? Summary of ERM paper How does this process changeFirst, circumstances. TheConditional distribution alapse CTE estimator depends on a wide range of variables such as How does change when we start estimating Tailof trying to quantify the precision of this theprocess estimated sitivities, such as a change inExpectations? assumptions John Manistre sample actualthat distribution youtoare sampling itself. First, the bad news. There is no generalsize, set ofthe formulas are guaranteed work in all from 1 and the method of estimation John Manistre, FSA, FCIA, is risk measure. for set example, react appropriately. parFirst, the bad news.The There is no general ofestimator formulasand that are guaranteed to work in all Insuch circumstances. distribution of a CTE depends on a wide range of variables as circumstances. The distribution of a CTE estimator depends on a wide range of variables such as vice president, group risk, sample size, the actual distributionticular, you are sampling from judge and the method of estimation itself. users can whether the precision The purpose of this article is to high level summary of and thethepaper of the CTE sample size,provide the actual a distribution you are sampling from method“Variance of estimation itself. 1 at AEGON NV in Baltimore, The Value-at-Risk (VaR) estimator per- H. of the estimate is high enough the intended Estimator” by B. John Manistre and(i.e., Geoffrey Hancock that appeared in thefor North American a high level summary of the paper “Variance of the CTE 1 Actuarial Journal invalue) 2005.isWe also expand on some of the results. 1 centile or quantile still often used application. eoffrey H. Hancock that appeared in the North American as a risk measure. pand on some of the results. However, the Conditional An actuary who is responsible for estimating reserves and capital using stochastic methods must Tail Expectation (CTE, also called Expected How does this process change when we start deal with a wide range of issues. Not only must the actuary worry about the underlying mating reserves andorcapital using stochasticincreasingly methods must estimating Conditional Tail Expectations? Shortfall Tail-VaR) is becoming stochastic model, it’s parameterization, data, assumptions and calculation formulae but he/she nly must prevalent thesummary actuary worry the underlying due topaper itsabout desirable properties e a high level ofshould the “Variance of theissue CTE ofand now has (or have) the new trying to quantify the precision of the estimated risk data, assumptions and calculation formulae but he/she Geoffrey H. ease Hancock that appeared in the North American measure. of interpretation. A tool that can quantify First, the bad news. There is no general set of trying to quantify the precision of the estimated risk xpand on some of the results. the statistical precision of an estimated CTE is of formulas that are guaranteed to work in all The Value-at-Risk (VaR) estimator (i.e., percentile or quantile value) is still oftenofused as a esrisk mating reserves and capital using stochastic methods must therefore important. That is, the sampling error circumstances. The distribution a CTE measure. However, theunderlying Conditional Tail Expectation (CTE, also called Expected Shortfall or only must the actuary worry about the in the estimate can be in used perspective ., data, percentile or quantile value) is placed still often as a risk timator depends on a wide range of variables n, assumptions and calculation formulae but he/she Tail-VaR) is becoming increasingly prevalent due to its desirable properties and ease of Expectation (CTE, also called Expected Shortfall or other modeling issues, including paramsuchprecision as sample the actual distribution eilof trying with to quantify the precision of the estimated risk interpretation. A tool that can quantify the statistical of size, an estimated CTE is therefore evalent due to its desirable properties and ease of eter, model and assumption risk. error in the estimate you are sampling from and the with method important. That is, the sampling can be placed in perspective otherof the statistical precision of an estimated CTE is therefore modeling issues, including parameter, model and assumption risk. estimation itself. e., percentile or quantile is still often used as awith risk other n the estimate can bevalue) placed in perspective ail Expectation also called Expected Shortfall or mean Ifassumption all(CTE, we were interested in was a regular model and risk. revalent dueIftoall its we desirable andinease ofa regular mean then we would know what to do. We draw a wereproperties interested thenprecision we would whatCTE towas do.therefore We draw a Now, the (really) good news. y the statistical of anknow estimated is 1 arinmean then we would know what to do. We draw a the estimate can be in perspective with other sample of size n from our model xplaced , x ,..., x sample of size n from our model and calculate the sample average P̂ ¦i xi . 1 2 n n model andand assumption risk. 1 calculatethe the sample average average P̂ f the sample size being used is large enough, r model and calculate ¦ xi . 1. Ithen Statistical theorysample tells us three things there are approximate formulas analolar mean then we would know what to do. We draw a n i 1 Statistical theory tells us three things gous to apply an ordinary mean. to P̂ is an xi . ur model and calculate thesample sample average [ Pˆ ] that P the 1. The average estimator i.e. Ethose trueformean. We expect ¦unbiased n i It is therefore possible to quantify the statistianswer. ˆ ] right E[ Pthe P the sed estimator i.e. get true mean. We expect to 1. The sample average is an unbiased estimator cal precision of an estimated CTE. ˆ E [ P ] P the true mean. We expect to ased estimatori.e., i.e. the true mean. We expect to get 2. The variance of the sample average is VAR[ Pˆ ] V 2 / n where V 2 is the true variance. the right answer. 2. There is a practical process, called a “variage is VAR[ Pˆ ] V 2 / n where V 2 is the true variance. 3. If size is large enough, and a few otherverification” technical conditions satisfied, the ance exercise are in this article, 2 the sample 2 erage is VAR[ Pˆ ] V / n where V is the true variance. estimator P̂ has an approximate normal distribution. gh, and a few other technical conditions are satisfied, the gh, and a few other technical conditions are satisfied, the 1 normal distribution. on page 38 the actuary would The distribution’s variance can be estimated from Vˆ 2 ¦i ( xi Pˆ ) 2 andcontinued e normal distribution. n 1 1 2 2 Pˆ )theand mated fromVˆV2ˆ 1 ¦ (¦ the actuary xi ( xPˆi) 2and actuary would would mated from Vˆ n n1i1 then report thei result of the work as Pˆ r giving any potential user a sense of how large the n V̂ V̂ any potential user user a sense of howof large thelarge the rr n giving giving any potential a sense how sampling error (i.e., statistical uncertainty) might be. Users can then judge whether this is large n inty) might or be.small Users relative can then judge whether this sensitivities, is large to other model such as a change in lapse assumptions for example, nty) might be.a change Usersincan then judge whether this is large vities, suchand as lapse assumptions for example, react appropriately. In particular, users can judge whether the precision of the estimate is Md. He can be reached at [email protected]. Geoffrey Hancock, FSA, CERA, FCIA, is a director at Oliver Wyman in Toronto, ON. He can be reached at geoffrey. [email protected]. Page 37 w Risk Management The standard approach to this problem is to start with a ran w August 2008 Chairman’s Variance ofCorner the CTE Estimator from x CTE Estimator … w continued from page 37 (1) the model and then sort the sample in descending t x ( 2 ) t ... t x ( n ) . We can then calculate the plug-in es by looking at the observed empirical distribution. Let D the plug-in estimators are 1 k CTˆE n D ¦ x(i ) , In terms of this notation statistical theory hask i 1 the following three things to say VaˆRn D x(k ) . Now, the (really) good news. 1 that one can execute to test (and confirm) the validity of the approximate formulas for any 1.Now, If the sample size being used is large enough, there are approximate formulas application. Anthen example is given the (really) good news. particular In terms of this notation statistical theory has the following th to size quantify analogous to those that apply an ordinary at the for end this article. mean. It is therefore 1. Ipossible f the sample is large enough, and a few Now, the (really) good news. the statistical precision of an estimated CTE. 1. If the sample size being used is large enough, then there are approximate formulas very technical hold, If the sample size is conditions large enough, andthen a few other 1. other Now, the (really) good news.mean. It is therefore possible to quantify analogous to those that ordinary thevariance sample size being used large enough, formulas 1.forIfan 3. Sapply ome standard reduction toolsissuch the pair hasan an approximatemultivaria (CTˆare E n ,approximate VaˆRn ) has thethen pairthere approximate precision an estimated CTE. analogous to those that applyvariate for anexercise ordinary m It isnormal therefore possible to quantify 2. Therethe is statistical a practical process, called a “variance verification” inean. this that asof sampling and control multivariate distribution. If the sample size being used is large enough, then there are article, approximate formulas 1.importance the statistical precision of an estimated CTE. te formulas forpossible any to quantify one can execute to test (and confirm) validity of the analogous to those that apply forproblem anapproxima ordinary mean. It is therefore methods can bethe adapted to the CTE 2. Thereapplication. is a practical An process, called a given “variance verification” exercise in this article, that the statistical precision of an estimated CTE. 2. T particular issometimes at dramatically, the end this article. toexample improve, the heverification” estimator pair is asymptotically unbiis a practical called a “variance exercise in this article, that 2. There te formulas for any one can execute to test (and confirm) the validityprocess, of the approxima one can execute to test (and confirm) the validity of the approxima te formulas for any precision estimator acalled givenacomased. For exercise any finite sample sizethat the CTE Thereofnews. isa aCTE “variance verification” in this article, 2. example Now, the news. (really) good particular application. An ispractical given atprocess, theforend article. Now, the (really) good thenews. (really) good particular application. Anthis example is given at(really) the end this article. Now, theunbiased. good news. 3. Some standardNow, te for variance reduction tools such as(and importance sampling and control varia putational cost. These techniques will be the plug-in estimator is negatively biased one can execute to test confirm) validity of the approxima te formulas any 2. The estimator pair is asymptotically For any finite sample size the CTE i.e. plug-in Now, the (really) good news. ˆ particular application. An example is given at the end this article. can be adapted to the CTE problem to improve, sometimes dramatically, the the sample size being used is large enough, then there are approximate formulas 1. Ifmethods 2. The estimator pair is asymptotically unbiased. For any finite sample size the CTE plug-in n o f. estimator is negatively biased i.e. , but the bias goes to 0 as E [ C T E ] CTE topic of3.asize second article variance. but the and bias goes sample size being used is large enough, then there aresampling approximate formulas 1. If the If the sample being usedonas isCTE large enough, then there are approximate formulas 1.variance 3. Some te standard reduction tools such importance control Some standard variance reduction tools1. such as sampling varia teto Ifand theimportance sample varia size being usedcontrol is large enough, then t analogous to those that apply for an ordinary m ean. It is therefore possible to quantify ˆ analogous to estimator those that an ordinary m ean. It issuggests therefore possible too quantify analogous toi.e. those that apply for an is ordinary m ean. Itdramatically, isthere therefore tothan quantify ational cost. These will be precision of acan for afor given comput If the sample size large enough, approximate formulas 1.toapply nare fpossible Practical experience thethen bias istechniques usually much smaller the sampling error. estimator isCTE negatively biased ,used but the bias goes to 02. as E [C TE ]can being CTE .sometimes Practical experience methods beto adapted to the CTE problem to improve, dramatically, thesugmethods be adapted the CTE problem improve, sometimes the The estimator pair is asymptotically unbiased. analogous to those that apply for an ordinary m ean.F Some standard variance reduction tools such as importance sampling and control varia 3. te the statistical the precision of an estimated CTE. the statistical precision of an estimated analogous to those apply for2.CTE. an ordinary mean. Itisis therefore possible to quantify precision ofthe an estimated CTE. The pair asymptotically unbiased. Forsmaller anywill finite sample size the precision of athat CTE estimator for aestimator given comput ational cost. These techniques be the topic of statistical aexperience second article CTE variance. ational cost. These techniques will be precision of a CTE estimator for a is given comput So what are the approximate formulas? First, gests the bias is usually much than Practical suggests bias usually much smaller than the sampling error. ˆE ] C methods can be adapted to the CTE problem to improve, sometimes dramatically, the the statistical precision of an estimated CTE. 2. on The estimator pair is asymptotically unbiased. For any finite sample size the CTE plug-in estimator is negatively biased i.e. E [ C T the statistical precision of anarticle estimated CTE. the topic of afollowing second on CTE variance. is negatively biased i.e. E goes to [Calso TˆE ] asymptotically CTE , but thevalid The approximate formulas are 3.variance. topicprocess, of a second article on CTE precision of CTE estimator for want a estimator given comput ational cost. These techniques will beas n bias need some notation. Suppose we tovariance/covariance the sampling error. a practical awe “variance verification” exercise in thisbiased that 2. There is the ˆE ]in There is a estimator practical process, called aarticle, “variance verification” exercise in this article, that 2.called o f isa“variance negatively i.e. , but the bias goes to 0 E [ C T CTE . samp Practical experience suggests the bias is usually There is a practical process, called a verification” exercise this article, that 2. Practical experience suggests the bias is usually much smaller than the Now, the (really) good news.to test (and the topic a second article on CTE variance. oneare can execute confirm) the2.validity ofConditional the approxima teneed formulas for any There issome aexercise practical process, called awant “variance verif The approximate variance/covariance formulas are also asymptotically valid 3. estimate theconfirm) Tail Expectation of a First, one can execute test (and confirm) the of2. the approxima tewe formulas for There isto a of practical process, called avalidity “variance verification” in this article, that 2any So what thefollowing approximate formulas? First, we some notation. Suppose want to So what are the approximate formulas? we need notation. Suppose we to one can execute to test (and the validity of the approxima te formulas for any D ( | ) ( ) VAR X X VaR CTE VaR t Practical experience suggests the(Cbias is usually much smaller than the sampling error. particular application. An example is given at thecan endexecute this article. ) |end , VAR TˆEthe one can execute to test (and confirm) the validity of particular application. An example is given at this article. one to test (and confirm) the validity of the approxima te formulas for any So what are the approximate formulas? First, we need some notation. Suppose we want to n random variable X, with cumulative distribu3. T he following approximate variance/covariIf the sample size particular being usedapplication. is largeExpectation enough, there approximate formulas 1. itional Expectation ofwea random variable X, cumulative distrto estimate the Cond ibution An then example is at Tail the end article. The approximate also asymptoti 3.this D ) with n(1variance/covariance distr estimate the Cond ibution itional Tail ofapplication. aaregiven random variable X,following cumulative The following approximate variance/covariance 3. some So what are the approximate formulas? First, need notation. Suppose weformulas want 2with particular An example isVaR given at the end this article. particular application. An example isare given at the en analogous to those that apply for an ordinary m ean. It is therefore possible to quantify D ( | ) ( ) VAR X X VaR CTE t estimate the Cond ibution itional Tail Expectation of a random variable X, with cumulative distr Thus, formulas arecumulative also asymptotically valid Pr^|as X dCond x` itional F ( xsampling ) atTail the level DThus, we want toance the conditional expectation ˆstandard reduction tools such and control varia tewe (variable 1calculate )sampling D expectation Dformulas (tion , CThus, TCTE. E Expectation of aRrandom distribution estimate theimportance Some variance reduction toolsvariance/covariance such as importance andalso control varia 3. VAR te Pr^the X3.statistical dSome x` standard Fprecision )variance at of the D we want to calculate the following approximate are asymptotically valid 2 3. n ) The anlevel estimated (conditional ,X, with VAR Vaˆsampling n) | 2 D ( 1 ) n `( xbeSome variance reduction tools such as importance and control varia te Pmethods r^X d 3. xcan Fadapted ( x)standard at D Thus, we want to calculate the conditional expectation D ( | ) ( VAR X X VaR CTE t VaR ) tothe the level CTE problem to improve, sometimes dramatically, the want to calculate the conditional expectation (VaR )]standard n[importance f XSome ˆEn ) |dramatically, methods adapted tolevel the CTE problem to improve, sometimes the Some variance tools such and control variate tools Pr^X3. dcan x` beFstandard ( x) at the Dreduction Thus, we want to as calculate the expectation (Cconditional VAR Tsampling variance reduction such 3. ( X |as VAR X, imp tV ˆEDn ) | can beprecision to the totoin sometimes dramatically, the of amethods CTE estimator for aaadapted given comput cost. These techniques will be (that ( 1 n ( VAR C T (1can ) adapted D CTE estimator Dproblem CTE D )ational > X | X ! q @ is a practical process, called “variance verification” exercise thisCTE article, 2. Thereprecision of aational CTE for aimprove, given comput cost. These techniques will be methods be the problem toE improve, sometimes dramatically, the D )can be adapted to 2 the CTE problem to impr ( CTE VaR D methods ˆ cost. (the , V aˆRvalidity ()X| | Xational ) D (CTE ) VAR t VaR n ˆRnThese (CT . These CoV n) | theexecute topic of atoprecision second article CTE variance. 2 comput ofon aVAR CTE estimator aofgiven ational techniques will betechniques one can test (and confirm) of the approxima te!formulas for any n , Va ˆvariance. D(1VaR D ) ,will be the topic of an[second precision for a Egiven cost. CTE (for )(VaR E)]CTE XX estimator ||on X VAR E fDX(D )(Vaˆof VaR n )(| CTE ) aarticle E > X CTE !(CqqT @ CTE D ) Ecomput > X | X nf!precision @ , for a given computatio VAR Rn )a|CTE estimator Xq(D DD particular application. An example is given at the end this article. 2 the topic of a second article on CTE variance. D ( 1 ) n D (1 D ) the topic second article on CTE variance. )] on CTE variance. n[ f X (VaR where isof theaVaR D -quantile, defined as to thethe smallest value satisfying the topic of af Xsecond article (x ) of probability density the random X at the 2 , DThe ) fSuppose D (qCTE (VaˆRvariable VAR X (VaR ) refers n) | So what are the approximateSo formulas? First, we need somenotation notation. we want toneed ˆEare 2. The estimator pair is asymptotically unbiase ˆ ( , ) . CoV C T V a R | what the approximate formulas? First, we some notation. Suppose we want to ( ) CTE VaR D ( 1 ) D D [ )] n f X (VaR n n such as importance sampling and control varia 3. Some standard variance reduction tools te ˆ ˆ ( , ) . CoV C T E V a R | point x = VaR. (-quantile, )withvalue nf VaR ˆR where defined the smallqas isthe the DTail defined the value satisfying (Vas , X, with VAR aas |smallest n notation. n itional Tail Expectation of Cond awhat random variable X, cumulative distr estimate Cond ibution X D n )random So are the approximate formulas? First, we need Suppose we wantbiased to 2 some So what are the approximate formulas? First, we need some notation. Suppose we want to where q is the D -quantile, defined smallest satisfying itional Expectation of a variable cumulative distr estimate the ibution where qthe is the D -quantile, defined as the smallest value satisfying estimator is negatively i.e. E [ C TˆE)] methods can be adapted to the CTE problem to improve, sometimes dramatically, the ( ) nf VaR D D X [ f^So nPr ` are X(VaR !what qD )] 1 Dthe approximate VaR D (CTE formulas? First, weneed X Pprecision r^X dThe x` estimate Fnotation (ax)CTE at thethe level D Thus, we want to calculate the conditional expectation ˆ est value satisfying f X Cond (VaR ) f (x ) refers to the probability density of the random variable X at the itional Tail Expectation of a random variable X, with cumulative distr estimate the Cond ibution of estimator for a given comput ational cost. These techniques will be CoV (Cdensity TEn ,VaˆfRX n(x) )|of the random vari Pr^X d x` Tail F ( xExpectation ) at the levelofDaThus, want toX, calculate the conditional expectation X we itional random variable with cumulative distr ibution f X (VaR ) refers The notation to the probability Several comments are in order Practical suggestsofthe bias is usua )theCond D!(CTE 1VaR itionalexperience Tail Expectation anfrandom the topic of a second article on CTE variance. ^|xwant Pr Xconditional qestimate expectation Dthe X (VaR )var Pr^X D d Thus, x` F ( xwe ) atwant the level D we toD `calculate conditional expectation ˆEThus, pointPr^xX= dVaR. ˆ ( , ) . CoV C T V a R x` F ( x) at the level to calculate the point = VaR. nD n ^ ` Pr X ! q 1 is Doften1called Value-at-Risk isFused extensively in ^(qVaR `notation PXr(VaR) XThe d )xand ( x) at the level Dthe Thus, we want to calculat D CTE (D ) E>The X | XDPr !-quantile qDX @ ! (VaR ) refers notation tofinancial probability The refers the probability D (D )easier Ewe >toX work |want Xnf ! x qnotation. CTE isCTE clearly with VaR. If wef Xwere using VaRto as athe risk measure thendensit So what are the approximate formulas? First, we need some Suppose toD @than management of trading risk over a fixed (usually short)3.time horizon. The following approximate variance/covaria Dnotation -quantile isf Xoften Value-at-Risk inrandom the financial (Dprobability ) comments E >Xto| and Xestimate ! isqinDused @order Several are f X (VaR ) in order wecalled have to find a (VaR) way probability density to )would fVaR. toCTE density the random variable X the at the Several comments are in order density of the variable X at point xthe =extensively X (x ) of itional Tail Expectation ofThe a random variable with distr estimate the Cond ibution CTE (D(VaR ) X, Erisk >refers X over |cumulative X ! qthe D @ (usually short) time horizon. management of trading aused fixed CTE ( D ) E >nXtails | X ! qD apply the asymptotic formula. This can be very difficult in practice, especially in the Dat -quantile often called Value-at-Risk (VaR) and is extensively in the financial qFD (isx)the Dthe -quantile, defined as the smallest value satisfying Pr^X dwhere x` The level is D Thus, we want to calculate the conditional expectation where qD ispoint the D -quantile, defined as to the smallest value satisfying The isapproach often and called Value-at-Risk point VaR. x =standard VaR. (VaR) x1 , xIf2 ,..., xn of sizeVaR The this problem is istoclearly start with a=torandom sample The Dm -quantile isisoften called Value-at-Risk is used in xthe financial x extensively CTE easier work with than VaR. we were using ofIf the distribution. ( Xa |risk VARas X x CTE easier workaqwith than VaR. we were using VaR as a risk measure then anagement ofclearly trading riskto over fixed (usually short) time horizon. where is the D -quantile, defined as the smallest value satisfying ˆ Several comments are in order D (Csize VAR TEstatistics | f X (VaR) in (VaR) and is(usually used extensively in thethe financial we wouldin have torandom find a way to estimate the probability density thenproblem sort sample descending order toxobtain the order from the model and n) , x ,..., x of n The standard approach to this is to start with a sample of trading risk over a fixed short) time horizon. management where q is the D -quantile, defined as the smallest value satisfying 1 2 n ) inq order we would haveCTE to find a !way the probability density f X (VaR toD-quantile, defined as the smallest value sati D (Pr D ^)XSeveral Eq>DX` to | Xestimate t!Dxqof D@ Pra^can X !then q(usually 1 where the Dareasymptotic in.order apply formula. This can be very difficult in practice, the We calculate the plug-in estimators for the required param eters especia x (11comments t x... are tThe xrisk management over fixed Several comments are in order and CTE estimators positively correlated. This makes intuitive sense. D is the D` ) model 2 ) trading ( n )VaR then sort in descending order to obtain the statistics from apply the asymptotic formula. This (is can bestart very difficult insample practice, especially in the x , x ,..., xntails The standard approach to this problem toand with athe random sample of sizeeasier n order x CTE is clearly to work with thanDVaR. ^ ` Pr X q 1 D (1 DIf) of!the distribution. 1 2 D short) time horizon. k for the required (Vaˆeters VAR Rn ) | for . observed We can then calculate the plug-in estimators param x (1) tby x (looking ...to t at xstart of the distribution. 2) t ( nthe ) ! where q is the D -quantile, defined as the smallest value satisfying ^ ` Pr X q 1 D empirical distribution. Let so that practical expressions 1 D x , x ,..., x The standard approach to this problem is with a random sample of size n D x The variance of the CTE estimator has two terms. The first term is the “obvious” we would to find arisk way estimate The D-quantile is often calledand Value-at-Risk (VaR) is used extensively D in the xis CTE is clearly easier to financial work withand VaR. If1 we using as awith measure 2 were n have [then )] VaR from the model then sort theand sample inValue-at-Risk descending order tothan obtain the The D -quantile often called (VaR) is used extensively in the VaR financial ^toThis `f Xintuitive Pr XVaR. !makes qnDIf 1(the Dpro n0statistics •CTE C TE isorder to work than x The VaR andin estimators are positively correlated. kclearly extension of(usually what was happening the first ( D the )easier case. The origin ofin the second term management of trading risk over a fixed (usually short) time horizon. apply asymptotic formula. This can be very f (VaR ) we would have to find a way to estimate the probability density order to m anagement of trading risk over a fixed short) time horizon. The D -quantile is often called Value-at-Risk (VaR) and is used extensively in the financial by looking at the observed empirical distribution. Let so that practical expressions for 1 D . Pr We the plug-in estimators theobservation required param eters x (1) xmodel t xThe ... tand xthen (CTE Va D theare estimators from the and sort sample in can descending order then using order statistics ( 2 ) tVaR ( n )CTE estimators positively correlated. This makes intuitive sense. The toare thisby problem is to tofor we were as a Xrisk measure then we be seen conditioning onobtain the of theVaR estimated write aˆT R ^X can `standard ! the qDthen 1plug-in calculate D approach ˆE. We (Vterms. )then CoV C Vaˆcan Rfirst |term m anagement trading risk over aisfixed short) time of theestimator distribution. k horizon. n ,The nin The D-quantile is often called Value-at-Risk (VaR) andformula. usedx(usually extensively the financial applyofthe asymptotic This can bein very difficult in practice, especially the tails The variance of the CTE has two the(VaR “obv) 1 nfand .The We can then calculate the estimators the required param xThe t x ( 2 ) tapproach ... t xto( nthis start with a arandom sample size nnafor would have toxfind a wayeters to estimate prob-is x1, xis2time kof extension is often called ,..., xnstart size nThe problem is to start with random sample the plug-in estimators areplug-in X is us D-quantile CTˆthat E Drandom x(expressions (1) standard ) of xValue-at-Risk size nthe(VaR) standard approach to this problem to with sample i ) ,happening of what was in2 ,..., the first The origin of the D 0 ) case. 1, x n (of management trading risk over a fixed (usually short) of the distribution. by looking at the observed empirical distribution. Let so practical for 1Cfinancial Tˆhorizon. D ˆEk ¦ ˆ ˆ ˆ k The variance of the CTE estimator has two terms. The first term is the “obvious” VAR [ E ] E { VAR [ C T | V a R ]} VAR { E [ C T E | V a R ]}. The D-quantile isxoften called Value-at-Risk (VaR) and is used extensively in the n to start n i 1 notation n) refers n size mwith anagement ofn sample trading over fixed short) time 1 by from the and ability apply the xrisk ,CTE x2order ,..., xinnaorder fthe The toto(usually the den of n probability The approach tothe thissort problem isnstatistics aorder random the sample in standard descending order tothen obtain the the ordersample from the model and then sort seen the observation of the estimated c VaˆR . We ˆdescending The VaRonand estimators are positively corr 1(VaR X then sort todensity obtain statistics model andmodel Daˆbe in Tk Ecan xx(xconditioning management of trading risk overoffrom a what fixedthe (usually short) time horizon. extension was happening in the first ( D sample The of¦ the term 0 ) case.inC norigin isecond ) ,. V R D byxlooking at the observed empirical distribution. Let so that practical expressions for n (k ) 1 D k . We can then calculate the plug-in estimators for the required param eters t x t ... t x descending order to obtain the order statistics asymptotic formula. This can be very difficult i 1 x point (1) the ( 2 ) plug-in ( n ) estimators then sort the sample in descending order to obtain the order statistics from the model and are . We can then calculate the plug-in estimators for the required param eters x t x t ... t x x = VaR. , x ,..., x of size n Thecan standard approach to this problem is to start with a random sample x The VaR and CTE estimators are positively correlated. This makes intuitive sense. ( 2) ) Intuitively, can say that when wecan the CTE to we this are estimating both the be seen by(1)conditioning on( nthe observation of the we estimated . We write VaˆRThe 1 estimate n ˆ2then ˆECTE standard approach is with a ra ˆtwo [xC Tpractice, EThe E {VAR [CTˆrequired E Vsay aˆproblem RCTE tails VAR {Eto [the CTstart Rn ]}.term k statistical of has three things to|the n] n ]} n | Va size Vxtheory aˆn R Dthe xthe start We can calculate plug-in estimators for the param eters x (1) tIn x (terms t xthis variance ofnthe estimator has xthen x 2calculate ,..., nVAR The standard approach to this problem is to with random sample k( n ) .notation We can then in especially in the dis(k ) . following and VaR uncertainty in the VaR estimate increases uncertainty of the CTE 2 ) ta ... 1 , and n n ofthe 1the kthe then sort the sample in descending order to obtain order statistics from the model and by looking at the observed empirical distribution. Let so that practical expressions for 1 D ˆ by looking at the observed empirical distribution. Let so that practical expressions for 1 D C T E D x , the plug-in estimators are extension of what happening the first ( D estimate. This isparameters the origin offollowing the second then sort the sample inin descending from the model and nn for ( ithe ) theory inEterms descending order order from the model and then sortVAR kterm. The variance of the CTE estimator has terms. The first term iswas “obvious” plug-in estimators the required tribution. Several comments are inthe order In of notation statistical three things toeters say nrequired Cn )TˆE. n We ] by {xVAR [Catthis Tˆcalculate E VaˆRnto ]} VAR TˆE nhas |Intuitively, Vstatistics aˆthe R ]}.two kobtain x t x(are can then the plug-in estimators the param t ...the t x[sample can sayso that when we estimate the CTEfor we are estimating bot n |observed nfor ik 1 {E [C looking the empirical distribution. Let that practical expressions 1 tother Dwe 2 ) can ple size is l ns hold, then If the sam arge enough, and a few very technical conditio 1. x(1) t the then (calculate the plug-in estimators for the required param eters x ( 2 )plug-in t ... t estimators x ( n()1). We can be seen by conditioning on the observation . We can then calculate the plug-in ofe x t x ... t x extension of what was happening in the first ( ) case. The origin of the second term D 0 1 ( 1 ) ( 2 ) ( n ) n looking at the observed empirical distribuand the VaR and uncertainty in the VaR estimate increases the uncertainty the plug-inby estimators ˆare V aˆR D by xE(kn ,)V. aˆxR(ni )) ,has on E D k approximate kCT ˆconditioning n ( C T the pair an multivariate normal distribution. k n k ˆ 1 can be seen the observation of the estimated . We can then write V a R estimate. This is the origin of the second term. x CTE is clearly easier to work with than VaR. plug-in are If the sam size large enough, andthat a for few very technical conditio ns hold, then 1.when looking at the empirical distribution. practical expressions 1expressions estimating DThe 1so tion. that practical •other Thethe VaR and CTEfor estimators are positively by looking at theby observed empirical distribution. Let that expressions 1Simple pleso Dxestimators Intuitively, we observed can that we estimate CTE we are both CTE Example European ktheisLet Csay TˆEthe i practical 1C– ¦ n D (i ) , TˆE xPut by at the observed empirical distribution. Let D ˆ ˆ k, n ¦ n D n ( i )looking ˆ VAR [ C T E ] E { VAR [ C T E | V a R ]}1 ˆ k we would have to find a way to estimate 1 In termsand of this notation statistical theory has the following things to say ˆRare i the 1 pair (Cestimate TE ) has anthree n n nthe the approximate normal distribution. the VaR and uncertainty in VaR the ofcorrelated. the CTE kD i 1 multivariate ˆE uncertainty n , Va nincreases This makes intuitive sense. for the plug-in estimators C T x , ¦ the plug-in estimators are n ( i ) ˆROne term.to[xC(ktest the plug-in way above toVAR pick that simple enoughThis that we ˆThe ˆRnis]}. n DVAR Tˆ)E. the CTˆ.Ethe Rn ]}isEuropean estimators {an EPut [example C TˆEasymptotic Vare aˆorigin Rn D of x the .Va kn |iV estimate.estimators This is the second apply the formula. can be v 1–a nV]a n | Va ˆformulas {DVAR presented RnESimple x[Example plug-in are ) 1 (kk) k can get closed1 form expressions(kfor all the relevant risk measures. We can then perform 2 k CTˆE n D x , Intuitively, we can say that when we estimate th of the distribution. ¦ ˆ (i ) V a R D x . ˆ terms ofof thisthis notation statistical theoryof thetheory following three things to nfollowing (k ) to If the sam ple size ishasthis large a few other very technical conditio nsestimators then 1. Csimulations TEand xmodel InInterms notation statistical has the things k enough, •the Tsay he variance ofabove the CTE estimator has two1is simple onsay the to see how well variance perform. In the ˆformal paper In terms notation statistical theory has the three things tohold, say i 1 ¦ n D following ( i ) , three Cexample Tin Ethe DCTE x(i ) ,en One way to test the formulas presented to pick anboth that and the VaR and VaR estimate 1 isuncertainty k ¦ nthe Intuitively, we can say that when we estimate the CTE we are estimating i 1 2 ˆE n ,V chose statistical the example of annormal “in-the-money” European Putall option at the D“obvious” 0.95 confidence ˆˆRRn In thiswe notation theory has the following three things to say aa D terms x (kan . approximate k ithen V pair (C–TThe multivariate distribution. terms. The term therisk )of 1term. can get closed form expressions forfirst the relevant measures. We can perfc n ) hasPut Simplethe Example European This isisCTE the origin theextensecond xestimate. The VaR and estimators are positively and the VaR and in the VaR increases the uncertainty ofofthe CTE level. Here isuncertainty an editedsimulations excerpt from theestimate paper.to ˆ V a R D x . If the sam ple size is l arge enough, and a few other very technical conditio ns hold, then 1. n (k ) on the model see how well the variance estimators perform. In the fo If the sam ple size is l arge enough, and a few other very technical conditio ns hold, then 1. In terms of this notation statistical theory has the following three things to say ˆ sion of what was happening in the first ( = 0) V a R D x . n1 (k ) If the(Csam sizeanisapproximate large1.enough, and a few other very technical conditio ns hold, then 1. the estimate. This isisthe origin of the second term. ˆE n ,Vple we chose example of an “in-the-money” European Putdate option at the D 0.95 c Tto 1 ˆ ple sam size lapproximate arge enough, and athe few very technical conditio ns hold, then n ) has CIfTˆthe E ) has Onepair test the formulas presented above isnormal pick an example that isthesimple enough that we thestatistical pair (multivariate an multivariate normal distribution. Inway terms ofaˆRthis notation theory has thedistribution. following three things to other say AR European put option gives the holder right to sell the underlying asset on the maturity for the specified n , Va n to x The variance of the CTE estimator has two te level. Here is an edited excerpt from the paper. ˆ In terms of this notation statistical theory has the following strike E n ,V aˆRn ) has theget pair approximate distribution. sam ple size(C is T large enough, andfor aan few other technical nsnormal hold, then 1. If the can (very C TˆE nprice. ,V aˆmultivariate Rnmeasures. ) hasconditio the pairrelevant an approximate multivariate normal distribution. closed form expressions all the risk We can then perform Example Simple – The European Put extension of what was happening in the first ˆ 2 1 perform. In the formal paper (CTE n , VaˆR the pair hasmodel an approximate normal estimators distribution. simulations on)the to see howmultivariate well the variance European put option gives the holder the right to sell the underlying asset on the maturity 3 date fo is largeExample enough, –and a few other1Astrike very technical conditions thenby conditioning on the observati 1. If nthe sample sizeSimple canhold, be seen The European Putprice. we chose the example of an “in-the-money” European Put option at the D 2 0.1. confidence If way the sam plethe size is largepresented enough,2 and a few 95 One to test formulas above is toother pick (CTˆexcerpt E n , VaˆRfrom pair anpaper. approximate multivariate normal distribution. n ) has level. Herethe is an edited the ˆ ˆ ( C T E , V a R ) the pair has an approximate multivaria get closed form expressions for all relevant 2 thethat n n is simple One way to test the formulas presented above is2can to pick an example that enough we ˆ ˆ VAR[CTE n ] E{VAR[CTE n | Vrisk aˆR w Page 38 2 how simulations on the model see well the variance ne 1 can get closed form expressions for all the relevant risk measures. We cantothen perform A European put option gives the holder the right to sell the underlying asset on the maturity date for the specified > ^ @ ` ¦ ¦ ˆ ˆ ˆ ˆ ˆVaRCoˆv0.09% CTE ,VaR >PˆTX%) T ZFSE @FSE CVoˆavˆRCTE ,VaR FSE0.22 VaˆRf CTE (95%)f (x7.72 FSE SC TTE Cvariable TS1.01 E(eV 95 TaˆER 0.19 RatVapproxima 0 CV Minimum ) refers ) of the ion f X (VaR the probability the one cantoexecute to testdensity (and confirm) therandom validity of athe te formulas for any X VaR. fˆ VaR T in V TT Z @10 years with a strike price of X August 2008 w Risk Management Toisbe more at specific, assume the4.50 particular application. An example given the end this article. 13.70 1.63 2.42 n/a Average T n/a matures Soption S e >P1.91 13.80 n/a 13.80 4.39 4.39 2.37 n/a n/a Closed Form Closed Form >PT V T Z0.40% @ 2.37110 . S T S e The current stock price is S 100 and assumed to follow a lognormal return process with omments are in 13.67 Trial Maximum where Z order isFirst a standard Normal .54 13.67 5.09 1stock .54 1.65 1.40 1.65 0.49% First Trial variate with mean andis,sampling unit variance. Ptools 81 % and the1.40 pricecontrol at 5.09 maturity is te given by: 0.49% Vas importance 15zero % . That such and varia 3. Some standard variance reduction where Z is a standard Normal variate with mean zero and unit variance. 1.63 0.18 0.77 1.06 0.19% methods can be Deviation adapted to thewere CTE problem improve, the with3.07 Std where Zsometimes is1.76 a standard Normal variate mean0.22% zero and unit variance. 0.22% 14.93 1.95 14.93 3.33 then 3.07 1.95 dramatically, 3.33 4.91 4.91 Trial Last Trial TE is clearly easierLast to work with than VaR. If we using VaR as ato risk measure >PT will V T Z @ fcomput ) inational precision of atoCTE estimator for adensity given These beCTE we wish to e would have to a way estimate the probability order to cost. S techniques T S ewhose X (VaR , the random variable Using a find continuous discount rate of G 6 % 7.72 1.01 7.720 0.19 1.01 0 0.22 0.19 0.09% 0.22 0.09% Minimum Minimum ply the asymptotic formula. This canUsing be very difficult in practice, especially inrate the tails the topic of a second article CTE variance. the random variable CTE wevariable wish towhose CTE we wish to aoncontinuous discount G 6% , discount , the random Using aofcontinuous rate of G 6%whose key13.70 takeaways from this table are: thecalculate distribution. isAverage thenSome the present value payoff function: 1.63 13.70 4.50 1.91 1.63 4.50 2.42 1.91 0.40% 2.42 0.40% Average where Z ispresent a standard Normal variate with mean zero and unit variance. calculate is then the value payoff function: calculate is then the present value payoff function: SoVaR what are estimators the approximate formulas? First, need some notation. we13.05 want to 7.31 3.65% he and CTE are positively correlated. M This makes intuitive sense. 18.89 2.27 we 18.89 9.17 7.31 2.27Suppose 9.17 13.05 3.65% Maximum aximum origin Any trial provides abereasonable estimate of what happen if the simulation wereto case. The ofgiven the second term Using spreadsheet werandom can generate the variable whose CTE we wish Using arandom continuous discount G software, 6% , would GaTcan Prate T V ofTcumulative Z itional Tail Expectation of variable X, with distr estimate the Cond ibution V T Z C Deviation max 0,1.63 X S eeG T 0.77 V T ZG T 0.77 0.19% 1.76 0.18 0function: 1.76 0.19% Deviation he variance of Std the CTE estimator has two1.63 terms. Std The firsteterm is “obvious” eFrom max S e PT1.06 max ,X Svariable. e PCT1.06 calculate isthe then thenpresent value payoff but a0.18 sample of there isconsiderable variability, for nterm 1000 by) at conditioning thewith observation ofof size =C 1000 samples of this this 0, X especially Pr^X of d what x`seen F ( xhappening therepeated, Don Thus, we want to calculate the conditional expectation tension was inlevel the first ( D 0 ) case. The origin the second n be seen by conditioning on the of the estimated then writewe can calculate the plug-in estimathe estimated then write VaˆR . We cansample, Vobservation aˆR. We can PT V T Zgenerate T key takeaways from Some this key table takeaways are: from this table are: Using spreadsheet software, this variable. From thi 1000 samples Using Some spreadsheet software, we can generate samples this this n 1000 e G and max 0, X S we variable. esamples theCgenerate CTE VaR using thecan formulas spreadsheet software, wefor can ofFrom thisn variable. From of this nof 1000 CTE (D )TˆE E > X | X ! qtors ˆE |Using D@ ˆ ˆ VAR[CTˆE n ] E{VAR [ C T V a R ]} VAR { E [ C | V a R ]}. n CTE n n n The plug-in estimator is biased below the true closed form value of 13.80 (i.e., average sample, we can calculate the plug-in estimators for the CTE and VaR using the formula sample, we can calculate theweplug-in estimators for the CTEToand VaR the VaR formulas developed earlier. estimate theusing probability sample, can calculate the plug-in estimators for the CTE and using the formulas ˆ developed earlier. To estimate the probability density use the estimator: f (VaR ) Using spreadsheet software, we can generate samples of this variable. From this n 1000 T Edeveloped isprovides 13.70). However the bias much smaller than the error. given trial Any asmallest reasonable given trial estimate provides of what reasonable happen estimate ifofthe simulation would happen were if the simulation were tuitively, we can say when we estimate theCTE we areestimate estimating both theisCTE developed earlier. To estimate the probability density use estimator: f a(VaR )would earlier. Tovalue estimate the probability density use the estimator: fsampling (VaR ) what weC can say that we density use thethe estimator: where qD isIntuitively, theAny Dthat -quantile, defined aswhen the sample, we canofsatisfying calculate d the VaR and uncertainty in the VaR estimate increases the uncertainty the CTE the plug-in estimators for the CTE and VaR using the formulas with butn with aestimate sample of n density therethe issize considerable variability, there) isuseconsiderable especially forvariability, especially for 1000 1000 CTE arebut estimating bothrepeated, the size CTE of and timate. This isthe therepeated, originwe of second term.a sample developed earlier. the estimator: f (VaRquite theThe asymptotic variance formulaTofor the CTEprobability estimator performs well on [average (i.e., ˆ (VaR) the V VaR aˆR.and uncertainty in thePrVaR V^aˆXR.!estimate f [ [fˆ (VaR qD ` 1 D standard 1 1 ˆˆ(EVaRempirical ) average FSE CfT deviation of CTˆE ). Fˆn (D ) Fˆn (D [ ) ) ˆ 1 (D ) Fˆ 1[(D [ ) 1 1 xample – Theincreases European the Put uncertainty of the CTE estimate. F ˆ ˆ ˆ n (VaR Fnbelow (D )estimator theFntrue (Dfis [ )n) form 1 value ThisThe CTE plug-in estimator term. TheisCTE biased plug-in closed biased true 13.80 closed (i.e.,form average value of 13.80 (i.e., average ˆ 1of is the origin of the second Fˆnbelow (D ) Fthe n (D [ ) The D -quantile is often called Value-at-Risk (VaR) and is used extensively in the financial The VaR plug-in estimator is biased high (average is 4.50), but again the bias is much smaller to test the formulas presented above is to pick an example that is simple enough that we 1 ˆ ˆ CTE is 13.70). However Cthe TEWe bias iscan13.70). is much However smaller thanbias sampling much smaller error. than the sampling error. with .the is [ the osed form expressions for all the relevant risk measures. then perform 100 management of trading risk over a fixed (usually short) time horizon. 1 error. thanwith the sampling . [ 1In the. formal paper [estimators ns on the model to1see how. well the variance perform. with 100 with [Simple Example — 1 100 variance option The formula theconfidence variance CTE estimator formulaperforms for the CTE quiteestimator well on average performs (i.e., quite well on average (i.e., 100asymptotic the example of anThe “in-the-money” European Put atasymptotic the D for 0.95 We then calculate theFormula Standard x1 , ators x2the ,..., xisFormula of size n ˆisError The approach to this problem is to startfor with random sample The European Put sample covariance alla 1000 pairs ofcalculate estim theThe which higher (lower) than re is anstandard edited excerpt from paper. We can can then Standard (FSE) of each estimator as n 2.37, ˆ ˆ ˆ standard FSE CTdeviation E the empirical ofofC Tstandard E ). deviation of estimator CTE ). as average FSE CTE empirical average Error (FSE) each estimator asof We can thenfrom calculate Formula Standard Error (FSE) oftoeach We can then calculate the Formula Standard Error (FSE) each estimator as the estimated covariance the first (last) trials, but close the mean of all covariance then sort the sample in descending order to obtain the order statistics from the model and an put option gives theway holder right the to sellformulas the underlying asset on the maturity for the specified Onethen tothetest presented abovedate can calculate the Formula Standard Error (FSE) offor each estimatorparam as eters ce. x(We . estimators. We can then thethat plug-in estimators the required t ... t xan X ( k ) ) the X ( k ) ) 2smaller D bias ( X (but (CTˆEismuch 1) t x ( 2is ) to ( n ) example pick thatestimator is simple enough The VaR plug-in calculate The is VaR biased plug-in high (average estimatorisis 4.50), biased but high again (average the bias isVAR 4.50), is much smaller again 1) ,..., CTE , VAR ( X FSE D (CˆTˆE X ( k ) ) 22 3 (1) ,...,(X (k ) ) ) k (CTE VAR ) ( X (1) ,..., X ( k ) ) D (CTE X ( k ) ) , n (1 D ) weatcan closed form expressions allsampling the D 1error. than the sampling error. thanforthe by looking theget observed empirical distribution. Let so that practical expressions for FSE 2 n D ( 1 ) FSE CTE ( ) , ˆ VAR Xn ( kdocumented Xn( kpaper ( X (1) ,...,are ) D )using the same methodology A measures. number ofWemore examples relevant risk can realistic then perform ) ) D (CTEinthe )(1 D 1 D 1 FSE (CTE ) allsample ,) higher D Dis 1 1 The sample covariance for The 1000 pairs covariance of estim for all 1000 pairs of estim ators is 2.37, which is (lower) 2.37,examples than which, is were higher (lower) than The the plug-in estimators simulations the modelabove. to see how well the asonare described In each case the asymptotic as expected. , ators FSE FSE (VaR n ((VaR 1 )1Dtheory ) ˆ Dworked ˆ n f (VaR ) D 1 n k f VaR ( ) the estimated the from estimated the covariance (last) from butthe close to(last) the problems trials,ofbut all close covariance to the mean of all covariance 1 first )trials, ,mean FSE first variance estimators perform. the formal chosen tocovariance test a wide of possible behaviours practical facing insurers. CIn TˆErange ,(VaR ˆ (VaRDand ˆE X ) n D 1 ¦ x ( i )D ˆE nX ) D 1 ( C T D f ) ( C T (k ) (k ) k i 1 Coˆv(CTE ,VaR paperestimators. we chose the example ofestimators. ) an “in-the,) FSE (VaR . ) . Coˆˆv(CTE , VaR ˆ (VaˆR ) ) f (VaˆR n D (CTˆE n X n f (VaR money” European Put option1 at 0.95 ) Vaˆthe Rn Dfˆ= x(k )). (k ) , VaR ) examples . Cfrom oˆmore v(CTE Athis number of more realistic A number examples are realistic inn to the paper are documented using the same in the methodology paper using the same methodology confidence level. Here is antheory edited excerpt ˆsay ˆE documented ˆ In terms of notation statistical has theof following three things ( ) f V a R ( ) C T X D Table shows the results oflast) two trials (first (k ) 1 1 shows the results two trials and also theand results of repeating the entire the as paper. described above. In each as described case asymptotic above. Table In each theory case worked the asymptotic asand expected. theory Theworked examples expected. were examples were the entir 1 shows the(first results of two trials (first andaslast) and alsoThe the results of repeating , VaR )Tablethe .of Coˆv(CTE simulation 1000 times. The tablethe also shows exact values of the CTE and VaR for this ˆVariance last) and also results of the repeating the entire ˆ A More Practical Example – Verification ( ) n f V a R simulation 1000 times. The table also shows the exact values of the CTE wide range chosen ofproblem, possible to test abehaviours wide range and of practical possible behaviours problems facing and practical insurers. problems facing insurers. and VaR for thi the samto pletest sizea is large enough, and few other very technical conditio ns hold, 1. If chosen which can be calculated from closed form expressions are given the paper. Table 1 shows the aresults ofsimulation two trials (first and last) and alsothen thethat results of in repeating the entire 1000 times. The table also shows the problem, which can be calculated from closed form expressions that are given in the paper. Topair be more specific, assume the option matures ˆ (CTˆS Euppose , V a R ) the has an approximate multivariate normal distribution. n nsimulation values the CTE and VaR for this probtimes. Theexact table alsotoofshows exact values of the CTE VaR for this you have a1000 model that takes all night run n =the 1,000 scenarios. It would be and impractical in1 T = 10 years with a strike pricetrials of X 1: =(first 110. Table shows therepeat results of and last) and also the results ofReduction repeating the entire Monte Carlo Simulation without Variance ume theThe option matures in years with a strike price of . Ttwo 10Table X 110 to the run process hundreds oflem, time in order test the validity ofthat the variance formulas as which can beto calculated from closed form problem, which can be calculated from closed form expressions are given in the paper. current stock price is S = 100 and asCTE(95%) for a 10-year European Put Option (1000 Trials), X=$110, S=$100 Table 1:with Monte Carlo Simulation without Variance Reduction simulation 1000 times. The table also shows the exact values of the CTE and VaR for this To be more specific, assume the option matures in years a strike price of . T 10 X 110 described above. expressions that are given in the paper. sfic,Sassume and assumed to follow a lognormal return process with 100 CTE(95%) a110 10-year European Put Option (1000 Trials), X=$110, S=$100 sumed to follow aprice log normal return process with the option matures years with a strike price of for . 10 Xreturn The current stock is in and assumed to follow a expressions lognormal process with S T100 2 paper. problem, which can be calculated from closed form that are given in the More Practical Example Athe More –stock Variance Verification Example – by: Variance Verification That is,P the price at maturity isPractical given by: Table 1: Monte Carlo Simulation without Variance Reduction price is and assumed to follow a lognormal return process with SA 8100 %stock and . stock price at maturity is given That is, price at V 15 % A more practical process for confirming the asymptotic variance formula, which we call To be more specific, assume the option matures in T 10 Table years with 1: a strike price of X 110 . 4 CTE(95%) 10-year (1000 Trials), X=$110, S=$100 the stock price at isisagiven by: European 15% . That maturity is given by: Theis, current stock price isVerification, assumed to follow a lognormalPut returnOption process with S maturity 100 andfor Variance as follows: SPuppose you have aThat model Sthe uppose that takes you all have night to run that n = takes 1,000 all scenarios. night to It run would n = 1,000 be impractical scenarios. It would be aimpractical Monte Carlo Simulation without Variance Reduction CTE(95%) for 10-year > PT Va Tmodel Z @ 8 % and . is, stock price at maturity is given by: V 15 % T S e Variance Reduction Table 1: Monte Carlo Simulation >PT tohundreds V repeat T Z S @without European Put Option (1000 Trials), X=$110, S=$100 to repeatSthe run S process the of time run process in order hundreds to test the of validity time in of order the to variance test the formulas validity of as the variance formulas as T e >PT V T Z @ Put Option (1000 Trials), X=$110, S=$100 CTE(95%) for aS10-year T S eEuropean > PT V T Z @ described above. described above. 5 S T mean e zero and unit variance. where with whereZZisisaastandard standardNormal Normalvariate variate with Smean 4 FSE VaˆR CoˆvCTE ,VaR CTˆE (95%) FSE CTˆE fˆ VaR VaˆR zero and unit variance. A A more process for more confirming practical the process asymptotic for confirming variance formula, the asymptotic which variance we call formula, which we call dard Normal variate with mean zero and variance. where is practical a standard Normal variate with and unit variance. ormal variate mean zero and unit , the random variable whose CTE we wish to Using aZ with continuous discount rate ofunit Gmean 6variance. %zero n/a 4.39 n/a 2.37 n/a Variance Verification, isVariance as follows: Verification,Closed is as Form follows:13.80 4 calculate is then the present value payoff function: 13.67 1 .54 5.09 1.40 1.65 0.49% First Trial Using discountrate rateof of G 6% , the the random variable whose CTE we wish to Using aacontinuous variable whoseLast CTE we wish to ous discount ratecontinuous of G 6discount % , the random 14.93we wish 1.95 to 3.33 3.07 4.91 0.22% Trial CTE calculate is then the present valuewe payoff function: , CTE the random variable whose count rate ofvariable G 6 % random whose wish to calcu G T PT V T Z he present value payoff function: C e max 0, X S e Minimum 7.72 1.01 0 0.19 5 0.22 0.09% 5 late payoff is then thefunction: present value payoff function: ent value 18.89 2.27 9.17 7.31 > >@ @ > 13.05 3.65% > @ @ > > @ @ 13.70 1.63 4.50 C e max 0, X S e PT V T Z Average G T PT generate V T Z Using C spreadsheet software, we can samples of this variable. From this n 1000 18.89 2.27 9.17 M aximum e max 0, X S e sample, we can calculate estimators for the CTEofand usingFrom the 0.18 formulas Using spreadsheet software,theweplug-in can generate this VaR variable. this 1.63 1.76 Std Deviation G T T V T Zn 1000 samples developed To0estimate density useVaR the estimator: Csample, e weearlier. canmax , Xthe plug-in Sthe probability e Pestimators calculate for thef (VaR CTE) and using the formulas software, we can generate samples of this variable. From this n the 1000 developed earlier. To estimate probability density estimator: f (VaR ) use Some keythe takeaways from this table are: > G T > @ @ 1.91 2.42 0.40% 7.31 13.05 3.65% 0.77 1.06 0.19% et calculate the plug-in estimators forfˆ (VaR the) CTE and [VaR using the formulas continued on page 40 1 1 . To estimate the probability f (VaR )Fˆnuse ware, we can generate density of this variable. From this estimate of what would happen if the simulation were n 1000 (D the )[ Fˆestimator: [ ) given Any trial provides a reasonable n (D fˆ (VaR ) samples 1 1 Fˆn (D ) Fˆn (D [ ) repeated, but with adate of n strike there is considerable variability, especially for 1000price. European putestimators option gives the holder right toCTE sell the underlying asset onusing the maturity for thesize specified ate the Aplug-in forthethe and VaR thesample formulas ˆR V a . with [ ˆ1 . [ f (VaR) use the estimator: stimate the density f (1100 VaR.) with [probability 1 1 1 100 Fˆn (D ) Fˆn (D [ ) The CTE plug-in estimator is biased below the true closed form value of 13.80 (i.e., average CTˆE estimator is 13.70). as However the bias is much smaller than the sampling error. We can then calculate the Formula Standard Error (FSE) of each We can then calculate the Formula Standard Error (FSE) of each estimator as The asymptotic variance formula for the CTE estimator performs quite well on average (i.e., fˆ (VaR) Page 39 w ˆ FSE 1 1 VAR X X Tˆ(ECT X) 2( kC) )Tˆ2E empirical standard deviation of CTˆE ). ( ,..., ) average (1) X EX )( k) D D (C Fˆ ( ) Fˆ VAR ( (X ,..., ) [ n D (CTE )n FSE D [ (1) (k ) (k ) , CTE (95%) Closed wForm Risk Management August 2008 13.80 First Trial 13.67 Last Trial 14.93 Minimum 7.72 Average 13.70 Closed M Form aximum 13.80 18.89 First Trial Std Deviation 13.67 1.63 CTˆE (95%) FSE CTE VaˆR ˆ CTE (Closed 9513.80 %) Form FSE Closed Form n/a 4.39 n/a CTˆE13.80 n/a f VaR CoˆvCTEn/a FSE Vaˆn/a R4.39 ,VaR VaˆR4.39 n/a 2.37 2.37 n/a First Trial n/a113.67 13.80 4.395.091.54 n/a1.405.09 13.67 .54 Chairman’s Corner Variance of the CTE 5.09 1.40 Estimator1.65 0.49% Closed First FormTrial 1.54 3.33 LastEstimator Trial CTE 1.01Minimum 0… 1.65 n/a0.49% 4.91 0.49% 0.22% 14.93 5.093.331.95 1.403.073.33 Last Trial 1.541.95 13.67 14.93 3.07 4.91 0.22% 7.72 0 VaˆR4.91 0.22 Minimum 0.19 CoˆvCTE0.22% ,VaR CTˆE (951.95 %)1.01 FSE 3.33 CTˆE 01.01V3.07 aˆR 0.19 FSE 14.937.72 0.22 0.09% FirstLast TrialTrial 1.95 1.40 2.37 1.65 3.07 1.654.91 fˆ VaR2.37 n/a 0.19 0.22 0.09% ˆ VaR 4.50 V13.80 page CTE FSE aˆfrom R1.01 C oˆv39 ,VaR 13.70 1.91 FSE C TˆE Closed f0.19 Average VaˆRForm Minimum w continued 13.70 1.63 1.91 2.42 n/a 04.501.634.39 n/a 0.22 1.63 Average 4.507.72 1.91 2.42 0.40% 2.42 0.40% 2.370.09% 13.05 3.65% 1.650.40% 18.89 2.37 7.31 Maximum n/a 4.39 n/a n/a7.319.171.40 2.42 13.70 1.632.27 4.509.172.275.09 1.91 Average 13.05 Maximum 13.67 1.54 First Trial18.89 7.31 2.27 9.17 13.05 3.65% 1.63 1.06 0.18 0.19% 0.77 1.06 Std Deviation 18.89 2.27 9.17A1.76 7.310.771.76 13.05 3.65% 1Some .54 5.09 1.40 0.49% M aximum 1.63 1.06 4.91 0.19% Std Deviation 14.93 1.95 3.07 for 0.18 1.76 0.77 Last Trial key takeaways from this table are:0.18 1.65 more 3.33 practical process confirming the 1.63 3.077.72 0.18 1.76asymptotic0.22% 0.77 0.19% 14.93 1.95 4.91 Std Deviation 3.33 Last Trial which we0.22 call 1.01 0variance formula, 0.19 1.06 Minimum Some keya takeaways from this table are: Some key takeaways from this table are: Some key takeaways from this table are: • A ny given trial provides reasonable esVariance Verification, is as follows: 7.72 1.01 0 0.19 0.22 0.09% Minimum 13.70 1.63 4.50 1.91 2.42 Average timate of what would happen if the simuSome key takeaways from this table are: 13.70 1.63 4.50 1.91 2.42 0.40% Average 18.89trial provides 2.27 a reasonable 9.17 7.31 13.05 Maximum ofAny estimate of what ifwould happen if we th Any given trial provides a reasonable estimate what would the1. simulation were were repeated, butgiven with a happen sample if estimate As part model development, orthe in an off lation Any given trial provides a reasonable of of what would happen simulation 18.89 2.27 9.17 7.31 13.05 3.65% Maximum 1.63 0.18 1.76 0.77 1.06 Std Deviation repeated, but size with aestimate sample size of is thereaifis considerable nwould 1000 repeated, but with a sample sizerepeated, of there considerable variability, especially for n 1000 size isis varipeak time, generate (once) larger sampleespecially of variabil butprovides with a considerable sample of there considerable variability, n 1000 Any given trial a reasonable of what happen the simulation were f 1.63 0.18 1.76 0.77 1.06 0.19% Std Deviation ˆ VaˆR. V a R . ability, especially for say N = 5,000 scenarios. Let CTE be the esˆ V a R . N repeated, but with a sample size of n 1000 there is considerable variability, especially for Some keyestimator takeaways from this table are:timated CTE based on this large sample and • The plug-in is biased below ˆR.CTE The CTE plug-in estimator isVa biased belowplug-in the true closed form value of 13.80 (i.e., average value The CTE plug-in estimator isthe biased below standard the true closed form valueavera of 1 The CTE estimator is biased below true closed form value offor13.80 (i.e., the true closed form of 13.80 (i.e., averFSE the formula error a given Some key ˆtakeaways from this table are: N ˆ CTE is 13.70). However theage bias isˆ much smaller than the sampling error. C Ttrial E isthe 13.70). However the bias isthan much smaller than the sampling error. 13.70). However is much confidence level α. TE isis 13.70). However the bias is much smaller the sampling The CCTE plug-in estimator isbias biased the true closed form value oferror. 13.80 (i.e., average Any given provides abelow reasonable estimate of what would happen if the simula smaller the sampling 2.on Fnrom the large sample of size N, drawvariability, m=100 buterror. with a sample size of there is considerable 1000 CTˆE for isthan 13.70). the bias isvariance much smaller than sampling error. Any given trial provides a reasonable estimate of what would happen ifwell the simulation were The asymptotic variance formula therepeated, CTE estimator performs quite average (i.e., However The asymptotic formula forthe the CTE estimator performs quiteespe we TheThe asymptotic variance formula for the CTE estimator performs quite well on average (i. • asymptotic variance formula for the CTE random sub-samples of size n=1,000 without ˆE ). Vdeviation aˆR.there repeated, butFSE with size of is considerable variability, especially for n 1000 CTaˆEsample empirical standard of C T average ˆ ˆ FSE C T E empirical standard deviation of C T E ). average average FSEvariance CTˆE wellformula empirical standard deviation of performs CTˆE ). quite well on average (i.e., The asymptotic for the CTE estimator estimator performs quite on average (i.e., replacement. VaˆR. ˆ ˆ The plug-in estimator is biased below the true closed form value of 13.80 (i.e 3. or F ofEpeak the mtime, sub-samples calculate a The VaR plug-in estimator isaverage biasedhigh (average isempirical 4.50), butstandard again thedeviation bias isor much smaller FSE TAs E part of T ). (average 1.CCTE of model development, ineach an C off generate (once) a larger samp The VaR plug-in estimator is biased high is 4.50), butisagain the bia The VaR plug-in estimator is biased high (average is 4.50), but again the bias much smal ˆ say N = 5,000 scenarios. Let CTE be the estimated CTE based on this large samp deviation of ). CTE estimate and an FSE estimate. Also Thethan CTEthe plug-in estimator true closed form value of 13.80 average sampling error. is biased below CTthe E is 13.70). However the bias is much smaller than the sampling error. N(i.e., than the sampling error. than the sampling error. bias VaR plug-in estimator biased high (average 4.50), but again the is much formula standard errorcheck for is a given D. FSE • The Theis VaR plug-in estimator isissampling biased high to see confidence whether ourlevel bestbias answer CTEN smaller N thethe CTˆE is 13.70). However the much smaller than error. 1000 pairs of estim The sample covariance for all ators is 2.37, which is higher (lower) than The asymptotic variance formula for the CTE estimator performs quite well on ave than the is sampling (average 4.50),buterror. againsample the biascovariance is much for lies in thepairs approximate 95% is The 1000 of estim ators 2.37, which is th h Thefirst sample covariance for sample allto 1000 pairs ofall estim the ators is 2.37, which isconfidence higher (lower) the estimated covariance from (last) trials, but close the mean of all covariance ˆ 2. s of size n=1,000 From the large of size N, draw m=100 random sub-sample ˆ The asymptotic variance formula for the estimator quite well on average (i.e., FSEestimated Cperforms TE empirical standard deviation of Ctrials, TIfE it). does average smaller thanCTE the sampling error. interval CTE+2xFSE. we set to thethe the covariance from the first (last) but close mean estimated covariance from firstof(last) close to the all covarian without replacement. sample covariance for 1000the pairs estimtrials, Thethe ators isbut 2.37, which is mean higherof(lower) than estimators. FSE CTˆE empirical of E all ). pairs average • Tstandard he sampledeviation covarianceestimators. forC allTˆ1000 of esCI (Confidence Interval) count to 1 and 0 estimators. the estimated covariance from the first (last) trials, but close to the mean of all covariance The VaR plug-in estimator is biased high (average is 4.50), but again the bias is muc timators is 2.37, is higher (lower) than otherwise.a CTE estimate and an FSE estimate. Also 3. which amples calculate For each of the m sub-s A number of more realistic examples are documented in the paper using the same methodology estimators. than the sampling error. The VaR plug-in estimator is biased high (average is 4.50), but again thebest bias issemuch smaller the estimated covariance from first (last) 4. U the standard of the the CTE esticheck to see whether our answer CTE lies indeviation the approximate 95% confidence N documented A number ofthe more realistic examples are in paper using the s as described above. In each case the asymptotic theory worked as expected. The examples were A number of more realistic examples are documented the paper using the same count methodolo than the sampling error. wefrom setinthe Interval) 2 u FSE . If it does trials, but close to interval the meanCTE of allrcovariance mates StepCI3 (Confidence to check the validity of the to 1 a The sample covariance for all 1000 pairs of estim ators is 2.37, which is higher (lo as described above. In each case the asymptotic theory worked as expected. T chosen to test a wide range of behaviours andIn practical problems facing insurers. as described above. each case the theory worked examples we Apossible number of more realistic examples areasymptotic documented in the paper as using methodology otherwise. estimators. asymptotic formula. Asexpected. we the will same seeThe shortly a the estimated covariance from first(lower) (last) trials, but close to the mean ofwere all c chosen tocase test athe wide range of theory possible behaviours and practical problems facin alldescribed 1000 pairs of aestim The sample covariance for as ators is 2.37, which is the higher than chosen toabove. test wide range of possible behaviours and practical problems insurers. In each asymptotic worked as expected. The examples simple adjustment needs to befacing made to this the estimated covariance from thetofirst (last) but close the mean ofand all practical covariance 4. trials, dard to deviation of the CTE estimates from Step 3 to check the validity o Use the stan chosen test aestimators. wide range of possible behaviours insurers. A number of more realistic examples are docunumber before problems comparingfacing it to the formula asymptotic formula. As we will see shortly a simple adjustment needs to be made t estimators. mented in the papernumber using the samecomparing method- it to the estimates. estimates. before A number of more realistic examples are formula documented in the paper using the same met as described above. In each case the A More Practical Example – ology Variance Verification as described above. In each case the asymptotic theory worked The exam A number of more realistic examples are documented in the paper using the same methodology asymptotic theory worked asPractical expected. The of applying The table below shows thessresults of applying The tab toasanexpected. inforce portfolio of U leMore below shows theExample results the above proce A – Variance Verification A More Practical Example – Variance Verification to test a wide range of possible behaviours and practical problems facing insurer as described above. each case the asymptotic theory worked as expected. The examples were Suppose you haveIna model thatexamples takes allchosen night to run n = 1,000 scenarios. It would be impractical ariable annuities with GMDB, GMAB and GMWB features. The book is slightly out of v were chosen to test a wide range of the above process to an inforce portfolio of U.S. A More Practical Example – Variance Verification to repeat the run process hundreds of time in order to test the validity of the variance formulas as money. 5000 real world (P measure) scenarios were generated and for each scenario the pr chosen to test a wide range of possible behaviours Sand practical problems facing insurers. uppose youthat have a model that all=night to runGMDB, n = 1,000 scenarios. It wou possible and practical problems variable with GMAB andimpractic Supposebehaviours you have a model takes all night totakes runpresent nannuities 1,000 scenarios. Itfees would be described above. value of guarantee benefits (claims) less the value of related was captured. repeat the run process hundreds time inThe order to istest thevariance validity of the var facing insurers. GMWBof features. book slightly out of theformulas to repeat runatomodel process hundreds ofnight time in test the validity of the Suppose you the have that takes all to order run n to = 1,000 scenarios. It would be impractical described above. money. 5000 world (P the measure) described above. A more practical process forto confirming asymptotic variance formula, we call The claims have been normalized so the mean CTE(0) is 1,000. Thevariance feesscenarios were scaled soastha repeat thethe run process hundreds of time inwhich order to test thereal validity of formulas Variance Verification, is as follows: A More Practical Example— TE(75) of the net (PV claims – PV fees) is zero. The example itself is for CTE(90). C described above. A More Practical Example – Variance Verification Aprocess more practical process the for asymptotic confirming variance the asymptotic variance A more practical for confirming formula, which formula, we call whic Variance Verification A More Practical Example – Variance Verification Variance is as2:follows: Variancevariance Verification Variance Verification, is Verification, as follows: Table A more practical process for confirming asymptotic which we It call Syou uppose you have a model thatthe takes all night to run n5formula, = 1,000 scenarios. would be im Suppose have a model that Variance Verification, is as follows: to repeat then run process hundreds of time in order to test the validity of the takesall allnight night Suppose you have a model that takes to run run = 1,000 scenarios. It would be impractical N = 5,000 Samples variance for D 90% described above. CTE = FSE 2214 to repeat the run process hundredsscenarios. of time in orderbe toimpractical test the validity of the variance formulas as N N = 159 It would described above. m = 100 Random Sub Samples of Size n = 1,000 to repeat the run process hundreds 5 A more practical process for confirming the asymptotic variance formula, which we call of time in order to test the validity of Variance Verification, as follows: FSE CI Count A more practical process for confirming the asymptotic variance isformula, which we call CTE the variance formulas as described Mean 2,111 346 94% Variance Verification, is as follows: above. First Last Min Max Std Dev'n Adjusted Std Dev'n w Page 40 2,004 2,242 1,558 3,120 316 354 5 355 353 262 376 30 1 1 0 0 2% =Std Dev'n/ (1-n/N)^1/2 The first thing we note is that if the FSE based on the sample of size 5,000 is correct then w expect an FSE of about 159 5 | 355 when dealing with a sample size of 1,000. The FSE estimates are clearly consistent with this, but the standard deviation of 316 (of the 100 CTE estimates) is not. One possible explanation for this discrepancy is sampling error, but more The table below shows the results of applying the above process to an inforce portfolio of U.S. variable annuities with GMDB, GMAB and GMWB features. The book is slightly out of the money. 5000 real world (P measure) scenarios were generated and for each scenario the present value of guarantee benefits (claims) less the present value of related fees was captured. August 2008 w Risk Management The claims have been normalized so the mean CTE(0) is 1,000. The fees were scaled so that CTE(75) of the net (PV claims – PV fees) is zero. The example itself is for CTE(90). Table 2: Variance Verification D 90% N = 5,000 Samples were generated and for each scenario the presOur variance verification test is therefore to CTE N = 2214 FSE N = 159 ent value of guarantee benefits compare empirical errorofestimate 354 from universe N=5,000 sowith they are not independent. Because the various CTE m = 100(claims) Randomless Subthe Samples of Size the nthe = same 1,000 are using some346. of the data they are positively correlated and so the set of present value of related fees was captured. theestimates mean formula estimate Wesame conclude more tightly clustered than if they were truly independent. This intuitive result CTE FSEasymptotic CI is Count thatestimates the theory appears to be working Mean 2,111 346easy 94% very to understand as n o N . The claims have been normalized so the mean for this model.1That is, the asymptotic variance First 2,004 from355 the same universe of N=5,000 so they are not independent. Because the various CTE Lastwere scaled so 2,242 353 1 2 CTE(0) is 1,000. The fees that of the Estimator agrees with “experiment,” It 262 isCTE possible to use theof meth our paper topositively show that,correlated in the large sample in data estimates are using some theods same they are and so thelimit, set ofthe Min 1,558 0 CTE(75) of the net (PV claims – PV fees) is zero. after adjusting for the non-independence of the orrelation of two sub-sample estimates is just c U n / N . A better estimate for theresult sampli estimates is more tightly clustered than if they were truly independent. This intuitive is Max 3,120 376 0 from the same universe N=5,000 sounderstand they are not Because the various CTE The example itself is for CTE(90). sub-samples. Std Dev'n 316 of 30 to 2% very easy asindependent. noN . 316 316 estimates are using some of error the same they are positively and so the of but3 whendata using a sample size of correlated 1,000 is therefore notset 316, | 35 n 1000 Dev'n/ (1-n/N)^1/2 Adjusted Std Dev'n 354 =Std 2 This intuitive result is estimates is more tightly clustered than if they were truly independent. 1 the It is possible use theismeth limit, in our paper The first thing we note is that if the FSE based The CI Counttoresult alsoods consistent with to show that, in the large1sample N for the 5000 very easy to understand as n o N . of two sub-sample estimates is just U n / N . A better estimate orrelation c thething sample sizeis5,000 correct theOur idea that theverification standard errors areto compare the empirical error estimate sampling Theonfirst weofnote that ifisthe FSE then basedwe onexthe sample of size 5,000 isformula correct then weis therefore variance test 354 with 316 316 2count | 355 pectan anFSE FSEofofabout about 159 It5 is when dealing working. The actual of 94 very expect whento dealing with a sample size of 1,000. FSE mean estimate 346 . ofis We conclude that the asymptotic possible use the meth our toThe show that, in theclose sample limit, ods informula paper error when using a sample size 1,000 islarge therefore not 316,the but3 theory appears to be workin | 354 n with 1000 estimates are clearly with the standard ofvalue 316That the 100 CTE thisestimates model. is,U(i.e., the nasymptotic variance of theforCTE agrees “experiment of two sub-sample is(of just correlation N . confidence A better estimate the Estimator sampling1 with a sample sizeconsistent of 1,000. Thethis, FSEbut estimates todeviation expected of 95 a/95% 1 estimates) is not. One possible explanation for this discrepancy sampling error, moreuniverse of N=5,000 from thenon-independence same so they are316 not independent. Because CTE afterisadjusting for thebut of the sub-samples. N the various 5000 316 are clearly consistent with this, but the standard interval). 3 estimates are using some of the same data they are positively correlated and so the set of testing shows this is not the case. error when using a sample size of 1,000 is therefore not 316, but . | 354 Our variance verification test is therefore to compare empirical error estimate 354 with the deviation of 316 (of the 100 CTE estimates) is estimates is more tightly clustered than were truly independent. This intuitive result is n if theythe 1000 ulaappears standardtoerrors are worfo Theformula CI Count result is also consistent with thethe idea that the form 1 1 mean estimate 346 . We conclude that asymptotic theory be working very easy to understand as . n o N e universeThe of not. N=5,000 so they are not independent. Because the various CTEis a biased standard deviation of the m=10 0 sub-samples estimate of the sampling error Table and 2 that N One possible explanation for this discrepFinally, it might appear from 5000 T he actual count of 94 is very close to expected value of 95 (i.e., a 95% confidence interval this model. That is, the asymptotic variance of the CTE Estimator agrees with “experiment”, using some of not the same data are positively correlated and so the set of it is tothey understand why. The various sub-samples (each ofevidence size n=1,000) were all ancyhard isthansampling but more testing shows there istherefore of small sample It ismaterial possible todrawn use the meth Our variance verification test is compare the empirical error 354 that, withinthe the large sample limit, the ods in ourestimate paper 2 to show more tightly clustered if they wereerror, truly independent. This intuitive result is after adjusting fortothe non-independence of the sub-samples. of two sub-sample estimates is just correlation / N . A better for the sampling understand as this noN mean formula estimate 346 . We conclude that theory appears toUbe nworking for estimate is. not the case. bias since mean ofthe theasymptotic 100Table sub-sample Finally, itthe might appear from 2 that there is evidence of material small sample bias si 316 316 6 universe of N=5,000 so they are not independent. Because the various CTE 3 this model. That is, the asymptotic of agrees with “experiment”, 2 error when usingEstimator a sample size of is 1,000 is that therefore 316, butstandard | workin 354 . thCI e mean ofresult the 100 sub-sample estimates 2,111, with an apparent precision of estimates isvariance 2,111, with an CTE apparent precision ula errors are to usesome the meth The Count is the also consistent with the idea the not form that, incorrelated the large and sample limit, odssame in our paper sing of the data they to areshow positively so the setthe of n 1000 adjusting theis non-independence ofofthe sub-samples. 1 1 ore tightly clustered than if is they This intuitive result two sub-sample estimates justwere U truly n / N independent. . A after estimate for thefor sampling T actual count is very close less to expected 95 (i.e., a 95%from confidence interval). deviation ofbetter the m=100 sub-samofhe316 which isismuch than / 100 | 32 ,94 which much less than thevalue valueof2,214 obtained of size N the sample 5000 nderstand as nThe . o N standard 316 316 sing a sample size of 1,000 is therefore not 316, but3 Our variance verification test is therefore to compare the empirical error estimate 354 with | 354 . 5,000. However, this analysis misleading, againerrors due to the non-independence of thethesubples is a biased estimate The of large the sampling erroris also theconsistent value 2,214 from theisform sample of nCount 1000 ula standard are working. withobtained the idea thatTable the o use the methods in our paper 2 to show that, in the the mean formula estimate 346 . We conclude that the asymptotic theory appears to be working for 1CI sample 1limit, result Finally, it might appear from 2 that there is evidence of material small sample bias since samples. N for 5000 wo sub-sample estimates is just U nto/ Nunderstand . A better the sampling and it is not hard why. The varisize 5,000. However, analysis misleadthisthis model. That is,isthe asymptotic variance of the CTE Estimator agrees with “experiment”, The estimate actual count of 94 is th very close to expected value of 95 (i.e., a 95% confidence interval). e mean of the 100after sub-sample is 2,111,ofwith an apparent precision of verification test is therefore to compare the empirical error 354 with the 316 estimate 316 adjusting forestimates the non-independence the sub-samples. 3 ng a sample size ofsub-samples 1,000 is therefore not 316, |all 354 . ofbut size n=1,000) were a estimate 346ous . We conclude that the(each asymptotic theory appears to 1000 be working for ing, again due to the non-independence of the n 316 / 100 | 32 , which is much less than value 2,214 obtained theissample When sam ples are positively correlated the variance of the samplefrom mean larger of size 1 with 1“experiment”, That is, the asymptotic of thesame CTE Estimator agrees Finally, might Table 2 that2 there material small bias The is CI evidence Count resultof is also consistent with sample the idea that the since formula standard errors are working. drawn variance from the universe ofitN=5,000 so from sub-samples. N 5000appear gerification for the non-independence of the sub-samples. 5,000. However, this analysis is misleading, again due to the non-independence of interval). the precisio subVARestimates ( x ) V [ Uis (1actual Uwith )count / m]an would be.ofA95better estimate for test is therefore to compare the empirical estimate 354 100 with the T2,111, he ofthan 94 is it very close tootherwise expected (i.e., a 95% confidence thBecause e error mean of the sub-sample apparent precision of value are notthat independent. various estimate 346 .they We conclude the asymptotic theory appearsthe to be working for samples. ula tatresult is asymptotic also consistent with the ideaCTE thatEstimator the form 1000 1000 is, the variance of the agrees witherrors “experiment”, 316 /standard | 32aredata , working. which When is muchsamples less than are the value 2,214appear obtained from the sample of size material small sample bias since Finally, it might from CTE estimates are using some of100 the same the 2,111 number ispositively then which is not materially differ 346 correlated (Table 1 2 that)there / 100is evidence | 158 of unt the of 94 is very close to expected value of 95 (i.e., a 95% confidence interval). for non-independence of the sub-samples. thagain e meandue of theto100 sub-sample estimates is 2,111, with suban apparent precision of 5,000. However, this analysis is misleading, the non-independence o f the 5000 5000 they are positively correlated and so the set of the variance sample correlated mean is larger When samplesofarethe positively the variance of the sample mean is larger 316 / 100 | 32 , which is much less than the value 2,214 obtained from the sample of size ula standard working. result is also consistent idea that the form ght appear from Table 2with that the there is evidence ofsamples. material small errors sampleare bias since 2 2,214 value.otherwise Thus, while ther from theVsampling error in] the e is some evidence of small sam estimates is more tightly clustered than if they than would nt of 94 is very close to expected value of 95 (i.e., a 95% confidence interval). VAR ( x ) [ U ( 1 U ) / m it would be. A better estimate for the 5,000. However, this analysis is misleading, again due to the non-independence of theprecision subhe 100 sub-sample estimates is 2,111, with an apparent precision of bias in that 2,111 < 2,214, there is not enough data here to quantify it. The bias is lost in th samples. estimate for the , which iswere much than the value 2,214 of obtained from thesample sample of size truly independent. This intuitive result is otherwise be. A better ht32appear from Tableless 2 that there is evidence material small bias since 1000 1000is larger When sam are positively correlated the variance of the sample mean error. This usually thecase. analysis misleading, againwith due to non-independence of the sub-CTE thesampling 0ever, sothis they areis not independent. Because theples 100 sub-sample 2,111, anthe apparent ofvarious 2,111 number thenis which is not materially different 346 correlated) /the 100 | 158of the 2 veryestimates easy tois understand asVAR n →precision N. precision ofwould theisWhen 2,111 is (1then samnumber plesbe. are positively variance sampleofmean is larger 5000 5000 ( x ) V [ U ( 1 U ) / m ] than it otherwise A better estimate for the precision 2 , which is much less than the value 2,214 obtained from the sample of size me data they are positively correlated and so the set of 2 VAR( x ) V [ U (1 U ) / m] than it would otherwise be. A better estimate for the precision of er, this analysis is misleading, again due to the non-independence of the subn exercise fewetimes a practicing kn going through verificatio the 2,214 value. Thus, whilea ther from the sampling errorainvariance is some evidence actuary of smallwill samp if theycorrelated were truly independent. This intuitive result isAfter esthan are positively the variance of the sample mean is larger 1000 1000 1000herematerially 1000working 2 dependent. Because the various CTE the is then which different 346 (1asymptotic < the ) / 100 | 158 whether formulas described are forwhich his/her situation. possible to use be. theAmethods innumber ourthepaper in that the 2,111 2,214, there isis not enough data quantify it. The bias is lost in the 2,111 number then is notparticular materially different 346 is not (1 here )to / 100 | 158 [.U (1 U ) / It m] isthan it would otherwise better2,111 estimate for precision of bias 5000 5000 5000 5000 are positively the variance the sample mean is larger sitively correlated and soof the set of sampling error. This is usually theerror case. tocorrelated show 1000 that, 1000in the large sample limit, the the 2,214 value. Thus, while ther from thewhile sampling in some e is some evidence of small sample U (1isthen U) / m ] than it would mber whichestimate isthe not materially different 346 (1 otherwise ) / 100 be. | 158 A better for the precision of in the 2,214 value. Thus, ther from sampling error e is evidence of small sample uly independent. This intuitive result is limit, isthe which is not materially 2 5000 5000 bias different in that 2,111from < 2,214, there is not enough data here to quantify it. The bias is lost in the correlation of two sub-sample estimates the samo to show that, in the large sample ur paper 1000 1000 bias innot that 2,111 <sample 2,214,After there is notthrough enough data here quantify The is lost in the while therewhich pling in the 2,214value. is some evidence of small ber is error then 346 is materially different (1 Thus, ) / 100 | 158 sampling error. to This is usually case. bias Conclusions nit.the exercise a few times a practicing actuary will know going variance verificatio 5000 mates is just Unot enough n /5000 N data . AA better estimate for the just better estimate the sampling error in the 2,214a value. Thus, while there 111 < 2,214, there is here to quantify it. The for bias is lost in sampling the is usually sampling error. This the case. the 2,214 value. Thus, while there is some evidence of small sample whether the asymptotic formulas described here are working for his/her particular situation. ingThis errorisinusually or. the case. After going throughbias a variance verification exercise a few times a practicing actuary will know pling error when a3 sample size ofis1,000 is is some evidence small sample in have that 11 < 2,214, there is not enough data using here to quantify it. The bias lost 316 316in the experience, so far, we yet to see practical situation wheresituation. the theory authors’ of whether the asymptotic formulas described hereaare working for his/her particular therefore 316, but |variance 354In. the w1,000 in the sample limit, the . that, This is is usually thelarge case. not 3 n exercise a few times a practicing actuary that will practitioners know After going through a verificatio n exercise a few times a practicing actuary 1000 will know 2,111 hrough a variance verificatio therefore not 316, but < 2,214, enough data to we believe herethere failsisinnot a material way.here Hence, can use the asym outlined n /ough N .a variance Aformulas better estimate symptotic described here arefor for his/her particular 1asymptotic situation. n exercise aworking fewthe times a1 practicing actuary will know formulas described here are working for his/her particular situation. verificatio wsampling hether the quantify it. The bias is lost in the sampling error. C onclusions Nparticular situation. 5000 ymptotic formulas described here are working for his/her Conclusions 316 316 3 2 is usually the case. This erefore compare the empirical error ot 316,tobut | estimate 354 . 354 with the This result is not in the original paper. 3 experience, so far, wetohave to see a practical situation where the the theory In the authors’ sodistributed, far, we have seeyet a practical situation where theory In for the n 1000 appears to be working conclude that the1 asymptotic If authors’ we have experience, m identically butyet correlated, samples from a distribution with finite Csituation onclusions theory here fails in aon material way. Hence, we believe that practitioners canuse use the asymptotic outlined continued page 42 utlined here fails in a material way. Hence, we believe that practitioners can the asympto o a practical where the theory s’ experience, so far, we have yet to see 1 2 cexperience, variance of the CTE Estimator agrees with “experiment”, N 5000 (x x) so far,way. we have yet we to see a practical situation where thethe theory fails in a material Hence, believe that practitioners can use asymptotic mean and variance V 2 then E[¦ i ] V 2 (1 U ) where U is the correlation betwee ails in a material Hence, we believe that practitioners can use the asymptotic 2 dence of theway. sub-samples. the empirical error estimate 354 the experience, so far, we have yet This result ispractical not in the original paper. to see a situation where the theory In with the authors’ 1 m i 2 Ifthe webelieve have m identically distributed, but correlated, from a distribution with finite This result is not in3 we original paper. s not in the original paper. appears to be o utlined here way. Hence, that practitioners can use thesamples asymptotic 2 theory ymptotic for fails in a material not in the original samples. Thispaper. result is not in the originalworking paper. 3 m identically distributed, but correlated, samples from adistribution distribution with finite finite ( xi x ) 2 samples 2 2 If we have m identically distributed, but correlated, from with finite 3 ula standard errors are working. ent with the idea that the form identically distributed, but correlated, samples from a with then E[¦ ] V (1 U ) wherea Udistribution is the correlation between If we agrees have m identically distributed, but correlated, samples from a distribution with finitemean meanand andvariance variance V then 2 with “experiment”, TE Estimator m 1 2 2 2 i x( xixx) ) of where is the correlation between samples. U E[ [¦(value then Vexpected ] V95 V (1(i.e., UU) )where where is the correlation between eiance to a 95% confidence interval). ( ) x x U E ance then V ] ( 1 is the correlation between 2 ¦i mm11 mples. This result is not in the original paper. mean and variance samples. V 2 then E[¦ i ] V 2 (1 U ) where U is the correlation between 3 1 m i If we have m identically distributed, but correlated, samples from a distribution with finite 7 e 2 that there is evidence of material small sample bias 7since 7 samples. ( x x ) 2 ula standard errors are working. at the form 2 2 timates is 2,111, with an apparent precision of and variance V then E[¦ i ] V (1 U ) where U is the correlation between of 95 (i.e., a 95% confidence mean interval). m 1 i ess than the value 2,214 obtained from the sample of size samples. misleading, again due to the non-independence of the sub2 i 2 2 i dence of material small sample bias since ith an apparent precision of 7 Page 41 w Risk Management w August 2008 Chairman’s Variance ofCorner the CTE Estimator CTE Estimator … w continued from page 41 After going through a variance verification exercise a few times a practicing actuary will know whether the asymptotic formulas described here are working for his/her particular situation. If simply increasing the run size is impractical then there are other tools that can be used to improve the precision of the CTE estimator without significantly increasing the computational cost. That will be the subject of our next article. F Conclusions In the authors’ experience, so far, we have yet to see a practical situation where the theory outlined here fails in a material way. Hence, we believe that practitioners can use the asymptotic formulas presented here to understand the sampling error in a given CTE estimate. A practitioner can go through a variance verification exercise if they want to prove to themselves, or others, that the asymptotic formulas are working in their particular situation. formulas presented here to understand the sampling error in a given CTE estimate. A practitioner can go through a variance verification exercise if they want to prove to themselves, Finally, we note that if a practitioner were or others, that the asymptotic formulas are working in their particular situation. using this tool then, on the basis of the first run of 1,000 scenarios, theythat would a CTEwere using this tool then, on the basis of the first run of Finally, we note if a report practitioner estimate of 2,004+355. Is a relative sampling scenarios, they would report a CTE estimate of 2,004 r 355 . Is a relative sampling error 1,000 formulas presented here to understand the sampling error in a given CTE estimate. A errorpractitioner ofofroughly acceptroughly acceptable? The answer to that question depends on the 355 / 2 , 004 | 18 % can go through a variance verification exercise if they want to prove to themselves, or others, that the asymptotic formulas are working in their particular situation. circumstances. able? The answer to that question depends on the circumstances. firstoption run of is to increase the Finally, we is note practitionererror werethen usingwhat this tool on the basis of the arethen, the alternatives? One If 18% notthat anifaacceptable 1,000 scenarios, they would report a CTE estimate of 2,004 r 355 . Is a relative sampling error number of scenarios. If we increase the number of scenarios by a factor of K then the sampling If 18% is not an acceptable then whatThe are answer to that question depends on the of roughly 355 / 2,004 | 18error % acceptable? error scales by a factor of 1 / K . Cutting the sampling error by a factor of 4 would require us t circumstances. the alternatives? One option is to increase the run 16,000 scenarios. number of scenarios. If we increase thewhat number error then are the alternatives? One option is to increase the If 18% is not an acceptable of scenarios by a factor of K then the sampling K then number of scenarios. If we increase the number of scenarios a factor the sampling thenbythere areofother tools that can be used to If simply increasing the run size is impractical afactor factor of 1 / ofKthe . Cutting Cutting the error bysignificantly a factor of 4 would requirethe us to errorerror scales bybyathe of thesampling imscales prove precision CTE estimator without increasing computational run 16,000 scenarios. cost. That will be the subject of our next article. sampling error by a factor of 4 would require us to run 16,000increasing scenarios.the run size is impractical then there are other tools that can be used to If simply improve the precision of the CTE estimator without significantly increasing the computational cost. That will be the subject of our next article. w Page 42
© Copyright 2026 Paperzz