Uncertainty Specification and Propagation for Loss Estimation Using FOSM Methods

Size: px
Start display at page:

Download "Uncertainty Specification and Propagation for Loss Estimation Using FOSM Methods"

Transcription

1 Uncertanty Specfcaton and Propagaton for Loss Estmaton Usng FOSM Methods J.W. Baer and C.A. Corne Dept. of Cv and Envronmenta Engneerng, Stanford Unversty, Stanford, CA Keywords: Sesmc, oss estmaton, correaton, Frst-Order Second-Moment ABSTRACT: Probabstc predcton of structura and nonstructura damage costs due to future earthquaes s one component of oss estmaton currenty beng deveoped for use n performance-based earthquae engneerng. Sources of uncertanty n ths predcton ncude epstemc and aeatory uncertanty n the ste ground moton hazard, the budng response, the damage measures of each of the many budng eements, and repar cost of each of the eements. These are nter- and cross-correated random varabes. Two desred resuts are the tota uncertanty n annua osses, and the contrbuton of each uncertanty source to the tota uncertanty. Monte Caro smuaton s a smpe souton, but t can be computatonay expensve. Ths study proposes an aternatve approach usng Frst-Order Second-Moment (FOSM) methods for a but the (domnant) ground moton ntensty varabe. Suggestons for characterzaton of correatons are presented. A procedure for appyng FOSM methods n the cacuaton of tota uncertanty s outned. The proposed technque s very effcent, and easy used for senstvty studes. 1 INTRODUCTION Estmaton of annua osses n a budng due to earthquae damage s a quantty of nterest to decson maers, and s a current topc of study n performance-based earthquae engneerng. Among the quanttes to be determned are the uncertanty n the resut, and the contrbuton of each source of uncertanty to the tota uncertanty. The Pacfc Earthquae Engneerng Research Center (PEER) has proposed the foowng framng equaton for ths anayss: TC ( TC DVE DVE DM DM EDP EDP u v y x z ) f ( z, u) f ( u, v) f ( v, y) f ( y, x) d ( x) (1) wth terms defned n Secton 1.1. Ths equaton aows for moduar consderaton of the ground moton hazard, budng response, damage to budng eements, eement repar costs, and tota repar cost (Corne and Krawner, 000, Krawner, 00). One opton for cacuatng uncertanty n the resut s through Monte Caro smuaton (Porter, 00). Athough straghtforward, t can be very expensve computatonay, especay when mutpe runs are requred to cacuate senstvtes. The obectve of ths study s to propose an aternatve method of cacuatng uncertanty usng the Frst- Order Second-Moment (FOSM) method (e.g. Mechers, 1999). We sha use ths approxmate method to coapse out severa of the ntermedate condtona random varabes, eavng a mean and varance of Tota Cost (TC) condtoned on the Intensty Measure () of the ground moton. Ths nformaton can then be combned wth the ground moton hazard, d (x), to obtan the expected annua tota cost, varance n annua tota cost, and the annua rate of exceedng a gven tota cost. 1.1 Expanaton of the Framng Equaton The varabes n Equaton 1 are defned as foows: TC (z) s the annua rate of exceedng a tota repar cost of z, where tota repar cost, TC s the decson varabe under study. f TC DVE (z,u) s the PDF of TC, condtoned on the vector of damage vaues of each eement (DVE s the damage vaue of eement ). The assumpton n the framewor descrbed beow s that the tota cost of repar s the sum of a eement repar costs, but ths can be easy generazed. f DVE DM (u,v) s the PDF of the vector of damage vaues of each eement, gven the vector of damage states of each eement (DM s the damage state of

2 eement ). Mean repar costs can be estmated from sources such as R.S. Means Co. s pubshed materas on constructon cost estmatng (R.S. Means Co. 00). Addtona quantfcaton of repar costs s a topc of current research. f DM EDP (v,y) s the PDF of the vector of damage states, gven the vector of engneerng demand parameters. In current research, these damage states are typcay dscrete, and each state s descrbed by a fragty functon, whch returns the probabty of an eement exceedng the damage state at a gven EDP eve. See Asan and Mranda (00) for exampes. f EDP (y,x) s the PDF of the vector of engneerng demand parameters, gven the ntensty measure. For aeatory uncertanty, ths dstrbuton can be determned usng, for exampe, Incrementa Dynamc Anayss (Vamvatsos and Corne 00). d (x) s the absoute vaue of the dervatve of the annua rate of exceedng a gven vaue of the ntensty measure (the sesmc hazard curve). The absoute vaue s needed because the dervatve s negatve. See, for exampe, Kramer (1995) for bacground on hazard curves. ASSUMPTIONS A Marovan dependence s assumed for a dstrbutons n the framewor. For exampe, t s assumed that the dstrbuton of the DM vector can be condtoned soey on the EDP vector, and that nowedge of the provdes no addtona nformaton. In ths way, prevous condtonng nformaton does not need to be carred forward through a future dstrbutons, reducng compexty. A condtonng varabe that contans a necessary condtona nformaton s deemed a suffcent descrptor (Luco 00). A damage s assumed to occur on an eement eve. The tota cost of damage to the structure s then the sum of the damage cost of each eement n the structure. The excepton to ths assumpton s when coapse occurs, and repar costs w be a functon of the coapse, rather than ndvdua eement responses. The treatment of ths excepton s expaned beow. A reatons n the framewor are assumed to be scaar functons. For exampe, the condtona dstrbuton of the Damage Measure of eement s a functon of ony the th Engneerng Demand Parameter. Or aternatvey, f DM EDP ( v, y) f DM EDP ( v, y ) () Note aso that the functon s not condtoned on varabes from any prevous steps, because of the Marovan process assumpton descrbed earer. To cacuate tota uncertanty n our decson varabe, t w be necessary to account for both epstemc and aeatory uncertanty. The framewor outned here s approprate for ether source of uncertanty. These two sources of uncertanty are uncorreated, aowng ther contrbutons to be cacuated separatey for smpcty, and towards the end of the procedure. Ths s further dscussed n Secton 3.6. These assumptons are beeved to be consstent wth the most advanced current sesmc oss estmaton efforts. Most can be reaxed wthout forma dffcuty. 3 PROCEDURE The procedure outned maes use of FOSM approxmatons to cacuate the mean and varance of TC gven. Ths nformaton can then be combned wth the ground moton hazard, d (x), to obtan the expected annua tota cost, varance n annua tota cost, and (together wth a dstrbuton type assumpton), the mean annua rate of exceedng a gven tota cost. For ths fna combnaton wth the ground moton hazard, FOSM approxmatons are not used. The FOSM approxmatons are ustfed by the assumpton that the uncertanty n the hazard curve s the most sgnfcant contrbutor to varance of the tota oss. Therefore, we are retanng the fu dstrbuton for tsef, but usng the FOSM approxmatons for a moments condtoned on. In addton, we ey do not have nformaton about the fu dstrbutons of some varabes (for exampe, repar costs), and so negectng hgher moments of these dstrbutons does not resut n a sgnfcant oss of avaabe nformaton. Note that we are worng wth natura ogarthms of the varabes descrbed prevousy. Ths aows us to wor wth sums of terms, rather than products. We revert to a non-og form for the fna resut. The procedure s outned n the foowng sectons. 3.1 Specfy n EDP The proposed mode n ths study s EDP =H () (), where H () s the (determnstc) mean vaue of EDP gven, and () s a random varabe wth mean of one, and condtona varance adusted to mode the varance n EDP. (We ntroduce the random varabe notaton X Y, to denote that the mode of X s condtoned on Y.) Then when we use the og form of EDP, we have a random varabe of the form n(edp )=n(h ())+n( ()). Note that the expected vaue of n(edp ) s n(h ()), and the varance of n(edp ) s equa to the varance of n( ()). Both n(h ()) and n( ())], as we as the correatons between nedp s, can be determned from Incrementa Dynamc Anayss.

3 We w need the foowng nformaton for our cacuatons: n EDP ], denoted h () for a EDP (3) Fgure 1: Exampe Eement Fragty Functons n EDP ], denoted h* () for a EDP (4) ρ(n EDP, n EDP ), denoted ĥ () for a {EDP, EDP } (5) These functons w be used n Secton 3.3 beow. 3. Specfy DM n EDP and n DVE DM, and coapse to n DVE n EDP The dscrete states of the Damage Measure varabe found n current oss estmaton (Asan and Mranda, 00, Porter, 00) are not compatbe wth the FOSM approach, whch requres contnuous functons for the moments. To dea wth the dscrete states, we tae advantage of the fact that we can aways coapse the two dstrbutons DM n EDP and n DVE DM nto one contnuous dstrbuton n DVE n EDP by ntegratng over the approprate varabe: f DVE EDP ( u, y) fdve DM ( u, v) fdm EDP ( v, y) (6) v For a gven eement wth n possbe damage states, we use a set of eement fragty functons F 1, F F n, such that F (y)=p(dm>d EDP=y) (see Fgure 1). We aso defne F 0 1 (the probabty that each eement has at east zero damage s one). These functons w have a correspondng set of dstrbutons c 1, c c n of eement repar costs such that c (v) s a probabty dstrbuton of DVE, gven that the damage state equas d (see Fgure ). Wth ths nformaton, we can determne the frst two moments of the coapsed dstrbutons. For exampe, Asan and Mranda (00) document the deveopment of one set of these functons. Fgure : Eement Repar Costs From the tota probabty theorem, we now that n ths case, Equaton 6 can be wrtten n scaar form for each DVE as: f (7) f P DVE EDP DVE DM d DM d EDP DamageStates (reca our assumpton that each DVE s dependent on a snge EDP). For our FOSM purposes, furthermore, t s suffcent to fnd smpy the condtona means, varances, and covarances of the DVE s gven the EDP s. Thus, tang the mean of ths PDF, we have the resut: [ DVE EDP] F ( EDP) F 1 Damage States E ( EDP) (8) Appyng the same thnng to E [ DVE EDP], and recognzng that X]=X ]- X, we have the foowng resut: DVE EDP] DVE Damage States ( Damage States (( F ( EDP) F ) 1 ( F ( EDP) F ( EDP)) 1 ( EDP)) (9)

4 Fgure 3: Coapsed dstrbuton DVE EDP Fgure 3 shows an exampe of the mean and mean pus or mnus one sgma, as generated from the exampe dstrbutons shown n Fgure 1 and Fgure Quantfyng Correatons We now need to determne correatons among the DVE s of a eements n the structure. Note that e the mean and varance, these correatons are condtoned on the EDP s. Whe these cacuatons are straghtforward, estmaton of the necessary correaton nputs s a dffcut tas due to a ac of data. In the absence of addtona nformaton, t may be hepfu to use the foowng characterzaton scheme. Let us assume for ths purpose a mode of the form: n DVE n EDP g (n EDP ) n n n (10) m E where represents uncertanty common to the entre structure, represents uncertanty m common ony to eements of cass m (e.g. drywa parttons, moment connectons, etc.), and E represents uncertanty unque to eement. A of these s are assumed to be mutuay uncorreated. We then defne: n n ] (11) EDP n n ] for a m (1) m EDP E n n ] for a (13) E EDP Then the varance of ndve nedp s the sum of these varances. For ths speca case, a smpe cosed form souton exsts for the correaton coeffcent between two eement DVE s. If the two eements are of the same eement cass, then: (n DVE,n DVE n EDP,n EDP ) (14) If the eements are of dfferent eement casses, then ther correaton coeffcent s gven by: (n DVE,n DVE n EDP,n EDP ) E (15) Note that ths formuaton requres to be equa for a eement casses, and E to be equa for a eements. If ths s excessvey mtng, a cosed form souton aso exsts that aows to vary by cass, and E to vary by eement (Baer and Corne, 00). The prncpes used for ths souton are deveoped n Dtevsen (1981). The use of more than two uncertan terms, and the use of terms that vary by cass or eement are both generazatons of the basc equcorreated mode. Thus, we w refer to a mode ncorporatng any of these generazatons as a generazed equcorreated mode. The correaton matrx for a generazed equcorreated mode w have off-dagona terms that vary from term to term, as opposed to the strct equcorreated mode, where a off-dagona terms are dentca. We have now concuded the coapse of the dstrbuton ndve nedp. We have the condtona mean and varance functons of ndve nedp, obtaned by coapsng the dstrbutons provded (see Equatons 8 and 9), and correaton coeffcents determned usng the generazed equcorreated mode (see Equatons 14 and 15). We choose for future notatona carty to denote these resuts as: n DVE n EDP ], denoted g (n EDP ) for a DVE (16) n DVE n EDP ], denoted g* (n EDP ) for a DVE (17) ρ(n DVE, n DVE n EDP, n EDP ), denoted ĝ (n EDP, n EDP ) for a {DVE, DVE } (18) Wth ths nformaton quantfed, we can now use t aong wth the resuts from Secton 3.1 to cacuate n DVE. 3.3 Cacuate n DVE Usng nformaton from above, we can cacuate the frst and second moments of n DVE. Ths nvoves coapsng out the dependence on EDP, as suggested n Equaton 19 beow. f DVE E ( u, x) fdve EDP ( u, y) fedp ( y, x) (19) y To mantan tractabty, we sha use an FOSM approxmaton here. To remove dependence on EDP, we tae the expectaton of n DVE of wth respect to n EDP (gven ). We wrte ths as ndve ]=E EDP [ndve nedp ]], where E EDP [] denotes ths partcuar condtona expec

5 taton operator. Substtutng our notaton from Equatons 16 and 17, we have: n DVE ] g h ( ) (0) Usng a smar approach to condtona moments, and usng the resut from probabty theory: [ X Z] Var X ] X ] E Var Z (1) we can aso derve the varance and covarances of n DVE wth the usua FOSM approxmatons: g n DVE ] g * ( h ( )) n EDP h ( ) h * ( ) () Cov[n DVE ˆ, DVE ] g( h ( ), h ( )) g * ( h ( )) g * ( h ( )) (3) g g hˆ ( ) n EDP n EDP h ( ) h ( ) h * ( ) h * ( ) 3.4 Swtch to the non-og form DVE To swtch to the non-og form of DVE, we can use the frst-order approxmaton e X ] e X]. Then we have the foowng resuts: g ( h ( )) DVE ] e (4) DVE Cov[ DVE, DVE g ( h ( )) ] e n DVE ] (5) g ( h ( )) g ( h ( )) ] e Cov[n DVE,n DVE ] (6) These eement resuts can now be used to compute the moments of TC. 3.5 Compute moments of TC Under the assumpton that Tota Cost s the sum of eement costs, we can now aggregate the resuts from a ndvdua eements to compute an expectaton and varance for the tota cost of damage to the entre budng. The expected tota cost s the sum of expected eement costs: # eements TC ] DVE ] (7) 1 We denote the expected vaue computed n Equaton 7 as q(). The varance of tota cost can s the sum of eement varances, ncudng covarances between eement costs: TC ] # eements 1 DVE # eements # eements 1 1 ] Cov[ DVE, DVE ] (8) We denote the expected vaue computed n Equaton 8 as q*(). 3.6 Repeat Procedure to Cacuate Epstemc Uncertanty We assume a mode of the form TC =q(), R U where q() s the best estmate of the condtona mean as cacuated n Equaton 7, and R and U are uncorreated random varabes representng aeatory and epstemc uncertanty, respectvey. Then ntc n q( ) n n (9) Because R and U are uncorreated, we may dea wth them n separate steps. The above procedure usng aeatory uncertanty aone aowed us to fnd the varance due to R. We must now repeat the procedure to cacuate the varance due to U. The tota uncertanty can then be cacuated by combnng the two uncertantes as foows: ntc ] n ] n ] (30) We denote ths TC. Note that we have swtched to ogs agan to aow use of sums rather than products. The change can be made usng the foowng reatonshp: Var R[ TC ] n ] n 1 R (31) ER[ TC ] We denote n R ] as R. From ths vaue, we can then cacuate Var [ntc ] usng Equaton 30. Once we have cacuated Var [ntc ], accountng for both aeatory and epstemc uncertanty, we can denote ths vaue as q*(), and use t n the equatons to foow. When repeatng the procedure to fnd the varance n TC gven due to epstemc uncertanty, epstemc uncertantes for each condtona random varabe w need to be estmated, and characterzaton of correaton s agan a chaengng tas. It s suggested that the generazed equcorreated mode deveoped n Secton 3..1 may be used effectvey for ths probem. 4 ACCOUNTING FOR COLLAPSE CASES At hgh eves, the potenta exsts for a structure to experence coapse (defned here as extreme defectons at one or more story eves). In ths budng state, repar costs are more ey a functon of the coapse, rather than ndvdua eement damage. In fact, the structure s ey not to be repared at a. Thus, our predcted oss may not be accurate n these cases. In addton, the arge defectons predcted n a few cases w sew our expected vaues of some EDP s such as nterstory drfts, athough coapse s ony occurrng n a fracton of cases. To account for the possbty of coapse, we woud e to use the technque outned above for no-coapse R R U U

6 cases, and aow for an aternate oss estmate when coapse occurs. The foowng modfcaton s suggested. Note, n the foowng cacuatons, we are condtonng on a coapse ndcator varabe. To communcate ths, we have denoted the coapse and no coapse condton as C and NC respectvey. At each eve, compute the probabty of no coapse. Ths probabty, p ( NC ), s smpy the fracton of anayss runs where no coapse occurs Cacuate resuts usng the FOSM anayss as before, but usng ony the runs that resuted n no coapse. We now denote these resuts E [ TC, NC] and Var [ TC, NC]. Defne an expected vaue and varance of tota cost, gven that coapse has occurred, denoted E [ TC, C] and Var [ TC, C]. These vaues w ey not be functons of, but the condtonng on s st noted for consstency. The expected vaue of TC for a gven eve s now the average of the coapse and no coapse TC, weghted by ther respectve probabtes of occurrng: TC ] p( NC ) TC, NC] 1 p( NC ) TC, C] (3) The varance can be computed usng the property from Equaton 1: p( NC ) TC, NC] Var [ TC ] 1 p( NC ) TC, C] (33) p( NC ) TC ] TC, NC] 1 p( NC ) TC ] TC, C] The procedure can now be mpemented as before, usng these moments. Ths coapse-case modfcaton s probaby necessary for any mpementaton of the mode, as anayss of shang () eves suffcent to cause arge fnanca oss are ey aso to cause coapse n some representatve ground moton records. 5 INCORPORATE THE SITE HAZARD The expected vaue and varance of TC gven can now be ncorporated wth the ste hazard to compute the expected annua oss, and the rate of exceedng a gven Tota Cost. 5.1 Annua Loss Usng the functons q() and q*(), and the dervatve of the hazard curve, d ( ), the mean and varance of TC per annum can be cacuated by numerca ntegraton: E [ TC] q( ) d ( ) (34) TC] TC ] TC ]] q *( ) d( ) q ( ) d( ) TC] (35) Note that the frst term of Equaton 35 s the contrbuton from uncertanty n the cost functon gven, and the second two terms are the contrbuton from uncertanty n the. 5. Rate of Exceedance of a Gven TC The frst and second moment nformaton for TC can be combned wth a ste hazard to compute TC (z), the annua frequency of exceedng a gven Tota Cost z. For ths cacuaton, t s necessary to assume a probabty dstrbuton for TC that has condtona mean and varance equa to the vaues cacuated prevousy. The rate of exceedance of a gven TC s then gven by: TC ( z) F ( z, x) d ( x) (36) TC 5.3 Anaytc Souton Generay, the ntegra above w requre a numerca ntegraton. However, f the foowng smpfyng assumptons are made, an anaytc souton s avaabe: The dstrbuton of TC s ognorma TC s approxmated as constant for a ; we ca ths constant vaue * TC TC ] s approxmated by a functon of the form a b, where a and b are constants. Note that ths s consstent wth fttng the medan of TC wth a b, where 1 * TC a a' e (37) An approxmate hazard curve of the form (x)= 0 x - s ft to the true ste hazard curve Under these condtons, the annua rate of exceedng a gven Tota Cost s gven by: b z 1 TC ( z) 0 exp 1 * TC (38) a' b b We note that f the a from Equaton 37 s substtuted nto Equaton 38, then the resut becomes b z 1 TC ( z) 0 exp * TC (39) a b Ths equaton s usefu as an effcent estmate of TC (z), but t s aso very nformatve as a measure of the reatve mportance of uncertanty n the cacuaton. The term: b z 0 (40) a'

7 n the Equaton 38 woud be the resut f * TC were to equa zero that s f we made a cacuatons ony usng expected vaues and negected cost uncertanty gven. The term: 1 exp b 1 * TC (41) b n Equaton 38 s an ampfcaton factor that vares wth the uncertanty n TC present n the probem. Thus for ths speca case, t s smpe to cacuate the effect of uncertanty on the rate of exceedng a gven Tota Cost. As we sha show beow, t may not be unreasonabe for ths factor to ncrease TC (z) by a factor of 10, so the effect of uncertanty may very we be sgnfcant. However, even for arge vaues of * TC, the annua rate of exceedance s st domnated by the term from Equaton 40. It s for ths reason that t has been proposed here that the FOSM approxmatons of * TC performed above are suffcent to provde an accurate resut. For ustraton, et us assume that the expected TC as a functon of has been estmated usng the above technque as: TC ] 1 e (4) Ths functon coud be approxmated by the functon: 1.8 TC ] 1.4 (43) A pot of these two functons s shown n Fgure 4. Note that the anaytc functon s a good ft over the range 0<TC<0.5. Tota Cost Expected TC Gven Exampe vaue of TC ] Functon ft for Anaytca Souton Fgure 4: Expected TC, for ft TC ] 1 e and 5.4 Comparson of Resuts from Numerca Integraton and Anaytc Souton Both the numerca souton and the anaytca souton outned above n Sectons 5. and 5.3, respectvey, can be evauated for severa vaues of TC. The resuts can then be potted to generate a oss curve, as shown beow n Fgure 5. Ths fgure was generated usng the expected TC curves shown n Fgure 4 above. For both soutons, we have assumed a hazard curve of the form: ( x) 0 x (44) where 0 and are constants equa to 0.00 and 3, respectvey. We have aso assumed TC equa to 0.6 for both soutons. Mean Annua Frequency of Exceedance 1E-01 1E-0 1E-03 1E-04 Annua Frequency of Exceedance of a Gven Tota Cost Anaytca Mode Numerca Integraton 1E Tota Cost Fgure 5: λ TC (z): Comparson of numerca ntegraton and anaytc souton Fgure 5 aso aows us to compare the anaytca souton to the numerca one. For the functons gven n Equatons 4 and 43, the anaytc souton s a good approxmaton of the numerc souton over the TC range where Equaton 43 cosey ft Equaton 4 (0<TC<0.5). As we move to hgher TC eves, where the anaytca souton was not a good ft, the λ TC (z) resuts aso dverge. 6 EFFECT OF UNCERTAINTY ON LOSSES The varance n annua osses s the resut that most expcty shows the effects of uncertanty. However, uncertanty aso has an effect on the annua rate of exceedng a gven TC. Ths s most ceary seen n the anaytca souton of Equaton 38, where * TC appears n the equaton for TC (z). Ths s dscussed n Secton 5.3. The antcpated effect s for uncertanty to ncrease the rate of occurrence of a gven TC. However, dependng on the sopes of the hazard curve and mean of TC as a functon of, (defned by the parameters and b), ncreasng * TC can potentay decrease TC (z), or have no effect at a. Usng the more genera numerca ntegraton of Equaton 36, we fnd smar resuts, as shown n Fgure 6. In ths exampe, usng the functons assumed prevousy, we see that the shft n resuts s mnor for * TC <1, but for * TC =, the expected annua frequency of occurrence of arge costs has ncreased by approxmatey an order of magntude. Because our expected repar cost functon does not ever produce a oss greater than 1, the ncuson of

8 uncertanty s crtca for estmatng occurrence of tota costs greater than 1, as seen n the fgure. Mean Annua Frequency of Exceedance 1E-01 1E-0 1E-03 1E-04 Annua Frequency of Exceedance of a Gven Tota Cost β = β = 1 β = 0 1E Tota Cost Fgure 6: Effect of Uncertanty on Frequency of Exceedance of Tota Cost The mportant concuson to be drawn from ths resut s that usng expected vaues aone and gnorng uncertantes, athough temptng because of ts ease, can potentay ead to naccurate resuts. CONCLUSION A procedure for estmaton of uncertanty n repar costs due to earthquae damage has been proposed. Ths procedure wors wthn the framewor proposed by PEER for performance-based earthquae engneerng. Tota cost defned s a functon of repar costs for ndvdua budng eements, except n the coapse case, where a separate cost estmaton s used. Identfed aeatory and epstemc uncertanty n ground moton hazard, budng response, damage to budng eements and eement repar costs s combned to produce an uncertanty n tota repar cost. The proposed procedure uses the Frst-Order Second-Moment (FOSM) method to coapse severa condtona random varabes nto a snge random varabe. Numerca ntegraton s then used to ncorporate the ground moton hazard, where the uncertanty s most sgnfcant. The resutng nformaton s expected annua oss, varance n annua oss, and the annua rate of exceedng a gven cost. Baer, J.W. and Corne, C.A. (00). Technca Report on Uncertanty Specfcaton and Loss Estmaton Usng FOSM Methods. Report n Preparaton. Corne, C. A. and Krawner, H. (000). Progress and Chaenges n Sesmc Performance Assessment. Peer Center News, 3(). Dtevsen, Ove (1981). Uncertanty Modeng wth Appcatons to Mutdmensona Cv Engneerng Systems. McGraw-H Internatona, New Yor. Kramer, S.L. (1995). Geotechnca Earthquae Engneerng. Prentce Ha, New Jersey. Krawner, H. (00). A Genera Approach to Sesmc Performance Assessment. Proceedngs, Internatona Conference on Advances and New Chaenges n Earthquae Engneerng Research, ICANCEER 00, Hong Kong, August 19-0, 00. Luco, Ncoas and Corne, C. An (00). ture-specfc Scaar Intensty Measures for Near-Source and Ordnary Earthquae Ground Motons. Earthquae Spectra, (submtted). Mechers, Robert E (1999). tura Reabty Anayss and Predcton. John Wey and Sons, Chchester. Porter, K.A, Bec, J.L., and Shahutdnov, R.V. (00). Investgaton of Senstvty of Budng Loss Estmates to Maor Uncertan Varabes for the Van Nuys Testbed. Report to Pacfc Earthquae Engneerng Research Center, Bereey. 41pp. RS Means Co., Inc. (00). Means Assembes Cost Data. Kngston, MA. Vamvatsos, D. and Corne, C.A. (00). Incrementa Dynamc Anayss. Earthquae Engneerng and tura Dynamcs, 31(3): ACKNOWLEDGEMENTS Fnanca support by PEER (proect number ) and the Shah Famy Fund feowshp s greaty apprecated. REFERENCES Asan, H. and Mranda, E. (00). Fragty of Sab-Coumn Connectons n Renforced Concrete Budngs. ASCE Journa of tura Engneerng, (submtted).

Pacific Earthquake Engineering Research Center

Pacific Earthquake Engineering Research Center Pacfc Earthquae Engneerng Research Center Uncertanty Specfcaton and Propagaton for Loss Estmaton Usng FOSM Method Jac W. Baer Stanford Unversty and C. Alln Cornell Stanford Unversty PEER 003/07 SEPT. 003

More information

MARKOV CHAIN AND HIDDEN MARKOV MODEL

MARKOV CHAIN AND HIDDEN MARKOV MODEL MARKOV CHAIN AND HIDDEN MARKOV MODEL JIAN ZHANG JIANZHAN@STAT.PURDUE.EDU Markov chan and hdden Markov mode are probaby the smpest modes whch can be used to mode sequenta data,.e. data sampes whch are not

More information

A DIMENSION-REDUCTION METHOD FOR STOCHASTIC ANALYSIS SECOND-MOMENT ANALYSIS

A DIMENSION-REDUCTION METHOD FOR STOCHASTIC ANALYSIS SECOND-MOMENT ANALYSIS A DIMESIO-REDUCTIO METHOD FOR STOCHASTIC AALYSIS SECOD-MOMET AALYSIS S. Rahman Department of Mechanca Engneerng and Center for Computer-Aded Desgn The Unversty of Iowa Iowa Cty, IA 52245 June 2003 OUTLIE

More information

COXREG. Estimation (1)

COXREG. Estimation (1) COXREG Cox (972) frst suggested the modes n whch factors reated to fetme have a mutpcatve effect on the hazard functon. These modes are caed proportona hazards (PH) modes. Under the proportona hazards

More information

Nested case-control and case-cohort studies

Nested case-control and case-cohort studies Outne: Nested case-contro and case-cohort studes Ørnuf Borgan Department of Mathematcs Unversty of Oso NORBIS course Unversty of Oso 4-8 December 217 1 Radaton and breast cancer data Nested case contro

More information

Chapter 6. Rotations and Tensors

Chapter 6. Rotations and Tensors Vector Spaces n Physcs 8/6/5 Chapter 6. Rotatons and ensors here s a speca knd of near transformaton whch s used to transforms coordnates from one set of axes to another set of axes (wth the same orgn).

More information

COMPONENT-LEVEL AND SYSTEM-LEVEL SENSITIVITY STUDY FOR EARTHQUAKE LOSS ESTIMATION

COMPONENT-LEVEL AND SYSTEM-LEVEL SENSITIVITY STUDY FOR EARTHQUAKE LOSS ESTIMATION 13 th World Conference on Earthquae Engneerng Vancouver, B.C., Canada August 1-6, 2004 Paper No. 1070 COMPONEN-LEVEL AND SYSEM-LEVEL SENSIIVIY SUDY FOR EARHQUAKE LOSS ESIMAION Hesameddn ASLANI 1, Eduardo

More information

Neural network-based athletics performance prediction optimization model applied research

Neural network-based athletics performance prediction optimization model applied research Avaabe onne www.jocpr.com Journa of Chemca and Pharmaceutca Research, 04, 6(6):8-5 Research Artce ISSN : 0975-784 CODEN(USA) : JCPRC5 Neura networ-based athetcs performance predcton optmzaton mode apped

More information

Predicting Model of Traffic Volume Based on Grey-Markov

Predicting Model of Traffic Volume Based on Grey-Markov Vo. No. Modern Apped Scence Predctng Mode of Traffc Voume Based on Grey-Marov Ynpeng Zhang Zhengzhou Muncpa Engneerng Desgn & Research Insttute Zhengzhou 5005 Chna Abstract Grey-marov forecastng mode of

More information

Research on Complex Networks Control Based on Fuzzy Integral Sliding Theory

Research on Complex Networks Control Based on Fuzzy Integral Sliding Theory Advanced Scence and Technoogy Letters Vo.83 (ISA 205), pp.60-65 http://dx.do.org/0.4257/ast.205.83.2 Research on Compex etworks Contro Based on Fuzzy Integra Sdng Theory Dongsheng Yang, Bngqng L, 2, He

More information

QUARTERLY OF APPLIED MATHEMATICS

QUARTERLY OF APPLIED MATHEMATICS QUARTERLY OF APPLIED MATHEMATICS Voume XLI October 983 Number 3 DIAKOPTICS OR TEARING-A MATHEMATICAL APPROACH* By P. W. AITCHISON Unversty of Mantoba Abstract. The method of dakoptcs or tearng was ntroduced

More information

COMBINING SPATIAL COMPONENTS IN SEISMIC DESIGN

COMBINING SPATIAL COMPONENTS IN SEISMIC DESIGN Transactons, SMRT- COMBINING SPATIAL COMPONENTS IN SEISMIC DESIGN Mchae O Leary, PhD, PE and Kevn Huberty, PE, SE Nucear Power Technooges Dvson, Sargent & Lundy, Chcago, IL 6060 ABSTRACT Accordng to Reguatory

More information

Lower bounds for the Crossing Number of the Cartesian Product of a Vertex-transitive Graph with a Cycle

Lower bounds for the Crossing Number of the Cartesian Product of a Vertex-transitive Graph with a Cycle Lower bounds for the Crossng Number of the Cartesan Product of a Vertex-transtve Graph wth a Cyce Junho Won MIT-PRIMES December 4, 013 Abstract. The mnmum number of crossngs for a drawngs of a gven graph

More information

G : Statistical Mechanics

G : Statistical Mechanics G25.2651: Statstca Mechancs Notes for Lecture 11 I. PRINCIPLES OF QUANTUM STATISTICAL MECHANICS The probem of quantum statstca mechancs s the quantum mechanca treatment of an N-partce system. Suppose the

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Image Classification Using EM And JE algorithms

Image Classification Using EM And JE algorithms Machne earnng project report Fa, 2 Xaojn Sh, jennfer@soe Image Cassfcaton Usng EM And JE agorthms Xaojn Sh Department of Computer Engneerng, Unversty of Caforna, Santa Cruz, CA, 9564 jennfer@soe.ucsc.edu

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Note 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2

Note 2. Ling fong Li. 1 Klein Gordon Equation Probablity interpretation Solutions to Klein-Gordon Equation... 2 Note 2 Lng fong L Contents Ken Gordon Equaton. Probabty nterpretaton......................................2 Soutons to Ken-Gordon Equaton............................... 2 2 Drac Equaton 3 2. Probabty nterpretaton.....................................

More information

Correspondence. Performance Evaluation for MAP State Estimate Fusion I. INTRODUCTION

Correspondence. Performance Evaluation for MAP State Estimate Fusion I. INTRODUCTION Correspondence Performance Evauaton for MAP State Estmate Fuson Ths paper presents a quanttatve performance evauaton method for the maxmum a posteror (MAP) state estmate fuson agorthm. Under dea condtons

More information

Supplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks

Supplementary Material: Learning Structured Weight Uncertainty in Bayesian Neural Networks Shengyang Sun, Changyou Chen, Lawrence Carn Suppementary Matera: Learnng Structured Weght Uncertanty n Bayesan Neura Networks Shengyang Sun Changyou Chen Lawrence Carn Tsnghua Unversty Duke Unversty Duke

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

On the Power Function of the Likelihood Ratio Test for MANOVA

On the Power Function of the Likelihood Ratio Test for MANOVA Journa of Mutvarate Anayss 8, 416 41 (00) do:10.1006/jmva.001.036 On the Power Functon of the Lkehood Rato Test for MANOVA Dua Kumar Bhaumk Unversty of South Aabama and Unversty of Inos at Chcago and Sanat

More information

Reliability Sensitivity Algorithm Based on Stratified Importance Sampling Method for Multiple Failure Modes Systems

Reliability Sensitivity Algorithm Based on Stratified Importance Sampling Method for Multiple Failure Modes Systems Chnese Journa o Aeronautcs 3(010) 660-669 Chnese Journa o Aeronautcs www.esever.com/ocate/ca Reabty Senstvty Agorthm Based on Strated Importance Sampng Method or Mutpe aure Modes Systems Zhang eng a, u

More information

L-Edge Chromatic Number Of A Graph

L-Edge Chromatic Number Of A Graph IJISET - Internatona Journa of Innovatve Scence Engneerng & Technoogy Vo. 3 Issue 3 March 06. ISSN 348 7968 L-Edge Chromatc Number Of A Graph Dr.R.B.Gnana Joth Assocate Professor of Mathematcs V.V.Vannaperuma

More information

3. Stress-strain relationships of a composite layer

3. Stress-strain relationships of a composite layer OM PO I O U P U N I V I Y O F W N ompostes ourse 8-9 Unversty of wente ng. &ech... tress-stran reatonshps of a composte ayer - Laurent Warnet & emo Aerman.. tress-stran reatonshps of a composte ayer Introducton

More information

The Application of BP Neural Network principal component analysis in the Forecasting the Road Traffic Accident

The Application of BP Neural Network principal component analysis in the Forecasting the Road Traffic Accident ICTCT Extra Workshop, Bejng Proceedngs The Appcaton of BP Neura Network prncpa component anayss n Forecastng Road Traffc Accdent He Mng, GuoXucheng &LuGuangmng Transportaton Coege of Souast Unversty 07

More information

The line method combined with spectral chebyshev for space-time fractional diffusion equation

The line method combined with spectral chebyshev for space-time fractional diffusion equation Apped and Computatona Mathematcs 014; 3(6): 330-336 Pubshed onne December 31, 014 (http://www.scencepubshnggroup.com/j/acm) do: 10.1164/j.acm.0140306.17 ISS: 3-5605 (Prnt); ISS: 3-5613 (Onne) The ne method

More information

Uncertainty as the Overlap of Alternate Conditional Distributions

Uncertainty as the Overlap of Alternate Conditional Distributions Uncertanty as the Overlap of Alternate Condtonal Dstrbutons Olena Babak and Clayton V. Deutsch Centre for Computatonal Geostatstcs Department of Cvl & Envronmental Engneerng Unversty of Alberta An mportant

More information

Delay tomography for large scale networks

Delay tomography for large scale networks Deay tomography for arge scae networks MENG-FU SHIH ALFRED O. HERO III Communcatons and Sgna Processng Laboratory Eectrca Engneerng and Computer Scence Department Unversty of Mchgan, 30 Bea. Ave., Ann

More information

Associative Memories

Associative Memories Assocatve Memores We consder now modes for unsupervsed earnng probems, caed auto-assocaton probems. Assocaton s the task of mappng patterns to patterns. In an assocatve memory the stmuus of an ncompete

More information

Distributed Moving Horizon State Estimation of Nonlinear Systems. Jing Zhang

Distributed Moving Horizon State Estimation of Nonlinear Systems. Jing Zhang Dstrbuted Movng Horzon State Estmaton of Nonnear Systems by Jng Zhang A thess submtted n parta fufment of the requrements for the degree of Master of Scence n Chemca Engneerng Department of Chemca and

More information

Andre Schneider P622

Andre Schneider P622 Andre Schneder P6 Probem Set #0 March, 00 Srednc 7. Suppose that we have a theory wth Negectng the hgher order terms, show that Souton Knowng β(α and γ m (α we can wrte β(α =b α O(α 3 (. γ m (α =c α O(α

More information

A finite difference method for heat equation in the unbounded domain

A finite difference method for heat equation in the unbounded domain Internatona Conerence on Advanced ectronc Scence and Technoogy (AST 6) A nte derence method or heat equaton n the unbounded doman a Quan Zheng and Xn Zhao Coege o Scence North Chna nversty o Technoogy

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

ON AUTOMATIC CONTINUITY OF DERIVATIONS FOR BANACH ALGEBRAS WITH INVOLUTION

ON AUTOMATIC CONTINUITY OF DERIVATIONS FOR BANACH ALGEBRAS WITH INVOLUTION European Journa of Mathematcs and Computer Scence Vo. No. 1, 2017 ON AUTOMATC CONTNUTY OF DERVATONS FOR BANACH ALGEBRAS WTH NVOLUTON Mohamed BELAM & Youssef T DL MATC Laboratory Hassan Unversty MORO CCO

More information

Multispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory

Multispectral Remote Sensing Image Classification Algorithm Based on Rough Set Theory Proceedngs of the 2009 IEEE Internatona Conference on Systems Man and Cybernetcs San Antono TX USA - October 2009 Mutspectra Remote Sensng Image Cassfcaton Agorthm Based on Rough Set Theory Yng Wang Xaoyun

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

LECTURE 21 Mohr s Method for Calculation of General Displacements. 1 The Reciprocal Theorem

LECTURE 21 Mohr s Method for Calculation of General Displacements. 1 The Reciprocal Theorem V. DEMENKO MECHANICS OF MATERIALS 05 LECTURE Mohr s Method for Cacuaton of Genera Dspacements The Recproca Theorem The recproca theorem s one of the genera theorems of strength of materas. It foows drect

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

GENERATIVE AND DISCRIMINATIVE CLASSIFIERS: NAIVE BAYES AND LOGISTIC REGRESSION. Machine Learning

GENERATIVE AND DISCRIMINATIVE CLASSIFIERS: NAIVE BAYES AND LOGISTIC REGRESSION. Machine Learning CHAPTER 3 GENERATIVE AND DISCRIMINATIVE CLASSIFIERS: NAIVE BAYES AND LOGISTIC REGRESSION Machne Learnng Copyrght c 205. Tom M. Mtche. A rghts reserved. *DRAFT OF September 23, 207* *PLEASE DO NOT DISTRIBUTE

More information

Prediction Error of the Multivariate Additive Loss Reserving Method for Dependent Lines of Business

Prediction Error of the Multivariate Additive Loss Reserving Method for Dependent Lines of Business Predcton Error of the Mutvarate Addtve Loss Reservng Method for Dependent Lnes of Busness by Mchae Merz and Maro V Wüthrch ABSTRACT Often n non-fe nsurance, cams reserves are the argest poston on the abty

More information

[WAVES] 1. Waves and wave forces. Definition of waves

[WAVES] 1. Waves and wave forces. Definition of waves 1. Waves and forces Defnton of s In the smuatons on ong-crested s are consdered. The drecton of these s (μ) s defned as sketched beow n the goba co-ordnate sstem: North West East South The eevaton can

More information

Numerical integration in more dimensions part 2. Remo Minero

Numerical integration in more dimensions part 2. Remo Minero Numerca ntegraton n more dmensons part Remo Mnero Outne The roe of a mappng functon n mutdmensona ntegraton Gauss approach n more dmensons and quadrature rues Crtca anass of acceptabt of a gven quadrature

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

Research Article H Estimates for Discrete-Time Markovian Jump Linear Systems

Research Article H Estimates for Discrete-Time Markovian Jump Linear Systems Mathematca Probems n Engneerng Voume 213 Artce ID 945342 7 pages http://dxdoorg/11155/213/945342 Research Artce H Estmates for Dscrete-Tme Markovan Jump Lnear Systems Marco H Terra 1 Gdson Jesus 2 and

More information

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

More information

EXPERIMENT AND THEORISATION: AN APPLICATION OF THE HYDROSTATIC EQUATION AND ARCHIMEDES THEOREM

EXPERIMENT AND THEORISATION: AN APPLICATION OF THE HYDROSTATIC EQUATION AND ARCHIMEDES THEOREM EXPERIMENT AND THEORISATION: AN APPLICATION OF THE HYDROSTATIC EQUATION AND ARCHIMEDES THEOREM Santos Lucía, Taaa Máro, Departamento de Físca, Unversdade de Avero, Avero, Portuga 1. Introducton Today s

More information

Optimum Selection Combining for M-QAM on Fading Channels

Optimum Selection Combining for M-QAM on Fading Channels Optmum Seecton Combnng for M-QAM on Fadng Channes M. Surendra Raju, Ramesh Annavajjaa and A. Chockangam Insca Semconductors Inda Pvt. Ltd, Bangaore-56000, Inda Department of ECE, Unversty of Caforna, San

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

Numerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes

Numerical Investigation of Power Tunability in Two-Section QD Superluminescent Diodes Numerca Investgaton of Power Tunabty n Two-Secton QD Superumnescent Dodes Matta Rossett Paoo Bardea Ivo Montrosset POLITECNICO DI TORINO DELEN Summary 1. A smpfed mode for QD Super Lumnescent Dodes (SLD)

More information

Analysis of Bipartite Graph Codes on the Binary Erasure Channel

Analysis of Bipartite Graph Codes on the Binary Erasure Channel Anayss of Bpartte Graph Codes on the Bnary Erasure Channe Arya Mazumdar Department of ECE Unversty of Maryand, Coege Par ema: arya@umdedu Abstract We derve densty evouton equatons for codes on bpartte

More information

Numerical Heat and Mass Transfer

Numerical Heat and Mass Transfer Master degree n Mechancal Engneerng Numercal Heat and Mass Transfer 06-Fnte-Dfference Method (One-dmensonal, steady state heat conducton) Fausto Arpno f.arpno@uncas.t Introducton Why we use models and

More information

A Dissimilarity Measure Based on Singular Value and Its Application in Incremental Discounting

A Dissimilarity Measure Based on Singular Value and Its Application in Incremental Discounting A Dssmarty Measure Based on Snguar Vaue and Its Appcaton n Incrementa Dscountng KE Xaou Department of Automaton, Unversty of Scence and Technoogy of Chna, Hefe, Chna Ema: kxu@ma.ustc.edu.cn MA Lyao Department

More information

Cyclic Codes BCH Codes

Cyclic Codes BCH Codes Cycc Codes BCH Codes Gaos Feds GF m A Gaos fed of m eements can be obtaned usng the symbos 0,, á, and the eements beng 0,, á, á, á 3 m,... so that fed F* s cosed under mutpcaton wth m eements. The operator

More information

D hh ν. Four-body charm semileptonic decay. Jim Wiss University of Illinois

D hh ν. Four-body charm semileptonic decay. Jim Wiss University of Illinois Four-body charm semeptonc decay Jm Wss Unversty of Inos D hh ν 1 1. ector domnance. Expected decay ntensty 3. SU(3) apped to D s φν 4. Anaytc forms for form factors 5. Non-parametrc form factors 6. Future

More information

Inthem-machine flow shop problem, a set of jobs, each

Inthem-machine flow shop problem, a set of jobs, each THE ASYMPTOTIC OPTIMALITY OF THE SPT RULE FOR THE FLOW SHOP MEAN COMPLETION TIME PROBLEM PHILIP KAMINSKY Industra Engneerng and Operatons Research, Unversty of Caforna, Bereey, Caforna 9470, amnsy@eor.bereey.edu

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Boundary Value Problems. Lecture Objectives. Ch. 27

Boundary Value Problems. Lecture Objectives. Ch. 27 Boundar Vaue Probes Ch. 7 Lecture Obectves o understand the dfference between an nta vaue and boundar vaue ODE o be abe to understand when and how to app the shootng ethod and FD ethod. o understand what

More information

Item calibration in incomplete testing designs

Item calibration in incomplete testing designs Pscoógca (20), 32, 07-32 Item cabraton n ncompete testng desgns Theo JHM Eggen * & Norman D Verhest** *Cto/Unversty of Twente, The etherands **Cto, The etherands Ths study dscusses the justfabty of tem

More information

ACCURATE COMPUTATION OF CRITICAL RESPONSE QUANTITES FOR LAMINATED COMPOSITE STRUCTURES

ACCURATE COMPUTATION OF CRITICAL RESPONSE QUANTITES FOR LAMINATED COMPOSITE STRUCTURES ACCURATE COMPUTATION OF CRITICAL RESPONSE QUANTITES FOR LAMINATED COMPOSITE STRUCTURES C.S. UPADHYAY, P.M. MOHITE and A. K. ONKAR Department of Aerospace Engneerng Indan Insttute of Technoogy Kanpur Kanpur

More information

Globally Optimal Multisensor Distributed Random Parameter Matrices Kalman Filtering Fusion with Applications

Globally Optimal Multisensor Distributed Random Parameter Matrices Kalman Filtering Fusion with Applications Sensors 2008, 8, 8086-8103; DOI: 10.3390/s8128086 OPEN ACCESS sensors ISSN 1424-8220 www.mdp.com/journa/sensors Artce Gobay Optma Mutsensor Dstrbuted Random Parameter Matrces Kaman Fterng Fuson wth Appcatons

More information

Comparison of Regression Lines

Comparison of Regression Lines STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence

More information

Journal of Multivariate Analysis

Journal of Multivariate Analysis Journa of Mutvarate Anayss 3 (04) 74 96 Contents sts avaabe at ScenceDrect Journa of Mutvarate Anayss journa homepage: www.esever.com/ocate/jmva Hgh-dmensona sparse MANOVA T. Tony Ca a, Yn Xa b, a Department

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Example: Suppose we want to build a classifier that recognizes WebPages of graduate students.

Example: Suppose we want to build a classifier that recognizes WebPages of graduate students. Exampe: Suppose we want to bud a cassfer that recognzes WebPages of graduate students. How can we fnd tranng data? We can browse the web and coect a sampe of WebPages of graduate students of varous unverstes.

More information

Sampling-based Approach for Design Optimization in the Presence of Interval Variables

Sampling-based Approach for Design Optimization in the Presence of Interval Variables 0 th Word Congress on Structura and Mutdscpnary Optmzaton May 9-4, 03, Orando, orda, USA Sampng-based Approach for Desgn Optmzaton n the Presence of nterva Varabes Davd Yoo and kn Lee Unversty of Connectcut,

More information

Module III, Lecture 02: Correlation Functions and Spectral Densities

Module III, Lecture 02: Correlation Functions and Spectral Densities Modue III, Lecture : orreaton Functons and Spectra enstes The ony part of Boch Redfed Wangsness spn reaxaton theory that we coud say nothng about n the prevous ecture were the correaton functons. They

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models Computaton of Hgher Order Moments from Two Multnomal Overdsperson Lkelhood Models BY J. T. NEWCOMER, N. K. NEERCHAL Department of Mathematcs and Statstcs, Unversty of Maryland, Baltmore County, Baltmore,

More information

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS

NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS IJRRAS 8 (3 September 011 www.arpapress.com/volumes/vol8issue3/ijrras_8_3_08.pdf NON-CENTRAL 7-POINT FORMULA IN THE METHOD OF LINES FOR PARABOLIC AND BURGERS' EQUATIONS H.O. Bakodah Dept. of Mathematc

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Polite Water-filling for Weighted Sum-rate Maximization in MIMO B-MAC Networks under. Multiple Linear Constraints

Polite Water-filling for Weighted Sum-rate Maximization in MIMO B-MAC Networks under. Multiple Linear Constraints 2011 IEEE Internatona Symposum on Informaton Theory Proceedngs Pote Water-fng for Weghted Sum-rate Maxmzaton n MIMO B-MAC Networks under Mutpe near Constrants An u 1, Youjan u 2, Vncent K. N. au 3, Hage

More information

Application of support vector machine in health monitoring of plate structures

Application of support vector machine in health monitoring of plate structures Appcaton of support vector machne n heath montorng of pate structures *Satsh Satpa 1), Yogesh Khandare ), Sauvk Banerjee 3) and Anrban Guha 4) 1), ), 4) Department of Mechanca Engneerng, Indan Insttute

More information

Optimization of JK Flip Flop Layout with Minimal Average Power of Consumption based on ACOR, Fuzzy-ACOR, GA, and Fuzzy-GA

Optimization of JK Flip Flop Layout with Minimal Average Power of Consumption based on ACOR, Fuzzy-ACOR, GA, and Fuzzy-GA Journa of mathematcs and computer Scence 4 (05) - 5 Optmzaton of JK Fp Fop Layout wth Mnma Average Power of Consumpton based on ACOR, Fuzzy-ACOR, GA, and Fuzzy-GA Farshd Kevanan *,, A Yekta *,, Nasser

More information

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University

PHYS 450 Spring semester Lecture 02: Dealing with Experimental Uncertainties. Ron Reifenberger Birck Nanotechnology Center Purdue University PHYS 45 Sprng semester 7 Lecture : Dealng wth Expermental Uncertantes Ron Refenberger Brck anotechnology Center Purdue Unversty Lecture Introductory Comments Expermental errors (really expermental uncertantes)

More information

Decentralized Adaptive Control for a Class of Large-Scale Nonlinear Systems with Unknown Interactions

Decentralized Adaptive Control for a Class of Large-Scale Nonlinear Systems with Unknown Interactions Decentrazed Adaptve Contro for a Cass of Large-Scae onnear Systems wth Unknown Interactons Bahram Karm 1, Fatemeh Jahangr, Mohammad B. Menhaj 3, Iman Saboor 4 1. Center of Advanced Computatona Integence,

More information

DISTRIBUTED PROCESSING OVER ADAPTIVE NETWORKS. Cassio G. Lopes and Ali H. Sayed

DISTRIBUTED PROCESSING OVER ADAPTIVE NETWORKS. Cassio G. Lopes and Ali H. Sayed DISTRIBUTED PROCESSIG OVER ADAPTIVE ETWORKS Casso G Lopes and A H Sayed Department of Eectrca Engneerng Unversty of Caforna Los Angees, CA, 995 Ema: {casso, sayed@eeucaedu ABSTRACT Dstrbuted adaptve agorthms

More information

Quantum Runge-Lenz Vector and the Hydrogen Atom, the hidden SO(4) symmetry

Quantum Runge-Lenz Vector and the Hydrogen Atom, the hidden SO(4) symmetry Quantum Runge-Lenz ector and the Hydrogen Atom, the hdden SO(4) symmetry Pasca Szrftgser and Edgardo S. Cheb-Terrab () Laboratore PhLAM, UMR CNRS 85, Unversté Le, F-59655, France () Mapesoft Let's consder

More information

CS 2750 Machine Learning. Lecture 5. Density estimation. CS 2750 Machine Learning. Announcements

CS 2750 Machine Learning. Lecture 5. Density estimation. CS 2750 Machine Learning. Announcements CS 750 Machne Learnng Lecture 5 Densty estmaton Mlos Hauskrecht mlos@cs.ptt.edu 539 Sennott Square CS 750 Machne Learnng Announcements Homework Due on Wednesday before the class Reports: hand n before

More information

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

Lower Bounding Procedures for the Single Allocation Hub Location Problem

Lower Bounding Procedures for the Single Allocation Hub Location Problem Lower Boundng Procedures for the Snge Aocaton Hub Locaton Probem Borzou Rostam 1,2 Chrstoph Buchhem 1,4 Fautät für Mathemat, TU Dortmund, Germany J. Faban Meer 1,3 Uwe Causen 1 Insttute of Transport Logstcs,

More information

Laboratory 1c: Method of Least Squares

Laboratory 1c: Method of Least Squares Lab 1c, Least Squares Laboratory 1c: Method of Least Squares Introducton Consder the graph of expermental data n Fgure 1. In ths experment x s the ndependent varable and y the dependent varable. Clearly

More information

Differentiating Gaussian Processes

Differentiating Gaussian Processes Dfferentatng Gaussan Processes Andrew McHutchon Aprl 17, 013 1 Frst Order Dervatve of the Posteror Mean The posteror mean of a GP s gven by, f = x, X KX, X 1 y x, X α 1 Only the x, X term depends on the

More information

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models ECO 452 -- OE 4: Probt and Logt Models ECO 452 -- OE 4 Maxmum Lkelhood Estmaton of Bnary Dependent Varables Models: Probt and Logt hs note demonstrates how to formulate bnary dependent varables models

More information

Lowest-Order e + e l + l Processes in Quantum Electrodynamics. Sanha Cheong

Lowest-Order e + e l + l Processes in Quantum Electrodynamics. Sanha Cheong Lowest-Order e + e + Processes n Quantum Eectrodynamcs Sanha Cheong Introducton In ths short paper, we w demonstrate some o the smpest cacuatons n quantum eectrodynamcs (QED), eadng to the owest-order

More information

RELIABILITY ASSESSMENT

RELIABILITY ASSESSMENT CHAPTER Rsk Analyss n Engneerng and Economcs RELIABILITY ASSESSMENT A. J. Clark School of Engneerng Department of Cvl and Envronmental Engneerng 4a CHAPMAN HALL/CRC Rsk Analyss for Engneerng Department

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family

Using T.O.M to Estimate Parameter of distributions that have not Single Exponential Family IOSR Journal of Mathematcs IOSR-JM) ISSN: 2278-5728. Volume 3, Issue 3 Sep-Oct. 202), PP 44-48 www.osrjournals.org Usng T.O.M to Estmate Parameter of dstrbutons that have not Sngle Exponental Famly Jubran

More information

arxiv: v1 [physics.comp-ph] 17 Dec 2018

arxiv: v1 [physics.comp-ph] 17 Dec 2018 Pressures nsde a nano-porous medum. The case of a snge phase fud arxv:1812.06656v1 [physcs.comp-ph] 17 Dec 2018 Oav Gateand, Dck Bedeaux, and Sgne Kjestrup PoreLab, Department of Chemstry, Norwegan Unversty

More information

6.3.4 Modified Euler s method of integration

6.3.4 Modified Euler s method of integration 6.3.4 Modfed Euler s method of ntegraton Before dscussng the applcaton of Euler s method for solvng the swng equatons, let us frst revew the basc Euler s method of numercal ntegraton. Let the general from

More information

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise. Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where y + = β + β e for =,..., y and are observable varables e s a random error How can an estmaton rule be constructed for the

More information

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation Statstcs for Managers Usng Mcrosoft Excel/SPSS Chapter 13 The Smple Lnear Regresson Model and Correlaton 1999 Prentce-Hall, Inc. Chap. 13-1 Chapter Topcs Types of Regresson Models Determnng the Smple Lnear

More information

Drift Design Method for High-rise Buildings using Dynamic Resizing Algorithm

Drift Design Method for High-rise Buildings using Dynamic Resizing Algorithm ctbuh.org/papers Tte: Authors: Subject: Keywords: Drft Desgn Method for Hgh-rse Budngs usng Dynac Reszng Agorth Hyo Seon Park, Ph.D Canddate, Yonse Unversty Seo J Hyun, Assocate Professor, Yonse Unversty

More information

SPATIAL KINEMATICS OF GEARS IN ABSOLUTE COORDINATES

SPATIAL KINEMATICS OF GEARS IN ABSOLUTE COORDINATES SATIAL KINEMATICS OF GEARS IN ABSOLUTE COORDINATES Dmtry Vasenko and Roand Kasper Insttute of Mobe Systems (IMS) Otto-von-Guercke-Unversty Magdeburg D-39016, Magdeburg, Germany E-ma: Dmtr.Vasenko@ovgu.de

More information

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression 11 MACHINE APPLIED MACHINE LEARNING LEARNING MACHINE LEARNING Gaussan Mture Regresson 22 MACHINE APPLIED MACHINE LEARNING LEARNING Bref summary of last week s lecture 33 MACHINE APPLIED MACHINE LEARNING

More information

22.51 Quantum Theory of Radiation Interactions

22.51 Quantum Theory of Radiation Interactions .51 Quantum Theory of Radaton Interactons Fna Exam - Soutons Tuesday December 15, 009 Probem 1 Harmonc oscator 0 ponts Consder an harmonc oscator descrbed by the Hamtonan H = ω(nˆ + ). Cacuate the evouton

More information

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure

CinChE Problem-Solving Strategy Chapter 4 Development of a Mathematical Model. formulation. procedure nhe roblem-solvng Strategy hapter 4 Transformaton rocess onceptual Model formulaton procedure Mathematcal Model The mathematcal model s an abstracton that represents the engneerng phenomena occurrng n

More information

18.1 Introduction and Recap

18.1 Introduction and Recap CS787: Advanced Algorthms Scrbe: Pryananda Shenoy and Shjn Kong Lecturer: Shuch Chawla Topc: Streamng Algorthmscontnued) Date: 0/26/2007 We contnue talng about streamng algorthms n ths lecture, ncludng

More information

e i is a random error

e i is a random error Chapter - The Smple Lnear Regresson Model The lnear regresson equaton s: where + β + β e for,..., and are observable varables e s a random error How can an estmaton rule be constructed for the unknown

More information