ON THE WORTH OF GROUNDWATER CONTAMINANT SAMPLING

Size: px
Start display at page:

Download "ON THE WORTH OF GROUNDWATER CONTAMINANT SAMPLING"

Transcription

1 ON THE WORTH OF GROUNDWATER CONTAMINANT SAMPLING 1. No-spatial case Lars Rosé Departmet of Geology Chalmers Uiversity of Techology SE Gotheburg, SWEDEN Tommy Norberg Departmet of Mathematical Statistics Chalmers Uiversity of Techology SE Gotheburg, SWEDEN

2 2 Abstract This is the first i a series of three papers studyig the ecoomical worth of groudwater cotamiat samplig i applicatios of icreased complexity. This paper presets the theoretical basis for Bayesia data-worth aalysis ad discusses the o-spatial case with oe hydraulically isolated cell of uiform ad isotropic coditios. The two subsequet papers will discuss data worth i two oe-dimesioal models ad i a two-dimesioal spatial settig. The mai cotributio of this paper is that it provides a ew theory for estimatig the worth of simultaeous samplig i more tha oe locatio, which is a especially importat aspect i cases of large samplig errors. Data worth depeds directly o the prior probability of cotamiatio ad we compare a case i which this probability is completely kow with a case i which it is kow with give ucertaity. I the latter case, we show that the data worth is sesitive with respect to the selected desig criteria. Itroductio Cotamiated soil ad groudwater is a problem of growig cocer i our society ad is today a major issue i lad use plaig ad maagemet, real estate assessmet, ad property sellig. Ivestigatio ad remediatio of cotamiated areas are ofte associated with high costs. As a example, the Swedish Evirometal Protectio Agecy (EPA) provides more tha the equivalet of $50 millio aually of govermetal resources to the Swedish couty authorities for ivestigatio ad remediatio of cotamiated areas where o resposible part ca be foud. Still, the bulk of ecessary ivestmets for ivestigatio ad remediatio come from resposible ladowers ad operators. The Swedish EPA (1999)

3 3 estimates that there are cotamiated sites i Swede, of which approximately 4000 are i eed of remediatio. Similar, or eve more severe, situatios are preset throughout Europe ad North America. The Polluter Pays Priciple, which is the regulatory philosophy i the Europea Uio, puts large fiacial pressure o resposible compaies ad authorities. The regulatory framework is therefore a strog icetive for both the public ad private compaies for applyig costefficiet ivestigatio ad remediatio strategies. I additio, the evirometal legislatio i may coutries, e.g. i Swede, state that the evirometal value of remediatio must be i proportio to the remediatio costs i order to provide proper prioritizatio of resources ad to achieve sustaiable use of lad ad water. Due to a combiatio of complex geological, hydrogeological, ad geochemical coditios at cotamiated sites ad high ivestigatio costs, it is usually ot possible to obtai complete iformatio or characterize a site with a high degree of certaity. Ivestigatios of cotamiated areas are therefore typically associated with large ucertaities regardig e.g. type ad extet of cotamiatio ad possible future cotamiat spreadig. These ucertaities trasform ito substatial ecoomical risks i the remediatio process. To hadle these risks properly ad to achieve a cost-efficiet maagemet of cotamiated sites, it is therefore of primary importace to desig ivestigatio strategies cosiderig both the samplig costs ad the ecoomical risks associated with makig icorrect decisios regardig remediatio.

4 4 The umber ad type of ew data i cotamiatio studies are ofte estimated with respect to the potetial for reductio of ucertaity. However, i order to be cost-effective the worth of ew data should be evaluated with respect to the reductio of the risk-cost of makig icorrect decisios, rather tha ucertaity. I a cost-beefit perspective, the worth of samplig depeds o the samplig cost ad the risk-cost reductio it provides with respect to failure of successfully solvig the specific problem at had. The risk-cost must the icorporate all expected costs for the decisio-maker related to improper evaluatio of the situatio, e.g. chage of remediatio strategy, prologed clea-up times, restrictios o lad-use, ad loss of evirometal values. A tool for achievig cost-effective ivestigatio is Bayesia data-worth aalysis, evolvig from a ecoomical decisio aalysis framework. Bayesia data-worth aalysis i hydrogeology has bee described by e.g. Freeze et al (1990; 1992), James & Freeze (1993) ad James & Gorelick (1994). The basic philosophy behid this approach i remediatio projects is to itegrate the cost of samplig, the reductio of the risk-cost of failure, ad the cost ad efficiecy of the remedial desig i order to fid the best samplig strategy amog a give set of alterative strategies. The data-worth aalysis is made iteratively accordig to Figure 1 util o further samplig is cosidered justifiable. Our work will be preseted i a series of three papers of which this is the first. The mai purpose is two-fold: (1) to aalyze ad describe how ecoomical factors, samplig ucertaity, ad prior assumptios regardig cotamiatio affect cost-efficiecy of the samplig strategy ad (2) to preset a spatial statistical model for Bayesia data-worth aalysis. We focus o the o-spatial aspects of data-worth i the first paper, 1-dimesioal aspects i the secod paper, ad 2-dimesioal aspects i the third paper. The fist two papers i this series are directed at aalyzig ecoomical ad statistical implicatios of data-worth

5 5 aalysis i a sigle cell 1 ad alog a lie of several cells. We preset coclusios ad give recommedatios regardig a umber of strategic decisios to be take whe applyig Bayesia data-worth aalysis. I the third paper we preset a 2D spatial statistical data-worth model, based o Mote Carlo Markov Chai (MCMC) aalysis. To fully uderstad ad appreciate the capabilities ad limitatios of Bayesia data-worth aalysis, implicatios o lower dimesioalities eed to be studied ad evaluated before full 2D ad/or 3D applicatios ca be made correctly. We believe that the work preseted i these three papers provide importat kowledge for a icreased cost-efficiecy i projects o ivestigatio ad remediatio of cotamiated soil ad groudwater. Geeral descriptio of Bayesia data-worth aalysis We state that the worth of samplig data ca be best evaluated by the use of Bayesia decisio aalysis, which has bee applied for cotamiated soil ad groudwater by several authors, e.g. Massma & Freeze (1987), Mari et al (1989), Freeze et al (1990, 1992), Massma et al (1991) James & Freeze (1993), James & Gorelick (1994), James et al (1996), Rosé & Wladis (1998), Wladis et al (1999), Norrma (2001), Back (2001) ad Rosé (2001). A special applicatio of Bayesia decisio aalysis is referred to as Bayesia dataworth aalysis, which is described below, primarily based o the works by Freeze et al (1992) ad James & Freeze (1993). We describe the methodology, assumig a simple case ivolvig a decisio of whether to remediate or ot, that this decisio depeds o the cotamiat levels with respect to a defied actio level, ad that failure occurs if o actio is take whe cotamiatio occurs. 1 A cell is a homogeeous part of the cotamiated site

6 6 Let C deote the evet of cotamiatio above a specific actio level i a hypothetical cell, which is hydraulically isolated, ad let C c deote the complemetary evet that cotamiatio is ot preset. The probability Pr(C) is the the prior probability of cotamiatio for this cell. Typically, Pr(C) is ot kow with certaity. I early stages, a estimate may therefore be based o professioal judgmet usig existig site iformatio ad experieces from other sites. Give that remediatio has to be made if cotamiatio occurs, the prior risk-cost is defied as: R = k Pr( ) (1) 0 F C where k F deotes the moetary cosequeces of failure. Thus, R 0 is the expected cost of o actio. Let k R deote the remediatio cost. Aother alterative is to take a sample to see if ay cotamiat is preset ad to act accordigly. Let D ad D c deote the evets that the cotamiat is detected ad ot detected, respectively, give that the cell is sampled. Bayes theorem gives us the pre-posterior probability of failure to remediate, give that o cotamiat is detected: Pr ( C c D Pr ( C) Pr ( D c C) ) = (2) c c c c Pr ( C) Pr ( D C) + Pr ( C ) Pr ( D C ) ad the pre-posterior risk-cost is: c R = k Pr( C D ) (3) 1 F

7 7 Let k M deote the samplig cost. The expected risk-cost of takig a sample is: c k R Pr( D) + mi[ k R, R1 ]Pr( D ) (4) sice the cell is remediated if cotamiatio is detected, ad if cotamiatio is ot detected the cell is remediated oly if the remediatio cost is less tha the pre-posterior risk-cost. Followig Freeze et al (1992) we defie the data worth as: c W = mi[ k R, R0 ] { k R Pr( D) + mi[ k R, R1 ]Pr( D )} (5) Sice the quatities ivolved i the calculatio of the worth W typically are kow oly approximately, oe may require that W be sigificatly larger tha the samplig cost i order to justify samplig. This aalysis applies to all cells of the regio of iterest ad the most elemetary ad aive way to proceed would be to say that the cells that should be aalyzed further are those with W sufficietly larger tha the samplig cost. However, such a aalysis does ot take ito accout the obvious spatial depedece that exists i ay remediatio site. The effects of spatial depedece are further described i the secod ad third papers i this series. Below, we preset a aalysis of ecoomical ad statistical implicatios o the dataworth without spatial cosideratios.

8 8 The worth of takig additioal samples Now cosider a cell to be ivestigated with respect to groudwater cotamiatio. Prior to aalyzig the value of additioal samples, we have two alterative optios: (1) to remediate the cell to the cost of k R ; (2) to do othig with a remaiig risk-cost of R 0 =k F Pr(C). The decisio rule is to do othig if Pr(C) k R /k F ad remediate if Pr(C) > k R /k F. A third alterative is to take ew samples ad act accordig to the outcome of these (see also Figure 1). A objective fuctio, Φ(), is defied for each samplig alterative i which the umber of ew samples is =1,2. The objective fuctio of the optimal actio is a priori, before cosiderig ay ew samples: Φ ( 0) = mi( k, k Pr( C)) (6) R F We will ow cosider the possibilities of takig oe or several imperfect samples with the coditioal detectio probabilities Pr(D C) (which may be much less tha 1) ad Pr(D C c ) (which may be cosiderably larger tha 0). Let k M deote the cost of oe ew sample ad k M () deote the cost of simultaeously takig ew samples. A sample cost fuctio that we have used i our work is: 1 k M ( ) = k M r, = 1, 2,... (7) where r 1 is a discout factor. The objective fuctio of takig ew samples is:

9 9 Φ( ) = k = k M M ( ) + k = 0 mi( k ( ) + E[mi( k R R, k, k F F Pr( C D Pr( C D ))] = k)) Pr( D = k) (8) where D (=1, 2, ) is the total umber of cotamiat detectios i idepedet samples. The coditioal probabilities Pr(C D = k) ca be show, by meas of Bayes rule, to be: Pr( C D Pr( C) Pr( D = k C) = k) = (9a) Pr( D = k) where c c Pr( D = k) = Pr( C) Pr( D = k C) + Pr( C ) Pr( D = k C ) (9b) Pr(D = k C) is the biomial probability of gettig k successes i idepedet trials with success probability Pr(D C) (a aalogous statemet holds for Pr(D = k C c )). The optimal umber of samples is: opt = arg mi Φ( ) 0 (10) Note that Φ = Φ ) is the expected cost of the optimal actio, which is do othig if opt ( opt opt = 0 ad Pr(C) < k R /k F, take opt samples ad act accordigly if opt 1, ad remediate if opt = 0 ad Pr(C) k R /k F.

10 10 The Expected Net Value (ENV) of takig samples is: ENV ( ) = Φ(0) Φ( ), = 1, 2,... (11) This is related to the Expected Value of Sample Iformatio (EVSI) discussed i Bruce & Freeze (1993) ad Bruce & Gorelick (1994). It is ecessary to iclude the sample cost i ENV, sice it depeds o the umber of samples, while it is ot icluded i the EVSI, which therefore is positive for all values of Pr(C). Aother differece of our approach to evaluatig the worth of additioal samples is that Bruce & Freeze (1993) ad Bruce & Gorelick (1994) discuss the value of takig oe additioal sample, whereas we are ope to the possibility that it may be cost-effective to simultaeously take several additioal samples. Obviously, it is cost-effective to take samples oly if: ENV = max ENV( ) > 0 1 Examples igorig ucertaity i Pr(C) estimates Havig described the geeral theory of Bayesia data-worth aalysis, we ow preset two hypothetical cases o usig the described methodology for evaluatio of the data-worth of additioal samples. I this first presetatio of cases, the purpose is to illustrate the applicatio of the data-worth assessmet approach described i this paper. I the subsequet

11 11 sectio, we itroduce ucertaity assessmet of iput variables ad ivestigate the impact of this o decisios regardig the desig of the samplig program. Case 1 Cosider a hydraulically isolated cell of uiform ad isotropic coditios. The cost fuctios ad the coditioal detectio probabilities, which reflect the errors of the chose samplig techique, are show i Table 1. As see from the table, we assume relatively large samplig errors. We further assume that the samplig cost is small compared to the remediatio cost, ad that the failure cost is sigificatly larger tha the remediatio cost. Note that we used o discout factor i the sample cost fuctio, k M, i this example. Usig the procedure for calculatig the optimal umber of samples, opt, described above, the optimum umber of samples for differet Pr(C) was calculated ad is show i Figure 2. Figure 3 shows the prior expected cost (Φ(0)) ad the optimal cost (Φ opt ) for differet Pr(C). Figure 4 shows the ENV for differet values of Pr(C). Cosiderig the results show i these figures, we take Pr(C) = k R /k F = 0.4 prior to ay samples, sice this represets the maximum decisio ucertaity. As ca be see from Figure 2, the optimal umber of samples for Pr(C) = 0.4 is 3. The four possible outcomes of takig three samples, D 3, are show i Table 2. Give the iformatio about the cell ad the data-worth aalysis, we thus coclude that three additioal samples represet the most cost-effective sample alterative. The example also illustrates that eve if cotamiatio was to be detected i oe sample (D 3 = 1), the decisio would still be to do othig, due to the large samplig error associated with the sample techique that we have chose. I real world applicatios the samplig error is ofte igored.

12 12 Case 2 The coditios of the secod hypothetical case, also assumig a cell of uiform ad isotropic coditios, are show i Table 3. Also here, we let Pr(C) = k R /k F = 0.4 prior to ay samples. I this case we estimate the samplig error to be smaller tha i the first example ad we use a discout factor for the k M cost fuctio. Figure 5 shows the prior expected cost (Φ(0)) ad the optimal cost (Φ opt ) for differet values of Pr(C). Figure 6 shows the ENV for differet Pr(C). Calculatig the optimum umber of additioal samples yields opt = 1. The two possible outcomes of takig oe sample, D 1, are show i Table 4. If the sigle ew sample is made the decisio would be do othig if cotamiatio was ot detected ad remediate if cotamiatio was detected. This example illustrates the situatio of a expesive, but also reasoably accurate samplig techique, which makes remediatio cost-effective whe cotamiatio is detected i a sigle sample. Examples icludig ucertaity i P(C) estimates I the two hypothetical cases just studied, Pr(C) is take to equal k R /k F sice for this value, the expected et value (ENV) is maximal. However, it is typical that oe has some idea about the true Pr(C), but that the estimate is associated with some ucertaity. A evaluator may therefore wat to express the estimatio as Pr(C) 0.65 ad 0.5 Pr(C) 0.8. This iformatio may be tured ito a prior distributio of Pr(C). The most obvious prior is the triagular distributio with mode 0.65 ad support (0.5, 0.8). This prior distributio gives mass 1 to the iterval (0.5, 0.8). Thus, it does ot accout for the fact that the expert specificatio 0.5 Pr(C) 0.8 might be wrog.

13 13 A more realistic prior is therefore the beta distributio, havig desity fuctio: Γ( α + β ) Γ( α) Γ( β ) α 1 β 1 f ( θ ) = θ (1 θ ) (12) for 0 < θ < 1 (i order to avoid a cumbersome otatio, we ow write θ istead of Pr(C)). The parameters α ad β cocretize the expert opiio ad Γ deotes the Euler gamma fuctio. A beta prior allows for wrog expert estimates. For example, ote that Pr{θ < 0.5} = Pr{θ > 0.8} = 0.05 if α = 16.4, β =8.6. This beta distributio is plotted i Figure 7. This is oe strategy to cocretize ad quatify the reliability of the experts. Aother strategy would be to specify the ucertaity usig a estimate of the most likely value, θ, ad a umber of equivalet samples represetig the reliability of this estimate. The appropriate beta distributio the is the oe with α = θ + 1, β = (1 θ ) + 1. This follows from the wellkow fact that the beta distributio is cojugate to the biomial (Bickel & Doksum, 1977). Providig θ (=Pr(C)) with a beta prior forces us to redefie opt. Firstly, we specify a desig criterio such as the expected value or a suitable percetile i the distributio of Φ() for various umber of samples. Secodly, if the expected value is take, the opt = arg mi E[ Φ( )] 0 Similarly, if the 100p th percetile is chose, the

14 14 opt = arg mi G 0 1 ( p) where G 1 ( p) deotes the percetile fuctio of the distributio fuctio G ( x) = Pr{ Φ( ) x}. We ow retur to the two hypothetical cases. Case 1 The coditios for this case are idetical to the first case i the previous sectio (see Table 1), except that we ow itroduce ucertaity i the θ (=Pr(C)) estimate. Let us assume that the experts agree o the followig: θ 0.45 ad the ucertaity (or reliability) of this statemet is equivalet to the oe achieved after takig 8 idepedet samples. A appropriate parameterizatio of the prior beta distributio is the α = 4.6 ad β = 5.4. This beta distributio is plotted i Figure 8. From Figure 2 it is see that the optimal umber of samples is opt = 3 if θ = However, if the ucertaity is take ito accout, the optimal umber of samples is 2 if the desig criterio is the expected value, i.e. = arg mi E[ Φ( )]. I Figure 9 the desity of Φ(2) is plotted. Notice that it has a discrete compoet: Pr{Φ(2) < 11} = ad Pr{Φ(2) = 11} = Note also that if the ucertaity i θ istead correspods to a larger umber of idepedet samples, e.g. 50, the opt =3 if the desig criterio is the expected value. opt 0 Case 2 I the secod hypothetical case, opt = 1 for Pr(C) 0.627, 2 for Pr(C) ad 0 otherwise. This meas that if it is estimated that θ 0.4, the opt = 1 also if ay

15 15 reasoable ucertaity is take ito accout. The situatio chages, however, ad becomes similar to what we observed i Case 1, if θ Assume that the expert opiio is 0.50 θ Assume also that the experts are regarded as highly qualified ad reliable by the resposible decisio-maker. A quatificatio of the reliability of the expert opiios may therefore be expressed as: Pr{θ < 0.5} = Pr{θ > 0.75} = This holds for α = ad β = I this case mi E[Φ()] is obtaied for both = 1 ad = 2 if the desig criterio is expectatio. This holds true also if the media is chose. However, if the desig criterio is the 75 th percetile, opt = 2. The outcome would differ if the decisio-maker would be less willig to deped o the expert opiio. If the reliability of the expert opiio was istead expressed as: Pr{θ < 0.5} = Pr{θ > 0.75} = 0.3, this would correspod to α = 2.70 ad β = I this case, the optimal umber of samples would be 1 if the desig criterio is expectatio, 2 if the media is the desig criterio ad 1 if the criterio is the 75 th percetile. Coclusios The followig coclusios were made with respect to the o-spatial case of Bayesia dataworth aalysis: 1. Give the expected costs for samplig, remediatio ad failure, the probability of cotamiatio, ad the sample error, the optimal umber of samples to be take simultaeously i the ext ivestigatio phase ca be calculated.

16 16 2. Bayesia data-worth aalysis is directed at miimizig the sum of ivestmet ad risk cost. The samplig cost is compared to the reductio of the ecoomical risk of failig to meet set up maagemet goals of the site, implyig that the risk reductio should exceed the samplig cost i order to make samplig cost-effective. This is i cotrast to more commoly applied variace reductio approaches, which are directed at reachig a level of acceptable ucertaity, rather tha lookig at the cost-efficiecy of the samplig program. 3. Bayesia data-worth aalysis provides a possibility for formal icorporatio of the samplig error i the valuatio of the cost-efficiecy of the samplig program. As show i the 1 st hypothetical case, that eve if cotamiatio is detected remediatio may ot be ecoomically motivated if the samplig error is large. 4. Differet results ca be obtaied, depedig o whether ucertaity i the prior estimate of the probability of cotamiatio is icluded or ot. 5. Differet results ca be obtaied, depedig o the strategy for hadlig the ucertaity of the prior estimatio of the probability of cotamiatio. If the ucertaity of the prior estimate is take ito accout, the decisio is sesitive to the selected desig criteria, e.g. expected value, media or aother percetile. We preset a approach to model this type of ucertaity usig a beta-distributio. Also remediatio costs ad failure costs may be associated with substatial ucertaities, which may have large impact o the fial decisio. We do ot preset specific approaches for hadlig these ucertaities, but recommed that they are icluded i the same maer as the cotamiatio ucertaity. Fially, the Bayesia data-worth aalysis approach preseted i this paper icorporates the possibility of takig several ew samples simultaeously, ad therefore differs to some extet

17 17 to the previous works by e.g. James & Freeze (1993) ad James & Gorelick (1994). We also illustrate the importace for hadlig ucertaities i the estimates of the key variables of the data-worth aalysis. We do this without the complexities of spatial statistical models. This will be further explored i the 2 d ad 3 rd papers i this series. Ackowledgmets Refereces Back, P-E, Samplig strategies ad data-worth aalysis for cotamiated lad a literature review. Departmet of Geology, Chalmers Uiversity of Techology. Publ B 486 ad SGI Varia 500. Bickel, P.J. & K.A. Doksum, Mathematical Statistics. Basic Ideas ad Selected Topics. Holde-Day, Ic. Freeze, R. A., Massma, J., Smith, J. L., Sperlig, T. & B. R. James, Hydrogeological decisio aalysis: 1. A Framework. Groud Water. Vol. 28, o. 5, pp Freeze, R. A., Massma, J., Smith, J. L., Sperlig, T. & B. R. James, Hydrogeological decisio aalysis: 4. The cocept of data worth ad its use i the developmet of site ivestigatio strategies. Groud Water. Vol. 28, o. 5, pp James, B. R. & S. M. Gorelick, Whe eough is eough: The worth of moitorig data i aquifer remediatio desig. Water Resources Research. Vol. 30, o. 12, pp

18 18 James, B. R., Huff, D. D., Trabalka, J. R., Ketelle, R. H. & C. T. Rightmire, Allocatio of evirometal fuds usig ecoomic risk-cost-beefit aalysis: A case study. Groud Water Remedaitio ad Moitorig. Fall 1996, pp James, B.R. & Freeze, R.A., 1993: The worth of data i predictig aquitard cotiuity i hydrogeological desig. Water Resources Research, Vol 29, No. 7, pp Mari, C. M., Media, M. A. & J. B. Butcher, Mote Carlo aalysis ad Bayesia decisio theory for assessig the effects of waste sites o groudwater, I: Theory. Joural of Cotamiat Hydrology. No. 5, pp Massma, J. & R. A. Freeze, Groudwater cotamiatio from waste sites: The iteractio betwee risk-based egieerig desig ad regulatory policy, 1: Methodology. Water Resources Research. Vol. 23, o. 2, pp Massma, J., Freeze, R. A., Smith, L., Sperlig, T. & B. James, Hydrogeological Decisio Aalysis: 2.Applicatios to Groud-Water Cotamiatio. Groud Water. Vol. 29, o. 4, pp Norrma, J, Decisio Aalysis uder Risk ad Ucertaity at Cotamiated Sites A literature review. Departmet of Geology, Chalmers Uiversity of Techology. Publ B 486 ad SGI Varia 500. Rosé, L. ad Wladis, D., 1998: Hydrogeological Decisio Aalysis for Moetary Evaluatio of Nitrate Reductio Measures i Europe. Proceedigs of the Iteratioal Associatio of Hydrogeologists XXVII Cogress. Wladis, D., Rosé, L. ad Kros, J., 1999: Risk-Based Decisio Aalysis of Emissio Alteratives to Reduce Groudwater Acidificatio o the Europea Scale. Groud Water, Vol. 37, No. 6, pp

19 19 Prior aalysis Use kow iformatio to choose best alterative Data-worth aalysis Is ew iformatio likely to be cost-effective? Yes Take sample No Stop Posterior aalysis Use best ew iformatio to choose best alterative Figure 1. Schematic descriptio of Bayesia data-worth aalysis (after James et al, 1996).

20 20 opt P C Figure 2. Plot of opt vs Pr(C) for the 1 st hypothetical case. 0, opt P C Figure 3. Plot of Φ(0) (dashed) ad Φ opt vs Pr(C) for the 1 st hypothetical case.

21 21 0, opt P C Figure 4. Plot of Φ(0) (dashed) ad Φ opt vs Pr(C) for the 2 d hypothetical case. ENV P C Figure 5. Plot of ENV vs Pr(C) for the 1 st hypothetical case.

22 22 ENV P C -1-2 Figure 6. Plot of ENV vs Pr(C) for the 2 d hypothetical case Figure 7. The beta(16.4, 8.6) desity fuctio.

23 Figure 8. The beta(4.6, 5.4) desity fuctio Figure 9. The distributio of Φ(2) for the 1 st hypothetical case for a beta(4.6, 5.4) prior o θ = Pr(C).

24 24 Table 1. Costs ad probabilities for 1 st hypothetical case. k M () k R k F Pr(D C) Pr(D C c ) Table 2. The four possible outcomes of takig three additioal samples, D 3. k Pr(D 3 = k) Pr(C D 3 = k) k F Pr(C D 3 = k) Decisio do othig do othig remediate remediate Table 3. Costs ad probabilities for 2 d hypothetical case. k M () k R k F Pr(D C) Pr(D C c ) Table 4. The two possible outcomes of takig oe additioal sample, D 1. k Pr(D 1 = k) Pr(C D 1 = k) k F Pr(C D 1 = k) Decisio do othig remediate

A statistical method to determine sample size to estimate characteristic value of soil parameters

A statistical method to determine sample size to estimate characteristic value of soil parameters A statistical method to determie sample size to estimate characteristic value of soil parameters Y. Hojo, B. Setiawa 2 ad M. Suzuki 3 Abstract Sample size is a importat factor to be cosidered i determiig

More information

On an Application of Bayesian Estimation

On an Application of Bayesian Estimation O a Applicatio of ayesia Estimatio KIYOHARU TANAKA School of Sciece ad Egieerig, Kiki Uiversity, Kowakae, Higashi-Osaka, JAPAN Email: ktaaka@ifokidaiacjp EVGENIY GRECHNIKOV Departmet of Mathematics, auma

More information

Bayesian Methods: Introduction to Multi-parameter Models

Bayesian Methods: Introduction to Multi-parameter Models Bayesia Methods: Itroductio to Multi-parameter Models Parameter: θ = ( θ, θ) Give Likelihood p(y θ) ad prior p(θ ), the posterior p proportioal to p(y θ) x p(θ ) Margial posterior ( θ, θ y) is Iterested

More information

AAEC/ECON 5126 FINAL EXAM: SOLUTIONS

AAEC/ECON 5126 FINAL EXAM: SOLUTIONS AAEC/ECON 5126 FINAL EXAM: SOLUTIONS SPRING 2015 / INSTRUCTOR: KLAUS MOELTNER This exam is ope-book, ope-otes, but please work strictly o your ow. Please make sure your ame is o every sheet you re hadig

More information

GUIDELINES ON REPRESENTATIVE SAMPLING

GUIDELINES ON REPRESENTATIVE SAMPLING DRUGS WORKING GROUP VALIDATION OF THE GUIDELINES ON REPRESENTATIVE SAMPLING DOCUMENT TYPE : REF. CODE: ISSUE NO: ISSUE DATE: VALIDATION REPORT DWG-SGL-001 002 08 DECEMBER 2012 Ref code: DWG-SGL-001 Issue

More information

The target reliability and design working life

The target reliability and design working life Safety ad Security Egieerig IV 161 The target reliability ad desig workig life M. Holický Kloker Istitute, CTU i Prague, Czech Republic Abstract Desig workig life ad target reliability levels recommeded

More information

An Introduction to Randomized Algorithms

An Introduction to Randomized Algorithms A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis

More information

This is an introductory course in Analysis of Variance and Design of Experiments.

This is an introductory course in Analysis of Variance and Design of Experiments. 1 Notes for M 384E, Wedesday, Jauary 21, 2009 (Please ote: I will ot pass out hard-copy class otes i future classes. If there are writte class otes, they will be posted o the web by the ight before class

More information

Statistics 511 Additional Materials

Statistics 511 Additional Materials Cofidece Itervals o mu Statistics 511 Additioal Materials This topic officially moves us from probability to statistics. We begi to discuss makig ifereces about the populatio. Oe way to differetiate probability

More information

The Random Walk For Dummies

The Random Walk For Dummies The Radom Walk For Dummies Richard A Mote Abstract We look at the priciples goverig the oe-dimesioal discrete radom walk First we review five basic cocepts of probability theory The we cosider the Beroulli

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

Double Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution

Double Stage Shrinkage Estimator of Two Parameters. Generalized Exponential Distribution Iteratioal Mathematical Forum, Vol., 3, o. 3, 3-53 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/.9/imf.3.335 Double Stage Shrikage Estimator of Two Parameters Geeralized Expoetial Distributio Alaa M.

More information

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY

EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY EXAMINATIONS OF THE ROYAL STATISTICAL SOCIETY GRADUATE DIPLOMA, 016 MODULE : Statistical Iferece Time allowed: Three hours Cadidates should aswer FIVE questios. All questios carry equal marks. The umber

More information

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures

FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals

More information

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering

CEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio

More information

Lecture 2: Monte Carlo Simulation

Lecture 2: Monte Carlo Simulation STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?

More information

Estimation for Complete Data

Estimation for Complete Data Estimatio for Complete Data complete data: there is o loss of iformatio durig study. complete idividual complete data= grouped data A complete idividual data is the oe i which the complete iformatio of

More information

Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample.

Statistical Inference (Chapter 10) Statistical inference = learn about a population based on the information provided by a sample. Statistical Iferece (Chapter 10) Statistical iferece = lear about a populatio based o the iformatio provided by a sample. Populatio: The set of all values of a radom variable X of iterest. Characterized

More information

Element sampling: Part 2

Element sampling: Part 2 Chapter 4 Elemet samplig: Part 2 4.1 Itroductio We ow cosider uequal probability samplig desigs which is very popular i practice. I the uequal probability samplig, we ca improve the efficiecy of the resultig

More information

Output Analysis and Run-Length Control

Output Analysis and Run-Length Control IEOR E4703: Mote Carlo Simulatio Columbia Uiversity c 2017 by Marti Haugh Output Aalysis ad Ru-Legth Cotrol I these otes we describe how the Cetral Limit Theorem ca be used to costruct approximate (1 α%

More information

Information-based Feature Selection

Information-based Feature Selection Iformatio-based Feature Selectio Farza Faria, Abbas Kazeroui, Afshi Babveyh Email: {faria,abbask,afshib}@staford.edu 1 Itroductio Feature selectio is a topic of great iterest i applicatios dealig with

More information

6.3 Testing Series With Positive Terms

6.3 Testing Series With Positive Terms 6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial

More information

The standard deviation of the mean

The standard deviation of the mean Physics 6C Fall 20 The stadard deviatio of the mea These otes provide some clarificatio o the distictio betwee the stadard deviatio ad the stadard deviatio of the mea.. The sample mea ad variace Cosider

More information

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5 CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio

More information

1 Review and Overview

1 Review and Overview DRAFT a fial versio will be posted shortly CS229T/STATS231: Statistical Learig Theory Lecturer: Tegyu Ma Lecture #3 Scribe: Migda Qiao October 1, 2013 1 Review ad Overview I the first half of this course,

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

Chapter 6 Sampling Distributions

Chapter 6 Sampling Distributions Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to

More information

1 Introduction to reducing variance in Monte Carlo simulations

1 Introduction to reducing variance in Monte Carlo simulations Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by

More information

1 Inferential Methods for Correlation and Regression Analysis

1 Inferential Methods for Correlation and Regression Analysis 1 Iferetial Methods for Correlatio ad Regressio Aalysis I the chapter o Correlatio ad Regressio Aalysis tools for describig bivariate cotiuous data were itroduced. The sample Pearso Correlatio Coefficiet

More information

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1

EECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1 EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum

More information

Mixtures of Gaussians and the EM Algorithm

Mixtures of Gaussians and the EM Algorithm Mixtures of Gaussias ad the EM Algorithm CSE 6363 Machie Learig Vassilis Athitsos Computer Sciece ad Egieerig Departmet Uiversity of Texas at Arligto 1 Gaussias A popular way to estimate probability desity

More information

Application Bayesian Approach in Tasks Decision-Making Support

Application Bayesian Approach in Tasks Decision-Making Support America Joural of Applied Mathematics ad Statistics, 05, Vol 3, No 6, 57-6 Available olie at http://pubssciepubcom/ajams/3/6/7 Sciece ad Educatio Publishig DOI:069/ajams-3-6-7 Applicatio Bayesia Approach

More information

UNIT 2 DIFFERENT APPROACHES TO PROBABILITY THEORY

UNIT 2 DIFFERENT APPROACHES TO PROBABILITY THEORY UNIT 2 DIFFERENT APPROACHES TO PROBABILITY THEORY Structure 2.1 Itroductio Objectives 2.2 Relative Frequecy Approach ad Statistical Probability 2. Problems Based o Relative Frequecy 2.4 Subjective Approach

More information

Properties and Hypothesis Testing

Properties and Hypothesis Testing Chapter 3 Properties ad Hypothesis Testig 3.1 Types of data The regressio techiques developed i previous chapters ca be applied to three differet kids of data. 1. Cross-sectioal data. 2. Time series data.

More information

Lecture 19: Convergence

Lecture 19: Convergence Lecture 19: Covergece Asymptotic approach I statistical aalysis or iferece, a key to the success of fidig a good procedure is beig able to fid some momets ad/or distributios of various statistics. I may

More information

Chapter 8: STATISTICAL INTERVALS FOR A SINGLE SAMPLE. Part 3: Summary of CI for µ Confidence Interval for a Population Proportion p

Chapter 8: STATISTICAL INTERVALS FOR A SINGLE SAMPLE. Part 3: Summary of CI for µ Confidence Interval for a Population Proportion p Chapter 8: STATISTICAL INTERVALS FOR A SINGLE SAMPLE Part 3: Summary of CI for µ Cofidece Iterval for a Populatio Proportio p Sectio 8-4 Summary for creatig a 100(1-α)% CI for µ: Whe σ 2 is kow ad paret

More information

Bayesian and E- Bayesian Method of Estimation of Parameter of Rayleigh Distribution- A Bayesian Approach under Linex Loss Function

Bayesian and E- Bayesian Method of Estimation of Parameter of Rayleigh Distribution- A Bayesian Approach under Linex Loss Function Iteratioal Joural of Statistics ad Systems ISSN 973-2675 Volume 12, Number 4 (217), pp. 791-796 Research Idia Publicatios http://www.ripublicatio.com Bayesia ad E- Bayesia Method of Estimatio of Parameter

More information

IP Reference guide for integer programming formulations.

IP Reference guide for integer programming formulations. IP Referece guide for iteger programmig formulatios. by James B. Orli for 15.053 ad 15.058 This documet is iteded as a compact (or relatively compact) guide to the formulatio of iteger programs. For more

More information

Minimax Estimation of the Parameter of Maxwell Distribution Under Different Loss Functions

Minimax Estimation of the Parameter of Maxwell Distribution Under Different Loss Functions America Joural of heoretical ad Applied Statistics 6; 5(4): -7 http://www.sciecepublishiggroup.com/j/ajtas doi:.648/j.ajtas.654.6 ISSN: 6-8999 (Prit); ISSN: 6-96 (Olie) Miimax Estimatio of the Parameter

More information

CS284A: Representations and Algorithms in Molecular Biology

CS284A: Representations and Algorithms in Molecular Biology CS284A: Represetatios ad Algorithms i Molecular Biology Scribe Notes o Lectures 3 & 4: Motif Discovery via Eumeratio & Motif Represetatio Usig Positio Weight Matrix Joshua Gervi Based o presetatios by

More information

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution

Mathematical Modeling of Optimum 3 Step Stress Accelerated Life Testing for Generalized Pareto Distribution America Joural of Theoretical ad Applied Statistics 05; 4(: 6-69 Published olie May 8, 05 (http://www.sciecepublishiggroup.com/j/ajtas doi: 0.648/j.ajtas.05040. ISSN: 6-8999 (Prit; ISSN: 6-9006 (Olie Mathematical

More information

Lecture 2: April 3, 2013

Lecture 2: April 3, 2013 TTIC/CMSC 350 Mathematical Toolkit Sprig 203 Madhur Tulsiai Lecture 2: April 3, 203 Scribe: Shubhedu Trivedi Coi tosses cotiued We retur to the coi tossig example from the last lecture agai: Example. Give,

More information

Output Analysis (2, Chapters 10 &11 Law)

Output Analysis (2, Chapters 10 &11 Law) B. Maddah ENMG 6 Simulatio Output Aalysis (, Chapters 10 &11 Law) Comparig alterative system cofiguratio Sice the output of a simulatio is radom, the comparig differet systems via simulatio should be doe

More information

x a x a Lecture 2 Series (See Chapter 1 in Boas)

x a x a Lecture 2 Series (See Chapter 1 in Boas) Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio

More information

Markov Decision Processes

Markov Decision Processes Markov Decisio Processes Defiitios; Statioary policies; Value improvemet algorithm, Policy improvemet algorithm, ad liear programmig for discouted cost ad average cost criteria. Markov Decisio Processes

More information

Topic 9: Sampling Distributions of Estimators

Topic 9: Sampling Distributions of Estimators Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be

More information

Discrete Mathematics for CS Spring 2008 David Wagner Note 22

Discrete Mathematics for CS Spring 2008 David Wagner Note 22 CS 70 Discrete Mathematics for CS Sprig 2008 David Wager Note 22 I.I.D. Radom Variables Estimatig the bias of a coi Questio: We wat to estimate the proportio p of Democrats i the US populatio, by takig

More information

n n i=1 Often we also need to estimate the variance. Below are three estimators each of which is optimal in some sense: n 1 i=1 k=1 i=1 k=1 i=1 k=1

n n i=1 Often we also need to estimate the variance. Below are three estimators each of which is optimal in some sense: n 1 i=1 k=1 i=1 k=1 i=1 k=1 MATH88T Maria Camero Cotets Basic cocepts of statistics Estimators, estimates ad samplig distributios 2 Ordiary least squares estimate 3 3 Maximum lielihood estimator 3 4 Bayesia estimatio Refereces 9

More information

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis

Recursive Algorithms. Recurrences. Recursive Algorithms Analysis Recursive Algorithms Recurreces Computer Sciece & Egieerig 35: Discrete Mathematics Christopher M Bourke cbourke@cseuledu A recursive algorithm is oe i which objects are defied i terms of other objects

More information

Let us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f.

Let us give one more example of MLE. Example 3. The uniform distribution U[0, θ] on the interval [0, θ] has p.d.f. Lecture 5 Let us give oe more example of MLE. Example 3. The uiform distributio U[0, ] o the iterval [0, ] has p.d.f. { 1 f(x =, 0 x, 0, otherwise The likelihood fuctio ϕ( = f(x i = 1 I(X 1,..., X [0,

More information

L = n i, i=1. dp p n 1

L = n i, i=1. dp p n 1 Exchageable sequeces ad probabilities for probabilities 1996; modified 98 5 21 to add material o mutual iformatio; modified 98 7 21 to add Heath-Sudderth proof of de Fietti represetatio; modified 99 11

More information

GG313 GEOLOGICAL DATA ANALYSIS

GG313 GEOLOGICAL DATA ANALYSIS GG313 GEOLOGICAL DATA ANALYSIS 1 Testig Hypothesis GG313 GEOLOGICAL DATA ANALYSIS LECTURE NOTES PAUL WESSEL SECTION TESTING OF HYPOTHESES Much of statistics is cocered with testig hypothesis agaist data

More information

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE

Provläsningsexemplar / Preview TECHNICAL REPORT INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE TECHNICAL REPORT CISPR 16-4-3 2004 AMENDMENT 1 2006-10 INTERNATIONAL SPECIAL COMMITTEE ON RADIO INTERFERENCE Amedmet 1 Specificatio for radio disturbace ad immuity measurig apparatus ad methods Part 4-3:

More information

Random Variables, Sampling and Estimation

Random Variables, Sampling and Estimation Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig

More information

Intro to Learning Theory

Intro to Learning Theory Lecture 1, October 18, 2016 Itro to Learig Theory Ruth Urer 1 Machie Learig ad Learig Theory Comig soo 2 Formal Framework 21 Basic otios I our formal model for machie learig, the istaces to be classified

More information

Read through these prior to coming to the test and follow them when you take your test.

Read through these prior to coming to the test and follow them when you take your test. Math 143 Sprig 2012 Test 2 Iformatio 1 Test 2 will be give i class o Thursday April 5. Material Covered The test is cummulative, but will emphasize the recet material (Chapters 6 8, 10 11, ad Sectios 12.1

More information

Infinite Sequences and Series

Infinite Sequences and Series Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet

More information

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence

Chapter 3. Strong convergence. 3.1 Definition of almost sure convergence Chapter 3 Strog covergece As poited out i the Chapter 2, there are multiple ways to defie the otio of covergece of a sequece of radom variables. That chapter defied covergece i probability, covergece i

More information

Robust Default Correlation for Cost Risk Analysis

Robust Default Correlation for Cost Risk Analysis Robust Default Correlatio for Cost Risk Aalysis Christia Smart, Ph.D., CCEA Director, Cost Estimatig ad Aalysis Missile Defese Agecy Preseted at the 03 ICEAA Professioal Developmet ad Traiig Workshop Jue,

More information

If, for instance, we were required to test whether the population mean μ could be equal to a certain value μ

If, for instance, we were required to test whether the population mean μ could be equal to a certain value μ STATISTICAL INFERENCE INTRODUCTION Statistical iferece is that brach of Statistics i which oe typically makes a statemet about a populatio based upo the results of a sample. I oesample testig, we essetially

More information

6 Sample Size Calculations

6 Sample Size Calculations 6 Sample Size Calculatios Oe of the major resposibilities of a cliical trial statisticia is to aid the ivestigators i determiig the sample size required to coduct a study The most commo procedure for determiig

More information

Posted-Price, Sealed-Bid Auctions

Posted-Price, Sealed-Bid Auctions Posted-Price, Sealed-Bid Auctios Professors Greewald ad Oyakawa 207-02-08 We itroduce the posted-price, sealed-bid auctio. This auctio format itroduces the idea of approximatios. We describe how well this

More information

Basics of Probability Theory (for Theory of Computation courses)

Basics of Probability Theory (for Theory of Computation courses) Basics of Probability Theory (for Theory of Computatio courses) Oded Goldreich Departmet of Computer Sciece Weizma Istitute of Sciece Rehovot, Israel. oded.goldreich@weizma.ac.il November 24, 2008 Preface.

More information

Goodness-of-Fit Tests and Categorical Data Analysis (Devore Chapter Fourteen)

Goodness-of-Fit Tests and Categorical Data Analysis (Devore Chapter Fourteen) Goodess-of-Fit Tests ad Categorical Data Aalysis (Devore Chapter Fourtee) MATH-252-01: Probability ad Statistics II Sprig 2019 Cotets 1 Chi-Squared Tests with Kow Probabilities 1 1.1 Chi-Squared Testig................

More information

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015

ECE 8527: Introduction to Machine Learning and Pattern Recognition Midterm # 1. Vaishali Amin Fall, 2015 ECE 8527: Itroductio to Machie Learig ad Patter Recogitio Midterm # 1 Vaishali Ami Fall, 2015 tue39624@temple.edu Problem No. 1: Cosider a two-class discrete distributio problem: ω 1 :{[0,0], [2,0], [2,2],

More information

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4

MATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4 MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.

More information

1.010 Uncertainty in Engineering Fall 2008

1.010 Uncertainty in Engineering Fall 2008 MIT OpeCourseWare http://ocw.mit.edu.00 Ucertaity i Egieerig Fall 2008 For iformatio about citig these materials or our Terms of Use, visit: http://ocw.mit.edu.terms. .00 - Brief Notes # 9 Poit ad Iterval

More information

Assessment of extreme discharges of the Vltava River in Prague

Assessment of extreme discharges of the Vltava River in Prague Flood Recovery, Iovatio ad Respose I 05 Assessmet of extreme discharges of the Vltava River i Prague M. Holický, K. Jug & M. Sýkora Kloker Istitute, Czech Techical Uiversity i Prague, Czech Republic Abstract

More information

PROBABILITY LOGIC: Part 2

PROBABILITY LOGIC: Part 2 James L Bec 2 July 2005 PROBABILITY LOGIC: Part 2 Axioms for Probability Logic Based o geeral cosideratios, we derived axioms for: Pb ( a ) = measure of the plausibility of propositio b coditioal o the

More information

MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND.

MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND. XI-1 (1074) MOST PEOPLE WOULD RATHER LIVE WITH A PROBLEM THEY CAN'T SOLVE, THAN ACCEPT A SOLUTION THEY CAN'T UNDERSTAND. R. E. D. WOOLSEY AND H. S. SWANSON XI-2 (1075) STATISTICAL DECISION MAKING Advaced

More information

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series

Comparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series

More information

Frequentist Inference

Frequentist Inference Frequetist Iferece The topics of the ext three sectios are useful applicatios of the Cetral Limit Theorem. Without kowig aythig about the uderlyig distributio of a sequece of radom variables {X i }, for

More information

4.3 Growth Rates of Solutions to Recurrences

4.3 Growth Rates of Solutions to Recurrences 4.3. GROWTH RATES OF SOLUTIONS TO RECURRENCES 81 4.3 Growth Rates of Solutios to Recurreces 4.3.1 Divide ad Coquer Algorithms Oe of the most basic ad powerful algorithmic techiques is divide ad coquer.

More information

Agreement of CI and HT. Lecture 13 - Tests of Proportions. Example - Waiting Times

Agreement of CI and HT. Lecture 13 - Tests of Proportions. Example - Waiting Times Sigificace level vs. cofidece level Agreemet of CI ad HT Lecture 13 - Tests of Proportios Sta102 / BME102 Coli Rudel October 15, 2014 Cofidece itervals ad hypothesis tests (almost) always agree, as log

More information

Statistical Pattern Recognition

Statistical Pattern Recognition Statistical Patter Recogitio Classificatio: No-Parametric Modelig Hamid R. Rabiee Jafar Muhammadi Sprig 2014 http://ce.sharif.edu/courses/92-93/2/ce725-2/ Ageda Parametric Modelig No-Parametric Modelig

More information

Measurement uncertainty of the sound absorption

Measurement uncertainty of the sound absorption Measuremet ucertaity of the soud absorptio coefficiet Aa Izewska Buildig Research Istitute, Filtrowa Str., 00-6 Warsaw, Polad a.izewska@itb.pl 6887 The stadard ISO/IEC 705:005 o the competece of testig

More information

Exponential Families and Bayesian Inference

Exponential Families and Bayesian Inference Computer Visio Expoetial Families ad Bayesia Iferece Lecture Expoetial Families A expoetial family of distributios is a d-parameter family f(x; havig the followig form: f(x; = h(xe g(t T (x B(, (. where

More information

Sample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D.

Sample Size Estimation in the Proportional Hazards Model for K-sample or Regression Settings Scott S. Emerson, M.D., Ph.D. ample ie Estimatio i the Proportioal Haards Model for K-sample or Regressio ettigs cott. Emerso, M.D., Ph.D. ample ie Formula for a Normally Distributed tatistic uppose a statistic is kow to be ormally

More information

CSE 527, Additional notes on MLE & EM

CSE 527, Additional notes on MLE & EM CSE 57 Lecture Notes: MLE & EM CSE 57, Additioal otes o MLE & EM Based o earlier otes by C. Grat & M. Narasimha Itroductio Last lecture we bega a examiatio of model based clusterig. This lecture will be

More information

There is no straightforward approach for choosing the warmup period l.

There is no straightforward approach for choosing the warmup period l. B. Maddah INDE 504 Discrete-Evet Simulatio Output Aalysis () Statistical Aalysis for Steady-State Parameters I a otermiatig simulatio, the iterest is i estimatig the log ru steady state measures of performace.

More information

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 +

62. Power series Definition 16. (Power series) Given a sequence {c n }, the series. c n x n = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + 62. Power series Defiitio 16. (Power series) Give a sequece {c }, the series c x = c 0 + c 1 x + c 2 x 2 + c 3 x 3 + is called a power series i the variable x. The umbers c are called the coefficiets of

More information

PH 425 Quantum Measurement and Spin Winter SPINS Lab 1

PH 425 Quantum Measurement and Spin Winter SPINS Lab 1 PH 425 Quatum Measuremet ad Spi Witer 23 SPIS Lab Measure the spi projectio S z alog the z-axis This is the experimet that is ready to go whe you start the program, as show below Each atom is measured

More information

Hashing and Amortization

Hashing and Amortization Lecture Hashig ad Amortizatio Supplemetal readig i CLRS: Chapter ; Chapter 7 itro; Sectio 7.. Arrays ad Hashig Arrays are very useful. The items i a array are statically addressed, so that isertig, deletig,

More information

Math 2784 (or 2794W) University of Connecticut

Math 2784 (or 2794W) University of Connecticut ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really

More information

ECE 901 Lecture 12: Complexity Regularization and the Squared Loss

ECE 901 Lecture 12: Complexity Regularization and the Squared Loss ECE 90 Lecture : Complexity Regularizatio ad the Squared Loss R. Nowak 5/7/009 I the previous lectures we made use of the Cheroff/Hoeffdig bouds for our aalysis of classifier errors. Hoeffdig s iequality

More information

Discrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Note 19

Discrete Mathematics and Probability Theory Spring 2016 Rao and Walrand Note 19 CS 70 Discrete Mathematics ad Probability Theory Sprig 2016 Rao ad Walrad Note 19 Some Importat Distributios Recall our basic probabilistic experimet of tossig a biased coi times. This is a very simple

More information

STAT Homework 1 - Solutions

STAT Homework 1 - Solutions STAT-36700 Homework 1 - Solutios Fall 018 September 11, 018 This cotais solutios for Homework 1. Please ote that we have icluded several additioal commets ad approaches to the problems to give you better

More information

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract

Goodness-Of-Fit For The Generalized Exponential Distribution. Abstract Goodess-Of-Fit For The Geeralized Expoetial Distributio By Amal S. Hassa stitute of Statistical Studies & Research Cairo Uiversity Abstract Recetly a ew distributio called geeralized expoetial or expoetiated

More information

DS 100: Principles and Techniques of Data Science Date: April 13, Discussion #10

DS 100: Principles and Techniques of Data Science Date: April 13, Discussion #10 DS 00: Priciples ad Techiques of Data Sciece Date: April 3, 208 Name: Hypothesis Testig Discussio #0. Defie these terms below as they relate to hypothesis testig. a) Data Geeratio Model: Solutio: A set

More information

Number of fatalities X Sunday 4 Monday 6 Tuesday 2 Wednesday 0 Thursday 3 Friday 5 Saturday 8 Total 28. Day

Number of fatalities X Sunday 4 Monday 6 Tuesday 2 Wednesday 0 Thursday 3 Friday 5 Saturday 8 Total 28. Day LECTURE # 8 Mea Deviatio, Stadard Deviatio ad Variace & Coefficiet of variatio Mea Deviatio Stadard Deviatio ad Variace Coefficiet of variatio First, we will discuss it for the case of raw data, ad the

More information

ANALYSIS OF EXPERIMENTAL ERRORS

ANALYSIS OF EXPERIMENTAL ERRORS ANALYSIS OF EXPERIMENTAL ERRORS All physical measuremets ecoutered i the verificatio of physics theories ad cocepts are subject to ucertaities that deped o the measurig istrumets used ad the coditios uder

More information

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall

Similarity Solutions to Unsteady Pseudoplastic. Flow Near a Moving Wall Iteratioal Mathematical Forum, Vol. 9, 04, o. 3, 465-475 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/0.988/imf.04.48 Similarity Solutios to Usteady Pseudoplastic Flow Near a Movig Wall W. Robi Egieerig

More information

Interval Intuitionistic Trapezoidal Fuzzy Prioritized Aggregating Operators and their Application to Multiple Attribute Decision Making

Interval Intuitionistic Trapezoidal Fuzzy Prioritized Aggregating Operators and their Application to Multiple Attribute Decision Making Iterval Ituitioistic Trapezoidal Fuzzy Prioritized Aggregatig Operators ad their Applicatio to Multiple Attribute Decisio Makig Xia-Pig Jiag Chogqig Uiversity of Arts ad Scieces Chia cqmaagemet@163.com

More information

Principle Of Superposition

Principle Of Superposition ecture 5: PREIMINRY CONCEP O RUCUR NYI Priciple Of uperpositio Mathematically, the priciple of superpositio is stated as ( a ) G( a ) G( ) G a a or for a liear structural system, the respose at a give

More information

Teaching Mathematics Concepts via Computer Algebra Systems

Teaching Mathematics Concepts via Computer Algebra Systems Iteratioal Joural of Mathematics ad Statistics Ivetio (IJMSI) E-ISSN: 4767 P-ISSN: - 4759 Volume 4 Issue 7 September. 6 PP-- Teachig Mathematics Cocepts via Computer Algebra Systems Osama Ajami Rashaw,

More information

The Sample Variance Formula: A Detailed Study of an Old Controversy

The Sample Variance Formula: A Detailed Study of an Old Controversy The Sample Variace Formula: A Detailed Study of a Old Cotroversy Ky M. Vu PhD. AuLac Techologies Ic. c 00 Email: kymvu@aulactechologies.com Abstract The two biased ad ubiased formulae for the sample variace

More information

MBACATÓLICA. Quantitative Methods. Faculdade de Ciências Económicas e Empresariais UNIVERSIDADE CATÓLICA PORTUGUESA 9. SAMPLING DISTRIBUTIONS

MBACATÓLICA. Quantitative Methods. Faculdade de Ciências Económicas e Empresariais UNIVERSIDADE CATÓLICA PORTUGUESA 9. SAMPLING DISTRIBUTIONS MBACATÓLICA Quatitative Methods Miguel Gouveia Mauel Leite Moteiro Faculdade de Ciêcias Ecoómicas e Empresariais UNIVERSIDADE CATÓLICA PORTUGUESA 9. SAMPLING DISTRIBUTIONS MBACatólica 006/07 Métodos Quatitativos

More information

Table 12.1: Contingency table. Feature b. 1 N 11 N 12 N 1b 2 N 21 N 22 N 2b. ... a N a1 N a2 N ab

Table 12.1: Contingency table. Feature b. 1 N 11 N 12 N 1b 2 N 21 N 22 N 2b. ... a N a1 N a2 N ab Sectio 12 Tests of idepedece ad homogeeity I this lecture we will cosider a situatio whe our observatios are classified by two differet features ad we would like to test if these features are idepedet

More information

6. Sufficient, Complete, and Ancillary Statistics

6. Sufficient, Complete, and Ancillary Statistics Sufficiet, Complete ad Acillary Statistics http://www.math.uah.edu/stat/poit/sufficiet.xhtml 1 of 7 7/16/2009 6:13 AM Virtual Laboratories > 7. Poit Estimatio > 1 2 3 4 5 6 6. Sufficiet, Complete, ad Acillary

More information

Confidence Intervals for the Population Proportion p

Confidence Intervals for the Population Proportion p Cofidece Itervals for the Populatio Proportio p The cocept of cofidece itervals for the populatio proportio p is the same as the oe for, the samplig distributio of the mea, x. The structure is idetical:

More information