Product Integration. Richard D. Gill. Mathematical Institute, University of Utrecht, Netherlands EURANDOM, Eindhoven, Netherlands August 9, 2001
|
|
- Chastity Ball
- 6 years ago
- Views:
Transcription
1 Produc Inegraion Richard D. Gill Mahemaical Insiue, Universiy of Urech, Neherlands EURANDOM, Eindhoven, Neherlands Augus 9, 21 Absrac This is a brief survey of produc-inegraion for biosaisicians. 1 Produc-Inegraion Produc-inegraion was inroduced more han 11 years ago by he Ialian mahemaician Vio Volerra, as a ool in he soluion of a cerain class of differenial equaions. I was sudied inensively by mahemaicians for half a cenury, bu finally he subjec became unfashionable and lapsed ino obscuriy. Tha is a piy, since ideas of produc-inegraion make a very naural appearance in survival analysis, and he developmen of his subjec (in paricular, of he Kaplan-Meier esimaor) could have been a lo smooher if produc-inegraion had been a familiar opic from he sar. The Kaplan-Meier esimaor is he produc-inegral of he Nelson-Aalen esimaor of he cumulaive hazard funcion; hese wo esimaors bear he same relaion o one anoher as he acual survival funcion and he acual cumulaive hazard funcion. There are many oher applicaions of produc-inegraion in survival analysis, for insance in he sudy of muli-sae processes (conneced o he heory of Markov processes), and in he heory of parial likelihood. Ordinary inegraion is a generalisaion of summaion, and properies of inegrals are ofen easily guessed by hinking of hem as sums of very, very many erms (all or mos of hem being very small). Similarly, produc-inegraion generalises he aking of producs; a produc inegral is a produc of many, many erms (all or mos of hem being very close o he number 1). Thinking of produc-inegrals in his simplisic way is acually very helpful. Properies of produc-inegrals are easy o guess and o undersand. The heory of producinegraion can be a grea help in sudying he saisical properies of saisical quaniies which explicily or implicily are defined in erms of produc-inegrals. Before defining produc-inegrals in general and exhibiing some of heir properies, we will discuss he relaion, in survival analysis, beween survival 1
2 funcion and hazard funcion. This will lead us naurally o he noion of produc-inegraion in he mos simple possible of conexs. Consider a survival ime T wih survival funcion S() = Pr(T > ), ; S() = 1. Suppose T is coninuously disribued wih a densiy f() and a hazard rae α(). These wo funcions have inuiive probabilisic meanings: for a small ime inerval, + h, he uncondiional probabiliy Pr( T + h) f().h; while he condiional probabiliy Pr( T + h T ) α().h. In fac, he probabiliy densiy f() = (d/d)s() while he hazard rae α() = f()/s(). One can mahemaically recover he disribuion funcion F () = 1 S() from he densiy by inegraion; F () = f(s)ds. Also one can recover he survival funcion from he hazard rae. Noing ha α() = (d/d) log S() one finds by inegraion (and using ha S() = 1, hence log S() = ), ha log S() = α(s)ds hence S() = exp( α(s)ds). This is simple enough, bu neiher he resul nor is derivaion have a probabilisic inerpreaion. If he survival ime T had a discree disribuion, one would inroduce he discree densiy f() = Pr(T = ) and he discree hazard α() = Pr(T = T ) = f()/s( ). Sill he survival funcion can be recovered from boh densiy and hazard, bu he formula in he laer case now seems quie differen: S() = (1 α(s)). The coninuous case formula S() = exp( α(s)ds) has herefore wo major defecs: firsly, i does no have any inuiive inerpreaion; and secondly; i gives he wrong generalisaion o he discree case. Here is how boh formulas can be unified and made inuiively inerpreable. Define he cumulaive hazard A() by, in he coninuous case, A() = α(s)ds, and in he discree case, A() = α(s). (These wo formulas are special cases of he compleely general expression A() = ds(s)/s(s ).) Now we can wrie, boh in he coninuous and he discree case, S() = (1 da(s)) (1) which can be inerpreed as he produc over many small ime inervals s, s + ds making up he inerval [,], of he probabiliy (1 da(s)). Since he hazard da(s) can be hough of as he probabiliy of dying in he inerval from s o s+ds given survival up o he beginning of ha ime inerval, 1 minus he hazard is he probabiliy of surviving hrough he small ime inerval given survival up o is sar. Muliplying over he small ime inervals making up [,] yields he uncondiional probabiliy of surviving pas ; in oher words, equaion (1) is jus he limiing form of he equaliy Pr(T > ) = k Pr(T > i T > i 1 )= i=1 k ( ) 1 Pr(T i T > i 1 ) where = < 1 < < k = is a pariion of he ime inerval [,]. Consider now he saisical problem of esimaing he survival curve S() given a sample of independenly censored survival imes. Le 1 < 2 <... i=1 2
3 denoe he disinc imes when deahs are observed; le r j denoe he number of individuals a risk jus before ime j and le d j denoe he number of observed deahs a ime j. We esimae he cumulaive hazard funcion A corresponding o S wih he Nelson-Aalen esimaor Â() = j This is a discree cumulaive hazard funcion, corresponding o he discree esimaed hazard α( j )=d j /r j, α() zero for no an observed deah ime. The produc-inegral of  is hen d j r j. Ŝ() = (1 dâ) = (1 d j ), r j j which is nohing else han he Kaplan-Meier esimaor. The acual definiion of he produc-inegral in (1) is he following: (1 da(s)) = lim (1 (A(i ) A( i 1 ))) max i i 1 where he limi is aken over a sequence of ever finer pariions = < 1 < < k = of he ime inerval [,]. From his poin we can choose eiher o sudy properies of he producinegral or define i in greaer generaliy. Boh aspecs are imporan in applicaions. Le us firs give a more general definiion. The imporan generalisaion is ha we will define produc-inegrals of marix-valued funcions, raher han jus scalar valued funcions. The concep now really comes ino is own, because when we muliply a sequence of marices ogeher he resul will generally depend on he order in which he marices are aken. Even in he coninuous case here will no be a simple exponenial formula expressing he resul in erms of an ordinary inegral. Muliplying producs of marices urns up in he heory of Markov processes, and his connecs direcly o he saisical analysis of muli-sae models in survival analysis. Suppose X() is a p p marix-valued funcion of ime. Suppose also X (or if you like, each componen of X) is righ coninuous wih lef hand limis. Le I denoe he ideniy marix. The produc-inegral of X over he inerval [,] is now defined as (I + dx(s)) = lim (I +(X(i ) X( i 1 ))) max i i 1 where as always he limi is aken over a sequence of ever finer pariions = < 1 < < k = of he ime inerval [,]. For he limi o exis, X has o be of bounded variaion; equivalenly, each componen of X is he difference of wo increasing funcions. 3
4 2 Applicaion o Markov processes We briefly skech he applicaion of produc-inegraion o Markov processes. Suppose an individual moves beween p differen saes as ime proceeds, saying in each sae for some random lengh of ime and hen jumping o anoher. Suppose he individual has inensiy α ij () of jumping from sae i o sae j a ime, given he whole pas hisory (in oher words, he process is Markov: he inensiy only depends on he presen ime and he presen sae). Define cumulaive inensiies A ij () = α ij(s)ds and negaive oal cumulaive inensiy of leaving a sae A ii = j i A ij. Collec hese ino a square marix valued funcion of ime A. Then one can show ha he marix of ransiion probabiliies P (,), whose ij componen is he probabiliy of being in sae j a ime given he individual sared a ime in sae i, is given by a produc-inegral of A: P (,)= (I + da(s)). This formula generalises he usual formula for ransiion probabiliies of a discree ime Markov chain, since he marix (I + da(s)) can be hough of as he ransiion probabiliy marix for he small ime inerval s, s + ds. Given possibly censored observaions from a Markov process, one can esimae he elemens of he marix of cumulaive inensiies A by Nelson-Aalen esimaors. The corresponding esimae of he ransiion probabiliies, he Aalen- Johansen esimaor, is found by aking he produc-inegral of Â. 3 Mahemaical properies A very obvious propery of produc-inegraion is is muliplicaiviy. Defining he produc-inegral over an arbirary ime inerval in he naural way, we have for < s < s (I + dx) = (I + dx) (I + dx). s We can guess many oher useful properies of produc-inegrals by looking a various simple ideniies for finie producs. For insance, i is ofen imporan o sudy he difference beween wo produc-inegrals. Now if a 1,...,a k and b 1,...,b k are wo sequences of numbers, we have he ideniy: (1 + ai ) (1 + b i )= + a i )(a j b j ) j i<j(1 (1 + b i ). i>j This can be easily proved by replacing he middle erm on he righ, (a j b j ), by (1 + a i ) (1 + b i ). Expanding abou his difference, he righ hand side becomes ( (1 + a i ) (1 + b i ) (1 + a i ) (1 + b i )). j i j i>j i>j 1 4
5 This is a elescoping sum; wriing ou he erms one by one he whole expression collapses o he wo ouside producs, giving he lef hand side of he ideniy. The same manipulaions work for marices. In general i is herefore no surprise, replacing sums by inegrals and producs by produc-inegrals, ha (I + dx) (I + dy )= s= s (I + dx)(dx(s) dy (s)) (I + dy ). s+ This valuable ideniy is called he Duhamel equaion. As an example, consider he scalar case, le A be a cumulaive hazard funcion and  he Nelson-Aalen esimaor based on a sample of censored survival imes. Le S be he corresponding survival fucnion and Ŝ he Kaplan-Meier esimaor. The Duhamel equaion hen becomes he ideniy Ŝ() S() = s= Ŝ(s )(dâ(s) da(s))(s()/s(s)) which can be exploied o ge boh small sample and asympoic resuls for he Kaplan-Meier esimaor. We illusrae one oher imporan ideniy in a similar manner. Noe ha (1 + a i ) (1 + a i )= (1 + a i )a j. i j Adding over j from 1 o k gives us i k(1 + a i ) 1= j (1 + a i )a j. Now we can guess he ideniy (I + dx) I = s= s (I + dx)dx(s). This is essenially Kolmogorov s forward equaion from he heory of Markov processes, and i is he ype of equaion solve Y () =I + Y (s )dx(s) for unknown Y, given X which originally moivaed Volerra o inven producinegraion. Y () = (I + dx) is he unique soluion of his equaion. (I is also jus a special case of he Duhamel equaion when we ake he second inegrand Y idenically equal o zero). 4 Concluding remarks The produc-inegral seems firs o have been used as a fundamenal ool in modern survival analysis by Aalen and Johansen (1978), hough i also appears in a more informal conex in Cox (1972) and in Kalbfleisch and Prenice 5
6 (198). Surveys of he heory of produc-inegraion are given by Gill and Johansen (199), Gill (1994). The former paper also pays aenion o he earlier hisory of he subjec. In paricular, i is worh menioning ha a large variey of noaions has been used for he produc-inegral, including large curly P s, produc-symbols, and he ordinary inegral sign embellished wih a half circle over he op. As well as playing a role in he heory of he Kaplan-Meier and he Aalen- Johansen esimaors, he produc-inegral is also a useful way o wrie likelihoods and parial-likelihoods in survival analysis, since hese can be usefully hough of as coninuous producs of condiional likelihoods for he daa in each new infiniesimal ime inerval given he pas. The produc-inegral is also useful in mulivariae survival analysis. In paricular, Dabrowska s 1988 mulivariae produc-limi esimaor is based on a represenaion of a mulivariae survival funcion in erms of produc-inegrals of a collecion of higher-dimensional join and condiional hazard funcions. The book Andersen e al. (1993) gives a brief survey of he heory and many deailed applicaions, covering all he opics menioned above. References Aalen, O. and S. Johansen (1978). An empirical ransiion marix for nonhomogenous Markov chains based on censored observaions. Scandinavian Journal of Saisics 5, Andersen, P., Ø. Borgan, R. Gill, and N. Keiding (1993). Saisical Models Based on Couning Processes. New York: Springer-Verlag. Cox, D. (1972). Regression models and life ables (wih discussion). Journal of he Royal Saisical Sociey, Series B 34, Dabrowska, D. (1988). Kaplan-meier esimae on he plane. Annals of Saisics 16, Gill, R. (1994). Lecures on survival analysis. In P. Bernard (Ed.), Lecures on Probabiliy Theory (Ecole d Éé de Probabiliés de Sain Flour XXII ), Berlin, pp Springer-Verlag (SLNM 1581). Gill, R. and S. Johansen (199). A survey of produc-inegraion wih a view owards applicaion in survival analysis. Annals of Saisics 18, Kalbfleisch, J. and R. Prenice (198). The Saisical Analysis of Failure Time Daa. New York: Wiley. 6
Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationChapter 2. Models, Censoring, and Likelihood for Failure-Time Data
Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationLecture Notes 2. The Hilbert Space Approach to Time Series
Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship
More informationPhysics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution
Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationLecture 10: The Poincaré Inequality in Euclidean space
Deparmens of Mahemaics Monana Sae Universiy Fall 215 Prof. Kevin Wildrick n inroducion o non-smooh analysis and geomery Lecure 1: The Poincaré Inequaliy in Euclidean space 1. Wha is he Poincaré inequaliy?
More informationLECTURE 1: GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS
LECTURE : GENERALIZED RAY KNIGHT THEOREM FOR FINITE MARKOV CHAINS We will work wih a coninuous ime reversible Markov chain X on a finie conneced sae space, wih generaor Lf(x = y q x,yf(y. (Recall ha q
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
More informationStochastic Model for Cancer Cell Growth through Single Forward Mutation
Journal of Modern Applied Saisical Mehods Volume 16 Issue 1 Aricle 31 5-1-2017 Sochasic Model for Cancer Cell Growh hrough Single Forward Muaion Jayabharahiraj Jayabalan Pondicherry Universiy, jayabharahi8@gmail.com
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationComparison between the Discrete and Continuous Time Models
Comparison beween e Discree and Coninuous Time Models D. Sulsky June 21, 2012 1 Discree o Coninuous Recall e discree ime model Î = AIS Ŝ = S Î. Tese equaions ell us ow e populaion canges from one day o
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationGMM - Generalized Method of Moments
GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................
More information) were both constant and we brought them from under the integral.
YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationThe Arcsine Distribution
The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationInternational Journal of Scientific & Engineering Research, Volume 4, Issue 10, October ISSN
Inernaional Journal of Scienific & Engineering Research, Volume 4, Issue 10, Ocober-2013 900 FUZZY MEAN RESIDUAL LIFE ORDERING OF FUZZY RANDOM VARIABLES J. EARNEST LAZARUS PIRIYAKUMAR 1, A. YAMUNA 2 1.
More information23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes
Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationEchocardiography Project and Finite Fourier Series
Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every
More informationFinish reading Chapter 2 of Spivak, rereading earlier sections as necessary. handout and fill in some missing details!
MAT 257, Handou 6: Ocober 7-2, 20. I. Assignmen. Finish reading Chaper 2 of Spiva, rereading earlier secions as necessary. handou and fill in some missing deails! II. Higher derivaives. Also, read his
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More information6. Stochastic calculus with jump processes
A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio
More informationWavelet Methods for Time Series Analysis. What is a Wavelet? Part I: Introduction to Wavelets and Wavelet Transforms. sines & cosines are big waves
Wavele Mehods for Time Series Analysis Par I: Inroducion o Waveles and Wavele Transforms waveles are analysis ools for ime series and images as a subjec, waveles are relaively new (983 o presen) a synhesis
More informationACE 562 Fall Lecture 4: Simple Linear Regression Model: Specification and Estimation. by Professor Scott H. Irwin
ACE 56 Fall 005 Lecure 4: Simple Linear Regression Model: Specificaion and Esimaion by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Simple Regression: Economic and Saisical Model
More informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More informationDISCRETE GRONWALL LEMMA AND APPLICATIONS
DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE MAA NORTH CENTRAL SECTION MEETING AT UND 24 OCTOBER 29 Gronwall s lemma saes an inequaliy ha is useful in he heory of differenial equaions. Here is
More informationKEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationThe Strong Law of Large Numbers
Lecure 9 The Srong Law of Large Numbers Reading: Grimme-Sirzaker 7.2; David Williams Probabiliy wih Maringales 7.2 Furher reading: Grimme-Sirzaker 7.1, 7.3-7.5 Wih he Convergence Theorem (Theorem 54) and
More informationLecture 33: November 29
36-705: Inermediae Saisics Fall 2017 Lecurer: Siva Balakrishnan Lecure 33: November 29 Today we will coninue discussing he boosrap, and hen ry o undersand why i works in a simple case. In he las lecure
More informationInventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions
Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationInstructor: Barry McQuarrie Page 1 of 5
Procedure for Solving radical equaions 1. Algebraically isolae one radical by iself on one side of equal sign. 2. Raise each side of he equaion o an appropriae power o remove he radical. 3. Simplify. 4.
More informationCOALESCENT RANDOM WALKS ON GRAPHS
COALESCENT RANDOM WALKS ON GRAPHS JIANJUN PAUL TIAN AND ZHENQIU LIU Absrac. Inspired by coalescen heory in biology, we inroduce a sochasic model called muli-person simple random walks or coalescen random
More informationarxiv: v1 [math.ca] 15 Nov 2016
arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy
More informationMODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE
Topics MODULE 3 FUNCTION OF A RANDOM VARIABLE AND ITS DISTRIBUTION LECTURES 2-6 3. FUNCTION OF A RANDOM VARIABLE 3.2 PROBABILITY DISTRIBUTION OF A FUNCTION OF A RANDOM VARIABLE 3.3 EXPECTATION AND MOMENTS
More information= ( ) ) or a system of differential equations with continuous parametrization (T = R
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationBasic definitions and relations
Basic definiions and relaions Lecurer: Dmiri A. Molchanov E-mail: molchan@cs.u.fi hp://www.cs.u.fi/kurssi/tlt-2716/ Kendall s noaion for queuing sysems: Arrival processes; Service ime disribuions; Examples.
More informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
More information2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes
Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion
More informationin Engineering Prof. Dr. Michael Havbro Faber ETH Zurich, Switzerland Swiss Federal Institute of Technology
Risk and Saey in Engineering Pro. Dr. Michael Havbro Faber ETH Zurich, Swizerland Conens o Today's Lecure Inroducion o ime varian reliabiliy analysis The Poisson process The ormal process Assessmen o he
More information13.3 Term structure models
13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More information( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
More informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
More informationES.1803 Topic 22 Notes Jeremy Orloff
ES.83 Topic Noes Jeremy Orloff Fourier series inroducion: coninued. Goals. Be able o compue he Fourier coefficiens of even or odd periodic funcion using he simplified formulas.. Be able o wrie and graph
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationBernoulli numbers. Francesco Chiatti, Matteo Pintonello. December 5, 2016
UNIVERSITÁ DEGLI STUDI DI PADOVA, DIPARTIMENTO DI MATEMATICA TULLIO LEVI-CIVITA Bernoulli numbers Francesco Chiai, Maeo Pinonello December 5, 206 During las lessons we have proved he Las Ferma Theorem
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationEffects of Coordinate Curvature on Integration
Effecs of Coordinae Curvaure on Inegraion Chrisopher A. Lafore clafore@gmail.com Absrac In his paper, he inegraion of a funcion over a curved manifold is examined in he case where he curvaure of he manifold
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
More informationCash Flow Valuation Mode Lin Discrete Time
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728,p-ISSN: 2319-765X, 6, Issue 6 (May. - Jun. 2013), PP 35-41 Cash Flow Valuaion Mode Lin Discree Time Olayiwola. M. A. and Oni, N. O. Deparmen of Mahemaics
More information20. Applications of the Genetic-Drift Model
0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0
More informationRANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY
ECO 504 Spring 2006 Chris Sims RANDOM LAGRANGE MULTIPLIERS AND TRANSVERSALITY 1. INTRODUCTION Lagrange muliplier mehods are sandard fare in elemenary calculus courses, and hey play a cenral role in economic
More informationMATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence
MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,
More informationOn Multicomponent System Reliability with Microshocks - Microdamages Type of Components Interaction
On Mulicomponen Sysem Reliabiliy wih Microshocks - Microdamages Type of Componens Ineracion Jerzy K. Filus, and Lidia Z. Filus Absrac Consider a wo componen parallel sysem. The defined new sochasic dependences
More informationA Uniform Asymptotic Estimate for Discounted Aggregate Claims with Subexponential Tails
A Uniform Asympoic Esimae for Discouned Aggregae Claims wih Subeponenial Tails Xuemiao Hao and Qihe Tang Deparmen of Saisics and Acuarial Science The Universiy of Iowa 241 Schae er Hall, Iowa Ciy, IA 52242,
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More informationExponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits
DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,
More informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
More informationSemi-Competing Risks on A Trivariate Weibull Survival Model
Semi-Compeing Risks on A Trivariae Weibull Survival Model Jenq-Daw Lee Graduae Insiue of Poliical Economy Naional Cheng Kung Universiy Tainan Taiwan 70101 ROC Cheng K. Lee Loss Forecasing Home Loans &
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More informationGeneralized Chebyshev polynomials
Generalized Chebyshev polynomials Clemene Cesarano Faculy of Engineering, Inernaional Telemaic Universiy UNINETTUNO Corso Viorio Emanuele II, 39 86 Roma, Ialy email: c.cesarano@unineunouniversiy.ne ABSTRACT
More informationFrom Complex Fourier Series to Fourier Transforms
Topic From Complex Fourier Series o Fourier Transforms. Inroducion In he previous lecure you saw ha complex Fourier Series and is coeciens were dened by as f ( = n= C ne in! where C n = T T = T = f (e
More informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
More informationBiol. 356 Lab 8. Mortality, Recruitment, and Migration Rates
Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationThis document was generated at 7:34 PM, 07/27/09 Copyright 2009 Richard T. Woodward
his documen was generaed a 7:34 PM, 07/27/09 Copyrigh 2009 Richard. Woodward 15. Bang-bang and mos rapid approach problems AGEC 637 - Summer 2009 here are some problems for which he opimal pah does no
More informationIMPLICIT AND INVERSE FUNCTION THEOREMS PAUL SCHRIMPF 1 OCTOBER 25, 2013
IMPLICI AND INVERSE FUNCION HEOREMS PAUL SCHRIMPF 1 OCOBER 25, 213 UNIVERSIY OF BRIISH COLUMBIA ECONOMICS 526 We have exensively sudied how o solve sysems of linear equaions. We know how o check wheher
More informationA Uniform Asymptotic Estimate for Discounted Aggregate Claims with Subexponential Tails
A Uniform Asympoic Esimae for Discouned Aggregae Claims wih Subeponenial Tails Xuemiao Hao and Qihe Tang Deparmen of Saisics and Acuarial Science The Universiy of Iowa 241 Schaeffer Hall, Iowa Ciy, IA
More information4.6 One Dimensional Kinematics and Integration
4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More information5. Stochastic processes (1)
Lec05.pp S-38.45 - Inroducion o Teleraffic Theory Spring 2005 Conens Basic conceps Poisson process 2 Sochasic processes () Consider some quaniy in a eleraffic (or any) sysem I ypically evolves in ime randomly
More informationACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.
ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple
More informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationWeyl sequences: Asymptotic distributions of the partition lengths
ACTA ARITHMETICA LXXXVIII.4 (999 Weyl sequences: Asympoic disribuions of he pariion lenghs by Anaoly Zhigljavsky (Cardiff and Iskander Aliev (Warszawa. Inroducion: Saemen of he problem and formulaion of
More informationGENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT
Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received
More informationEE 315 Notes. Gürdal Arslan CLASS 1. (Sections ) What is a signal?
EE 35 Noes Gürdal Arslan CLASS (Secions.-.2) Wha is a signal? In his class, a signal is some funcion of ime and i represens how some physical quaniy changes over some window of ime. Examples: velociy of
More informationInequality measures for intersecting Lorenz curves: an alternative weak ordering
h Inernaional Scienific Conference Financial managemen of Firms and Financial Insiuions Osrava VŠB-TU of Osrava, Faculy of Economics, Deparmen of Finance 7 h 8 h Sepember 25 Absrac Inequaliy measures for
More informationarxiv:math/ v1 [math.nt] 3 Nov 2005
arxiv:mah/0511092v1 [mah.nt] 3 Nov 2005 A NOTE ON S AND THE ZEROS OF THE RIEMANN ZETA-FUNCTION D. A. GOLDSTON AND S. M. GONEK Absrac. Le πs denoe he argumen of he Riemann zea-funcion a he poin 1 + i. Assuming
More informationAverage Number of Lattice Points in a Disk
Average Number of Laice Poins in a Disk Sujay Jayakar Rober S. Sricharz Absrac The difference beween he number of laice poins in a disk of radius /π and he area of he disk /4π is equal o he error in he
More informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More information