BAYESIAN PREMIUM AND ASYMPTOTIC APPROXIMATION FOR A CLASS OF LOSS FUNCTIONS

Size: px
Start display at page:

Download "BAYESIAN PREMIUM AND ASYMPTOTIC APPROXIMATION FOR A CLASS OF LOSS FUNCTIONS"

Transcription

1 TH USHNG HOUS ROCDNGS OF TH ROMANAN ACADMY Series A OF TH ROMANAN ACADMY Volume 6 Number 3/005 pp AYSAN RMUM AND ASYMTOTC AROMATON FOR A CASS OF OSS FUNCTONS Roxana CUMARA Department of Mathematics Academy of conomic Studies ucharest; Roxana_Ciumara@hotmailcom We derie ayesian premium for a class of loss functions An exact calculation of ayesian premium is possible under restrictie circumstances only therefore we apply a theorem of Gatto [5] in order to derie an approximation for the premium This approximation is based on aplace's method is simply computed and allows for analytical interpretations Different forms of weighted quadratic loss functions are considered We proe that the aplace approximation of ayesian premium does not depend on the loss function inoled in exact ealuation of the ayesian premium if this loss function is of weighted quadratic type or exponentially scaled Key words: ayesian premium aplace approximation weighted quadratic loss function exponentially scaled loss function NTRODUCTON The determination of the premium for an insured risk is one of the important problems of risk theory y premium calculation it is meant the determination of the adequate premium for the risk assumed by an insured indiidual within a collectiity Moreoer if some claim experience for the gien indiidual is aailable we seek to determine the premium which exploits these actual claim amounts Therefore we compute ayesian premium which is a premium based on claim experience and it is defined as the alue which minimizes an expected loss with respect to a posterior distribution An exact calculation for the ayesian premium is possible under restrictie circumstances only regarding the prior distribution the likelihood function and the loss function That is why it would be helpful if we had a good approximation of the ayesian premium simply to compute and which allows for analytical interpretations n order to derie an approximation for the ayesian premium we employ aplace method by using a theorem of Gatto [5] We do this considering a weighted quadratic loss function mportant research concerning ayesian analysis in actuarial science is due to ailey [] uhlmann [4] Klugman [7] Mako [8] Schnieper [0] n what concerns the use of aplace approximation in actuarial science we hae to mention the work of leistein and Handelsman [3] Tierney and Kadane [] Gatto [5] AYSAN RMUM AND TS AAC AROMATON We consider random ariables n of claims in n consecutie periods for a gien indiidual Reccomended by Marius OSFSCU member of the Romanian Academy on a probability field { Ω K r} i takes positie real alues representing the size

2 et { ρ} Roxana CUMARA S be a measurable space (with measure ρ ) and Θ : Ω S a random ariable taking the alues S We now consider a probability measure Π on { S } which is absolutely continuous with respect to ρ and π is its corresponding density We say that π is the prior density of Θ The parameter (alues of Θ ) is the indiidual risk within a collectiity We assume that gien Θ the random ariables n are independently and identically distributed with common distribution function F Θ( x ) and density f Θ( x ) (for a fixed indiidual the claim amounts of consecutie periods are independently and identically distributed) et y ( y y y n ) be an obsered claim experience during the last n periods ie a realization of Θ Denote by the conditional random ector n l ( y y yn ) f Θ( yi ) the likelihood function Then the posterior distribution is gien by and the posterior density (from ayes' theorem) by Next define the loss function where ( x ) π n i ( y) r( Θ y) Π Θ ( y) l( y y yn ) π() l( y y yn ) π() S : R + R + is the loss incurred by a decision maker taking the action and facing the outcome x ( - the premium paid by the insured and x - the claim size) Definition The (indiidual) risk premium is a function : S R + which is measurable on S where { } or equialently () arg min ( ) d ( Θ ) f Θ () arg min ( x ) f Θ( x ) dx Ω if the expectation inoled in the aboe expression exists Similarly we define ayesian premium Definition The ayesian premium is the measurable function R n : gien by or equialently ( y) ( y) arg min ( ) π Θ min ( Θ ) S π ( ) y arg y d if the expectation inoled in the aboe expression exists

3 3 ayesian premium and asymptote approximation for a class of loss functions Note that the ayesian premium is the real number ( y) ( ) y π oer where Θ is the indiidual risk premium Remark n fact the ayesian premium is where ( y) min ( ) minimizing the posterior expected loss ( y) arg π Θ ( Θ ) () arg min ( ) and : are loss functions not necessarily identical f nothing else is mentioned we will consider that n order to compute the ayesian premium we need to ealuate the integral () y d R S This could be soled analitically under rather restrictie circumstances only which concerns the choice of the loss function and the prior density and the likelihood function giing the posterior density Anyway this choice should be based on objectie criteria possibly related to the real enironment of the insurer rather than on the feasibility of the ensuing calculations Of course numerical integration (wheneer it is possible) proides the most general solutions but they do not yield an analytical result which would allow for interpretations n this context the aplace method gies an asymptotic approximation of the integral aboe with a small asymptotic relati error and this allows to obtain an accurate asymptotic approximation of the ayesian premium (Gatto [5]) The results obtained in this paper are based on the result below alue of Theorem (Gatto [5]) et ( y y ) y y n denote a realization of n gien a fixed Θ ( n Θ ) and suppose that Π Θ ( y) is an absolutely continuous π Θ y which is positie oer an interal ( a b) f π Θ ( y) is differentiable distribution with density oer ( a b) has a unique maximum at ( ab) ( y) π Θ exists and is strictly negatie then (ie is the only point satisfying ( π ) 0 arg min ( ( ) ) y as n y and ( y) ( () ) ) min n Theorem (leistein Handelsman 986) et g :( a b) R + be a differentiable function oer ( ab) arg such that it has a unique minimum at an interior point which is the only point satisfying g () 0 and () g exists and is strictly positie f h ( a b ) R argument on ( a b) then : is differentiable with respect to the first ng() ng ( arg min h e d arg min h e ) as n

4 Roxana CUMARA 4 ng() ng n fact these authors proe that ( ) arg min h e d arg min h e + O t should be n mentioned that the error term contains otherwise the approximation would be an exact one The aboe result is based on aplace's method of approximation riefly this means that g and h are approximated by a Taylor polynomial of order around and then their product is integrated on a small neighborhood of Remark 3 f the loss function inoled in the calculation of ayesian premium is ( x ) different than the one inoled in indiidual risk premium ( x ) then under the assumptions of Theorem we get y arg min as n ( ) Remark 4 from Theorem is the posterior mode of Θ gien y Gatto [5] proposed an alternatie approximation of ayesian premium as y arg min ˆ as n ( ) where ˆ is the maximum likelihood estimate of the risk parameter (ie the solution of the likelihood equation l ( y y y n ) 0 ) The main disadantage of this approximation is that the prior information is lost Denote ( ( y arg min ) ) or ( y ) arg min ( ( ) ) (if ) and ( y ) arg min ( ( ˆ ) ) or arg min ( ( ˆ ) ) y (if ) n what follows we will consider a loss function of weighted quadratic type and suppose that the conditions of Theorem are fulfilled Theorem et : R+ R + be the loss function gien by ( x ) w( x)( x ) where w : R is a measurable function Then the indiidual risk premium is while the ayesian premium is gien by ( y) ( ) ( w Θ ) ( w Θ ) ( ( ) y) π w( ) y π w Furthermore an approximation of the ayesian premium is (! ) y where is the posterior mode of Θ assuming that all expectations aboe exist roof: Apply Theorem for the weighted quadratic loss function roposition Under the assumptions of Theorems and and assuming that w is differentiable of order k at 0 the indiidual risk premium is

5 5 ayesian premium and asymptote approximation for a class of loss functions () while the ayesian premium is gien by which could be approximated by ( y) k 0 k 0 k 0 k 0 ( k ) k+ ( Θ ) w ( k ) k ( Θ ) w ( k w ) π ( k w ) k 0 k + ( y) k ( y) ( k ) k+ f ( Θ ) w Θ k 0 y k w k f ( Θ ) Θ assuming that all expections inoled exist roof: Apply Theorem We will next analyze what happens if the loss functions inoled in calculations of risk premium and ayesian premium are both of weighted quadratic type but they are not identical Theorem 3 et R : be the loss functions gien by w ()( x x ) i i i where w i : R are measurable functions Suppose that the risk premium is computed from and the ayesian premium is based on loss function Then the indiidual risk premium is the ayesian premium is gien by () and an approximation of the ayesian premium is ( w ( ) Θ ) w ( ) Θ f ( ) Θ ( w ( ( Θ ) y) w ( ) y ( y) ( w ( ) Θ ) w ( ) Θ y () where is the posterior mode of Θ assuming that all expectations aboe exist roof: Apply Theorem and use the fact that w ()( x x ) i i for i Remark 5 f the loss function inoled in the calculation of the ayesian premium is of weighted quadratic type and different from the loss function inoled in ealuating the risk premium the approximation of ayesian premium does not depend on The weighted quadratic loss function inoled in the calculation of ayesian premium is not the only loss function that leads to an approximation of ayesian premium that depends only on the loss function

6 Roxana CUMARA 6 inoled in ealuating the risk premium n the next theorem we consider an exponentially scaled loss function Theorem 4 et R x : be the exponentially scaled loss function gien by ( x ) ( e e ) where > 0 Suppose that all conditions of Theorem are fulfilled Then the indiidual risk premium is while the ayesian premium is gien by () ln ( e Θ ) π ( y) ln e ( y) Furthermore an approximation of the ayesian premium is y where is the posterior mode of Θ assuming that all expectations aboe exist x roof: Apply Theorem for ( x ) ( e e ) The next propositions analyze the cases when the loss functions inoled in ealuating risk premium and ayesian premium are different and at least one of them is exponentially scaled roposition et i : x ( x ) ( e e ) R x w x x i be loss functions gien by where w : R is a measurable function and > 0 Suppose that the risk premium is computed from and the ayesian premium is based on loss function Then the indiidual risk premium is the ayesian premium is gien by () ( w ( ) Θ ) w ( ) Θ f ( ) Θ π ( y) ln e ( y) and an approximation of the ayesian premium is y where is the posterior mode of Θ assuming that all expectations aboe exist roof: Apply Theorem for the two loss functions mentioned aboe roposition 3 et R i : ( x ) w ( x)( x ) x i be loss functions gien by x e e where w : R is a measurable function and > 0 Suppose that the risk premium is computed from and the ayesian premium is based on loss function Then the indiidual risk premium is the ayesian premium is gien by ( y) () ln ( e Θ ) ( w ( ( Θ ) y) w ( ) y

7 7 ayesian premium and asymptote approximation for a class of loss functions and an approximation of the ayesian premium is ( ) y where is the posterior mode of Θ assuming that all expectations aboe exist roof: Apply Theorem for the two loss functions mentioned aboe roposition 4 et R i : βx β ( x ) ( e e ) x i be loss functions gien by x e e where β > 0 Suppose that the risk premium is computed from and the ayesian premium is based on loss function Then the indiidual risk premium is the ayesian premium is gien by and an approximation of the ayesian premium is () ln ( e Θ ) π β ( y) ln e β ( ) y ( y) where is the posterior mode of Θ assuming that all expectations aboe exist roof: Apply Theorem for the two loss functions mentioned aboe Remark 6 n all the cases presented in ropositions through 4 the approximation of the ayesian premium does not depend on the loss function inoled in the calculation of ayesian premium RFRNCS AY A Credibility procedures aplace's generalizations of ayes rule and the combination of collateral knowledge with obsered data roceedings of the Casualty Actuarial Society 37 pp RGR J O Statistical Decision Theory and ayesian Analysis New York Springer STN N HANDSMAN RA Asymptotic xpansions of ntegrals nd dition Doer UHMANN H Mathematical Methods in Risk Theory Springer GATTO R An accurate asymptotic approximation for experience rated premiums Astin ulletin 34 pp OSFSCU M MHOC G THODORSCU R Teoria robabilităţilor şi Statistica Matematică ucureşti ditura Tehnică KUGMAN S ayesian Statistics in Actuarial Science Kluwer 99 8 MAKOV U oss robustness ia Fisher-weighted squared-error loss function nsurance: Mathematics and conomics 6 pp RDA V Statistical Decision Theory ucharest Romanian Academy ubl House 99 0 SCHNR R Robust ayesian experience rating Astin ulletin 34 pp TRNY KADAN J Accurate approximations for posterior moments and marginal densities Journal of the American Statistical Association 8 pp ZĂGANU G Metode matematice în teoria riscului şi actuariat ucureşti ditura Uniersităţii 004 Receied h September 9 005

Final Exam (Solution) Economics 501b Microeconomic Theory

Final Exam (Solution) Economics 501b Microeconomic Theory Dirk Bergemann and Johannes Hoerner Department of Economics Yale Uniersity Final Exam (Solution) Economics 5b Microeconomic Theory May This is a closed-book exam. The exam lasts for 8 minutes. Please write

More information

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS AP PHYSICS

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS AP PHYSICS DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS AP PHYSICS LSN 7-4: CONSERVATION OF ENERGY AND MOMENTUM IN COLLISIONS LSN 7-5: ELASTIC COLLISIONS IN ONE DIMENSION LSN 7-6: INELASTIC COLLISIONS Questions From

More information

Online Companion to Pricing Services Subject to Congestion: Charge Per-Use Fees or Sell Subscriptions?

Online Companion to Pricing Services Subject to Congestion: Charge Per-Use Fees or Sell Subscriptions? Online Companion to Pricing Serices Subject to Congestion: Charge Per-Use Fees or Sell Subscriptions? Gérard P. Cachon Pnina Feldman Operations and Information Management, The Wharton School, Uniersity

More information

Algorithms and Data Structures 2014 Exercises and Solutions Week 14

Algorithms and Data Structures 2014 Exercises and Solutions Week 14 lgorithms and ata tructures 0 xercises and s Week Linear programming. onsider the following linear program. maximize x y 0 x + y 7 x x 0 y 0 x + 3y Plot the feasible region and identify the optimal solution.

More information

THE FIFTH DIMENSION EQUATIONS

THE FIFTH DIMENSION EQUATIONS JP Journal of Mathematical Sciences Volume 7 Issues 1 & 013 Pages 41-46 013 Ishaan Publishing House This paper is aailable online at http://www.iphsci.com THE FIFTH DIMENSION EQUATIONS Niittytie 1B16 03100

More information

Econometrics II - EXAM Outline Solutions All questions have 25pts Answer each question in separate sheets

Econometrics II - EXAM Outline Solutions All questions have 25pts Answer each question in separate sheets Econometrics II - EXAM Outline Solutions All questions hae 5pts Answer each question in separate sheets. Consider the two linear simultaneous equations G with two exogeneous ariables K, y γ + y γ + x δ

More information

Collective Risk Models with Dependence Uncertainty

Collective Risk Models with Dependence Uncertainty Collectie Risk Models with Dependence Uncertainty Haiyan Liu and Ruodu Wang ebruary 27, 2017 Abstract We bring the recently deeloped framework of dependence uncertainty into collectie risk models, one

More information

Would you risk your live driving drunk? Intro

Would you risk your live driving drunk? Intro Martha Casquete Would you risk your lie driing drunk? Intro Motion Position and displacement Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Constant acceleration: A special

More information

Chapter 4: Techniques of Circuit Analysis

Chapter 4: Techniques of Circuit Analysis Chapter 4: Techniques of Circuit Analysis This chapter gies us many useful tools for soling and simplifying circuits. We saw a few simple tools in the last chapter (reduction of circuits ia series and

More information

Decision Making in Committees: Transparency, Reputation and Voting Rules

Decision Making in Committees: Transparency, Reputation and Voting Rules Decision Making in Committees: Transparency, Reputation and Voting Rules Gilat Ley, Department of Economics, LSE Noember 16, 2005 Abstract In this paper I analyze the e ect of the transparency of the decision

More information

Collective circular motion of multi-vehicle systems with sensory limitations

Collective circular motion of multi-vehicle systems with sensory limitations Collectie circular motion of multi-ehicle systems with sensory limitations Nicola Ceccarelli, Mauro Di Marco, Andrea Garulli, Antonello Giannitrapani Dipartimento di Ingegneria dell Informazione Uniersità

More information

On the Linear Threshold Model for Diffusion of Innovations in Multiplex Social Networks

On the Linear Threshold Model for Diffusion of Innovations in Multiplex Social Networks On the Linear Threshold Model for Diffusion of Innoations in Multiplex Social Networks Yaofeng Desmond Zhong 1, Vaibha Sriastaa 2 and Naomi Ehrich Leonard 1 Abstract Diffusion of innoations in social networks

More information

The Inverse Function Theorem

The Inverse Function Theorem Inerse Function Theorem April, 3 The Inerse Function Theorem Suppose and Y are Banach spaces and f : Y is C (continuousl differentiable) Under what circumstances does a (local) inerse f exist? In the simplest

More information

Motion in Two and Three Dimensions

Motion in Two and Three Dimensions PH 1-1D Spring 013 Motion in Two and Three Dimensions Lectures 5,6,7 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter

More information

(a) Taking the derivative of the position vector with respect to time, we have, in SI units (m/s),

(a) Taking the derivative of the position vector with respect to time, we have, in SI units (m/s), Chapter 4 Student Solutions Manual. We apply Eq. 4- and Eq. 4-6. (a) Taking the deriatie of the position ector with respect to time, we hae, in SI units (m/s), d ˆ = (i + 4t ˆj + tk) ˆ = 8tˆj + k ˆ. dt

More information

ADVANCED General Certificate of Education Mathematics Assessment Unit M2. assessing. Module M2: Mechanics 2 [AMM21] THURSDAY 2 JUNE, AFTERNOON

ADVANCED General Certificate of Education Mathematics Assessment Unit M2. assessing. Module M2: Mechanics 2 [AMM21] THURSDAY 2 JUNE, AFTERNOON ADVANCED General Certificate of Education 206 Mathematics Assessment Unit M2 assessing Module M2: Mechanics 2 [AMM2] THURSDAY 2 JUNE, AFTERNOON MARK SCHEME 9846.0 F GCE ADVANCED/ADVANCED SUBSIDIARY (AS)

More information

Section 6: PRISMATIC BEAMS. Beam Theory

Section 6: PRISMATIC BEAMS. Beam Theory Beam Theory There are two types of beam theory aailable to craft beam element formulations from. They are Bernoulli-Euler beam theory Timoshenko beam theory One learns the details of Bernoulli-Euler beam

More information

Network Flow Problems Luis Goddyn, Math 408

Network Flow Problems Luis Goddyn, Math 408 Network Flow Problems Luis Goddyn, Math 48 Let D = (V, A) be a directed graph, and let s, t V (D). For S V we write δ + (S) = {u A : u S, S} and δ (S) = {u A : u S, S} for the in-arcs and out-arcs of S

More information

Motion in Two and Three Dimensions

Motion in Two and Three Dimensions PH 1-A Fall 014 Motion in Two and Three Dimensions Lectures 4,5 Chapter 4 (Halliday/Resnick/Walker, Fundamentals of Physics 9 th edition) 1 Chapter 4 Motion in Two and Three Dimensions In this chapter

More information

Two-sided bounds for L p -norms of combinations of products of independent random variables

Two-sided bounds for L p -norms of combinations of products of independent random variables Two-sided bounds for L p -norms of combinations of products of independent random ariables Ewa Damek (based on the joint work with Rafał Latała, Piotr Nayar and Tomasz Tkocz) Wrocław Uniersity, Uniersity

More information

Math Mathematical Notation

Math Mathematical Notation Math 160 - Mathematical Notation Purpose One goal in any course is to properly use the language of that subject. Finite Mathematics is no different and may often seem like a foreign language. These notations

More information

Review of Matrices and Vectors 1/45

Review of Matrices and Vectors 1/45 Reiew of Matrices and Vectors /45 /45 Definition of Vector: A collection of comple or real numbers, generally put in a column [ ] T "! Transpose + + + b a b a b b a a " " " b a b a Definition of Vector

More information

Some Lecture Notes on Auctions

Some Lecture Notes on Auctions Some Lecture Notes on Auctions John Morgan February 25 1 Introduction These notes sketch a 1 hour and 5 minute lecture on auctions. The main points are: 1. Set up the basic auction problem. 2. Derie the

More information

Collective circular motion of multi-vehicle systems with sensory limitations

Collective circular motion of multi-vehicle systems with sensory limitations Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 2005 Seille, Spain, December 12-15, 2005 MoB032 Collectie circular motion of multi-ehicle systems with

More information

The Land Acquisition Problem: a Mechanism Design Approach

The Land Acquisition Problem: a Mechanism Design Approach The Land Acquisition Problem: a Mechanism Design Approach Soumendu Sarkar February 2012 (Preliminary and Incomplete. Please do not quote.) Abstract The land acquisition problem is an exchange problem with

More information

4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion.

4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion. 4. A Physical Model for an Electron with Angular Momentum. An Electron in a Bohr Orbit. The Quantum Magnet Resulting from Orbital Motion. We now hae deeloped a ector model that allows the ready isualization

More information

On resilience of distributed routing in networks under cascade dynamics

On resilience of distributed routing in networks under cascade dynamics On resilience of distributed routing in networks under cascade dynamics Ketan Sala Giacomo Como Munther A. Dahleh Emilio Frazzoli Abstract We consider network flow oer graphs between a single origin-destination

More information

Parameters Identification of Equivalent Circuit Diagrams for Li-Ion Batteries

Parameters Identification of Equivalent Circuit Diagrams for Li-Ion Batteries Parameters Identification of Equialent Circuit Diagrams for Li-Ion eries Ahmad ahmoun, Helmuth Biechl Uniersity of Applied ciences Kempten Ahmad.ahmoun@stud.fh-empten.de, biechl@fh-empten.de Abstract-eries

More information

Lect-19. In this lecture...

Lect-19. In this lecture... 19 1 In this lecture... Helmholtz and Gibb s functions Legendre transformations Thermodynamic potentials The Maxwell relations The ideal gas equation of state Compressibility factor Other equations of

More information

LECTURE 2: CROSS PRODUCTS, MULTILINEARITY, AND AREAS OF PARALLELOGRAMS

LECTURE 2: CROSS PRODUCTS, MULTILINEARITY, AND AREAS OF PARALLELOGRAMS LECTURE : CROSS PRODUCTS, MULTILINEARITY, AND AREAS OF PARALLELOGRAMS MA1111: LINEAR ALGEBRA I, MICHAELMAS 016 1. Finishing up dot products Last time we stated the following theorem, for which I owe you

More information

UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS

UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS SUBAREA I. COMPETENCY 1.0 UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS MECHANICS Skill 1.1 Calculating displacement, aerage elocity, instantaneous elocity, and acceleration in a gien frame of reference

More information

Towards Universal Cover Decoding

Towards Universal Cover Decoding International Symposium on Information Theory and its Applications, ISITA2008 Auckland, New Zealand, 7-10, December, 2008 Towards Uniersal Coer Decoding Nathan Axig, Deanna Dreher, Katherine Morrison,

More information

Optimal Joint Detection and Estimation in Linear Models

Optimal Joint Detection and Estimation in Linear Models Optimal Joint Detection and Estimation in Linear Models Jianshu Chen, Yue Zhao, Andrea Goldsmith, and H. Vincent Poor Abstract The problem of optimal joint detection and estimation in linear models with

More information

Patterns of Non-Simple Continued Fractions

Patterns of Non-Simple Continued Fractions Patterns of Non-Simple Continued Fractions Jesse Schmieg A final report written for the Uniersity of Minnesota Undergraduate Research Opportunities Program Adisor: Professor John Greene March 01 Contents

More information

An Alternative Characterization of Hidden Regular Variation in Joint Tail Modeling

An Alternative Characterization of Hidden Regular Variation in Joint Tail Modeling An Alternatie Characterization of Hidden Regular Variation in Joint Tail Modeling Grant B. Weller Daniel S. Cooley Department of Statistics, Colorado State Uniersity, Fort Collins, CO USA May 22, 212 Abstract

More information

Chapter 11 Collision Theory

Chapter 11 Collision Theory Chapter Collision Theory Introduction. Center o Mass Reerence Frame Consider two particles o masses m and m interacting ia some orce. Figure. Center o Mass o a system o two interacting particles Choose

More information

A Geometric Review of Linear Algebra

A Geometric Review of Linear Algebra A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. I assume the reader is familiar with basic (i.e., high school) algebra and trigonometry.

More information

different formulas, depending on whether or not the vector is in two dimensions or three dimensions.

different formulas, depending on whether or not the vector is in two dimensions or three dimensions. ectors The word ector comes from the Latin word ectus which means carried. It is best to think of a ector as the displacement from an initial point P to a terminal point Q. Such a ector is expressed as

More information

CHAPTER (i) No. For each coefficient, the usual standard errors and the heteroskedasticity-robust ones are practically very similar.

CHAPTER (i) No. For each coefficient, the usual standard errors and the heteroskedasticity-robust ones are practically very similar. SOLUTIONS TO PROBLEMS CHAPTER 8 8.1 Parts (ii) and (iii). The homoskedasticity assumption played no role in Chapter 5 in showing that OLS is consistent. But we know that heteroskedasticity causes statistical

More information

General Lorentz Boost Transformations, Acting on Some Important Physical Quantities

General Lorentz Boost Transformations, Acting on Some Important Physical Quantities General Lorentz Boost Transformations, Acting on Some Important Physical Quantities We are interested in transforming measurements made in a reference frame O into measurements of the same quantities as

More information

PRODUCTS IN CONDITIONAL EXTREME VALUE MODEL

PRODUCTS IN CONDITIONAL EXTREME VALUE MODEL PRODUCTS IN CONDITIONAL ETREME VALUE MODEL RAJAT SUBHRA HAZRA AND KRISHANU MAULIK Abstract. The classical multiariate extreme alue theory tries to capture the extremal dependence between the components

More information

Insights into Cross-validation

Insights into Cross-validation Noname manuscript No. (will be inserted by the editor) Insights into Cross-alidation Amit Dhurandhar Alin Dobra Receied: date / Accepted: date Abstract Cross-alidation is one of the most widely used techniques,

More information

On computing Gaussian curvature of some well known distribution

On computing Gaussian curvature of some well known distribution Theoretical Mathematics & Applications, ol.3, no.4, 03, 85-04 ISSN: 79-9687 (print), 79-9709 (online) Scienpress Ltd, 03 On computing Gaussian curature of some well known distribution William W.S. Chen

More information

Bullwhip and Reverse Bullwhip Effects Under the Rationing Game

Bullwhip and Reverse Bullwhip Effects Under the Rationing Game Bullwhip and Reerse Bullwhip Effects Under the Rationing Game Ying Rong Antai College of Economics and Management, Shanghai Jiao Tong Uniersity, Shanghai, China 200052, yrong@sjtu.edu.cn Lawrence V. Snyder

More information

Kinematics on oblique axes

Kinematics on oblique axes Bolina 1 Kinematics on oblique axes arxi:physics/01111951 [physics.ed-ph] 27 No 2001 Oscar Bolina Departamento de Física-Matemática Uniersidade de São Paulo Caixa Postal 66318 São Paulo 05315-970 Brasil

More information

Mean-variance receding horizon control for discrete time linear stochastic systems

Mean-variance receding horizon control for discrete time linear stochastic systems Proceedings o the 17th World Congress The International Federation o Automatic Control Seoul, Korea, July 6 11, 008 Mean-ariance receding horizon control or discrete time linear stochastic systems Mar

More information

An Optimal Split-Plot Design for Performing a Mixture-Process Experiment

An Optimal Split-Plot Design for Performing a Mixture-Process Experiment Science Journal of Applied Mathematics and Statistics 217; 5(1): 15-23 http://www.sciencepublishinggroup.com/j/sjams doi: 1.11648/j.sjams.21751.13 ISSN: 2376-9491 (Print); ISSN: 2376-9513 (Online) An Optimal

More information

Chapter 2: 1D Kinematics Tuesday January 13th

Chapter 2: 1D Kinematics Tuesday January 13th Chapter : D Kinematics Tuesday January 3th Motion in a straight line (D Kinematics) Aerage elocity and aerage speed Instantaneous elocity and speed Acceleration Short summary Constant acceleration a special

More information

Astrometric Errors Correlated Strongly Across Multiple SIRTF Images

Astrometric Errors Correlated Strongly Across Multiple SIRTF Images Astrometric Errors Correlated Strongly Across Multiple SIRTF Images John Fowler 28 March 23 The possibility exists that after pointing transfer has been performed for each BCD (i.e. a calibrated image

More information

Some Lecture Notes on Auctions

Some Lecture Notes on Auctions Some Lecture Notes on Auctions John Morgan Haas School of Business and Department of Economics Uniersity of California, Berkeley Preliminaries Perhaps the most fruitful area for the application of optimal

More information

(a) During the first part of the motion, the displacement is x 1 = 40 km and the time interval is t 1 (30 km / h) (80 km) 40 km/h. t. (2.

(a) During the first part of the motion, the displacement is x 1 = 40 km and the time interval is t 1 (30 km / h) (80 km) 40 km/h. t. (2. Chapter 3. Since the trip consists of two parts, let the displacements during first and second parts of the motion be x and x, and the corresponding time interals be t and t, respectiely. Now, because

More information

A Geometric Review of Linear Algebra

A Geometric Review of Linear Algebra A Geometric Reiew of Linear Algebra The following is a compact reiew of the primary concepts of linear algebra. The order of presentation is unconentional, with emphasis on geometric intuition rather than

More information

PHYS 1443 Section 004 Lecture #4 Thursday, Sept. 4, 2014

PHYS 1443 Section 004 Lecture #4 Thursday, Sept. 4, 2014 PHYS 1443 Section 004 Lecture #4 Thursday, Sept. 4, 014 One Dimensional Motion Motion under constant acceleration One dimensional Kinematic Equations How do we sole kinematic problems? Falling motions

More information

Magnetic Fields Part 3: Electromagnetic Induction

Magnetic Fields Part 3: Electromagnetic Induction Magnetic Fields Part 3: Electromagnetic Induction Last modified: 15/12/2017 Contents Links Electromagnetic Induction Induced EMF Induced Current Induction & Magnetic Flux Magnetic Flux Change in Flux Faraday

More information

Dynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates

Dynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates ACC 9, Submitted St. Louis, MO Dynamic Vehicle Routing with Moing Demands Part II: High speed demands or low arrial rates Stephen L. Smith Shaunak D. Bopardikar Francesco Bullo João P. Hespanha Abstract

More information

arxiv: v1 [math.co] 25 Apr 2016

arxiv: v1 [math.co] 25 Apr 2016 On the zone complexity of a ertex Shira Zerbib arxi:604.07268 [math.co] 25 Apr 206 April 26, 206 Abstract Let L be a set of n lines in the real projectie plane in general position. We show that there exists

More information

Geostrophy & Thermal wind

Geostrophy & Thermal wind Lecture 10 Geostrophy & Thermal wind 10.1 f and β planes These are planes that are tangent to the earth (taken to be spherical) at a point of interest. The z ais is perpendicular to the plane (anti-parallel

More information

NEW FRONTIERS IN APPLIED PROBABILITY

NEW FRONTIERS IN APPLIED PROBABILITY J. Appl. Prob. Spec. Vol. 48A, 183 194 (2011) Applied Probability Trust 2011 NEW FRONTIERS IN APPLIED PROBABILITY A Festschrift for SØREN ASMUSSEN Edited by P. GLYNN, T. MIKOSCH and T. ROLSKI Part 4. Simulation

More information

Should Social Network Structure Be Taken into Account in Elections?

Should Social Network Structure Be Taken into Account in Elections? Should Social Network Structure Be Taken into Account in Elections? Vincent Conitzer Duke Uniersity April 2, 2011 Abstract If the social network structure among the oters in an election is known, how should

More information

VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION

VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION VISUAL PHYSICS ONLINE RECTLINEAR MOTION: UNIFORM ACCELERATION Predict Obsere Explain Exercise 1 Take an A4 sheet of paper and a heay object (cricket ball, basketball, brick, book, etc). Predict what will

More information

A nonstandard cubic equation

A nonstandard cubic equation MATH-Jan-05-0 A nonstandard cubic euation J S Markoitch PO Box West Brattleboro, VT 050 Dated: January, 05 A nonstandard cubic euation is shown to hae an unusually econoical solution, this solution incorporates

More information

Efficient solution of interval optimization problem

Efficient solution of interval optimization problem Math Meth Oper Res (212) 76:273 288 DOI 1.17/s186-12-399- ORIGINAL ARTICLE Efficient solution of interal optimization problem A. K. Bhurjee G. Panda Receied: 6 June 212 / Accepted: 4 August 212 / Published

More information

OPTIMIZATION OF FLOWS AND ANALYSIS OF EQUILIBRIA IN TELECOMMUNICATION NETWORKS

OPTIMIZATION OF FLOWS AND ANALYSIS OF EQUILIBRIA IN TELECOMMUNICATION NETWORKS OPTIMIZATION OF FLOWS AND ANALYSIS OF EQUILIBRIA IN TELECOMMUNICATION NETWORKS Ciro D Apice (a, Rosanna Manzo (b, Luigi Rarità (c (a, (b, (c Department of Information Engineering and Applied Mathematics,

More information

Unit 11: Vectors in the Plane

Unit 11: Vectors in the Plane 135 Unit 11: Vectors in the Plane Vectors in the Plane The term ector is used to indicate a quantity (such as force or elocity) that has both length and direction. For instance, suppose a particle moes

More information

DESIGN METHOD BASED ON THE CONCEPT OF EMERGENCE AND ITS APPLICATION

DESIGN METHOD BASED ON THE CONCEPT OF EMERGENCE AND ITS APPLICATION DESIGN METHOD BASED ON THE CONCEPT OF EMERGENCE AND ITS APPLICATION Koichiro Sato¹ and Yoshiyuki Matsuoka 2 ¹ Graduate School of Science and Technology, Keio Uniersity, Yokohama, Japan Koichiro_Sato@a7.keio.jp,

More information

Dynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates

Dynamic Vehicle Routing with Moving Demands Part II: High speed demands or low arrival rates Dynamic Vehicle Routing with Moing Demands Part II: High speed demands or low arrial rates Stephen L. Smith Shaunak D. Bopardikar Francesco Bullo João P. Hespanha Abstract In the companion paper we introduced

More information

On general error distributions

On general error distributions ProbStat Forum, Volume 6, October 3, Pages 89 95 ISSN 974-335 ProbStat Forum is an e-journal. For details please isit www.probstat.org.in On general error distributions R. Vasudea, J. Vasantha Kumari Department

More information

1 :: Mathematical notation

1 :: Mathematical notation 1 :: Mathematical notation x A means x is a member of the set A. A B means the set A is contained in the set B. {a 1,..., a n } means the set hose elements are a 1,..., a n. {x A : P } means the set of

More information

Physics 4A Solutions to Chapter 4 Homework

Physics 4A Solutions to Chapter 4 Homework Physics 4A Solutions to Chapter 4 Homework Chapter 4 Questions: 4, 1, 1 Exercises & Problems: 5, 11, 3, 7, 8, 58, 67, 77, 87, 11 Answers to Questions: Q 4-4 (a) all tie (b) 1 and tie (the rocket is shot

More information

Scalar multiplication and algebraic direction of a vector

Scalar multiplication and algebraic direction of a vector Roberto s Notes on Linear Algebra Chapter 1: Geometric ectors Section 5 Scalar multiplication and algebraic direction of a ector What you need to know already: of a geometric ectors. Length and geometric

More information

Position in the xy plane y position x position

Position in the xy plane y position x position Robust Control of an Underactuated Surface Vessel with Thruster Dynamics K. Y. Pettersen and O. Egeland Department of Engineering Cybernetics Norwegian Uniersity of Science and Technology N- Trondheim,

More information

Towards Green Distributed Storage Systems

Towards Green Distributed Storage Systems Towards Green Distributed Storage Systems Abdelrahman M. Ibrahim, Ahmed A. Zewail, and Aylin Yener Wireless Communications and Networking Laboratory (WCAN) Electrical Engineering Department The Pennsylania

More information

Dynamic Vehicle Routing for Translating Demands: Stability Analysis and Receding-Horizon Policies

Dynamic Vehicle Routing for Translating Demands: Stability Analysis and Receding-Horizon Policies 1 Dynamic Vehicle Routing for Translating Demands: Stability Analysis and Receding-Horizon Policies Shaunak D. Bopardikar Stephen L. Smith Francesco Bullo João P. Hespanha Abstract We introduce a problem

More information

LECTURE 3 3.1Rules of Vector Differentiation

LECTURE 3 3.1Rules of Vector Differentiation LETURE 3 3.1Rules of Vector Differentiation We hae defined three kinds of deriaties inoling the operator grad( ) i j k, x y z 1 3 di(., x y z curl( i x 1 j y k z 3 d The good news is that you can apply

More information

Displacement, Time, Velocity

Displacement, Time, Velocity Lecture. Chapter : Motion along a Straight Line Displacement, Time, Velocity 3/6/05 One-Dimensional Motion The area of physics that we focus on is called mechanics: the study of the relationships between

More information

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems

COE-RES Discussion Paper Series Center of Excellence Project The Normative Evaluation and Social Choice of Contemporary Economic Systems COE-RES Discussion Paper Series Center of Excellence Project The Normatie Ealuation and Social Choice of Contemporary Economic Systems Graduate School of Economics and Institute of Economic Research Hitotsubashi

More information

LESSON 4: INTEGRATION BY PARTS (I) MATH FALL 2018

LESSON 4: INTEGRATION BY PARTS (I) MATH FALL 2018 LESSON 4: INTEGRATION BY PARTS (I) MATH 6 FALL 8 ELLEN WELD. Integration by Parts We introduce another method for ealuating integrals called integration by parts. The key is the following : () u d = u

More information

Probabilistic Engineering Design

Probabilistic Engineering Design Probabilistic Engineering Design Chapter One Introduction 1 Introduction This chapter discusses the basics of probabilistic engineering design. Its tutorial-style is designed to allow the reader to quickly

More information

Note: the net distance along the path is a scalar quantity its direction is not important so the average speed is also a scalar.

Note: the net distance along the path is a scalar quantity its direction is not important so the average speed is also a scalar. PHY 309 K. Solutions for the first mid-term test /13/014). Problem #1: By definition, aerage speed net distance along the path of motion time. 1) ote: the net distance along the path is a scalar quantity

More information

Dynamic Vehicle Routing with Moving Demands Part I: Low speed demands and high arrival rates

Dynamic Vehicle Routing with Moving Demands Part I: Low speed demands and high arrival rates Dynamic Vehicle Routing with Moing Demands Part I: Low speed demands and high arrial rates Shaunak D. Bopardikar Stephen L. Smith Francesco Bullo João P. Hespanha Abstract We introduce a dynamic ehicle

More information

PHYS1169: Tutorial 8 Solutions

PHYS1169: Tutorial 8 Solutions PHY69: Tutorial 8 olutions Wae Motion ) Let us consier a point P on the wae with a phase φ, so y cosϕ cos( x ± ωt) At t0, this point has position x0, so ϕ x0 ± ωt0 Now, at some later time t, the position

More information

CHARACTERIZATIONS OF THE PARETO DISTRIBUTION BY THE INDEPENDENCE OF RECORD VALUES. Se-Kyung Chang* 1. Introduction

CHARACTERIZATIONS OF THE PARETO DISTRIBUTION BY THE INDEPENDENCE OF RECORD VALUES. Se-Kyung Chang* 1. Introduction JOURNAL OF THE CHUNGCHEONG MATHEMATICAL SOCIETY Volume 20, No., March 2007 CHARACTERIZATIONS OF THE PARETO DISTRIBUTION BY THE INDEPENDENCE OF RECORD VALUES Se-Kyung Chang* Abstract. In this paper, we

More information

Reaction/Diffusion at Electrode/Solution Interfaces: The EC 2 Reaction

Reaction/Diffusion at Electrode/Solution Interfaces: The EC 2 Reaction Int. J. Electrochem. Sci., 4(9) 1116-117 International Journal of ELECTROCHEMICAL SCIENCE www.electrochemsci.org Reaction/Diffusion at Electrode/Solution Interfaces: The EC Reaction Michael E G Lyons *

More information

OPTIMAL RESOLVABLE DESIGNS WITH MINIMUM PV ABERRATION

OPTIMAL RESOLVABLE DESIGNS WITH MINIMUM PV ABERRATION Statistica Sinica 0 (010), 715-73 OPTIMAL RESOLVABLE DESIGNS WITH MINIMUM PV ABERRATION J. P. Morgan Virginia Tech Abstract: Amongst resolable incomplete block designs, affine resolable designs are optimal

More information

Mon Apr dot product, length, orthogonality, projection onto the span of a single vector. Announcements: Warm-up Exercise:

Mon Apr dot product, length, orthogonality, projection onto the span of a single vector. Announcements: Warm-up Exercise: Math 2270-004 Week 2 notes We will not necessarily finish the material from a gien day's notes on that day. We may also add or subtract some material as the week progresses, but these notes represent an

More information

Dynamic Lot Size Problems with One-way Product Substitution 1

Dynamic Lot Size Problems with One-way Product Substitution 1 This is the Pre-Published Version. 03/2002/100057 ynamic Lot Size Problems with One-way Product Substitution 1 Vernon Ning Hsu School of Management George Mason Uniersity Fairfax, VA 22030, USA email hsu@gmu.edu

More information

Transmission lines using a distributed equivalent circuit

Transmission lines using a distributed equivalent circuit Cambridge Uniersity Press 978-1-107-02600-1 - Transmission Lines Equialent Circuits, Electromagnetic Theory, and Photons Part 1 Transmission lines using a distributed equialent circuit in this web serice

More information

Chromatic characterization of biclique cover. Denis Cornaz. Jean Fonlupt. Equipe Combinatoire, UFR 921, Universite Pierre et Marie Curie, 4 place

Chromatic characterization of biclique cover. Denis Cornaz. Jean Fonlupt. Equipe Combinatoire, UFR 921, Universite Pierre et Marie Curie, 4 place Chromatic characterization of biclique coer Denis Cornaz France Telecom R&D, 38-40 rue du General Leclerc, Issy-les-Moulineaux 92794 Jean Fonlupt Equipe Combinatoire, UFR 921, Uniersite Pierre et Marie

More information

Dynamics ( 동역학 ) Ch.2 Motion of Translating Bodies (2.1 & 2.2)

Dynamics ( 동역학 ) Ch.2 Motion of Translating Bodies (2.1 & 2.2) Dynamics ( 동역학 ) Ch. Motion of Translating Bodies (. &.) Motion of Translating Bodies This chapter is usually referred to as Kinematics of Particles. Particles: In dynamics, a particle is a body without

More information

Dynamic Vehicle Routing with Moving Demands Part I: Low speed demands and high arrival rates

Dynamic Vehicle Routing with Moving Demands Part I: Low speed demands and high arrival rates ACC 2009, Submitted St. Louis, MO Dynamic Vehicle Routing with Moing Demands Part I: Low speed demands and high arrial rates Shaunak D. Bopardikar Stephen L. Smith Francesco Bullo João P. Hespanha Abstract

More information

NP and NP-Completeness

NP and NP-Completeness 0/2/206 Algorithms NP-Completeness 7- Algorithms NP-Completeness 7-2 Efficient Certification NP and NP-Completeness By a solution of a decision problem X we understand a certificate witnessing that an

More information

Balanced Partitions of Vector Sequences

Balanced Partitions of Vector Sequences Balanced Partitions of Vector Sequences Imre Bárány Benjamin Doerr December 20, 2004 Abstract Let d,r N and be any norm on R d. Let B denote the unit ball with respect to this norm. We show that any sequence

More information

A possible mechanism to explain wave-particle duality L D HOWE No current affiliation PACS Numbers: r, w, k

A possible mechanism to explain wave-particle duality L D HOWE No current affiliation PACS Numbers: r, w, k A possible mechanism to explain wae-particle duality L D HOWE No current affiliation PACS Numbers: 0.50.-r, 03.65.-w, 05.60.-k Abstract The relationship between light speed energy and the kinetic energy

More information

The Riemann-Roch Theorem: a Proof, an Extension, and an Application MATH 780, Spring 2019

The Riemann-Roch Theorem: a Proof, an Extension, and an Application MATH 780, Spring 2019 The Riemann-Roch Theorem: a Proof, an Extension, and an Application MATH 780, Spring 2019 This handout continues the notational conentions of the preious one on the Riemann-Roch Theorem, with one slight

More information

Advances in Environmental Biology

Advances in Environmental Biology Adances in Enironmental Biology, 9() Special 5, Pages: 6- AENSI Journals Adances in Enironmental Biology ISSN-995-756 EISSN-998-66 Journal home page: http://www.aensiweb.com/aeb/ Mathematical Model for

More information

Notes on Linear Minimum Mean Square Error Estimators

Notes on Linear Minimum Mean Square Error Estimators Notes on Linear Minimum Mean Square Error Estimators Ça gatay Candan January, 0 Abstract Some connections between linear minimum mean square error estimators, maximum output SNR filters and the least square

More information

Stable Network Flow with Piece-wise Linear Constraints

Stable Network Flow with Piece-wise Linear Constraints Stable Network Flow with Piece-wise Linear Constraints Young-San Lin Thành Nguyen Abstract We consider a general stable flow problem in a directed and capacitated network, where each ertex has a strict

More information

NOTES ON THE REGULAR E-OPTIMAL SPRING BALANCE WEIGHING DESIGNS WITH CORRE- LATED ERRORS

NOTES ON THE REGULAR E-OPTIMAL SPRING BALANCE WEIGHING DESIGNS WITH CORRE- LATED ERRORS REVSTAT Statistical Journal Volume 3, Number 2, June 205, 9 29 NOTES ON THE REGULAR E-OPTIMAL SPRING BALANCE WEIGHING DESIGNS WITH CORRE- LATED ERRORS Authors: Bronis law Ceranka Department of Mathematical

More information

0 a 3 a 2 a 3 0 a 1 a 2 a 1 0

0 a 3 a 2 a 3 0 a 1 a 2 a 1 0 Chapter Flow kinematics Vector and tensor formulae This introductory section presents a brief account of different definitions of ector and tensor analysis that will be used in the following chapters.

More information

P(c y) = W n (y c)p(c)

P(c y) = W n (y c)p(c) ENEE 739C: Adanced Topics in Signal Processing: Coding Theory Instructor: Alexander Barg Lecture 9 (draft; /7/3). Codes defined on graphs. The problem of attaining capacity. Iteratie decoding. http://www.enee.umd.edu/

More information