A NOTE ON BAYESIAN ANALYSIS OF THE. University of Oxford. and. A. C. Davison. Swiss Federal Institute of Technology. March 10, 1998.

Size: px
Start display at page:

Download "A NOTE ON BAYESIAN ANALYSIS OF THE. University of Oxford. and. A. C. Davison. Swiss Federal Institute of Technology. March 10, 1998."

Transcription

1 A NOTE ON BAYESIAN ANALYSIS OF THE POLY-WEIBULL MODEL F. Luzada-Net University f Oxfrd and A. C. Davisn Swiss Federal Institute f Technlgy March 10, 1998 Summary We cnsider apprximate Bayesian analysis fr the ply-weibull mdel, which arises in cmpetitive risk scenaris when the risks are independent and there is n knwledge abut which factr was respnsible fr failure. Real and generated datasets illustrate the methdlgy, which is based n Laplace's methd fr integrals. Sme key wrds: Apprximate Bayesian analysis, bi-weibull hazard mdel; cmpeting risks; Laplace's methd fr integrals; Weibull distributin. Authr's ftnte: F. Luzada-Net is a research student at the University f Oxfrd, Department f Statistics, 1 Suth Parks Rad, Oxfrd OX1 3TG, UK, and A. C. Davisn is Prfessr f Statistics, Swiss Federal Institute f Technlgy, Department f Mathematics, 1015 Lausanne, Switzerland. F. Luzada-Net is grateful t the Brazilian institutins UFSCar and CNPq fr their nancial supprt. 1

2 1 INTRODUCTION In a recent paper Berger and Sun (1993) describe Bayesian analysis via the Gibbs sampler fr a ply-weibull mdel, which arises in scenaris f cmpeting risks (Cx and Oakes, 1984, Ch. 7) when nly the minimum f several Weibull failure times is bserved. The advantage f the ply-weibull mdel in relatin t the standard Weibull mdel is that it allws nt nly increasing, cnstant r decreasing hazard functins, but als nnmntne nes. This is an attractive prperty since such hazards are nt uncmmn in practice. The purpse f this paper is t pint ut that iterative simulatin is nt necessary fr Bayesian analysis f ply-weibull mdels, which can be dealt with mre simply using Laplace's methd (Tierney and Kadane, 1986) r the Bayesian btstrap (Smith and Gelfand, 1992). This reduces cnsiderably the cmputatinal burden in btaining psterir prbabilities, densities, r mments, and makes pssible the use f a wider class f prir densities because there is n need t cast the prblem in the Gibbs sampler framewrk. By the same tken, and despite the cmments in Berger and Sun (1993, Sectin 1.3), standard likelihd analysis f the mdel is entirely straightfrward, the nly diculty arising because f the pssible nn-identiability f the mdel, which we discuss in the next sectin. Mdel frmulatin is discussed in Sectin 2, where we briey discuss the identiability f the parameters, and the use f Laplace's methd is utlined in Sectin 3. In Sectin 4, the methdlgy is illustrated with datasets, tw generated and ne real, and an alternative analysis via the Bayesian btstrap is briey utlined. 2 MODEL FORMULATION Suppse that an individual r system is cmpsed f m 2 elements, and that the lifetime X f the th element has hazard functin h (t). The quantity bserved is T = min(x1; : : : ; X m ), which is said t have a ply-weibull mdel if its hazard functin is given by mx mx t?1 h(t) = h (t) = =1 =1 ; t > 0; (1) where ; > 0 are parameters assciated t the th element. This arises if the individual elements X have independent Weibull distributins with parameters and. An identiability prblem arises if the shape parameters 1; : : : ; m are equal, fr then (1) can be written as t?1, which is the hazard functin f a single Weibull randm P variable; in this case?1 =?. The cause f the prblem is that a redundant vectr f parameters has been kept in the parametrizatin (Bernard and Smith, 1994, p. 239). In practice the parameters are nly identiable when the risk factrs have suciently separated hazard functins. We have fund that graphical gdness-f-t methds such as hazard plts (Lawless, 1982, pp. 88, 278) are very useful fr checking when ply-hazard mdels are likely 2

3 lg H(t) lg H(t) lg t lg t Figure 1: Hazard plts fr simulated bi-weibull data. Left panel: data generated with 1 = 1000, 1 = 0:5, 2 = 3200 and 2 = 5. Right panel: data generated with 1 = 1000, 1 = 0:6, 2 = 3200 and 2 = 1:6 (Berger and Sun, 1993, p. 1416). The dtted lines represent maximum likelihd ts f a single Weibull distributin. t be identiable. T illustrate this, Figure 1 shws hazard plts fr tw datasets generated frm a ply-weibull mdel with m = 2 elements; we call this a bi-weibull mdel. The left panel f the gure shws a clear departure frm the straight line that crrespnds t a single Weibull distributin, but this is nt the case fr the sample in the right panel, which was simulated by Berger and Sun (1993) frm a bi-weibull mdel. This suggests nnidentiable parameters r parameters that are nly identiable thrugh functins q(; ) (Basu and Klein, 1982), where = (1; : : : ; m ) T, and = (1; : : : ; m ) T. 3 ESTIMATION Cnsider a sample f independent psitive randm variables T1; : : : ; T n with cmmn hazard functin h(t) and such that T i has assciated an indicatr variable dened by i = 1 if T i is an bserved failure time and i = 0 if T i is a right-censred bservatin. Up t an additive cnstant, the lg likelihd functin fr the parameters f any set f randmly-censred survival data frm this mdel, (t1; 1); : : : ; (t n ; n ), may be expressed as lg L(; ) = nx i lg h(t i )? H(t i ); (2) i=1 where H(t) is the cumulative hazard functin H(t) = R t 0 h(u) du. On substituting the bi- Weibull hazard (1) int (2) we see that the lg likelihd functin fr the parameter vectrs 3

4 , is lg L(; ) = ( nx X i lg A i? m i=1 =1 ) ti ; (3) P m where A i = =1 t?1 i?. Cntrary t the discussin in Sectins 1.3 and 1.4 f Berger and Sun (1993) there is n diculty in direct calculatin f (3), which simply invlves tw nested d-lps and can be perfrmed in a single S-Plus statement, and cnsequently there is n bstacle t standard frequentist likelihd analysis f the mdel, prvided its parameters are identiable. Tests fr nn-identiability culd be develped, but seem unlikely t be necessary in mst applicatins, where plts like thse in Figure 1 will suce t detect ptential prblems with likelihd analysis. The maximum likelihd estimates f and can be btained by slving the system f 2m nnlinear equatins nx i=1 where g ir = r?r r nx i=1 t i? i A i r i g ir A i?1 r + lg(t i = r )? t i r g ir lg(t i = r ) r g ir = 0; r = 1; : : : ; m; (4) = 0; r = 1; : : : ; m; (5) t r?1 i. Frm a frequentist viewpint interval inference fr the parameters is prbably best cnducted using their prle lg likelihds, as experience suggests that the asympttic apprximatins t the distributins f likelihd rati statistics are likely t be mre accurate in small samples than are asympttic apprximatins t maximum likelihd estimatrs. The parametric btstrap can be applied t avid dependence n asympttic apprximatins when samples are small; see Davisn and Hinkley (1997, p. 346) fr a related example. As the methds invlved are entirely standard, we turn t Bayesian inference rather than cnsider them further. When there is n strng prir infrmatin, the simplest pssibility is t assume a cnstant prir density fr (lg ; lg ), in which case the int lg psterir density fr and is f curse (3), apart frm the nrmalizing cnstant. This als crrespnds t a Jereys prir fr the parameters f the separate Weibull distributins, because the lg Weibull distributin is a lcatin-scale family with lcatin parameter lg and scale parameter. Unfrtunately this results in an imprper psterir distributin, and less vague prir infrmatin must be emplyed. In the examples belw we used independent nrmal distributins with means at the maximum likelihd estimates and standard deviatins 100 as prirs fr the parameters lg and lg. This is purely intended t represent very weak but prper prir infrmatin; f curse in an applicatin ne wuld nt use a data-based prir such as this. With the understanding that L 0 (; ) hencefrth dentes the sum f (2) and the lg prir density, and using Laplace's methd fr apprximatin f the tw integrals invlved in the marginal density f ( r ; r ) (Tierney and Kadane, 1986; Kass, Tierney and Kadane, 1990), 4

5 the apprximate lg int psterir density fr ( r ; r ), fr r = 1; 2; : : : ; m, is lg L 0 ( ; )? 1 2 lg D?r? lg L 0 (^; ^)? 1 2 lg ^D? lg(2): (6) Here = ( r ; ^?r ) and = ( r ; ^?r ), are the slutins f the analgues f (4) and (5) with ( r ; r ) held xed, D?r is the determinant f the Hessian matrix f minus the lg int psterir density fr?r and?r evaluated at and, and the crrespnding quantities maximizing ver all the parameters are ^, ^ and ^D. Expressin (6) must be evaluated ver a grid f values f ( r ; r ), and then cntured, but the ert invlved is much less than using Markv chain Mnte Carl methds, and there is n restrictin t particular classes f prirs in rder that it is pssible easily t simulate frm the necessary cnditinal distributins. Laplace's methd may fail if the psterir density is seriusly multi-mdal, but we have nt encuntered this prblem in any f the examples t which we have applied the apprximatin. If genuine prir knwledge abut the parameters is available in the frm f a density (; ), it can be incrprated int the analysis in the bvius way. Imprper psterirs can be bserved when ratis f the 's are near ne, a phenmenn related t the near-nn-identiability f the parameters in that situatin. This can be alleviated if there is strng prir infrmatin n the parameters (Bernard and Smith, 1994, p. 239). Berger and Sun (1993) utline methds f eliciting the parameters linked with the prir. 4 EXAMPLES We nw test the methdlgy prpsed in this paper n sme datasets. In each case we use the parametrizatin (lg ; lg ) fr the Laplace apprximatin. This has the benets that the ranges f the lg-transfrmed parameters are unbunded, and that the Hessian matrix f the lg psterir is rughly independent f the parameters. Fr the data in the left panel f Figure 1 the int psterir mdes fr 1 and 1, and fr 2 and 2, with the prir utlined in the previus sectin, are ( ~ 1; ~ 1) = (1304; 0:48) and ( ~ 2; ~ 2) = (3236; 4:88). The apprximate 95% highest psterir density (HPD) intervals fr 1, 1, 2 and 2 are (657; 3191), (0:34; 0:69), (2400; 4023), and (1:58; 13:39), respectively. The Schwarz criterin fr cmparisn f the tw mdels, 3.33, gives an apprximate Bayes factr f 27.9, which accrding t the rugh classicatin in Sectin 3.2 f Kass and Raftery (1995) is `strng' evidence in favur f the bi-weibull mdel. Fr the data in the right panel f Figure 1, the HPD interval fr 2 is (1636; 50000) using ur chice f prir; here we give as a likely upper bund fr an exact value that is dicult t establish. In practice this wuld be regarded as an uninfrmative psterir inference fr 2, and this is related t the identiability prblem pinted ut in Sectin 2. Hwever, if we use the same prirs as Berger and Sun (1993), the int psterir mdes using 5

6 Laplace apprximatin are ( ~ 1; ~ 1) = (866; 0:50) and ( ~ 2; ~ 2) = (3464; 1:92), and apprximate 95% HPD intervals fr 1, 1, 2, and 2 are (508; 1894), (0:44; 0:56), (2513; 7090), and (1:76; 2:03), respectively, similar t the results f Berger and Sun (1993). This suggests that their results depend crucially n the chice f prir, because the parameters cannt be identied frm the data, as Figure 1 shws. Further evidence f this is given by the Schwarz criterin fr cmparisn f the mdels, which is {2.85, crrespnding t a rugh Bayes factr f 17.3 in favur f the single Weibull mdel; this is `psitive' evidence in terms f the classicatin mentined abve. Of curse this evidence is crude, but mre detailed calculatins depend crucially n the chice f prir. It is n accident that the single Weibull mdel is adequate in this case: data sets simulated using the parameter values f Berger and Sun are eectively samples frm a single Weibull distributin, s that any infrmatin that they are in fact bi-weibull must cme frm the prir. The prir used by Berger and Sun eectively rules ut the single Weibull mdel, thereby inecting strng infrmatin that is nt crrbrated by the data. Fr a real example, we cnsider data frm Table 1 f Lagaks and Luis (1988) n the survival f 50 male rats receiving 60 mg/kg f tluene diiscyanate. The rats may have died frm side-eects f the drug, as well as frm the tumurs f the testis it was administered t prevent, but the cause f death is unknwn. The hazard plt in the left panel f Figure 2 suggests that a single Weibull mdel will be inadequate, and that dierent risk factrs may be at wrk. The lwer panels shw the apprximate cntur plts fr the lg int psterir density (6) fr (lg 1; lg 1), and fr (lg 2; lg 2), assuming n prir infrmatin abut the parameters. The int psterir mdes are ( ~ 1; ~ 1) = (122; 0:84) and ( ~ 2; ~ 2) = (107:6; 6:5). Apprximate 95% HPD intervals fr 1 and 1, and 2 and 2 are (78:0; 256:7) and (0:55; 1:17), and (98:8; 122:6) and (3:28; 12:1), respectively. In this case the Schwarz criterin equals 3.45, giving `strng' evidence fr the bi-weibull mdel. This is brne ut by cmparisn f the Kaplan{Meier estimate and the mdal tted mdels, shwn in the tp right panel f Figure 2. Fr a limited assessment f the numerical accuracy f Laplace apprximatin in this cntext, we estimated the survivr functin fr the rat data using sampling-imprtance resampling. We generated 10,000 values f the parameters frm a 4-variate nrmal distributin estimated frm the Laplace apprximatin t the psterir density, and then used a Bayesian btstrap (Smith and Gelfand, 1992) t btain a sample f size 1000 frm the true psterir distributin. This is easily achieved because the prir and likelihd are readily calculated. We then used this sample t give the expected marginal psterir survivr functin fr the rat data, which is shwn as a dashed line in the upper right panel f Figure 2. The curve is very clse t the result f the Laplace apprximatin, which seems adequate fr practical use in this setting. We wuld nt expect these curves t be equal, as sampling-imprtance resampling gives the expected psterir survivr functin, whereas 6

7 ur applicatin f Laplace's methd gives the mdal psterir survivr functin. A large number f applicatins f Laplace's methd wuld be required t apprximate the expected psterir survivr functin, which shuld nt dier frm the mdal value by much. The diference illustrates an interesting facet f the tw appraches, that mdal quantities are dicult t btain by simulatin methds, but relatively easily btained using asympttic apprximatin. REFERENCES Basu, A. and Klein, J. (1982). Sme recent develpments in cmpeting risks thery, in Crwley, J. and Jhnsn, R. A. (eds), Survival Analysis. Hayward, CA: Institute f Mathematical Statistics, 216{229. Berger, J. O. and Sun, D. (1993). Bayesian analysis fr the ply-weibull distributin. Jurnal f the American Statistical Assciatin 88, 1412{1417. Bernard, J. M. and Smith, A. F. M. (1994). Bayesian Thery. New Yrk: Wiley. Cx, D. R. and Oakes, D. (1984). Analysis f Survival Data. Lndn: Chapman & Hall. Davisn, A. C. and Hinkley, D. V. (1997). Cambridge: Cambridge University Press. Btstrap Methds and their Applicatin. Kass, R. E. and Raftery, A. E. (1995). Bayes factrs. Jurnal f the American Statistical Assciatin 90, 773{795. Kass, R. E., Tierney, L. and Kadane, J. B. (1990). The validity f psterir expansins based n Laplace's methd, in S. Geisser, J. S. Hdges, S. J. Press and A. Zellner (eds), Bayesian and Likelihd Methds in Statistics and Ecnmetrics: Essays in Hnr f Gerge A. Barnard, Amsterdam: Nrth-Hlland, 473{488. Lagaks, S. W. and Luis, T. A. (1988). Use f tumur lethality t interpret tumrigenicity experiments lacking cause-f-death data. Applied Statistics 37, 169{179. Lawless, J. F. (1982). Statistical Mdels and Methds fr Lifetime Data. New Yrk: Wiley. Smith, A. F. M. and Gelfand, A. E. (1992). Bayesian statistics withut tears: A samplingresampling perspective. The American Statistician 46, 84{88. Tierney, L. and Kadane, J. B. (1986). Accurate apprximatins fr psterir mments and marginal densities. Jurnal f the American Statistical Assciatin 81, 82{86. 7

8 lg H(t) S(t) lg t t lg beta lg beta lg theta1 lg theta2 Figure 2: Results fr ts t rat data (Lagaks and Luis, 1988). Upper left: hazard plt with tted single Weibull mdel shwn by the dtted line. Upper right: Kaplan{Meier estimate f survivr functin, (agged slid), with psterir mdal bi-weibull survivr functin btained using Laplace apprximatin (smth slid) and psterir expected bi-weibull survivr functin btained using sampling-imprtance resampling (dashes), and t f single Weibull (dts). Lwer panels: apprximate 50, 70, 90 and 95% cnturs fr the lg int psterir densities fr (lg 1; lg 1) (left) and fr (lg 2; lg 2) (right). 8

SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical model for microarray data analysis

SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical model for microarray data analysis SUPPLEMENTARY MATERIAL GaGa: a simple and flexible hierarchical mdel fr micrarray data analysis David Rssell Department f Bistatistics M.D. Andersn Cancer Center, Hustn, TX 77030, USA rsselldavid@gmail.cm

More information

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) >

Bootstrap Method > # Purpose: understand how bootstrap method works > obs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(obs) > Btstrap Methd > # Purpse: understand hw btstrap methd wrks > bs=c(11.96, 5.03, 67.40, 16.07, 31.50, 7.73, 11.10, 22.38) > n=length(bs) > mean(bs) [1] 21.64625 > # estimate f lambda > lambda = 1/mean(bs);

More information

Math Foundations 20 Work Plan

Math Foundations 20 Work Plan Math Fundatins 20 Wrk Plan Units / Tpics 20.8 Demnstrate understanding f systems f linear inequalities in tw variables. Time Frame December 1-3 weeks 6-10 Majr Learning Indicatrs Identify situatins relevant

More information

On Huntsberger Type Shrinkage Estimator for the Mean of Normal Distribution ABSTRACT INTRODUCTION

On Huntsberger Type Shrinkage Estimator for the Mean of Normal Distribution ABSTRACT INTRODUCTION Malaysian Jurnal f Mathematical Sciences 4(): 7-4 () On Huntsberger Type Shrinkage Estimatr fr the Mean f Nrmal Distributin Department f Mathematical and Physical Sciences, University f Nizwa, Sultanate

More information

, which yields. where z1. and z2

, which yields. where z1. and z2 The Gaussian r Nrmal PDF, Page 1 The Gaussian r Nrmal Prbability Density Functin Authr: Jhn M Cimbala, Penn State University Latest revisin: 11 September 13 The Gaussian r Nrmal Prbability Density Functin

More information

Differentiation Applications 1: Related Rates

Differentiation Applications 1: Related Rates Differentiatin Applicatins 1: Related Rates 151 Differentiatin Applicatins 1: Related Rates Mdel 1: Sliding Ladder 10 ladder y 10 ladder 10 ladder A 10 ft ladder is leaning against a wall when the bttm

More information

CAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank

CAUSAL INFERENCE. Technical Track Session I. Phillippe Leite. The World Bank CAUSAL INFERENCE Technical Track Sessin I Phillippe Leite The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Phillippe Leite fr the purpse f this wrkshp Plicy questins are causal

More information

CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS

CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS 1 Influential bservatins are bservatins whse presence in the data can have a distrting effect n the parameter estimates and pssibly the entire analysis,

More information

UNIV1"'RSITY OF NORTH CAROLINA Department of Statistics Chapel Hill, N. C. CUMULATIVE SUM CONTROL CHARTS FOR THE FOLDED NORMAL DISTRIBUTION

UNIV1'RSITY OF NORTH CAROLINA Department of Statistics Chapel Hill, N. C. CUMULATIVE SUM CONTROL CHARTS FOR THE FOLDED NORMAL DISTRIBUTION UNIV1"'RSITY OF NORTH CAROLINA Department f Statistics Chapel Hill, N. C. CUMULATIVE SUM CONTROL CHARTS FOR THE FOLDED NORMAL DISTRIBUTION by N. L. Jlmsn December 1962 Grant N. AFOSR -62..148 Methds f

More information

Admissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs

Admissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs Admissibility Cnditins and Asympttic Behavir f Strngly Regular Graphs VASCO MOÇO MANO Department f Mathematics University f Prt Oprt PORTUGAL vascmcman@gmailcm LUÍS ANTÓNIO DE ALMEIDA VIEIRA Department

More information

Eric Klein and Ning Sa

Eric Klein and Ning Sa Week 12. Statistical Appraches t Netwrks: p1 and p* Wasserman and Faust Chapter 15: Statistical Analysis f Single Relatinal Netwrks There are fur tasks in psitinal analysis: 1) Define Equivalence 2) Measure

More information

A Matrix Representation of Panel Data

A Matrix Representation of Panel Data web Extensin 6 Appendix 6.A A Matrix Representatin f Panel Data Panel data mdels cme in tw brad varieties, distinct intercept DGPs and errr cmpnent DGPs. his appendix presents matrix algebra representatins

More information

Resampling Methods. Chapter 5. Chapter 5 1 / 52

Resampling Methods. Chapter 5. Chapter 5 1 / 52 Resampling Methds Chapter 5 Chapter 5 1 / 52 1 51 Validatin set apprach 2 52 Crss validatin 3 53 Btstrap Chapter 5 2 / 52 Abut Resampling An imprtant statistical tl Pretending the data as ppulatin and

More information

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA

February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA February 28, 2013 COMMENTS ON DIFFUSION, DIFFUSIVITY AND DERIVATION OF HYPERBOLIC EQUATIONS DESCRIBING THE DIFFUSION PHENOMENA Mental Experiment regarding 1D randm walk Cnsider a cntainer f gas in thermal

More information

Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart

Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Key Wrds: Autregressive, Mving Average, Runs Tests, Shewhart Cntrl Chart Perfrmance f Sensitizing Rules n Shewhart Cntrl Charts with Autcrrelated Data Sandy D. Balkin Dennis K. J. Lin y Pennsylvania State University, University Park, PA 16802 Sandy Balkin is a graduate student

More information

Math Foundations 10 Work Plan

Math Foundations 10 Work Plan Math Fundatins 10 Wrk Plan Units / Tpics 10.1 Demnstrate understanding f factrs f whle numbers by: Prime factrs Greatest Cmmn Factrs (GCF) Least Cmmn Multiple (LCM) Principal square rt Cube rt Time Frame

More information

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion

22.54 Neutron Interactions and Applications (Spring 2004) Chapter 11 (3/11/04) Neutron Diffusion .54 Neutrn Interactins and Applicatins (Spring 004) Chapter (3//04) Neutrn Diffusin References -- J. R. Lamarsh, Intrductin t Nuclear Reactr Thery (Addisn-Wesley, Reading, 966) T study neutrn diffusin

More information

Lecture 13: Markov Chain Monte Carlo. Gibbs sampling

Lecture 13: Markov Chain Monte Carlo. Gibbs sampling Lecture 13: Markv hain Mnte arl Gibbs sampling Gibbs sampling Markv chains 1 Recall: Apprximate inference using samples Main idea: we generate samples frm ur Bayes net, then cmpute prbabilities using (weighted)

More information

AP Statistics Practice Test Unit Three Exploring Relationships Between Variables. Name Period Date

AP Statistics Practice Test Unit Three Exploring Relationships Between Variables. Name Period Date AP Statistics Practice Test Unit Three Explring Relatinships Between Variables Name Perid Date True r False: 1. Crrelatin and regressin require explanatry and respnse variables. 1. 2. Every least squares

More information

Distributions, spatial statistics and a Bayesian perspective

Distributions, spatial statistics and a Bayesian perspective Distributins, spatial statistics and a Bayesian perspective Dug Nychka Natinal Center fr Atmspheric Research Distributins and densities Cnditinal distributins and Bayes Thm Bivariate nrmal Spatial statistics

More information

BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS. Christopher Costello, Andrew Solow, Michael Neubert, and Stephen Polasky

BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS. Christopher Costello, Andrew Solow, Michael Neubert, and Stephen Polasky BOUNDED UNCERTAINTY AND CLIMATE CHANGE ECONOMICS Christpher Cstell, Andrew Slw, Michael Neubert, and Stephen Plasky Intrductin The central questin in the ecnmic analysis f climate change plicy cncerns

More information

GMM with Latent Variables

GMM with Latent Variables GMM with Latent Variables A. Rnald Gallant Penn State University Raffaella Giacmini University Cllege Lndn Giuseppe Ragusa Luiss University Cntributin The cntributin f GMM (Hansen and Singletn, 1982) was

More information

NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION

NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION NUROP Chinese Pinyin T Chinese Character Cnversin NUROP CONGRESS PAPER CHINESE PINYIN TO CHINESE CHARACTER CONVERSION CHIA LI SHI 1 AND LUA KIM TENG 2 Schl f Cmputing, Natinal University f Singapre 3 Science

More information

initially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur

initially lcated away frm the data set never win the cmpetitin, resulting in a nnptimal nal cdebk, [2] [3] [4] and [5]. Khnen's Self Organizing Featur Cdewrd Distributin fr Frequency Sensitive Cmpetitive Learning with One Dimensinal Input Data Aristides S. Galanpuls and Stanley C. Ahalt Department f Electrical Engineering The Ohi State University Abstract

More information

MATHEMATICS SYLLABUS SECONDARY 5th YEAR

MATHEMATICS SYLLABUS SECONDARY 5th YEAR Eurpean Schls Office f the Secretary-General Pedaggical Develpment Unit Ref. : 011-01-D-8-en- Orig. : EN MATHEMATICS SYLLABUS SECONDARY 5th YEAR 6 perid/week curse APPROVED BY THE JOINT TEACHING COMMITTEE

More information

I. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is

I. Analytical Potential and Field of a Uniform Rod. V E d. The definition of electric potential difference is Length L>>a,b,c Phys 232 Lab 4 Ch 17 Electric Ptential Difference Materials: whitebards & pens, cmputers with VPythn, pwer supply & cables, multimeter, crkbard, thumbtacks, individual prbes and jined prbes,

More information

Lead/Lag Compensator Frequency Domain Properties and Design Methods

Lead/Lag Compensator Frequency Domain Properties and Design Methods Lectures 6 and 7 Lead/Lag Cmpensatr Frequency Dmain Prperties and Design Methds Definitin Cnsider the cmpensatr (ie cntrller Fr, it is called a lag cmpensatr s K Fr s, it is called a lead cmpensatr Ntatin

More information

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION U. S. FOREST SERVICE RESEARCH PAPER FPL 50 DECEMBER U. S. DEPARTMENT OF AGRICULTURE FOREST SERVICE FOREST PRODUCTS LABORATORY OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

More information

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA

Modelling of Clock Behaviour. Don Percival. Applied Physics Laboratory University of Washington Seattle, Washington, USA Mdelling f Clck Behaviur Dn Percival Applied Physics Labratry University f Washingtn Seattle, Washingtn, USA verheads and paper fr talk available at http://faculty.washingtn.edu/dbp/talks.html 1 Overview

More information

ON-LINE PROCEDURE FOR TERMINATING AN ACCELERATED DEGRADATION TEST

ON-LINE PROCEDURE FOR TERMINATING AN ACCELERATED DEGRADATION TEST Statistica Sinica 8(1998), 207-220 ON-LINE PROCEDURE FOR TERMINATING AN ACCELERATED DEGRADATION TEST Hng-Fwu Yu and Sheng-Tsaing Tseng Natinal Taiwan University f Science and Technlgy and Natinal Tsing-Hua

More information

CS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007

CS 477/677 Analysis of Algorithms Fall 2007 Dr. George Bebis Course Project Due Date: 11/29/2007 CS 477/677 Analysis f Algrithms Fall 2007 Dr. Gerge Bebis Curse Prject Due Date: 11/29/2007 Part1: Cmparisn f Srting Algrithms (70% f the prject grade) The bjective f the first part f the assignment is

More information

IN a recent article, Geary [1972] discussed the merit of taking first differences

IN a recent article, Geary [1972] discussed the merit of taking first differences The Efficiency f Taking First Differences in Regressin Analysis: A Nte J. A. TILLMAN IN a recent article, Geary [1972] discussed the merit f taking first differences t deal with the prblems that trends

More information

Lecture 17: Free Energy of Multi-phase Solutions at Equilibrium

Lecture 17: Free Energy of Multi-phase Solutions at Equilibrium Lecture 17: 11.07.05 Free Energy f Multi-phase Slutins at Equilibrium Tday: LAST TIME...2 FREE ENERGY DIAGRAMS OF MULTI-PHASE SOLUTIONS 1...3 The cmmn tangent cnstructin and the lever rule...3 Practical

More information

Determining the Accuracy of Modal Parameter Estimation Methods

Determining the Accuracy of Modal Parameter Estimation Methods Determining the Accuracy f Mdal Parameter Estimatin Methds by Michael Lee Ph.D., P.E. & Mar Richardsn Ph.D. Structural Measurement Systems Milpitas, CA Abstract The mst cmmn type f mdal testing system

More information

COMP 551 Applied Machine Learning Lecture 5: Generative models for linear classification

COMP 551 Applied Machine Learning Lecture 5: Generative models for linear classification COMP 551 Applied Machine Learning Lecture 5: Generative mdels fr linear classificatin Instructr: Herke van Hf (herke.vanhf@mail.mcgill.ca) Slides mstly by: Jelle Pineau Class web page: www.cs.mcgill.ca/~hvanh2/cmp551

More information

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India

CHAPTER 3 INEQUALITIES. Copyright -The Institute of Chartered Accountants of India CHAPTER 3 INEQUALITIES Cpyright -The Institute f Chartered Accuntants f India INEQUALITIES LEARNING OBJECTIVES One f the widely used decisin making prblems, nwadays, is t decide n the ptimal mix f scarce

More information

AP Statistics Notes Unit Two: The Normal Distributions

AP Statistics Notes Unit Two: The Normal Distributions AP Statistics Ntes Unit Tw: The Nrmal Distributins Syllabus Objectives: 1.5 The student will summarize distributins f data measuring the psitin using quartiles, percentiles, and standardized scres (z-scres).

More information

Part 3 Introduction to statistical classification techniques

Part 3 Introduction to statistical classification techniques Part 3 Intrductin t statistical classificatin techniques Machine Learning, Part 3, March 07 Fabi Rli Preamble ØIn Part we have seen that if we knw: Psterir prbabilities P(ω i / ) Or the equivalent terms

More information

Pressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects

Pressure And Entropy Variations Across The Weak Shock Wave Due To Viscosity Effects Pressure And Entrpy Variatins Acrss The Weak Shck Wave Due T Viscsity Effects OSTAFA A. A. AHOUD Department f athematics Faculty f Science Benha University 13518 Benha EGYPT Abstract:-The nnlinear differential

More information

4th Indian Institute of Astrophysics - PennState Astrostatistics School July, 2013 Vainu Bappu Observatory, Kavalur. Correlation and Regression

4th Indian Institute of Astrophysics - PennState Astrostatistics School July, 2013 Vainu Bappu Observatory, Kavalur. Correlation and Regression 4th Indian Institute f Astrphysics - PennState Astrstatistics Schl July, 2013 Vainu Bappu Observatry, Kavalur Crrelatin and Regressin Rahul Ry Indian Statistical Institute, Delhi. Crrelatin Cnsider a tw

More information

Bayesian Analysis of Normal Scalar Populations with a Common Mean

Bayesian Analysis of Normal Scalar Populations with a Common Mean Pertania (),1 - (1) Bayesian Analysis f Nrmal Scalar Ppulatins with a Cmmn Mean MAT YUSOFF ABDULLAH Mathematics Department, Faculty f Science and Envirnmental Studies, Universiti Pertanian Malaysia, 00

More information

Bayesian nonparametric modeling approaches for quantile regression

Bayesian nonparametric modeling approaches for quantile regression Bayesian nnparametric mdeling appraches fr quantile regressin Athanasis Kttas Department f Applied Mathematics and Statistics University f Califrnia, Santa Cruz Department f Statistics Athens University

More information

Optimization Programming Problems For Control And Management Of Bacterial Disease With Two Stage Growth/Spread Among Plants

Optimization Programming Problems For Control And Management Of Bacterial Disease With Two Stage Growth/Spread Among Plants Internatinal Jurnal f Engineering Science Inventin ISSN (Online): 9 67, ISSN (Print): 9 676 www.ijesi.rg Vlume 5 Issue 8 ugust 06 PP.0-07 Optimizatin Prgramming Prblems Fr Cntrl nd Management Of Bacterial

More information

Sequential Allocation with Minimal Switching

Sequential Allocation with Minimal Switching In Cmputing Science and Statistics 28 (1996), pp. 567 572 Sequential Allcatin with Minimal Switching Quentin F. Stut 1 Janis Hardwick 1 EECS Dept., University f Michigan Statistics Dept., Purdue University

More information

ENSC Discrete Time Systems. Project Outline. Semester

ENSC Discrete Time Systems. Project Outline. Semester ENSC 49 - iscrete Time Systems Prject Outline Semester 006-1. Objectives The gal f the prject is t design a channel fading simulatr. Upn successful cmpletin f the prject, yu will reinfrce yur understanding

More information

Resampling Methods. Cross-validation, Bootstrapping. Marek Petrik 2/21/2017

Resampling Methods. Cross-validation, Bootstrapping. Marek Petrik 2/21/2017 Resampling Methds Crss-validatin, Btstrapping Marek Petrik 2/21/2017 Sme f the figures in this presentatin are taken frm An Intrductin t Statistical Learning, with applicatins in R (Springer, 2013) with

More information

arxiv:hep-ph/ v1 2 Jun 1995

arxiv:hep-ph/ v1 2 Jun 1995 WIS-95//May-PH The rati F n /F p frm the analysis f data using a new scaling variable S. A. Gurvitz arxiv:hep-ph/95063v1 Jun 1995 Department f Particle Physics, Weizmann Institute f Science, Rehvt 76100,

More information

AMERICAN PETROLEUM INSTITUTE API RP 581 RISK BASED INSPECTION BASE RESOURCE DOCUMENT BALLOT COVER PAGE

AMERICAN PETROLEUM INSTITUTE API RP 581 RISK BASED INSPECTION BASE RESOURCE DOCUMENT BALLOT COVER PAGE Ballt ID: Title: USING LIFE EXTENSION FACTOR (LEF) TO INCREASE BUNDLE INSPECTION INTERVAL Purpse: 1. Prvides a methd t increase a bundle s inspectin interval t accunt fr LEF. 2. Clarifies Table 8.6.5 Als

More information

BASD HIGH SCHOOL FORMAL LAB REPORT

BASD HIGH SCHOOL FORMAL LAB REPORT BASD HIGH SCHOOL FORMAL LAB REPORT *WARNING: After an explanatin f what t include in each sectin, there is an example f hw the sectin might lk using a sample experiment Keep in mind, the sample lab used

More information

The Kullback-Leibler Kernel as a Framework for Discriminant and Localized Representations for Visual Recognition

The Kullback-Leibler Kernel as a Framework for Discriminant and Localized Representations for Visual Recognition The Kullback-Leibler Kernel as a Framewrk fr Discriminant and Lcalized Representatins fr Visual Recgnitin Nun Vascncels Purdy H Pedr Mren ECE Department University f Califrnia, San Dieg HP Labs Cambridge

More information

Dr M. BROUARD. 5. Thermodynamic formulation of Transition State Theory Entropy of activation. Thermochemical kinetics. CHEMICAL REACTION RATES

Dr M. BROUARD. 5. Thermodynamic formulation of Transition State Theory Entropy of activation. Thermochemical kinetics. CHEMICAL REACTION RATES CHEMICAL REACTION RATES Dr M. BROUARD Trinity Term 2003 A. Bimlecular Reactins 5 Lectures 1. Intrductin Simple cllisin thery. Ptential energy curves and surfaces. The reactin crdinate and barriers t reactin.

More information

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax

Revision: August 19, E Main Suite D Pullman, WA (509) Voice and Fax .7.4: Direct frequency dmain circuit analysis Revisin: August 9, 00 5 E Main Suite D Pullman, WA 9963 (509) 334 6306 ice and Fax Overview n chapter.7., we determined the steadystate respnse f electrical

More information

ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES

ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES ANALYTICAL SOLUTIONS TO THE PROBLEM OF EDDY CURRENT PROBES CONSISTING OF LONG PARALLEL CONDUCTORS B. de Halleux, O. Lesage, C. Mertes and A. Ptchelintsev Mechanical Engineering Department Cathlic University

More information

The blessing of dimensionality for kernel methods

The blessing of dimensionality for kernel methods fr kernel methds Building classifiers in high dimensinal space Pierre Dupnt Pierre.Dupnt@ucluvain.be Classifiers define decisin surfaces in sme feature space where the data is either initially represented

More information

MATCHING TECHNIQUES. Technical Track Session VI. Emanuela Galasso. The World Bank

MATCHING TECHNIQUES. Technical Track Session VI. Emanuela Galasso. The World Bank MATCHING TECHNIQUES Technical Track Sessin VI Emanuela Galass The Wrld Bank These slides were develped by Christel Vermeersch and mdified by Emanuela Galass fr the purpse f this wrkshp When can we use

More information

Modelling Record Times in Sport with Extreme Value Methods

Modelling Record Times in Sport with Extreme Value Methods utput 2017/5/9 5:13 page 1 #1 Malaysian Jurnal f Mathematical Sciences *(*): ** ** (201*) MALAYSIAN JOURNAL OF MATHEMATICAL SCIENCES Jurnal hmepage: http://einspem.upm.edu.my/jurnal Mdelling Recrd Times

More information

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322

ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA. December 4, PLP No. 322 ELECTRON CYCLOTRON HEATING OF AN ANISOTROPIC PLASMA by J. C. SPROTT December 4, 1969 PLP N. 3 These PLP Reprts are infrmal and preliminary and as such may cntain errrs nt yet eliminated. They are fr private

More information

5.4 Measurement Sampling Rates for Daily Maximum and Minimum Temperatures

5.4 Measurement Sampling Rates for Daily Maximum and Minimum Temperatures 5.4 Measurement Sampling Rates fr Daily Maximum and Minimum Temperatures 1 1 2 X. Lin, K. G. Hubbard, and C. B. Baker University f Nebraska, Lincln, Nebraska 2 Natinal Climatic Data Center 1 1. INTRODUCTION

More information

SIZE BIAS IN LINE TRANSECT SAMPLING: A FIELD TEST. Mark C. Otto Statistics Research Division, Bureau of the Census Washington, D.C , U.S.A.

SIZE BIAS IN LINE TRANSECT SAMPLING: A FIELD TEST. Mark C. Otto Statistics Research Division, Bureau of the Census Washington, D.C , U.S.A. SIZE BIAS IN LINE TRANSECT SAMPLING: A FIELD TEST Mark C. Ott Statistics Research Divisin, Bureau f the Census Washingtn, D.C. 20233, U.S.A. and Kenneth H. Pllck Department f Statistics, Nrth Carlina State

More information

Kinetic Model Completeness

Kinetic Model Completeness 5.68J/10.652J Spring 2003 Lecture Ntes Tuesday April 15, 2003 Kinetic Mdel Cmpleteness We say a chemical kinetic mdel is cmplete fr a particular reactin cnditin when it cntains all the species and reactins

More information

SURVIVAL ANALYSIS WITH SUPPORT VECTOR MACHINES

SURVIVAL ANALYSIS WITH SUPPORT VECTOR MACHINES 1 SURVIVAL ANALYSIS WITH SUPPORT VECTOR MACHINES Wlfgang HÄRDLE Ruslan MORO Center fr Applied Statistics and Ecnmics (CASE), Humbldt-Universität zu Berlin Mtivatin 2 Applicatins in Medicine estimatin f

More information

Least Squares Optimal Filtering with Multirate Observations

Least Squares Optimal Filtering with Multirate Observations Prc. 36th Asilmar Cnf. n Signals, Systems, and Cmputers, Pacific Grve, CA, Nvember 2002 Least Squares Optimal Filtering with Multirate Observatins Charles W. herrien and Anthny H. Hawes Department f Electrical

More information

Theoretical study of third virial coefficient with Kihara potential

Theoretical study of third virial coefficient with Kihara potential Theretical study f third virial cefficient with Kihara ptential Jurnal: Manuscript ID cjp-017-0705.r Manuscript Type: Article Date Submitted by the Authr: 6-Dec-017 Cmplete List f Authrs: Smuncu E.; Giresun

More information

Five Whys How To Do It Better

Five Whys How To Do It Better Five Whys Definitin. As explained in the previus article, we define rt cause as simply the uncvering f hw the current prblem came int being. Fr a simple causal chain, it is the entire chain. Fr a cmplex

More information

Computational modeling techniques

Computational modeling techniques Cmputatinal mdeling techniques Lecture 4: Mdel checing fr ODE mdels In Petre Department f IT, Åb Aademi http://www.users.ab.fi/ipetre/cmpmd/ Cntent Stichimetric matrix Calculating the mass cnservatin relatins

More information

Verification of Quality Parameters of a Solar Panel and Modification in Formulae of its Series Resistance

Verification of Quality Parameters of a Solar Panel and Modification in Formulae of its Series Resistance Verificatin f Quality Parameters f a Slar Panel and Mdificatin in Frmulae f its Series Resistance Sanika Gawhane Pune-411037-India Onkar Hule Pune-411037- India Chinmy Kulkarni Pune-411037-India Ojas Pandav

More information

[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y )

[COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t o m a k e s u r e y o u a r e r e a d y ) (Abut the final) [COLLEGE ALGEBRA EXAM I REVIEW TOPICS] ( u s e t h i s t m a k e s u r e y u a r e r e a d y ) The department writes the final exam s I dn't really knw what's n it and I can't very well

More information

You need to be able to define the following terms and answer basic questions about them:

You need to be able to define the following terms and answer basic questions about them: CS440/ECE448 Sectin Q Fall 2017 Midterm Review Yu need t be able t define the fllwing terms and answer basic questins abut them: Intr t AI, agents and envirnments Pssible definitins f AI, prs and cns f

More information

Internal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9.

Internal vs. external validity. External validity. This section is based on Stock and Watson s Chapter 9. Sectin 7 Mdel Assessment This sectin is based n Stck and Watsn s Chapter 9. Internal vs. external validity Internal validity refers t whether the analysis is valid fr the ppulatin and sample being studied.

More information

A New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation

A New Evaluation Measure. J. Joiner and L. Werner. The problems of evaluation and the needed criteria of evaluation III-l III. A New Evaluatin Measure J. Jiner and L. Werner Abstract The prblems f evaluatin and the needed criteria f evaluatin measures in the SMART system f infrmatin retrieval are reviewed and discussed.

More information

Slide04 (supplemental) Haykin Chapter 4 (both 2nd and 3rd ed): Multi-Layer Perceptrons

Slide04 (supplemental) Haykin Chapter 4 (both 2nd and 3rd ed): Multi-Layer Perceptrons Slide04 supplemental) Haykin Chapter 4 bth 2nd and 3rd ed): Multi-Layer Perceptrns CPSC 636-600 Instructr: Ynsuck Che Heuristic fr Making Backprp Perfrm Better 1. Sequential vs. batch update: fr large

More information

Introduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem

Introduction: A Generalized approach for computing the trajectories associated with the Newtonian N Body Problem A Generalized apprach fr cmputing the trajectries assciated with the Newtnian N Bdy Prblem AbuBar Mehmd, Syed Umer Abbas Shah and Ghulam Shabbir Faculty f Engineering Sciences, GIK Institute f Engineering

More information

Maximum A Posteriori (MAP) CS 109 Lecture 22 May 16th, 2016

Maximum A Posteriori (MAP) CS 109 Lecture 22 May 16th, 2016 Maximum A Psteriri (MAP) CS 109 Lecture 22 May 16th, 2016 Previusly in CS109 Game f Estimatrs Maximum Likelihd Nn spiler: this didn t happen Side Plt argmax argmax f lg Mther f ptimizatins? Reviving an

More information

1996 Engineering Systems Design and Analysis Conference, Montpellier, France, July 1-4, 1996, Vol. 7, pp

1996 Engineering Systems Design and Analysis Conference, Montpellier, France, July 1-4, 1996, Vol. 7, pp THE POWER AND LIMIT OF NEURAL NETWORKS T. Y. Lin Department f Mathematics and Cmputer Science San Jse State University San Jse, Califrnia 959-003 tylin@cs.ssu.edu and Bereley Initiative in Sft Cmputing*

More information

5 th grade Common Core Standards

5 th grade Common Core Standards 5 th grade Cmmn Cre Standards In Grade 5, instructinal time shuld fcus n three critical areas: (1) develping fluency with additin and subtractin f fractins, and develping understanding f the multiplicatin

More information

Hierarchical priors and mixture models, with applications in regression and density estimation

Hierarchical priors and mixture models, with applications in regression and density estimation Hierarchical prirs and mixture mdels, with applicatins in regressin and density estimatin By Mike West, Peter Müller & Michael D Escbar In: Aspects f Uncertainty: A Tribute t D.V. Lindley (eds: P Freeman

More information

Fall 2013 Physics 172 Recitation 3 Momentum and Springs

Fall 2013 Physics 172 Recitation 3 Momentum and Springs Fall 03 Physics 7 Recitatin 3 Mmentum and Springs Purpse: The purpse f this recitatin is t give yu experience wrking with mmentum and the mmentum update frmula. Readings: Chapter.3-.5 Learning Objectives:.3.

More information

ALE 21. Gibbs Free Energy. At what temperature does the spontaneity of a reaction change?

ALE 21. Gibbs Free Energy. At what temperature does the spontaneity of a reaction change? Name Chem 163 Sectin: Team Number: ALE 21. Gibbs Free Energy (Reference: 20.3 Silberberg 5 th editin) At what temperature des the spntaneity f a reactin change? The Mdel: The Definitin f Free Energy S

More information

Computational modeling techniques

Computational modeling techniques Cmputatinal mdeling techniques Lecture 2: Mdeling change. In Petre Department f IT, Åb Akademi http://users.ab.fi/ipetre/cmpmd/ Cntent f the lecture Basic paradigm f mdeling change Examples Linear dynamical

More information

SAMPLING DYNAMICAL SYSTEMS

SAMPLING DYNAMICAL SYSTEMS SAMPLING DYNAMICAL SYSTEMS Melvin J. Hinich Applied Research Labratries The University f Texas at Austin Austin, TX 78713-8029, USA (512) 835-3278 (Vice) 835-3259 (Fax) hinich@mail.la.utexas.edu ABSTRACT

More information

NUMBERS, MATHEMATICS AND EQUATIONS

NUMBERS, MATHEMATICS AND EQUATIONS AUSTRALIAN CURRICULUM PHYSICS GETTING STARTED WITH PHYSICS NUMBERS, MATHEMATICS AND EQUATIONS An integral part t the understanding f ur physical wrld is the use f mathematical mdels which can be used t

More information

11. DUAL NATURE OF RADIATION AND MATTER

11. DUAL NATURE OF RADIATION AND MATTER 11. DUAL NATURE OF RADIATION AND MATTER Very shrt answer and shrt answer questins 1. Define wrk functin f a metal? The minimum energy required fr an electrn t escape frm the metal surface is called the

More information

On Boussinesq's problem

On Boussinesq's problem Internatinal Jurnal f Engineering Science 39 (2001) 317±322 www.elsevier.cm/lcate/ijengsci On Bussinesq's prblem A.P.S. Selvadurai * Department f Civil Engineering and Applied Mechanics, McGill University,

More information

Rangely RE 4 Curriculum Development 5 th Grade Mathematics

Rangely RE 4 Curriculum Development 5 th Grade Mathematics Unit Title Dctr We Still Need t Operate... Length f Unit 12 weeks Fcusing Lens(es) Inquiry Questins (Engaging Debatable): Structure Systems Standards and Grade Level Expectatins Addressed in this Unit

More information

Methods for Determination of Mean Speckle Size in Simulated Speckle Pattern

Methods for Determination of Mean Speckle Size in Simulated Speckle Pattern 0.478/msr-04-004 MEASUREMENT SCENCE REVEW, Vlume 4, N. 3, 04 Methds fr Determinatin f Mean Speckle Size in Simulated Speckle Pattern. Hamarvá, P. Šmíd, P. Hrváth, M. Hrabvský nstitute f Physics f the Academy

More information

EDA Engineering Design & Analysis Ltd

EDA Engineering Design & Analysis Ltd EDA Engineering Design & Analysis Ltd THE FINITE ELEMENT METHOD A shrt tutrial giving an verview f the histry, thery and applicatin f the finite element methd. Intrductin Value f FEM Applicatins Elements

More information

Chemistry 20 Lesson 11 Electronegativity, Polarity and Shapes

Chemistry 20 Lesson 11 Electronegativity, Polarity and Shapes Chemistry 20 Lessn 11 Electrnegativity, Plarity and Shapes In ur previus wrk we learned why atms frm cvalent bnds and hw t draw the resulting rganizatin f atms. In this lessn we will learn (a) hw the cmbinatin

More information

CHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came.

CHAPTER 24: INFERENCE IN REGRESSION. Chapter 24: Make inferences about the population from which the sample data came. MATH 1342 Ch. 24 April 25 and 27, 2013 Page 1 f 5 CHAPTER 24: INFERENCE IN REGRESSION Chapters 4 and 5: Relatinships between tw quantitative variables. Be able t Make a graph (scatterplt) Summarize the

More information

( ) kt. Solution. From kinetic theory (visualized in Figure 1Q9-1), 1 2 rms = 2. = 1368 m/s

( ) kt. Solution. From kinetic theory (visualized in Figure 1Q9-1), 1 2 rms = 2. = 1368 m/s .9 Kinetic Mlecular Thery Calculate the effective (rms) speeds f the He and Ne atms in the He-Ne gas laser tube at rm temperature (300 K). Slutin T find the rt mean square velcity (v rms ) f He atms at

More information

The Law of Total Probability, Bayes Rule, and Random Variables (Oh My!)

The Law of Total Probability, Bayes Rule, and Random Variables (Oh My!) The Law f Ttal Prbability, Bayes Rule, and Randm Variables (Oh My!) Administrivia Hmewrk 2 is psted and is due tw Friday s frm nw If yu didn t start early last time, please d s this time. Gd Milestnes:

More information

THE LIFE OF AN OBJECT IT SYSTEMS

THE LIFE OF AN OBJECT IT SYSTEMS THE LIFE OF AN OBJECT IT SYSTEMS Persns, bjects, r cncepts frm the real wrld, which we mdel as bjects in the IT system, have "lives". Actually, they have tw lives; the riginal in the real wrld has a life,

More information

Checking the resolved resonance region in EXFOR database

Checking the resolved resonance region in EXFOR database Checking the reslved resnance regin in EXFOR database Gttfried Bertn Sciété de Calcul Mathématique (SCM) Oscar Cabells OECD/NEA Data Bank JEFF Meetings - Sessin JEFF Experiments Nvember 0-4, 017 Bulgne-Billancurt,

More information

A mathematical model for complete stress-strain curve prediction of permeable concrete

A mathematical model for complete stress-strain curve prediction of permeable concrete A mathematical mdel fr cmplete stress-strain curve predictin f permeable cncrete M. K. Hussin Y. Zhuge F. Bullen W. P. Lkuge Faculty f Engineering and Surveying, University f Suthern Queensland, Twmba,

More information

Lyapunov Stability Stability of Equilibrium Points

Lyapunov Stability Stability of Equilibrium Points Lyapunv Stability Stability f Equilibrium Pints 1. Stability f Equilibrium Pints - Definitins In this sectin we cnsider n-th rder nnlinear time varying cntinuus time (C) systems f the frm x = f ( t, x),

More information

Dead-beat controller design

Dead-beat controller design J. Hetthéssy, A. Barta, R. Bars: Dead beat cntrller design Nvember, 4 Dead-beat cntrller design In sampled data cntrl systems the cntrller is realised by an intelligent device, typically by a PLC (Prgrammable

More information

MATCHING TECHNIQUES Technical Track Session VI Céline Ferré The World Bank

MATCHING TECHNIQUES Technical Track Session VI Céline Ferré The World Bank MATCHING TECHNIQUES Technical Track Sessin VI Céline Ferré The Wrld Bank When can we use matching? What if the assignment t the treatment is nt dne randmly r based n an eligibility index, but n the basis

More information

A RATIONALE FOR AN ASYMPTOTIC LOGNORMAL FORM OF WORD-FREQUENCY DISTRIBUTIONS. John B. Carroll

A RATIONALE FOR AN ASYMPTOTIC LOGNORMAL FORM OF WORD-FREQUENCY DISTRIBUTIONS. John B. Carroll R ~ ~8 R L ~ ~ TI A RATIONALE FOR AN ASYMPTOTIC LOGNORMAL FORM OF WORD-FREQUENCY DISTRIBUTIONS Jhn B. Carrll RB-69-9 N This Bulletin is a draft fr interffice circulatin. Crrectins and suggestins fr revisin

More information

Snow avalanche runout from two Canadian mountain ranges

Snow avalanche runout from two Canadian mountain ranges Annals f Glacilgy 18 1993 Internati n al Glaciigicai Sciety Snw avalanche runut frm tw Canadian muntain ranges D.J. NIXON AND D. M. MCCLUNG Department f Civil Engineering, University f British Clumbia,

More information

Drought damaged area

Drought damaged area ESTIMATE OF THE AMOUNT OF GRAVEL CO~TENT IN THE SOIL BY A I R B O'RN EMS S D A T A Y. GOMI, H. YAMAMOTO, AND S. SATO ASIA AIR SURVEY CO., l d. KANAGAWA,JAPAN S.ISHIGURO HOKKAIDO TOKACHI UBPREFECTRAl OffICE

More information

ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS

ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS Particle Acceleratrs, 1986, Vl. 19, pp. 99-105 0031-2460/86/1904-0099/$15.00/0 1986 Grdn and Breach, Science Publishers, S.A. Printed in the United States f America ROUNDING ERRORS IN BEAM-TRACKING CALCULATIONS

More information