Joe Hirschberg and Jenny Lye Economics, University of Melbourne. September 2017
|
|
- Solomon Hood
- 6 years ago
- Views:
Transcription
1 A graphc comparson of the Feller and Delta ntervals for ratos of parameter estmates. Joe Hrschberg and Jenny Lye Economcs, Unversty of Melbourne. September 017 1
2 Introducton Feller s 1954 proposal for the use of an nverse test to construct confdence ntervals (Cs) for the rato of normally dstrbuted statstcs has been shown to be superor to the applcaton of the Delta method n several applcatons. In ths presentaton, we demonstrate how a smple graphc exposton can be generated to llustrate the relatonshp between the Delta and the Feller Cs. The advantage of the graphcal presentaton over the numerc result avalable n the feller Stata command (Coveney 004) s that t may ndcate how the level of sgnfcance can be changed to result n fnte upper and lower bounds. And to be used n cases where the Feller results are not defned.
3 We defne two normally dstrbuted random varables as: ˆ σ {[ ] } 1 σ 1 ˆ, γβ ~ N γβ,, Ω, where Ω= σ1 σ Where the quantty of nterest s the rato defned as γ θ= β ˆ Whch has been estmated by θ= ˆ γ β ˆ The objectve of ths analyss s to determne the confdence nterval for ˆθ. 3
4 The estmaton of the confdence nterval can be done va the applcaton of the Delta approxmaton method (.e. the use of the Stata nlcom command) or alternatvely they may employ the Feller method for the approxmaton of the confdence nterval. Many papers have compared these approxmatons and found that they concde for cases when the denomnator varable s βˆ estmated wth a low relatve varance (.e. tβ = σˆ s large (> 3)). However, n other cases the Feller has been shown to provde superor coverage (See Hrschberg and Lye 010) for comparsons). Note that Feller ntervals are not forced to be symmetrc as are the Delta ntervals. 4
5 The Inverse Test for the rato (Feller) 1 One method for determnng the value of the rato s to rewrte ˆ the formula ˆ γ θ= as γˆ θβ ˆˆ= 0. β ˆ Thus we could use the values of γβ ˆ, ˆ and try dfferent values of θ=θ to search for the value of θ where γ θβ= ˆ *ˆ 0 s true then conclude that ˆθ=θ. * 1 See our paper Lye and Hrschberg (01) for more examples of smlar nverse tests. 5
6 The beneft of ths approach s that the lnear functon γ θβ ˆ ˆ of normally dstrbuted random varables s normally dstrbuted. ˆ Whle the rato of normally dstrbuted random varables ˆ γ θ=, s β ˆ dstrbuted va a Cauchy dstrbuton that does not have fnte moments. However, ths does not preclude the constructon of probablty statements concernng ˆθ. 6
7 by: We can test our canddate values of θ va the hypothess defned H H 0 1 : γ θβ= ˆ ˆ 0 : γ θβ ˆ ˆ 0 Under H 0 we can construct a test statstc defned by the rato of the lnear combnaton evaluated at θ : ( ) Var ( ) T ˆ = γ θβ ˆ ˆ γ θβ ˆ ˆ Would be dstrbuted as the square of a standard normal or a Chsquare wth 1 degree of freedom. 7
8 We nvert the problem by determnng the value(s) of θ where the test statstc s equal to the value of the square of standard z α : normal for the level of α ( ) ( ˆ ) Var ( ˆ ) z α = γ θβˆ γ θβ ˆ as the lmts of the confdence nterval. We desgnate these lmts as θ 0. Usng the estmate of ˆ σˆ1 σˆ 1 Ω=, ths s equvalent to the quadratc n θ 0: σˆ1 σˆ ( ˆ ) ( ) ( ) ˆ z ˆ α θ ˆ ˆ 0 + ˆ1zα θ 0 + ˆ1 zα β σ σ γβ γ σ = 0 The two roots of ths quadratc defne the potental solutons for θ 0 and provde the lmts of the confdence nterval. 8
9 The soluton to ths quadratc are defned by: θ θ = {, } U L ( ˆ ˆ ˆ) 1z σ βγ ± α α ( βˆ σˆ ) z 0 0 z βˆ σ ˆ +γˆ σ ˆ +σˆ 1 1 σˆ σˆ 1 1 α z α βγσ ˆ ˆ ˆ Computng the soluton to ths quadratc equaton s wdely avalable (.e. the polyval command). However, the quadratc may not have real roots. Furthermore, the numbers generated my not provde much ntuton as to how the nterval may change wth dfferng values of α. z α 9
10 The Lne-plot Verson of the Feller The alternatve approach to computng the roots of the mpled quadratc equaton s the lne-plot approach. Ths approach returns to the formulaton: ( ) γ θβ ˆ ˆ 0 z α = ( ) Var γ θβ ˆ ˆ 0 After takng the square-root of both sdes we have: ˆ γ θβ ˆ ± Var γ θβ ˆ = 0 ( ˆ ) ( ˆ ) z α 0 0 Where ( γ θβ ˆ ) ± Var ( ˆ ) z α of the lnear combnaton ( ˆ ) γ θβ ˆ defnes the confdence nterval γ θβ ˆ for dfferent θs. 10
11 Thus by plottng the lne defned by γ θβ ˆ ˆ over dfferent values ˆ of θ around the value of the estmated rato ˆ γ θ= we can fnd the β ˆ confdence bounds { θ } 0U, θ0l at the values of confdence nterval of the lne defned by: cuts the zero-axs lne. ˆ ( γ θβ ˆ ) ± Var ( γ θβ ˆ ) z α ˆ θ where the 11
12 An Example of the rato of the means of a common varable. The means of dfferent sze samples of a common varable s a specal case where the covarance between the means ( σ 1 ) s set to zero and the varances are equal ( σ 1 =σ ). Here we use the property of a least squares regresson wth two dummy varables and no ntercept. Y = γ D +β D +ε t 1t t t Ths specfcaton wll estmate the means of Y for two groups when we assume the varance of Y s the same for both groups. 1
13 The estmate of the rato of the means s defned by the rato of the ˆ two estmated parameters ˆ γ θ=. β ˆ By usng a regresson we can employ the margns command, as well as the nlcom command to provde a comparson to the Delta method. The Stata code below reads the smulated blood pressure data, creates two dummes by the varable sex and computes the two means va a regresson. 13
14 Ths data has an ndcator for sex and we frst create two dummes for each sex. Then we generate a new varable _xx_ that wll only be used here. The regresson estmates the two means (note that _xx_ s not used n the regresson but we need t for the margns command also we need to drop the constant n ths case) tabulate sex, generate(sex) gen _xx_ = 0 label var _xx_ "Rato" reg bp Sex1 Sex _xx_, noconstant See program Feller_for_nd_corr_means.do for detals.. 14
15 We frst use the nlcom command to determne the Delta method 95% confdence nterval of the rato. nlcom _b[sex1] / _b[sex], level(95) Ths results n: bp Coef. Std. Err. z P> z [95% Conf. Interval] _nl_ Thus the rato s and the Delta CI s to a symmetrc nterval. 15
16 To generate the lne-plot verson of the Feller we can use the margns command to estmate the value of the lnear combnaton γ θβ ˆ ˆ where _xx_ s the value of θ along wth the 95% confdence nterval. margns, expresson (_b[sex1]-_b[sex]*_xx_) at(_xx_=(-1(.1)5)) level(95) Note that we have set the lmts of the expresson to be between - 1 and 5 as the lmtng values over whch we evaluate the possble values of the confdence nterval. 16
17 Next we use the margnsplot command to plot the lnear combnaton along wth the 95% nterval. margnsplot, ylne(0) xlne(0) recast(lne) recastc(rlne) name(rato_nd, replace) ttle(feller CI- Rato of Independent Means) We add the zero-zero axs to evaluate our hypothess. 17
18 Feller CI - Rato of Independent Means _b[sex1]-_b[sex]*_xx_ Rato Ths s the Stata produced graph. 18
19 ( γ θβ ˆ ˆ) + z α Var ( γ θβ ˆ ˆ) γ θβ ˆ ˆ ( γ θβ ˆ ˆ) z α Var ( γ θβ ˆ ˆ) θ= ˆ γˆ β ˆ = θ The tems plotted. 19
20 Delta 95% CI Feller 95% CI The dentfcaton of the 95% Feller and Delta CIs. 0
21 Alternatvely, the 99% Delta nterval s lsted as: nlcom _b[sex1] / _b[sex], level(99) wll always result n a fnte symmetrc CI: bp Coef. Std. Err. z P> z [99% Conf. Interval] _nl_ The results of the feller command however: feller bp, by(sex) level(99) results n: Quotent not statstcally sgnfcantly dfferent from zero at the 99 level. Confdence nterval not estmable 1
22 The equvalent plot however does provde more nformaton: _b[sex1]-_b[sex]*_xx_ Feller CI - Rato of Independent Means Feller lower bound Rato From ths plot we can observe a fnte lower bound wth an nfnte upper bound
23 An Example of the rato of the means from the same number of observatons. Ths example data contans data on observatons that are desgnated as before and after for the same ndvdual. We can use a seemngly unrelated regresson whch relaxes the assumptons ( σ 1 = 0, and σ 1 =σ ). However, t does requre that the number of observatons be equal. To use ths same example data as above, we frst need to reshape the data from a long format to a wde format usng the reshape command. reshape wde bp, (patent) j(when) 3
24 Then we estmate two equatons usng a seemngly unrelated regressons (SUR) specfcaton wth the sureg command: sureg (bp1 _xx_) (bp _xx_ ), corr Agan, we nclude the _xx_ varable thus we are nterested n the constants n equaton. Note that although the SUR procedure s a form of GLS the estmates n ths case are the same as OLS because the regressors are the same n both equatons, however the resultng covarance between the ntercepts wll be estmated as non-zero and the errors for each equaton wll be used to estmate separate varances. 4
25 The Delta nterval can be determned for ths rato can be estmated va the nlcom command where we dentfy the estmates of the constants from each equaton. nlcom [bp1]_b[_cons] / [bp]_b[_cons] Resultng n Coef. Std. Err. z P> z [95% Conf. Interval] _nl_ The rato s estmated as and the CI s approxmated as to
26 As before we now use the margns and margnsplot commands to estmate and plot the confdence nterval for the lnear combnaton. The code for these steps are: margns, expresson ([bp1]_b[_cons]- [bp]_b[_cons]*_xx_) at(_xx_=(-1(.1)5)) level(95) margnsplot, ylne(0) xlne(0) recast(lne) recastc(rlne) name(rato_corr, replace) ttle(feller CI for Rato of Correlated Means) 6
27 Feller CI - Rato of Correlated Means [bp1]_b[_cons]-[bp]_b[_cons]*_xx_ Delta CI Feller CI Rato The comparson of the CIs for the correlated means. 7
28 The turnng pont of a quadratc. A common case of the rato of parameter estmates would be where one has ft a quadratc specfcaton and would lke to determne the locaton of the turnng pont. For example, the typcal specfcaton would be of the form: Y =α+γ X +β X +ε t t t t The turnng pont occurs when the margnal effect of the X varable s zero: Yt * γ =γ+ β Xt, thus X = X β t 8
29 * Where X s the value of X where the sgn of the margnal effect reverses n other words where the turnng pont occurs. For example, usng the Stata example data set auto we can estmate a quadratc between a car s prce and the weght of the car here we only consder the domestc cars n the data. 3 3 See Turnng_pont.do for the full set of code for ths example. 9
30 Frst we estmate a quadratc equaton where weght s the square of weght: reg prce weght weght Source SS df MS Number of obs = 5 F(, 49) = Model Prob > F = Resdual R-squared = Adj R-squared = Total Root MSE = prce Coef. Std. Err. t P> t [95% Conf. Interval] weght weght _cons
31 The ftted regresson lne and data are shown below: 0 5,000 10,000 15,000,000 3,000 4,000 5,000 Weght (lbs.) Prce to weght relatonshp 31
32 To fnd the turnng pont we can use the followng code: nlcom (_b[weght] /( -*_b[weght]) ), level(95) Where the nlcom command provdes the Delta estmaton of the CIs for the turnng pont as: prce Coef. Std. Err. z P> z [95% Conf. Interval] _nl_ The turnng pont s estmated at, pounds and the approxmate 95% CI s, to,
33 Alternatvely, we can plot the margnal effect of the change n weght as t nfluences the prce. Usng the Stata code: margns, expresson (_b[weght] + *_b[weght]*weght) at (weght = (1500(00)3000)) level(95) margnsplot, ylne(0) recast(lne) recastc(rlne) name(turnng_pont, replace) ttle(feller CI - Turnng pont of quadratc ) 33
34 Feller CI - Turnng pont of quadratc _b[weght] + *_b[weght]*weght Feller CI Delta CI 1,500 1,700 1,900,100,300,500,700,900 Weght (lbs.) Turnng Pont Note that the Feller CI s asymmetrc and has a lower lmt that s less than the Delta nterval. (see Lye and Hrschberg 01 for more examples). 34
35 An Example of the 50% dose and Wllngness to Pay. The 50% dose problem was the orgnal (1944) applcaton that Feller consdered when he frst proposed ths method. In that case, he was nterested n the dose of a medcaton/poson that would result n a greater than 50% chance that the subject would lve/de. In the wllngness to pay case the varable appled s the prce pad (e.g. Hole 007). How much would someone pay for a product/servce? Thus, the data would be a dchotomous response to a contnuous stmulus-drug-prce. 35
36 Typcally, ths problem s modelled by a probablty of occurrence usng a logt or probt estmator. PY ( = 1) =α+γ X + β W +ε t t t t = 1 Where X s the varable dose/prce and the Ws are other covarates. The queston of the 50% dose s posed as: What value of X (desgnated as X * ) makes ths relatonshp k true? ½ k * =α+γ X + β = 1 W t 36
37 Solvng for X * we have the rato gven by: X * = ½ k α+ βw = 1 γ t The rato can be evaluated at specfc values of W for cases of nterest or at the means. 37
38 As an example, we use data on the survvors of the Ttanc dsaster were the dose/prce s the fare pad for the trp (see Frey et al 010). We ft the model to the data wth covarates for sex and age. 4 probt survve fare age female Whch results n: Probt regresson Number of obs = 1,033 LR ch(3) = Prob > ch = Log lkelhood = Pseudo R = survved Coef. Std. Err. z P> z [95% Conf. Interval] fare age female _cons See the program Wllngness_to_pay.do for the full code for ths example. 38
39 Usng the Delta method for the 50% fare for the average age and assumng 50% gender splt: nlcom (-(_b[_cons]+ _b[age]*m_age + _b[female] *.5) / _b[fare] ), level(95) Results n the estmated CI on the 50% fare of: survved Coef. Std. Err. z P> z [95% Conf. Interval] _nl_ Thus, we fnd a fare of 48.7 pounds wth an approxmate 95% nterval of (3.8 to 64.54). 39
40 To estmate the Feller nterval we use the margns command to predct the probablty of survval: margns, expresson ( normal( _b[_cons]+_b[fare]*fare + _b[age]*age + _b[female] * female ) ) at(age (mean) female =.5 fare=(0(5)100)) level(95) Followed by the margnsplot command: margnsplot, ylne(.5) recast(lne) recastc(rlne) 40
41 Feller CI - lmt of 50% dose Probablty of survval % Fare Feller CI Delta CI Fare pad n pounds In ths case the Delta and Feller ntervals concde qute closely. 41
42 Conclusons In addton to the examples dscussed here we have also added the routne used our 017 paper that estmates the ndrect leastsquares estmate for the just dentfed IV estmator whch can be shown to be the rato of parameters estmated from two seemngly unrelated regressons. 5 The types of code for these cases vares dependng on the best way to proceed. Usng the margns and the margnsplot command s a very compact method for generatng these plots. The lmts of the margns may need some adjustment after the frst attempt so that the bounds are plotted wthn the plot wndow. In addton, the level opton allows the constructon of CIs wth dfferent values of α. 5 See the program Indrect_Least_squares.do for ths example. 4
43 All programs/data used n ths talk are avalable at: do fles Feller_for_nd_corr.do Indrect_Least_Squares.do Turnng_pont.do Wllngness_to_pay.do Data bplong_example.dta maketable4.dta auto_domestc.dta ttanc1.dta 43
44 References Coveney, Joseph, (004), FIELLER: Stata module to calculate confdence ntervals of quotents of normal varate, downloaded 5/9/17, from : Feller, E. C., (1944), A Fundamental Formula n the Statstcs of Bologcal Assay and Some Applcatons, Quarterly Journal of Pharmacy and Pharmacology, 17, ,(1954), Some Problems n Interval Estmaton, Journal of the Royal Statstcal Socety, Ser. B, 16, Frey, B. S., D. A. Savage, and B. Torgler, 010, Interacton of natural survval nstncts and nternalzed socal norms explorng the Ttanc and Lustana dsasters, Proceedngs of the Natonal Academy of Scence, 107, Hole, A. R., 007, A comparson of approaches to estmatng confdence ntervals for wllngness to pay measures, Health Economcs, 16, Hrschberg, Joseph G. and Jenny N. Lye, 010, A Geometrc Comparson of the Delta and Feller Confdence Intervals, The Amercan Statstcan, 64, , 017, Invertng the Indrect - The Ellpse and the Boomerang: Vsualsng the Confdence Intervals of the Structural Coeffcent from Two-Stage Least Squares, Journal of Econometrcs, 199, Lye, Jenny N. and Joseph G. Hrschberg, 010, Two Geometrc Representatons of Confdence Intervals for Ratos of Lnear Combnatons of Regresson Parameters: An Applcaton to the NAIRU, Economcs Letters, 108, , 01, Inverse Test Confdence Intervals for Turnng Ponts: A Demonstraton wth Hgher Order Polynomals, Dek Terrell and Danel Mllmet (ed.) 30 th Annversary Edton (Advances n Econometrcs, Volume 30) Emerald Group Publshng Lmted, pp
MgtOp 215 Chapter 13 Dr. Ahn
MgtOp 5 Chapter 3 Dr Ahn Consder two random varables X and Y wth,,, In order to study the relatonshp between the two random varables, we need a numercal measure that descrbes the relatonshp The covarance
More informationNotes are not permitted in this examination. Do not turn over until you are told to do so by the Invigilator.
UNIVERSITY OF EAST ANGLIA School of Economcs Man Seres PG Examnaton 2016-17 BANKING ECONOMETRICS ECO-7014A Tme allowed: 2 HOURS Answer ALL FOUR questons. Queston 1 carres a weght of 30%; queston 2 carres
More informationModule Contact: Dr P Moffatt, ECO Copyright of the University of East Anglia Version 2
UNIVERSITY OF EAST ANGLIA School of Economcs Man Seres PG Examnaton 2012-13 FINANCIAL ECONOMETRICS ECO-M017 Tme allowed: 2 hours Answer ALL FOUR questons. Queston 1 carres a weght of 25%; Queston 2 carres
More informationWhich of the following provides the most reasonable approximation to the least squares regression line? (a) y=50+10x (b) Y=50+x (d) Y=1+50x
Whch of the followng provdes the most reasonable approxmaton to the least squares regresson lne? (a) y=50+10x (b) Y=50+x (c) Y=10+50x (d) Y=1+50x (e) Y=10+x In smple lnear regresson the model that s begn
More information3/3/2014. CDS M Phil Econometrics. Vijayamohanan Pillai N. Truncated standard normal distribution for a = 0.5, 0, and 0.5. CDS Mphil Econometrics
Lmted Dependent Varable Models: Tobt an Plla N 1 CDS Mphl Econometrcs Introducton Lmted Dependent Varable Models: Truncaton and Censorng Maddala, G. 1983. Lmted Dependent and Qualtatve Varables n Econometrcs.
More informationECONOMETRICS - FINAL EXAM, 3rd YEAR (GECO & GADE)
ECONOMETRICS - FINAL EXAM, 3rd YEAR (GECO & GADE) May 17, 2016 15:30 Frst famly name: Name: DNI/ID: Moble: Second famly Name: GECO/GADE: Instructor: E-mal: Queston 1 A B C Blank Queston 2 A B C Blank Queston
More informationTests for Two Correlations
PASS Sample Sze Software Chapter 805 Tests for Two Correlatons Introducton The correlaton coeffcent (or correlaton), ρ, s a popular parameter for descrbng the strength of the assocaton between two varables.
More information/ Computational Genomics. Normalization
0-80 /02-70 Computatonal Genomcs Normalzaton Gene Expresson Analyss Model Computatonal nformaton fuson Bologcal regulatory networks Pattern Recognton Data Analyss clusterng, classfcaton normalzaton, mss.
More informationCapability Analysis. Chapter 255. Introduction. Capability Analysis
Chapter 55 Introducton Ths procedure summarzes the performance of a process based on user-specfed specfcaton lmts. The observed performance as well as the performance relatve to the Normal dstrbuton are
More informationII. Random Variables. Variable Types. Variables Map Outcomes to Numbers
II. Random Varables Random varables operate n much the same way as the outcomes or events n some arbtrary sample space the dstncton s that random varables are smply outcomes that are represented numercally.
More informationTests for Two Ordered Categorical Variables
Chapter 253 Tests for Two Ordered Categorcal Varables Introducton Ths module computes power and sample sze for tests of ordered categorcal data such as Lkert scale data. Assumng proportonal odds, such
More informationGraphical Methods for Survival Distribution Fitting
Graphcal Methods for Survval Dstrbuton Fttng In ths Chapter we dscuss the followng two graphcal methods for survval dstrbuton fttng: 1. Probablty Plot, 2. Cox-Snell Resdual Method. Probablty Plot: The
More informationCHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS
CHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS QUESTIONS 9.1. (a) In a log-log model the dependent and all explanatory varables are n the logarthmc form. (b) In the log-ln model the dependent varable
More informationAnalysis of Variance and Design of Experiments-II
Analyss of Varance and Desgn of Experments-II MODULE VI LECTURE - 4 SPLIT-PLOT AND STRIP-PLOT DESIGNS Dr. Shalabh Department of Mathematcs & Statstcs Indan Insttute of Technology Kanpur An example to motvate
More informationRandom Variables. b 2.
Random Varables Generally the object of an nvestgators nterest s not necessarly the acton n the sample space but rather some functon of t. Techncally a real valued functon or mappng whose doman s the sample
More informationChapter 3 Descriptive Statistics: Numerical Measures Part B
Sldes Prepared by JOHN S. LOUCKS St. Edward s Unversty Slde 1 Chapter 3 Descrptve Statstcs: Numercal Measures Part B Measures of Dstrbuton Shape, Relatve Locaton, and Detectng Outlers Eploratory Data Analyss
More informationSolutions to Odd-Numbered End-of-Chapter Exercises: Chapter 12
Introducton to Econometrcs (3 rd Updated Edton) by James H. Stock and Mark W. Watson Solutons to Odd-Numbered End-of-Chapter Exercses: Chapter 1 (Ths verson July 0, 014) Stock/Watson - Introducton to Econometrcs
More informationSpatial Variations in Covariates on Marriage and Marital Fertility: Geographically Weighted Regression Analyses in Japan
Spatal Varatons n Covarates on Marrage and Martal Fertlty: Geographcally Weghted Regresson Analyses n Japan Kenj Kamata (Natonal Insttute of Populaton and Socal Securty Research) Abstract (134) To understand
More informationPASS Sample Size Software. :log
PASS Sample Sze Software Chapter 70 Probt Analyss Introducton Probt and lot analyss may be used for comparatve LD 50 studes for testn the effcacy of drus desned to prevent lethalty. Ths proram module presents
More informationTesting for Omitted Variables
Testng for Omtted Varables Jeroen Weese Department of Socology Unversty of Utrecht The Netherlands emal J.weese@fss.uu.nl tel +31 30 2531922 fax+31 30 2534405 Prepared for North Amercan Stata users meetng
More informationChapter 5 Student Lecture Notes 5-1
Chapter 5 Student Lecture Notes 5-1 Basc Busness Statstcs (9 th Edton) Chapter 5 Some Important Dscrete Probablty Dstrbutons 004 Prentce-Hall, Inc. Chap 5-1 Chapter Topcs The Probablty Dstrbuton of a Dscrete
More informationData Mining Linear and Logistic Regression
07/02/207 Data Mnng Lnear and Logstc Regresson Mchael L of 26 Regresson In statstcal modellng, regresson analyss s a statstcal process for estmatng the relatonshps among varables. Regresson models are
More informationoccurrence of a larger storm than our culvert or bridge is barely capable of handling? (what is The main question is: What is the possibility of
Module 8: Probablty and Statstcal Methods n Water Resources Engneerng Bob Ptt Unversty of Alabama Tuscaloosa, AL Flow data are avalable from numerous USGS operated flow recordng statons. Data s usually
More information4. Greek Letters, Value-at-Risk
4 Greek Letters, Value-at-Rsk 4 Value-at-Rsk (Hull s, Chapter 8) Math443 W08, HM Zhu Outlne (Hull, Chap 8) What s Value at Rsk (VaR)? Hstorcal smulatons Monte Carlo smulatons Model based approach Varance-covarance
More informationISyE 512 Chapter 9. CUSUM and EWMA Control Charts. Instructor: Prof. Kaibo Liu. Department of Industrial and Systems Engineering UW-Madison
ISyE 512 hapter 9 USUM and EWMA ontrol harts Instructor: Prof. Kabo Lu Department of Industral and Systems Engneerng UW-Madson Emal: klu8@wsc.edu Offce: Room 317 (Mechancal Engneerng Buldng) ISyE 512 Instructor:
More informationThe Mack-Method and Analysis of Variability. Erasmus Gerigk
The Mac-Method and Analyss of Varablty Erasmus Gerg ontents/outlne Introducton Revew of two reservng recpes: Incremental Loss-Rato Method han-ladder Method Mac s model assumptons and estmatng varablty
More informationMeasures of Spread IQR and Deviation. For exam X, calculate the mean, median and mode. For exam Y, calculate the mean, median and mode.
Part 4 Measures of Spread IQR and Devaton In Part we learned how the three measures of center offer dfferent ways of provdng us wth a sngle representatve value for a data set. However, consder the followng
More informationRisk and Return: The Security Markets Line
FIN 614 Rsk and Return 3: Markets Professor Robert B.H. Hauswald Kogod School of Busness, AU 1/25/2011 Rsk and Return: Markets Robert B.H. Hauswald 1 Rsk and Return: The Securty Markets Lne From securtes
More informationAn Application of Alternative Weighting Matrix Collapsing Approaches for Improving Sample Estimates
Secton on Survey Research Methods An Applcaton of Alternatve Weghtng Matrx Collapsng Approaches for Improvng Sample Estmates Lnda Tompkns 1, Jay J. Km 2 1 Centers for Dsease Control and Preventon, atonal
More informationEvaluating Performance
5 Chapter Evaluatng Performance In Ths Chapter Dollar-Weghted Rate of Return Tme-Weghted Rate of Return Income Rate of Return Prncpal Rate of Return Daly Returns MPT Statstcs 5- Measurng Rates of Return
More informationInterval Estimation for a Linear Function of. Variances of Nonnormal Distributions. that Utilize the Kurtosis
Appled Mathematcal Scences, Vol. 7, 013, no. 99, 4909-4918 HIKARI Ltd, www.m-hkar.com http://dx.do.org/10.1988/ams.013.37366 Interval Estmaton for a Lnear Functon of Varances of Nonnormal Dstrbutons that
More informationLinear Combinations of Random Variables and Sampling (100 points)
Economcs 30330: Statstcs for Economcs Problem Set 6 Unversty of Notre Dame Instructor: Julo Garín Sprng 2012 Lnear Combnatons of Random Varables and Samplng 100 ponts 1. Four-part problem. Go get some
More informationSimple Regression Theory II 2010 Samuel L. Baker
SIMPLE REGRESSIO THEORY II Smple Regresson Theory II 00 Samuel L. Baker Assessng how good the regresson equaton s lkely to be Assgnment A gets nto drawng nferences about how close the regresson lne mght
More informationEXAMINATIONS OF THE HONG KONG STATISTICAL SOCIETY
EXAMINATIONS OF THE HONG KONG STATISTICAL SOCIETY HIGHER CERTIFICATE IN STATISTICS, 2013 MODULE 7 : Tme seres and ndex numbers Tme allowed: One and a half hours Canddates should answer THREE questons.
More informationFinal Exam. 7. (10 points) Please state whether each of the following statements is true or false. No explanation needed.
Fnal Exam Fall 4 Econ 8-67 Closed Book. Formula Sheet Provded. Calculators OK. Tme Allowed: hours Please wrte your answers on the page below each queston. (5 ponts) Assume that the rsk-free nterest rate
More informationSampling Distributions of OLS Estimators of β 0 and β 1. Monte Carlo Simulations
Addendum to NOTE 4 Samplng Dstrbutons of OLS Estmators of β and β Monte Carlo Smulatons The True Model: s gven by the populaton regresson equaton (PRE) Y = β + β X + u = 7. +.9X + u () where β = 7. and
More informationσ may be counterbalanced by a larger
Questons CHAPTER 5: TWO-VARIABLE REGRESSION: INTERVAL ESTIMATION AND HYPOTHESIS TESTING 5.1 (a) True. The t test s based on varables wth a normal dstrbuton. Snce the estmators of β 1 and β are lnear combnatons
More informationAppendix - Normally Distributed Admissible Choices are Optimal
Appendx - Normally Dstrbuted Admssble Choces are Optmal James N. Bodurtha, Jr. McDonough School of Busness Georgetown Unversty and Q Shen Stafford Partners Aprl 994 latest revson September 00 Abstract
More informationA Bootstrap Confidence Limit for Process Capability Indices
A ootstrap Confdence Lmt for Process Capablty Indces YANG Janfeng School of usness, Zhengzhou Unversty, P.R.Chna, 450001 Abstract The process capablty ndces are wdely used by qualty professonals as an
More informationEconomic Design of Short-Run CSP-1 Plan Under Linear Inspection Cost
Tamkang Journal of Scence and Engneerng, Vol. 9, No 1, pp. 19 23 (2006) 19 Economc Desgn of Short-Run CSP-1 Plan Under Lnear Inspecton Cost Chung-Ho Chen 1 * and Chao-Yu Chou 2 1 Department of Industral
More informationChapter 3 Student Lecture Notes 3-1
Chapter 3 Student Lecture otes 3-1 Busness Statstcs: A Decson-Makng Approach 6 th Edton Chapter 3 Descrbng Data Usng umercal Measures 005 Prentce-Hall, Inc. Chap 3-1 Chapter Goals After completng ths chapter,
More informationNumerical Analysis ECIV 3306 Chapter 6
The Islamc Unversty o Gaza Faculty o Engneerng Cvl Engneerng Department Numercal Analyss ECIV 3306 Chapter 6 Open Methods & System o Non-lnear Eqs Assocate Pro. Mazen Abualtaye Cvl Engneerng Department,
More informationCalibration Methods: Regression & Correlation. Calibration Methods: Regression & Correlation
Calbraton Methods: Regresson & Correlaton Calbraton A seres of standards run (n replcate fashon) over a gven concentraton range. Standards Comprsed of analte(s) of nterest n a gven matr composton. Matr
More informationMode is the value which occurs most frequency. The mode may not exist, and even if it does, it may not be unique.
1.7.4 Mode Mode s the value whch occurs most frequency. The mode may not exst, and even f t does, t may not be unque. For ungrouped data, we smply count the largest frequency of the gven value. If all
More information3: Central Limit Theorem, Systematic Errors
3: Central Lmt Theorem, Systematc Errors 1 Errors 1.1 Central Lmt Theorem Ths theorem s of prme mportance when measurng physcal quanttes because usually the mperfectons n the measurements are due to several
More informationEDC Introduction
.0 Introducton EDC3 In the last set of notes (EDC), we saw how to use penalty factors n solvng the EDC problem wth losses. In ths set of notes, we want to address two closely related ssues. What are, exactly,
More informationPrinciples of Finance
Prncples of Fnance Grzegorz Trojanowsk Lecture 6: Captal Asset Prcng Model Prncples of Fnance - Lecture 6 1 Lecture 6 materal Requred readng: Elton et al., Chapters 13, 14, and 15 Supplementary readng:
More informationInternational ejournals
Avalable onlne at www.nternatonalejournals.com ISSN 0976 1411 Internatonal ejournals Internatonal ejournal of Mathematcs and Engneerng 7 (010) 86-95 MODELING AND PREDICTING URBAN MALE POPULATION OF BANGLADESH:
More informationA Note on Robust Estimation of Repeat Sales Indexes with Serial Correlation in Asset Returns
A Note on Robust Estmaton of Repeat Sales Indexes wth Seral Correlaton n Asset Returns Kathryn Graddy Department of Economcs and Internatonal Busness School Brandes Unversty (kgraddy@brandes.edu) Jonathan
More informationElton, Gruber, Brown and Goetzmann. Modern Portfolio Theory and Investment Analysis, 7th Edition. Solutions to Text Problems: Chapter 4
Elton, Gruber, Brown and Goetzmann Modern ortfolo Theory and Investment Analyss, 7th Edton Solutons to Text roblems: Chapter 4 Chapter 4: roblem 1 A. Expected return s the sum of each outcome tmes ts assocated
More informationA Comparison of Statistical Methods in Interrupted Time Series Analysis to Estimate an Intervention Effect
Transport and Road Safety (TARS) Research Joanna Wang A Comparson of Statstcal Methods n Interrupted Tme Seres Analyss to Estmate an Interventon Effect Research Fellow at Transport & Road Safety (TARS)
More information1 Omitted Variable Bias: Part I. 2 Omitted Variable Bias: Part II. The Baseline: SLR.1-4 hold, and our estimates are unbiased
Introductory Appled Econometrcs EEP/IAS 118 Sprng 2014 Andrew Crane-Droesch Secton #5 Feb 26 2014 1 Omtted Varable Bas: Part I Remember that a key assumpton needed to get an unbased estmate of β 1 n the
More informationCHAPTER 3: BAYESIAN DECISION THEORY
CHATER 3: BAYESIAN DECISION THEORY Decson makng under uncertanty 3 rogrammng computers to make nference from data requres nterdscplnary knowledge from statstcs and computer scence Knowledge of statstcs
More informationUNIVERSITY OF VICTORIA Midterm June 6, 2018 Solutions
UIVERSITY OF VICTORIA Mdterm June 6, 08 Solutons Econ 45 Summer A0 08 age AME: STUDET UMBER: V00 Course ame & o. Descrptve Statstcs and robablty Economcs 45 Secton(s) A0 CR: 3067 Instructor: Betty Johnson
More informationUsing Conditional Heteroskedastic
ITRON S FORECASTING BROWN BAG SEMINAR Usng Condtonal Heteroskedastc Varance Models n Load Research Sample Desgn Dr. J. Stuart McMenamn March 6, 2012 Please Remember» Phones are Muted: In order to help
More informationChapter 10 Making Choices: The Method, MARR, and Multiple Attributes
Chapter 0 Makng Choces: The Method, MARR, and Multple Attrbutes INEN 303 Sergy Butenko Industral & Systems Engneerng Texas A&M Unversty Comparng Mutually Exclusve Alternatves by Dfferent Evaluaton Methods
More informationElements of Economic Analysis II Lecture VI: Industry Supply
Elements of Economc Analyss II Lecture VI: Industry Supply Ka Hao Yang 10/12/2017 In the prevous lecture, we analyzed the frm s supply decson usng a set of smple graphcal analyses. In fact, the dscusson
More informationMidterm Exam. Use the end of month price data for the S&P 500 index in the table below to answer the following questions.
Unversty of Washngton Summer 2001 Department of Economcs Erc Zvot Economcs 483 Mdterm Exam Ths s a closed book and closed note exam. However, you are allowed one page of handwrtten notes. Answer all questons
More informationSupplementary material for Non-conjugate Variational Message Passing for Multinomial and Binary Regression
Supplementary materal for Non-conjugate Varatonal Message Passng for Multnomal and Bnary Regresson October 9, 011 1 Alternatve dervaton We wll focus on a partcular factor f a and varable x, wth the am
More informationAlternatives to Shewhart Charts
Alternatves to Shewhart Charts CUSUM & EWMA S Wongsa Overvew Revstng Shewhart Control Charts Cumulatve Sum (CUSUM) Control Chart Eponentally Weghted Movng Average (EWMA) Control Chart 2 Revstng Shewhart
More informationTeaching Note on Factor Model with a View --- A tutorial. This version: May 15, Prepared by Zhi Da *
Copyrght by Zh Da and Rav Jagannathan Teachng Note on For Model th a Ve --- A tutoral Ths verson: May 5, 2005 Prepared by Zh Da * Ths tutoral demonstrates ho to ncorporate economc ves n optmal asset allocaton
More informationBasket options and implied correlations: a closed form approach
Basket optons and mpled correlatons: a closed form approach Svetlana Borovkova Free Unversty of Amsterdam CFC conference, London, January 7-8, 007 Basket opton: opton whose underlyng s a basket (.e. a
More informationUnderstanding price volatility in electricity markets
Proceedngs of the 33rd Hawa Internatonal Conference on System Scences - 2 Understandng prce volatlty n electrcty markets Fernando L. Alvarado, The Unversty of Wsconsn Rajesh Rajaraman, Chrstensen Assocates
More informationIntroduction. Why One-Pass Statistics?
BERKELE RESEARCH GROUP Ths manuscrpt s program documentaton for three ways to calculate the mean, varance, skewness, kurtoss, covarance, correlaton, regresson parameters and other regresson statstcs. Although
More informationMultifactor Term Structure Models
1 Multfactor Term Structure Models A. Lmtatons of One-Factor Models 1. Returns on bonds of all maturtes are perfectly correlated. 2. Term structure (and prces of every other dervatves) are unquely determned
More informationPrice Formation on Agricultural Land Markets A Microstructure Analysis
Prce Formaton on Agrcultural Land Markets A Mcrostructure Analyss Martn Odenng & Slke Hüttel Department of Agrcultural Economcs, Humboldt-Unverstät zu Berln Department of Agrcultural Economcs, Unversty
More informationAppendix for Solving Asset Pricing Models when the Price-Dividend Function is Analytic
Appendx for Solvng Asset Prcng Models when the Prce-Dvdend Functon s Analytc Ovdu L. Caln Yu Chen Thomas F. Cosmano and Alex A. Hmonas January 3, 5 Ths appendx provdes proofs of some results stated n our
More informationProblem Set 6 Finance 1,
Carnege Mellon Unversty Graduate School of Industral Admnstraton Chrs Telmer Wnter 2006 Problem Set 6 Fnance, 47-720. (representatve agent constructon) Consder the followng two-perod, two-agent economy.
More informationOCR Statistics 1 Working with data. Section 2: Measures of location
OCR Statstcs 1 Workng wth data Secton 2: Measures of locaton Notes and Examples These notes have sub-sectons on: The medan Estmatng the medan from grouped data The mean Estmatng the mean from grouped data
More informationNotes on experimental uncertainties and their propagation
Ed Eyler 003 otes on epermental uncertantes and ther propagaton These notes are not ntended as a complete set of lecture notes, but nstead as an enumeraton of some of the key statstcal deas needed to obtan
More informationEconomics 1410 Fall Section 7 Notes 1. Define the tax in a flexible way using T (z), where z is the income reported by the agent.
Economcs 1410 Fall 2017 Harvard Unversty Yaan Al-Karableh Secton 7 Notes 1 I. The ncome taxaton problem Defne the tax n a flexble way usng T (), where s the ncome reported by the agent. Retenton functon:
More informationA Utilitarian Approach of the Rawls s Difference Principle
1 A Utltaran Approach of the Rawls s Dfference Prncple Hyeok Yong Kwon a,1, Hang Keun Ryu b,2 a Department of Poltcal Scence, Korea Unversty, Seoul, Korea, 136-701 b Department of Economcs, Chung Ang Unversty,
More informationPrice and Quantity Competition Revisited. Abstract
rce and uantty Competton Revsted X. Henry Wang Unversty of Mssour - Columba Abstract By enlargng the parameter space orgnally consdered by Sngh and Vves (984 to allow for a wder range of cost asymmetry,
More informationApplications of Myerson s Lemma
Applcatons of Myerson s Lemma Professor Greenwald 28-2-7 We apply Myerson s lemma to solve the sngle-good aucton, and the generalzaton n whch there are k dentcal copes of the good. Our objectve s welfare
More informationIntroduction. Chapter 7 - An Introduction to Portfolio Management
Introducton In the next three chapters, we wll examne dfferent aspects of captal market theory, ncludng: Brngng rsk and return nto the pcture of nvestment management Markowtz optmzaton Modelng rsk and
More information- contrast so-called first-best outcome of Lindahl equilibrium with case of private provision through voluntary contributions of households
Prvate Provson - contrast so-called frst-best outcome of Lndahl equlbrum wth case of prvate provson through voluntary contrbutons of households - need to make an assumpton about how each household expects
More informationFinance 402: Problem Set 1 Solutions
Fnance 402: Problem Set 1 Solutons Note: Where approprate, the fnal answer for each problem s gven n bold talcs for those not nterested n the dscusson of the soluton. 1. The annual coupon rate s 6%. A
More informationMonetary Tightening Cycles and the Predictability of Economic Activity. by Tobias Adrian and Arturo Estrella * October 2006.
Monetary Tghtenng Cycles and the Predctablty of Economc Actvty by Tobas Adran and Arturo Estrella * October 2006 Abstract Ten out of thrteen monetary tghtenng cycles snce 1955 were followed by ncreases
More informationThe Integration of the Israel Labour Force Survey with the National Insurance File
The Integraton of the Israel Labour Force Survey wth the Natonal Insurance Fle Natale SHLOMO Central Bureau of Statstcs Kanfey Nesharm St. 66, corner of Bach Street, Jerusalem Natales@cbs.gov.l Abstact:
More informationLikelihood Fits. Craig Blocker Brandeis August 23, 2004
Lkelhood Fts Crag Blocker Brandes August 23, 2004 Outlne I. What s the queston? II. Lkelhood Bascs III. Mathematcal Propertes IV. Uncertantes on Parameters V. Mscellaneous VI. Goodness of Ft VII. Comparson
More informationASSESSING GOODNESS OF FIT OF GENERALIZED LINEAR MODELS TO SPARSE DATA USING HIGHER ORDER MOMENT CORRECTIONS
ASSESSING GOODNESS OF FIT OF GENERALIZED LINEAR MODELS TO SPARSE DATA USING HIGHER ORDER MOMENT CORRECTIONS S. R. PAUL Department of Mathematcs & Statstcs, Unversty of Wndsor, Wndsor, ON N9B 3P4, Canada
More informationA Set of new Stochastic Trend Models
A Set of new Stochastc Trend Models Johannes Schupp Longevty 13, Tape, 21 th -22 th September 2017 www.fa-ulm.de Introducton Uncertanty about the evoluton of mortalty Measure longevty rsk n penson or annuty
More informationUsing Cumulative Count of Conforming CCC-Chart to Study the Expansion of the Cement
IOSR Journal of Engneerng (IOSRJEN) e-issn: 225-32, p-issn: 2278-879, www.osrjen.org Volume 2, Issue (October 22), PP 5-6 Usng Cumulatve Count of Conformng CCC-Chart to Study the Expanson of the Cement
More informationSequential equilibria of asymmetric ascending auctions: the case of log-normal distributions 3
Sequental equlbra of asymmetrc ascendng auctons: the case of log-normal dstrbutons 3 Robert Wlson Busness School, Stanford Unversty, Stanford, CA 94305-505, USA Receved: ; revsed verson. Summary: The sequental
More informationCracking VAR with kernels
CUTTIG EDGE. PORTFOLIO RISK AALYSIS Crackng VAR wth kernels Value-at-rsk analyss has become a key measure of portfolo rsk n recent years, but how can we calculate the contrbuton of some portfolo component?
More informationComparison of Singular Spectrum Analysis and ARIMA
Int. Statstcal Inst.: Proc. 58th World Statstcal Congress, 0, Dubln (Sesson CPS009) p.99 Comparson of Sngular Spectrum Analss and ARIMA Models Zokae, Mohammad Shahd Behesht Unverst, Department of Statstcs
More informationIntroduction to PGMs: Discrete Variables. Sargur Srihari
Introducton to : Dscrete Varables Sargur srhar@cedar.buffalo.edu Topcs. What are graphcal models (or ) 2. Use of Engneerng and AI 3. Drectonalty n graphs 4. Bayesan Networks 5. Generatve Models and Samplng
More informationAsian basket options. in oil markets
Asan basket optons and mpled correlatons n ol markets Svetlana Borovkova Vre Unverstet Amsterdam, he etherlands Jont work wth Ferry Permana (Bandung) Basket opton: opton whose underlyng s a basket (e a
More informationCOS 511: Theoretical Machine Learning. Lecturer: Rob Schapire Lecture #21 Scribe: Lawrence Diao April 23, 2013
COS 511: Theoretcal Machne Learnng Lecturer: Rob Schapre Lecture #21 Scrbe: Lawrence Dao Aprl 23, 2013 1 On-Lne Log Loss To recap the end of the last lecture, we have the followng on-lne problem wth N
More informationEstimation of Wage Equations in Australia: Allowing for Censored Observations of Labour Supply *
Estmaton of Wage Equatons n Australa: Allowng for Censored Observatons of Labour Supply * Guyonne Kalb and Rosanna Scutella* Melbourne Insttute of Appled Economc and Socal Research The Unversty of Melbourne
More informationRisk Reduction and Real Estate Portfolio Size
Rsk Reducton and Real Estate Portfolo Sze Stephen L. Lee and Peter J. Byrne Department of Land Management and Development, The Unversty of Readng, Whteknghts, Readng, RG6 6AW, UK. A Paper Presented at
More informationHarmonised Labour Cost Index. Methodology
Harmonsed Labour Cost Index Methodology March 2013 Index 1 Introducton 3 2 Scope, coverage and reference perod 4 3 Defntons 5 4 Sources of nformaton 7 5 Formulae employed 9 6 Results obtaned 10 7 Seres
More informationFinal Exam - section 1. Thursday, December hours, 30 minutes
Econometrics, ECON312 San Francisco State University Michael Bar Fall 2013 Final Exam - section 1 Thursday, December 19 1 hours, 30 minutes Name: Instructions 1. This is closed book, closed notes exam.
More informationFinancial mathematics
Fnancal mathematcs Jean-Luc Bouchot jean-luc.bouchot@drexel.edu February 19, 2013 Warnng Ths s a work n progress. I can not ensure t to be mstake free at the moment. It s also lackng some nformaton. But
More informationMutual Funds and Management Styles. Active Portfolio Management
utual Funds and anagement Styles ctve Portfolo anagement ctve Portfolo anagement What s actve portfolo management? How can we measure the contrbuton of actve portfolo management? We start out wth the CP
More informationANOVA Procedures for Multiple Linear Regression Model with Non-normal Error Distribution: A Quantile Function Distribution Approach
Amercan Journal of Mathematcs and Statstcs 07, 7(4): 69-78 DOI: 0.593/j.ajms.070704.05 ANOVA Procedures for Multple Lnear Regresson Model wth Non-normal Error Dstrbuton: A Quantle Functon Dstrbuton Approach
More informationEU Banks Rating Assignments: Is There Heterogeneity Between New And Old Member Countries?
Department of Economcs and Fnance Workng Paper No. 09-39 Economcs and Fnance Workng Paper Seres Guglelmo Mara Caporale, Roman Matousek and Chrs Stewart EU Banks Ratng Assgnments: Is There Heterogenety
More informationMaximum Likelihood Estimation of Isotonic Normal Means with Unknown Variances*
Journal of Multvarate Analyss 64, 183195 (1998) Artcle No. MV971717 Maxmum Lelhood Estmaton of Isotonc Normal Means wth Unnown Varances* Nng-Zhong Sh and Hua Jang Northeast Normal Unversty, Changchun,Chna
More informationStochastic ALM models - General Methodology
Stochastc ALM models - General Methodology Stochastc ALM models are generally mplemented wthn separate modules: A stochastc scenaros generator (ESG) A cash-flow projecton tool (or ALM projecton) For projectng
More informationBayesian belief networks
CS 2750 achne Learnng Lecture 12 ayesan belef networks los Hauskrecht mlos@cs.ptt.edu 5329 Sennott Square CS 2750 achne Learnng Densty estmaton Data: D { D1 D2.. Dn} D x a vector of attrbute values ttrbutes:
More information