Three new ways to calculate average (U Th)/He ages

Size: px
Start display at page:

Download "Three new ways to calculate average (U Th)/He ages"

Transcription

1 Avalable onlne at Chemcal Geology 249 (2008) Three new ways to calculate average (U Th)/He ages eter ermeesch School of Earth Scences; Brkbeck, Unversty of London; London WC1E 7HX, Unted Kngdom Receved 15 June 2007; receved n revsed form 8 December 2007; accepted 20 January 2008 Edtor: S.L. Goldsten Abstract Tradtonally the average age of multple (U Th)/He analyses has been calculated as the arthmetc mean age. Ths paper presents three alternatve methods: (a) n analogy wth the fsson track method, the pooled age s calculated by addng the respectve U, Th and He abundances of several grans together, thereby generatng one synthetc mult-gran measurement; (b) the sochron age s the slope of helum concentraton versus present-day helum producton; (c) the central age s computed from the geometrc mean U Th He composton. Each of these methods s more approprate than the arthmetc mean age n certan applcatons. The pooled age s useful for comparng sngle-gran wth mult-gran analyses, whle the sochron age can be used to detect parentless helum. The central age s the most accurate and statstcally robust way to calculate a sample average of several sngle-gran analyses because U, Th and He form a ternary system and only the central age adequately captures the statstcs of ths compostonal data space. Fortunately, the expected dfference between the arthmetc mean age and the central age s relatvely small, less than 1% f the external age reproducblty s better than 15% (1σ). Fnally, the (U Th)/He age equaton s vsualzed on a ternary dagram to llustrate that the α-ejecton correcton should be appled before, and not after age calculaton, n order to avod a partal lnearzaton of the age equaton. Includng Sm as a fourth parent element precludes a straghtforward vsualzaton of the age equaton on a twodmensonal plot. Nevertheless, the pooled, sochron and central age methods can be easly generalzed to the case of (U Th Sm)/He datng. To facltate the calculaton of the central age, a web-based calculator s provded at Elsever B.. All rghts reserved. Keywords: (U TH)/He; Thermochronology; Compostonal data analyss 1. Introducton Radogenc helum-geochronology s based on a summed set of dfferental equatons: d½heš dt Xn d½ Š wth : d ½ Š dt dt k ½ Š ð1þ where t=tme, [He]=helum abundance, [ ]=abundance of the th parent nuclde and λ =decay constant of ths nuclde (for 1 n). Despte the smplcty of Eq. (1), there are several ways to solve t, three of whch wll be dscussed n Secton 2. A lnear approxmaton s accurate to better than 1% for ages up to 100 Ma, whch can be consdered satsfactory n comparson wth the external reproducblty of (U Th)/He datng (20 E-mal address: 30%; e.g., Stock et al., 2006). Nevertheless, most researchers rghtly decde to calculate an exact age by numercal teraton. Ths paper rases the pont that the accuracy ganed by dong so s easly lost by two common practces: (1) performng the α- ejecton after, rather than before age calculaton and (2) usng the arthmetc mean age to summarze a dataset of several sngle-gran measurements. After Secton 3 presents two smlarly based alternatves to the arthmetc mean age that are approprate for specfc applcatons, Secton 4 ntroduces the central age as the most accurate way to compute average (U Th)/He ages. The accuracy ganed by usng the central age nstead of the arthmetc mean age s comparable to that ganed by teratvely solvng the (U Th)/He age equaton nstead of usng the lnear approxmaton. The only cost of the new procedure s computatonal complexty. To facltate the calculatons, they are mplemented n an onlne calculator (http://pvermees.andropov.org/central) and llustrated on a publshed dataset of ncluson-bearng apattes. Fnally, Secton /$ - see front matter 2008 Elsever B.. All rghts reserved. do: /j.chemgeo

2 340. ermeesch / Chemcal Geology 249 (2008) presents a generalzed method to calculate central ages for datasets that also nclude a fourth radoactve parent, 147 Sm. 2. Calculatng sngle-gran ages: many ways to skn a cat Ten naturally-occurng long-lved α-emttng radonucldes exst on Earth: 144 Nd, 147 Sm, 148 Sm, 152 Gd, 174 Hf, 186 Os, 190 t, 232 Th, 235 U and 238 U. For the purpose of helum-thermochronology, all but the heavest three of these nucldes can often be neglected because of ther low abundance and low helumyeld. For example, only one α-partcle s produced per 147 Sm, whereas sx to eght are formed n the Th and U decay seres. Further smplfcaton s possble because the present-day 238 U/ 235 U-rato s constant n the solar system (=137.88; Steger and Jäger, 1977). Therefore, the ngrowth of helum wth tme (t) can be wrtten as a functon of the elemental U, Th and He abundances or concentratons: ½HeŠ 8 137:88 138:88 ek 238t 7 1 þ 138:88 ek 235t 1 ½UŠþ6 e k232t 1 ½ThŠ ð2þ the average has been estmated by the arthmetc mean of the sngle-gran ages. Ths secton wll ntroduce two alternatve methods for calculatng average ages, and the next secton wll add a thrd. Each of these new methods s more approprate than the arthmetc mean age n specfc applcatons The pooled age Helum can be extracted from the host gran ether n a resstance furnace or by laser-heatng n a mcro-oven (House et al., 2000). In the former, but sometmes also n the latter case, t may be necessary to analyze multple mneral grans together (e.g., ersano et al., 2007). oolng several grans boosts the sgnal strength and sometmes averages out α-ejecton correcton errors caused by zonng and mneral nclusons. ermeesch et al. (2007) ntroduced the pooled age as the best way to compare multple sngle-gran ages wth one or more multgran ages, or to compare two sets of mult-gran ages wth each other (Fg. 1). The pooled age s calculated by addng the respectve U, Th and He abundances (n moles) of several measurements together, thereby generatng one synthetc mult- wth λ 232, λ 235 and λ 238 the decay constants of 232 Th, 235 U and 238 U, respectvely. Eq. (2) has no analytcal soluton but s easy to solve teratvely. However, for young ages (t 1/λ 235 ), a reasonably accurate lnear approxmaton also exsts: t ½HeŠ ð3þ wth the present-day helum producton rate: 8 137:88 138:88 k 238 þ 7 138:88 k 235 ½UŠþ6k 232 ½ThŠ ð4þ The accuracy of ths soluton wll be dscussed n Secton 4. Besdes beng easy to mplement, the lnear age equaton s useful for llustratve purposes and opens up some new applcatons whch wll be dscussed n Secton 3. Meesters and Duna (2005) ntroduced an alternatve drect soluton to the (U Th)/He age equaton: t 1 ln 1 þ k wm k wm ½HeŠ ð5þ wth λ wm the weghted mean decay constant: k wm 8 137:88 138:88 k2 238 þ 7 138:88 k2 235 ½UŠþ6k ½ThŠ ð6þ As shown by Meesters and Duna (2005) and n Secton 4, ths soluton s remarkably accurate for all practcal applcatons. 3. Mult-gran ages Eqs. (2 6) can be used to calculate (U Th)/He ages from ndvdual U, Th and He measurements, but do not explan how to calculate the average value of multple analyses. Tradtonally, Fg. 1. Lnearzed (U Th)/He dagram wth a subset of the HF-treated nclusonbearng apatte data of ermeesch et al. (2007). 30% uncertanty (2σ) was added to the helum abundances to account for α-ejecton correcton nduced scatter (ermeesch et al., 2007). (a) Accordng to the lnear age equaton, (U Th)/He ages are gven by the slopes of lnes connectng each (,[He])-pont wth the orgn; (b) the pooled age s a synthetc mult-gran age calculated from the summed producton rates and helum abundances of all the measurements. The box n 1.b marks the outlne of 1.a. The pooled age of the sample s 11.28±0.14 Ma.

3 . ermeesch / Chemcal Geology 249 (2008) gran measurement. The age of the pooled measurement can then be calculated usng any of the equatons gven n Secton 2. An obvous dsadvantage of poolng compostonal data s that the resultng age s based to the hgh U, Th or He compostons. Because such bas can be assocated wth anomalous grans affected by radaton damage or mplanted helum, the pooled age may be wrong by effectvely gvng extra weght to outlers. However, these objectons are also true for standard mult-gran analyses, whch cannot be avoded when datng small, young, or U Th-poor grans (e.g., ersano et al., 2007). In short, the pooled age must be used for and only for averagng mult-gran alquots The (U Th)/He sochron The prevous secton showed that (U Th)/He data can be vsualzed on a two-dmensonal plot of helum abundance or concentraton versus-producton (Fg. 1). To calculate a pooled age, t s mportant that [U], [Th] and [He] are elemental abundances, expressed n moles. If the data are recast n unts of concentraton, some of the bas towards hgh U Th-grans dsappears and the He dagram can be used to defne a (U Th)/He sochron. Ths s an unconstraned lnear ft through a seres of sngle-gran (,[He]) measurements. For an applcaton of the sochron method, consder the U Th rch mneral nclusons n apatte whch are often held responsble for erroneously old (U Th)/He ages, because they produce parentless He. Ths problem can be detected wth the (U Th)/ He sochron. In the absence of mneral nclusons, the sochron goes through the orgn (=[He]=0). However, n the presence of α-emttng nclusons, the sochron s ether not defned or does not go through the orgn. For example, consder the worst-case scenaro of an α-emttng zrcon ncluson contaned n an apatte wthout U and Th. The ncluson ejects He nto the surroundng apatte that s measured followng degassng by heatng wth a laser or n a resstance furnace. However, the zrcon ncluson wll not dssolve n the concentrated HNO 3 that s commonly used to dgest apattes pror to U Th analyss. Therefore, the apparent (U Th)-producton of such a sample s zero, and ts sochron does not go through the orgn of the [He] dagram. ermeesch et al. (2007) solved the parentless helum problem by dssoluton of the apatte and ts nclusons n hot HF. The effectveness of ths technque s llustrated by comparng an ncluson-rch sample from Naxos usng the tradtonal HNO 3 method wth an HF-treated alquot of the same sample. The latter defnes a well-constraned (U Th)/He sochron wth zero ntercept, whereas the former does not (Fg. 2). Calculaton of the sochron age, ncludng error propagaton, can easly be done usng the Isoplot Excel add-n (Ludwg, 2003). form a ternary system, can be plotted on a ternary dagram, and are subject to the pecular mathematcs of the ternary dataspace. In a three-component system (A+B+C=1), ncreasng one component (e.g., A) causes a decrease n the two other components (B and C). Another consequence of so-called data closure s that the arthmetc mean of compostonal data has no physcal meanng (Weltje, 2002) lottng the (U Th)/He age equaton on ternary dagrams Followng the nomenclature of Atchson (1986), the ternary dagram s a 2-smplex (Δ 2 ). The very fact that t s possble to plot ternary data on a two-dmensonal sheet of paper tells us that the sample space really has only two, and not three dmensons. As a soluton to the compostonal data problem, Atchson (1986) suggested to transform the data from Δ 2 to R 2 usng the lograto transformaton. After performng the desred ( tradtonal ) statstcal analyss on the transformed data n R 2, the results can be transformed back to Δ 2 usng the nverse lograto transformaton (Fg. 3). Implementaton detals about the lograto transformaton wll be gven n Secton 4.3. Ternary dagrams and lograto plots are useful tools for vsualzng U Th He data and the (U Th)/He age equaton. 4. U Th He as a ternary system The age-eqs. (2), (3) and (5) do not specfy the measurement unts of [U], [Th] and [He]. These can be expressed n moles or moles/g, but they can also be non-dmensonalzed by normalzaton to a constant sum: [U ] [U]/([U]+[Th]+[He]), [Th ] [Th]/([U]+[Th]+[He]) and [He ] [He]/([U]+[Th] +[He]) so that [U ]+[Th ]+[He ]=1. Therefore, U, Th and He Fg. 2. (U Th)/He sochron plots for ncluson-bearng apattes from Naxos, Greece (ermeesch et al., 2007). (a) HNO 3 -treated apattes do not plot on a lne, ndcatng parentless helum caused by undssolved mneral nclusons contanng mssng U and Th; (b) HF-treated apattes from the same sample do form a well-defned sochron ntersectng the orgn, ndcatng that all parent and daughter nucldes are accounted for. The sochron age s 12.0±4.2 Ma.

4 342. ermeesch / Chemcal Geology 249 (2008) Fg. 3. Ternary U Th He data can be mapped to two-dmensonal lograto-space wthout loss of nformaton. Thus, t can be shown that the lnear age equaton s accurate to better than 1% for ages up to 100 Ma (Fg. 4.a) whereas the equaton of Meesters and Duna (2005) reaches the same accuracy at 1 Ga (Fg. 4.b). Fg. 4.c represents a warnng aganst applyng the α-ejecton correcton after, rather than before the age calculaton. Ths causes a partal lnearzaton of the age equaton and results n a loss of accuracy. For example, dvdng an uncorrected (U Th)/He age by an α-retenton factor F t of 0.7 results n a msft that s 30% of the lnear age equaton msft. To take full advantage of the accuracy of the exact age equaton, one must dvde [He] by F t before calculatng the (U Th)/He age The central age The lograto transformaton s useful for more than just the purpose of vsualzaton. It provdes a fourth and arguably best way to calculate the average age of a populaton of sngle-gran (U Th)/He measurements. The central age s calculated from the average U Th He composton of the dataset, where average s defned as the geometrc mean of the sngle-gran U, Th and He measurements. The geometrc mean of compostonal data equals the arthmetc mean of lograto transformed data. How mportant s the dfference between the arthmetc mean age and the central age? To smplfy ths queston, consder the specal case of a sample wth only one radoactve parent, say Th. Assume that W = ln([he]/[th]) s normally dstrbuted wth mean µ and standard devaton σ. Usng the lnearzed age equaton for clarty, the central age t c s gven by: t c Ce A ð7þ wth C=1/(6λ 232 ) for Th. Usng the frst raw moment of the lognormal dstrbuton (Atchson and Brown, 1957), the arthmetc mean age t m s: t m Ce Aþr2 =2 so that the relatve dfference between t m and t c s: t m t c t c e r2 =2 1 ð9þ ð8þ Usng the second central moment of the lognormal dstrbuton (Atchson and Brown, 1957), the varance of the sngle-gran ages s gven by: r 2 t C 2 e r2 1 e 2Aþr2 ð10þ lottng (t m t c )/t c versus σ t /t c reveals that the central age s systematcally younger than the mean age. Fortunately, the dfference s small. For example, for a typcal external reproducblty of ~25% (e.g. σ t /t =11% for Stock et al., 2006), the expected dfference s b1% (Fg. 5). Fnally, t s nterestng to note that the geometrc mean of the lognormal dstrbuton equals ts medan. Therefore, the central age asymptotcally converges to the medan age. However, typcal numbers of replcate analyses are not suffcent for ths approach to be truly benefcal Applcaton to HF-treated Naxos apattes We now return to the sample of HF-treated ncluson-bearng apattes from Naxos that was prevously used to llustrate the pooled and sochron age (Fgs. 1 and 2). The raw data and the dfferent steps of the central age calculaton are gven n Table 1. We wll now walk through the dfferent parts (labeled a, b and c) of ths table. (a) The upper left part of Table 1 lsts the U, Th and He abundances of 11 sngle-gran analyses. Ther respectve sngle-gran ages (t) were calculated usng the exact age equaton, even though the lnear age approxmaton (Eq. (3)) s accurate to better than 0.1% for such young ages (Fg. 4. a). The pooled U, Th and He abundances are obtaned by smple summaton of the consttuent grans. Note that the helum abundances are corrected for α-ejecton pror to beng pooled. A nomnal σ = 15% statstcal uncertanty s assocated wth F t, assumng randomly dstrbuted mneral nclusons (ermeesch et al., 2007). The pooled abundances were normalzed to unty to facltate comparson wth the geometrc mean composton (see below).

5 . ermeesch / Chemcal Geology 249 (2008) (b) To calculate the sochron age, the abundances are frst rescaled to unts of concentraton (e.g. n nmol/g). Ths removes the bas towards large grans, whch can domnate the pooled age calculaton. The α-producton rate s gven by Eq. (4). The lnear regresson (Fg. 2) was done Fg. 5. Expected dfference between the mean age t c and the central age t m plotted aganst the relatve spread of the sngle-gran ages. usng Isochron 3.0 (Ludwg, 2003), yeldng a slope of 12.0 ±4.2 Ma wth an ntercept of 0.05±0.45 nmol/g He. (c) Central ages are somewhat more complcated to calculate than arthmetc mean ages, pooled ages or sochron ages. Therefore, these calculatons wll be dscussed n more detal. Frst, transform each of the n sngle-gran analyses to lograto-space (Fg. 6): ½U Š ½Th Š ln ; W ln ð11þ ½He Š ½He Š Fg. 4. (a) Relatve msft of the lnear approxmaton (Eq. (3)) to the exact age equaton (Eq. (2)); (b) msft of the drect soluton of Meesters and Duna (2005) (Eq. (5)); (c) msft caused by makng the α-ejecton correcton after the age calculaton (wth F t =0.7). For = 1,...,n. Note that ths transformaton can be done rrespectve of whether the U, Th and He measurements are expressed n abundance unts or n unts of concentraton. Followng standard error propagaton, the (co)varances of these quanttes are estmated by: ru U 2 ; He W rth Th 2 ; cov ;W rhe 2 He He ð12þ Next, calculate the arthmetc mean of the lograto transformed data: 1 ; W 1 W ; n n Wth the followng (co)varances: 1 n 2 ; W 1 n 2 W ; cov ;W 1 n 2 cov ;W ð13þ ð14þ Note that Eq. (14) only propagates the nternal (.e. analytcal) uncertanty, and not the external error. Sngle-gran (U Th)/He ages tend to suffer from overdsperson wth respect to the formal analytcal precson for a number of reasons (Ftzgerald et al., 2006; ermeesch et al., 2007). Therefore, t may be better to use an alternatve equaton propagatng the external error: n 2 n ; W W W 2 ; ð15þ nn ð 1Þ nn ð 1Þ n cov ;W W W nn ð 1Þ

6 Table 1 Step-by-step data reducton of the Naxos dataset of 11 ncluson-bearng apattes (wth an arthmetc mean age of 11.58±0.20 Ma) 344. ermeesch / Chemcal Geology 249 (2008) Ths table should be read from a to c and from left to rght: (a) U, Th and He abundances and pooled age (11.28±0.14Ma).[He ]stheα-ejecton corrected He-abundance ([He ] [He]/F t ); (b) The U, Th and He concentratons, requred for the calculaton of a (U Th)/He sochron. The sochron age s 12.0±4.2 Ma. (c) lograto transformed data (,W) and the (dmensonless) geometrc mean composton, resultng n the central age of 11.38±0.20 Ma.

7 . ermeesch / Chemcal Geology 249 (2008) Error-weghtng can be done by trval generalzatons of Eqs. (13) (15), whch are mplemented n the web-calculator. The geometrc mean composton s gven by the nverse lograto transformaton (Atchson, 1986; Weltje, 2002): ½U Š e e þ e W þ 1 ; ½ThŠ e W e þ e W þ 1 ; ½HeŠ 1 e þ e W þ 1 ð16þ Wth varances: U a 2 b 2 A b 2 c 2 2bc d 2 e 2 2de He W cov ;W 1 A where e e W þ 1 e þ W a 2 b 2 e þ e W þ 1 e þ e W þ 1 e W e þ 1 e c 2 d 2 e þ e W þ 1 e þ e W þ 1 e W e 2 e þ e W þ 1 ð17þ The central age s then smply calculated by pluggng [U ], [Th ] and [He ] and ther uncertantes nto Eqs. (2), (3) or (5). As predcted (Fg. 5), the arthmetc mean age s older than the central age. There s less than 2% dsagreement between the arthmetc mean age (~11.58 Ma) and the central age (~ Ma), and 7% dfference between the pooled age (~11.28 Ma) and the sochron age (~12.0 Ma). 5. Generalzed equatons for (U Th Sm)/He datng For reasons gven n the Introducton, 147 Sm s often neglected n helum thermochronometry. However, n rare cases t does happen that apatte contans hgh abundances of Sm, affectng the helum age on the percent level. Ths secton wll explan how to add a fourth radoactve parent to the methods descrbed above. The exact age equaton (Eq. (2)) and the present-day helum producton rate (Eq. (4)) can easly be generalzed to nclude Sm: ½HeŠ 8 137:88 138:88 ek 238t 7 1 þ 138:88 ek 235t 1 ½UŠþ6 e k232t 1 ½ThŠþ0:1499 e k 147 t 1 ½SmŠ ð18þ and 8 137:88 138:88 k 238 þ 7 138:88 k 235 ½UŠþ6k 232 ½ThŠþ0:1499k 147 ½SmŠ ð19þ Wth λ 147 the decay constant of 147 Sm and all other parameters as n Eqs. (2) and (4). Usng Eq. (19), calculatng an sochron age for (U Th Sm)/He proceeds n exactly the same way as for the ordnary (U Th)/He method, and the same s true for the pooled age (Fg. 6). Calculatng (U Th Sm)/He central ages s also very smlar, although the equatons are a bt longer. In addton to and W (Eq. (11)), we defne a thrd lograto varable X (1 n): ½U Š ½Th Š ½Sm Š ln ; W ln ; X ln ½He Š ½He Š ½He Š ð20þ Because there are three nstead of two lograto varables, the (U Th Sm)/He age equaton cannot be vsualzed on a straghtforward bvarate dagram, but forms a set of hypersurfaces n trvarate lograto-space (Fg. 7). Lkewse, (U Th Sm)/He data do not form a ternary, but a tetrahedral system n Fg. 6. (a) Ternary dagram of the Naxos data (Table 1). (b) The same data plotted n lograto-space. Error ellpes are 2σ.

8 346. ermeesch / Chemcal Geology 249 (2008) Fg. 7. The lograto transformaton can easly be generalzed to the case of 3 radoactve parents but ths precludes a straghtforward vsualzaton on a two-dmensonal dagram. (a) the (U Th Sm)/He age equaton n compostonal dataspace. (b) two hypersurfaces representng two ages n lograto-space. compostonal dataspace (Fg. 7). Generalzng the (co)varances of Eq. (12): ru 2 ; rˆ W 2 U He rth 2 ; Th He X rsm 2 ð21þ Sm He cov ;W cov ;X cov W ;X rhe 2 He Calculatng the arthmetc lograto-means: 1 n ; W 1 n W ; X 1 n X ð22þ The (co-)varances of the lograto-means, propagatng only the nternal error: 1 n 2 cov ;W 1 n 2 cov W ;X 1 n 2 ; W 1 n 2 W ; X 1 n 2 cov ;W ; cov ;X 1 n 2 cov W ;X cov ;X ; X ð23þ The (co-)varances of the lograto-means, propagatng the external error: n 2 n ; rˆ2 W W W 2 nn ð 1Þ nn ð 1Þ n cov ;W W W ;cov nn ð 1Þ ;X n cov W ; X W W X X nn ð 1Þ n ; X X X 2 nn ð 1Þ n ð nn ð 1Þ Þ X X ð24þ The nverse lograto transformaton: e ½U Š e þ e W þ e X þ 1 ; ½ThŠ e þ e W þ e X þ 1 e ½SmŠ X e þ e W þ e X þ 1 ; ½HeŠ 1 e þ e W þ e X þ 1 ð25þ Fnally, calculatng the standard error propagaton of the geometrc mean compostons: rˆ U a rˆ 2 2 b 2 c 2 2ab 2ac 2bc Th C Sm A b 2 d 2 e 2 W 2bd 2be 2de B c 2 e 2 f 2 X 2ce 2cf 2ef A cov B ;W C rˆ 2 g 2 h 2 2 2gh 2g cov A He ;X cov W ;X wth e e W þ e X þ 1 e þ W a 2 ; b 2 e þ e W þ e X þ 1 e þ e W þ e X þ 1 e þ X e W e þ e X þ 1 c 2 ; d 2 e þ e W þ e X þ 1 e þ e W þ e X þ 1 e W þ X e X e þ e W þ 1 e 2 ; f 2 e þ e W þ e X þ 1 e þ e W þ e X þ 1 e g 2 ; h 2 e þ e W þ e X þ 1 e þ e W þ e X þ 1 e X 2 e þ e W þ e X þ 1 e W e W ð26þ An example of a well-behaved (U Th Sm)/He dataset from the Fsh Lake alley apatte standard (provded by rof. Danel Stockl, Unversty of Kansas) s gven n the web-calculator (http://pvermees.andropov.org/central). The arthmetc mean of 28 sngle-gran ages s 6.36±0.11 Ma, the pooled age 6.43±

9 . ermeesch / Chemcal Geology 249 (2008) Ma, the sochron ages 6.44±0.67 Ma (wth an ntercept of 0.005±0.056 fmol/µg, and the central age 6.41±0.14 Ma. Note that the central age s older and not younger than the arthmetc mean age. Ths ndcates that random varatons exceed the very small systematc dfference between the arthmetc and geometrc mean compostons. However, the central age probably stll s more accurate than the arthmetc mean age because t s less senstve to outlers. 6. Conclusons Ths paper compared three ways to calculate an age from a sngle set of U, Th and He measurements and four ways to calculate the average of several alquots of the same sample. U, Th and He form a ternary system, and the ternary dagram was ntroduced as an elegant way to make such a comparson. Ths reveals that the accuracy ganed by the exact soluton of the (U Th)/He equaton s easly lost f the average age of replcate measurements s calculated by the arthmetc mean. As a better alternatve, the central age s calculated from the geometrc mean composton of a dataset. In addton to the central age, the paper also ntroduced the pooled age and the sochron age as valuable alternatves to the arthmetc mean age n certan applcatons. The pooled age s calculated by addng the U, Th, (Sm) and He contents of several sngle-and/or mult-gran alquots of the same sample. ooled ages are based to hgh U Th-grans whch may be affected by radaton damage, but are the only sensble way to average mult-gran alquots. The sochron age s gven by the slope of a lnear ft of a dagram that plots helum content aganst present-day helum producton rate. In order to reduce the bas towards large grans, t s a good dea to transform the nput data to unts of concentraton, whch can be done by dvdng the atomc abundances by the estmated volume or mass of the component grans. Dong so wll translate the dataponts along a straght lne through the orgn of the sochron plot and mproves ts power for detectng parentless helum. If there s no parentless helum, data must plot on a sngle lne gong through the orgn. But f, nstead, the data do not defne a lne, or ths lne does not go through the orgn of the sochron dagram, parentless helum or a smlar problem may be present. The sochron age s less well suted for ages older than 100 Ma because t uses the lnearzed age equaton (Eq. (3)). Although most (U Th)/He geochronologsts are probably already aware that some accuracy s lost by calculatng the α- ejecton correcton by smply dvdng the uncorrected (U Th)/ He age by the α-retenton factor F t, t bears repeatng that nstead, the measured helum concentraton should be dvded by F t before the age calculaton. The effects dscussed n ths paper are relatvely mnor, affectng the calculated ages by at most a few percent. Nevertheless, the added computatonal cost of followng the above recommendatons pales n comparson wth the cost of collectng, separatng and analyzng samples. Therefore, there s no reason why not to gan the extra percent of accuracy. To facltate the calculaton of the central age, a webbased calculator s provded at central. It mplements the calculatons of central ages wth or wthout Sm, and offers several optons for propagatng nternal or external uncertantes. The web-calculator also allows the calculaton of error-weghted central ages, and ncludes two dataset for testng purposes: the ncluson-bearng (U Th)/He data from Naxos whch s also summarzed n Table 1, and (U Th Sm)/He data from a Fsh Lake alley apatte lab standard. Acknowledgments I would lke to thank Jeremy Hourgan and an anonymous revewer for the constructve comments, and Danny Stockl for sharng hs Fsh Lake alley (U Th Sm)/He data. Ths work was done whle the author was a Mare Cure postdoctoral fellow at ETH-Zürch n the framework of the CRONUS-EU ntatve (RTN project reference ). References Atchson, J., The Statstcal Analyss of Compostonal Data. Chapman and Hall, London. 416 pp. Atchson, J., Brown, J.A.C., The Lognormal Dstrbuton. Cambrdge Unversty ress, London. 176 pp. Ftzgerald,.G., Baldwn, S.L., Webb, L.E., O Sullvan,.B., Interpretaton of (U Th)/He sngle gran ages from slowly cooled crustal terranes: a case study from the Transantarctc Mountans of southern ctora Land. Chemcal Geology 225, House, M.A., Farley, K.A., Stockl, D., Helum chronometry of apatte and ttante usng Nd-YAG laser heatng. Earth and lanetary Scence Letters 183 (3 4), Ludwg, K.R., Isoplot 3.00 user manual, Berkeley Geochronology Center Specal ublcaton, vol. 4. Meesters, A.G.C.A., Duna, T.J., A nonteratve soluton of the (U Th)/ He age equaton. Geochemstry, Geophyscs, Geosystems 6 (4). ersano, C., Barfod, D.N., Stuart, F.M., Bshop,., Constrants on early Cenozoc underplatng-drven uplft and denudaton of western Scotland from low temperature thermochronometry. Earth and lanetary Scence Letters 263, Steger, R.H., Jäger, E., Subcommsson on geochronology: conventon on the use of decay constants n geo-and cosmochronology. Earth and lanetary Scence Letters 36, Stock, G.M., Ehlers, T.A., Farley, K.A., Where does sedment come from? Quantfyng catchment eroson wth detrtal apatte (U Th)/He thermochronometry. Geology 34, ermeesch,., Seward, D., Latkoczy, C., Wpf, M., Guenther, D., Baur, H., Alpha-emttng mneral nclusons n apatte, ther effect on (U Th)/ He ages, and how to reduce t. Geochmca et Cosmochmca Acta 71, Weltje, G.J., Quanttatve analyss of detrtal modes; statstcally rgorous confdence regons n ternary dagrams and ther use n sedmentary petrology. Earth-Scence Revews 57 (3 4),

NUMERICAL DIFFERENTIATION

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

More information

x i1 =1 for all i (the constant ).

x i1 =1 for all i (the constant ). Chapter 5 The Multple Regresson Model Consder an economc model where the dependent varable s a functon of K explanatory varables. The economc model has the form: y = f ( x,x,..., ) xk Approxmate ths by

More information

Analytical Chemistry Calibration Curve Handout

Analytical Chemistry Calibration Curve Handout I. Quck-and Drty Excel Tutoral Analytcal Chemstry Calbraton Curve Handout For those of you wth lttle experence wth Excel, I ve provded some key technques that should help you use the program both for problem

More information

Negative Binomial Regression

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

More information

AS-Level Maths: Statistics 1 for Edexcel

AS-Level Maths: Statistics 1 for Edexcel 1 of 6 AS-Level Maths: Statstcs 1 for Edecel S1. Calculatng means and standard devatons Ths con ndcates the slde contans actvtes created n Flash. These actvtes are not edtable. For more detaled nstructons,

More information

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

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

More information

Chapter 8 Indicator Variables

Chapter 8 Indicator Variables Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n

More information

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method

Comparison of the Population Variance Estimators. of 2-Parameter Exponential Distribution Based on. Multiple Criteria Decision Making Method Appled Mathematcal Scences, Vol. 7, 0, no. 47, 07-0 HIARI Ltd, www.m-hkar.com Comparson of the Populaton Varance Estmators of -Parameter Exponental Dstrbuton Based on Multple Crtera Decson Makng Method

More information

THE CURRENT BALANCE Physics 258/259

THE CURRENT BALANCE Physics 258/259 DSH 1988, 005 THE CURRENT BALANCE Physcs 58/59 The tme average force between two parallel conductors carryng an alternatng current s measured by balancng ths force aganst the gravtatonal force on a set

More information

4DVAR, according to the name, is a four-dimensional variational method.

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

More information

Bayesian predictive Configural Frequency Analysis

Bayesian predictive Configural Frequency Analysis Psychologcal Test and Assessment Modelng, Volume 54, 2012 (3), 285-292 Bayesan predctve Confgural Frequency Analyss Eduardo Gutérrez-Peña 1 Abstract Confgural Frequency Analyss s a method for cell-wse

More information

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

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

More information

F statistic = s2 1 s 2 ( F for Fisher )

F statistic = s2 1 s 2 ( F for Fisher ) Stat 4 ANOVA Analyss of Varance /6/04 Comparng Two varances: F dstrbuton Typcal Data Sets One way analyss of varance : example Notaton for one way ANOVA Comparng Two varances: F dstrbuton We saw that the

More information

RELIABILITY ASSESSMENT

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

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

Systematic Error Illustration of Bias. Sources of Systematic Errors. Effects of Systematic Errors 9/23/2009. Instrument Errors Method Errors Personal

Systematic Error Illustration of Bias. Sources of Systematic Errors. Effects of Systematic Errors 9/23/2009. Instrument Errors Method Errors Personal 9/3/009 Sstematc Error Illustraton of Bas Sources of Sstematc Errors Instrument Errors Method Errors Personal Prejudce Preconceved noton of true value umber bas Prefer 0/5 Small over large Even over odd

More information

Case A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k.

Case A. P k = Ni ( 2L i k 1 ) + (# big cells) 10d 2 P k. THE CELLULAR METHOD In ths lecture, we ntroduce the cellular method as an approach to ncdence geometry theorems lke the Szemeréd-Trotter theorem. The method was ntroduced n the paper Combnatoral complexty

More information

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

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

More information

III. Econometric Methodology Regression Analysis

III. Econometric Methodology Regression Analysis Page Econ07 Appled Econometrcs Topc : An Overvew of Regresson Analyss (Studenmund, Chapter ) I. The Nature and Scope of Econometrcs. Lot s of defntons of econometrcs. Nobel Prze Commttee Paul Samuelson,

More information

Cokriging Partial Grades - Application to Block Modeling of Copper Deposits

Cokriging Partial Grades - Application to Block Modeling of Copper Deposits Cokrgng Partal Grades - Applcaton to Block Modelng of Copper Deposts Serge Séguret 1, Julo Benscell 2 and Pablo Carrasco 2 Abstract Ths work concerns mneral deposts made of geologcal bodes such as breccas

More information

Markov Chain Monte Carlo Lecture 6

Markov Chain Monte Carlo Lecture 6 where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways

More information

Color Rendering Uncertainty

Color Rendering Uncertainty Australan Journal of Basc and Appled Scences 4(10): 4601-4608 010 ISSN 1991-8178 Color Renderng Uncertanty 1 A.el Bally M.M. El-Ganany 3 A. Al-amel 1 Physcs Department Photometry department- NIS Abstract:

More information

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding

Lecture 9: Linear regression: centering, hypothesis testing, multiple covariates, and confounding Recall: man dea of lnear regresson Lecture 9: Lnear regresson: centerng, hypothess testng, multple covarates, and confoundng Sandy Eckel seckel@jhsph.edu 6 May 8 Lnear regresson can be used to study an

More information

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for U Charts. Dr. Wayne A. Taylor

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for U Charts. Dr. Wayne A. Taylor Taylor Enterprses, Inc. Adjusted Control Lmts for U Charts Copyrght 207 by Taylor Enterprses, Inc., All Rghts Reserved. Adjusted Control Lmts for U Charts Dr. Wayne A. Taylor Abstract: U charts are used

More information

Chapter 4: Regression With One Regressor

Chapter 4: Regression With One Regressor Chapter 4: Regresson Wth One Regressor Copyrght 2011 Pearson Addson-Wesley. All rghts reserved. 1-1 Outlne 1. Fttng a lne to data 2. The ordnary least squares (OLS) lne/regresson 3. Measures of ft 4. Populaton

More information

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry Workshop: Approxmatng energes and wave functons Quantum aspects of physcal chemstry http://quantum.bu.edu/pltl/6/6.pdf Last updated Thursday, November 7, 25 7:9:5-5: Copyrght 25 Dan Dll (dan@bu.edu) Department

More information

x yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting.

x yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting. The Practce of Statstcs, nd ed. Chapter 14 Inference for Regresson Introducton In chapter 3 we used a least-squares regresson lne (LSRL) to represent a lnear relatonshp etween two quanttatve explanator

More information

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram

Adiabatic Sorption of Ammonia-Water System and Depicting in p-t-x Diagram Adabatc Sorpton of Ammona-Water System and Depctng n p-t-x Dagram J. POSPISIL, Z. SKALA Faculty of Mechancal Engneerng Brno Unversty of Technology Techncka 2, Brno 61669 CZECH REPUBLIC Abstract: - Absorpton

More information

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger JAB Chan Long-tal clams development ASTIN - September 2005 B.Verder A. Klnger Outlne Chan Ladder : comments A frst soluton: Munch Chan Ladder JAB Chan Chan Ladder: Comments Black lne: average pad to ncurred

More information

Singular Value Decomposition: Theory and Applications

Singular Value Decomposition: Theory and Applications Sngular Value Decomposton: Theory and Applcatons Danel Khashab Sprng 2015 Last Update: March 2, 2015 1 Introducton A = UDV where columns of U and V are orthonormal and matrx D s dagonal wth postve real

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem

Speeding up Computation of Scalar Multiplication in Elliptic Curve Cryptosystem H.K. Pathak et. al. / (IJCSE) Internatonal Journal on Computer Scence and Engneerng Speedng up Computaton of Scalar Multplcaton n Ellptc Curve Cryptosystem H. K. Pathak Manju Sangh S.o.S n Computer scence

More information

18. SIMPLE LINEAR REGRESSION III

18. SIMPLE LINEAR REGRESSION III 8. SIMPLE LINEAR REGRESSION III US Domestc Beers: Calores vs. % Alcohol Ftted Values and Resduals To each observed x, there corresponds a y-value on the ftted lne, y ˆ ˆ = α + x. The are called ftted values.

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Effective plots to assess bias and precision in method comparison studies

Effective plots to assess bias and precision in method comparison studies Effectve plots to assess bas and precson n method comparson studes Bern, November, 016 Patrck Taffé, PhD Insttute of Socal and Preventve Medcne () Unversty of Lausanne, Swtzerland Patrck.Taffe@chuv.ch

More information

Entropy generation in a chemical reaction

Entropy generation in a chemical reaction Entropy generaton n a chemcal reacton E Mranda Área de Cencas Exactas COICET CCT Mendoza 5500 Mendoza, rgentna and Departamento de Físca Unversdad aconal de San Lus 5700 San Lus, rgentna bstract: Entropy

More information

Electrical double layer: revisit based on boundary conditions

Electrical double layer: revisit based on boundary conditions Electrcal double layer: revst based on boundary condtons Jong U. Km Department of Electrcal and Computer Engneerng, Texas A&M Unversty College Staton, TX 77843-318, USA Abstract The electrcal double layer

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor

Copyright 2017 by Taylor Enterprises, Inc., All Rights Reserved. Adjusted Control Limits for P Charts. Dr. Wayne A. Taylor Taylor Enterprses, Inc. Control Lmts for P Charts Copyrght 2017 by Taylor Enterprses, Inc., All Rghts Reserved. Control Lmts for P Charts Dr. Wayne A. Taylor Abstract: P charts are used for count data

More information

e i is a random error

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

More information

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM

DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM Ganj, Z. Z., et al.: Determnaton of Temperature Dstrbuton for S111 DETERMINATION OF TEMPERATURE DISTRIBUTION FOR ANNULAR FINS WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY BY HPM by Davood Domr GANJI

More information

Meta-Analysis of Correlated Proportions

Meta-Analysis of Correlated Proportions NCSS Statstcal Softare Chapter 457 Meta-Analyss of Correlated Proportons Introducton Ths module performs a meta-analyss of a set of correlated, bnary-event studes. These studes usually come from a desgn

More information

STATISTICS QUESTIONS. Step by Step Solutions.

STATISTICS QUESTIONS. Step by Step Solutions. STATISTICS QUESTIONS Step by Step Solutons www.mathcracker.com 9//016 Problem 1: A researcher s nterested n the effects of famly sze on delnquency for a group of offenders and examnes famles wth one to

More information

SIMPLE REACTION TIME AS A FUNCTION OF TIME UNCERTAINTY 1

SIMPLE REACTION TIME AS A FUNCTION OF TIME UNCERTAINTY 1 Journal of Expermental Vol. 5, No. 3, 1957 Psychology SIMPLE REACTION TIME AS A FUNCTION OF TIME UNCERTAINTY 1 EDMUND T. KLEMMER Operatonal Applcatons Laboratory, Ar Force Cambrdge Research Center An earler

More information

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY.

A PROCEDURE FOR SIMULATING THE NONLINEAR CONDUCTION HEAT TRANSFER IN A BODY WITH TEMPERATURE DEPENDENT THERMAL CONDUCTIVITY. Proceedngs of the th Brazlan Congress of Thermal Scences and Engneerng -- ENCIT 006 Braz. Soc. of Mechancal Scences and Engneerng -- ABCM, Curtba, Brazl,- Dec. 5-8, 006 A PROCEDURE FOR SIMULATING THE NONLINEAR

More information

Lecture 12: Classification

Lecture 12: Classification Lecture : Classfcaton g Dscrmnant functons g The optmal Bayes classfer g Quadratc classfers g Eucldean and Mahalanobs metrcs g K Nearest Neghbor Classfers Intellgent Sensor Systems Rcardo Guterrez-Osuna

More information

experimenteel en correlationeel onderzoek

experimenteel en correlationeel onderzoek expermenteel en correlatoneel onderzoek lecture 6: one-way analyss of varance Leary. Introducton to Behavoral Research Methods. pages 246 271 (chapters 10 and 11): conceptual statstcs Moore, McCabe, and

More information

Thermodynamics General

Thermodynamics General Thermodynamcs General Lecture 1 Lecture 1 s devoted to establshng buldng blocks for dscussng thermodynamcs. In addton, the equaton of state wll be establshed. I. Buldng blocks for thermodynamcs A. Dmensons,

More information

Week 2. This week, we covered operations on sets and cardinality.

Week 2. This week, we covered operations on sets and cardinality. Week 2 Ths week, we covered operatons on sets and cardnalty. Defnton 0.1 (Correspondence). A correspondence between two sets A and B s a set S contaned n A B = {(a, b) a A, b B}. A correspondence from

More information

Topic 23 - Randomized Complete Block Designs (RCBD)

Topic 23 - Randomized Complete Block Designs (RCBD) Topc 3 ANOVA (III) 3-1 Topc 3 - Randomzed Complete Block Desgns (RCBD) Defn: A Randomzed Complete Block Desgn s a varant of the completely randomzed desgn (CRD) that we recently learned. In ths desgn,

More information

Fundamental loop-current method using virtual voltage sources technique for special cases

Fundamental loop-current method using virtual voltage sources technique for special cases Fundamental loop-current method usng vrtual voltage sources technque for specal cases George E. Chatzaraks, 1 Marna D. Tortorel 1 and Anastasos D. Tzolas 1 Electrcal and Electroncs Engneerng Departments,

More information

arxiv:quant-ph/ Jul 2002

arxiv:quant-ph/ Jul 2002 Lnear optcs mplementaton of general two-photon proectve measurement Andrze Grudka* and Anton Wóck** Faculty of Physcs, Adam Mckewcz Unversty, arxv:quant-ph/ 9 Jul PXOWRZVNDR]QDRODQG Abstract We wll present

More information

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

CHAPTER 8. Exercise Solutions

CHAPTER 8. Exercise Solutions CHAPTER 8 Exercse Solutons 77 Chapter 8, Exercse Solutons, Prncples of Econometrcs, 3e 78 EXERCISE 8. When = N N N ( x x) ( x x) ( x x) = = = N = = = N N N ( x ) ( ) ( ) ( x x ) x x x x x = = = = Chapter

More information

Min Cut, Fast Cut, Polynomial Identities

Min Cut, Fast Cut, Polynomial Identities Randomzed Algorthms, Summer 016 Mn Cut, Fast Cut, Polynomal Identtes Instructor: Thomas Kesselhem and Kurt Mehlhorn 1 Mn Cuts n Graphs Lecture (5 pages) Throughout ths secton, G = (V, E) s a mult-graph.

More information

Chemometrics. Unit 2: Regression Analysis

Chemometrics. Unit 2: Regression Analysis Chemometrcs Unt : Regresson Analyss The problem of predctng the average value of one varable n terms of the known values of other varables s the problem of regresson. In carryng out a regresson analyss

More information

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens THE CHINESE REMAINDER THEOREM KEITH CONRAD We should thank the Chnese for ther wonderful remander theorem. Glenn Stevens 1. Introducton The Chnese remander theorem says we can unquely solve any par of

More information

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva Econ 39 - Statstcal Propertes of the OLS estmator Sanjaya DeSlva September, 008 1 Overvew Recall that the true regresson model s Y = β 0 + β 1 X + u (1) Applyng the OLS method to a sample of data, we estmate

More information

2 Finite difference basics

2 Finite difference basics Numersche Methoden 1, WS 11/12 B.J.P. Kaus 2 Fnte dfference bascs Consder the one- The bascs of the fnte dfference method are best understood wth an example. dmensonal transent heat conducton equaton T

More information

F8: Heteroscedasticity

F8: Heteroscedasticity F8: Heteroscedastcty Feng L Department of Statstcs, Stockholm Unversty What s so-called heteroscedastcty In a lnear regresson model, we assume the error term has a normal dstrbuton wth mean zero and varance

More information

Introduction to Econometrics (3 rd Updated Edition, Global Edition) Solutions to Odd-Numbered End-of-Chapter Exercises: Chapter 13

Introduction to Econometrics (3 rd Updated Edition, Global Edition) Solutions to Odd-Numbered End-of-Chapter Exercises: Chapter 13 Introducton to Econometrcs (3 rd Updated Edton, Global Edton by James H. Stock and Mark W. Watson Solutons to Odd-Numbered End-of-Chapter Exercses: Chapter 13 (Ths verson August 17, 014 Stock/Watson -

More information

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing

Pop-Click Noise Detection Using Inter-Frame Correlation for Improved Portable Auditory Sensing Advanced Scence and Technology Letters, pp.164-168 http://dx.do.org/10.14257/astl.2013 Pop-Clc Nose Detecton Usng Inter-Frame Correlaton for Improved Portable Audtory Sensng Dong Yun Lee, Kwang Myung Jeon,

More information

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy

Comparative Studies of Law of Conservation of Energy. and Law Clusters of Conservation of Generalized Energy Comparatve Studes of Law of Conservaton of Energy and Law Clusters of Conservaton of Generalzed Energy No.3 of Comparatve Physcs Seres Papers Fu Yuhua (CNOOC Research Insttute, E-mal:fuyh1945@sna.com)

More information

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

More information

Chapter 4 High Breakdown Regression Procedures

Chapter 4 High Breakdown Regression Procedures Chapter 4 Hgh Breakdown Regresson Procedures Introducton Whle M and BI estmators provde an mprovement over OLS f the data has an outler or hgh nfluence pont, they cannot provde protecton aganst data wth

More information

Grey prediction model in world women s pentathlon performance prediction applied research

Grey prediction model in world women s pentathlon performance prediction applied research Avalable onlne www.jocpr.com Journal of Chemcal and Pharmaceutcal Research, 4, 6(6):36-4 Research Artcle ISSN : 975-7384 CODEN(USA) : JCPRC5 Grey predcton model n world women s pentathlon performance predcton

More information

Why Bayesian? 3. Bayes and Normal Models. State of nature: class. Decision rule. Rev. Thomas Bayes ( ) Bayes Theorem (yes, the famous one)

Why Bayesian? 3. Bayes and Normal Models. State of nature: class. Decision rule. Rev. Thomas Bayes ( ) Bayes Theorem (yes, the famous one) Why Bayesan? 3. Bayes and Normal Models Alex M. Martnez alex@ece.osu.edu Handouts Handoutsfor forece ECE874 874Sp Sp007 If all our research (n PR was to dsappear and you could only save one theory, whch

More information

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

Pulse Coded Modulation

Pulse Coded Modulation Pulse Coded Modulaton PCM (Pulse Coded Modulaton) s a voce codng technque defned by the ITU-T G.711 standard and t s used n dgtal telephony to encode the voce sgnal. The frst step n the analog to dgtal

More information

1 The Mistake Bound Model

1 The Mistake Bound Model 5-850: Advanced Algorthms CMU, Sprng 07 Lecture #: Onlne Learnng and Multplcatve Weghts February 7, 07 Lecturer: Anupam Gupta Scrbe: Bryan Lee,Albert Gu, Eugene Cho he Mstake Bound Model Suppose there

More information

Thermo-Calc Software. Modelling Multicomponent Precipitation Kinetics with CALPHAD-Based Tools. EUROMAT2013, September 8-13, 2013 Sevilla, Spain

Thermo-Calc Software. Modelling Multicomponent Precipitation Kinetics with CALPHAD-Based Tools. EUROMAT2013, September 8-13, 2013 Sevilla, Spain Modellng Multcomponent Precptaton Knetcs wth CALPHAD-Based Tools Kasheng Wu 1, Gustaf Sterner 2, Qng Chen 2, Åke Jansson 2, Paul Mason 1, Johan Bratberg 2 and Anders Engström 2 1 Inc., 2 AB EUROMAT2013,

More information

Module 3: Element Properties Lecture 1: Natural Coordinates

Module 3: Element Properties Lecture 1: Natural Coordinates Module 3: Element Propertes Lecture : Natural Coordnates Natural coordnate system s bascally a local coordnate system whch allows the specfcaton of a pont wthn the element by a set of dmensonless numbers

More information

Lab 4: Two-level Random Intercept Model

Lab 4: Two-level Random Intercept Model BIO 656 Lab4 009 Lab 4: Two-level Random Intercept Model Data: Peak expratory flow rate (pefr) measured twce, usng two dfferent nstruments, for 17 subjects. (from Chapter 1 of Multlevel and Longtudnal

More information

Singer & Willett, 2003 October 13, 2003

Singer & Willett, 2003 October 13, 2003 Snger & Wllett, October, Dong Data Analyss n n the the Multlevel Model for for Change Judy Snger & John Wllett Harvard Unversty Graduate School of Educaton What What we we wll wll cover? cover? Composte

More information

Code_Aster. Identification of the model of Weibull

Code_Aster. Identification of the model of Weibull Verson Ttre : Identfcaton du modèle de Webull Date : 2/09/2009 Page : /8 Responsable : PARROT Aurore Clé : R70209 Révson : Identfcaton of the model of Webull Summary One tackles here the problem of the

More information

Explicit constructions of all separable two-qubits density matrices and related problems for three-qubits systems

Explicit constructions of all separable two-qubits density matrices and related problems for three-qubits systems Explct constructons of all separable two-qubts densty matrces and related problems for three-qubts systems Y. en-ryeh and. Mann Physcs Department, Technon-Israel Insttute of Technology, Hafa 2000, Israel

More information

Energy Storage Elements: Capacitors and Inductors

Energy Storage Elements: Capacitors and Inductors CHAPTER 6 Energy Storage Elements: Capactors and Inductors To ths pont n our study of electronc crcuts, tme has not been mportant. The analyss and desgns we hae performed so far hae been statc, and all

More information

Problem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)

Problem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs) Problem adapted reduced models based on Reacton-Dffuson Manfolds (REDIMs) V Bykov, U Maas Thrty-Second Internatonal Symposum on ombuston, Montreal, anada, 3-8 August, 8 Problem Statement: Smulaton of reactng

More information

Andreas C. Drichoutis Agriculural University of Athens. Abstract

Andreas C. Drichoutis Agriculural University of Athens. Abstract Heteroskedastcty, the sngle crossng property and ordered response models Andreas C. Drchouts Agrculural Unversty of Athens Panagots Lazards Agrculural Unversty of Athens Rodolfo M. Nayga, Jr. Texas AMUnversty

More information

of Nebraska - Lincoln

of Nebraska - Lincoln Unversty of Nebraska - Lncoln DgtalCommons@Unversty of Nebraska - Lncoln MAT Exam Expostory Papers Math n the Mddle Insttute Partnershp 008 The Square Root of Tffany Lothrop Unversty of Nebraska-Lncoln

More information

The Gaussian classifier. Nuno Vasconcelos ECE Department, UCSD

The Gaussian classifier. Nuno Vasconcelos ECE Department, UCSD he Gaussan classfer Nuno Vasconcelos ECE Department, UCSD Bayesan decson theory recall that we have state of the world X observatons g decson functon L[g,y] loss of predctng y wth g Bayes decson rule s

More information

CHAPTER 4. Vector Spaces

CHAPTER 4. Vector Spaces man 2007/2/16 page 234 CHAPTER 4 Vector Spaces To crtcze mathematcs for ts abstracton s to mss the pont entrel. Abstracton s what makes mathematcs work. Ian Stewart The man am of ths tet s to stud lnear

More information

ARTICLE IN PRESS. Available online at Mathematics and Computers in Simulation xxx (2009) xxx xxx

ARTICLE IN PRESS. Available online at  Mathematics and Computers in Simulation xxx (2009) xxx xxx Avalable onlne at wwwscencedrectcom Mathematcs and Computers n Smulaton xxx 29 xxx xxx Dervatve based global senstvty measures and ther lnk wth global senstvty ndces IM Sobol a, S Kucherenko b, a Insttute

More information

Chapter 8 Multivariate Regression Analysis

Chapter 8 Multivariate Regression Analysis Chapter 8 Multvarate Regresson Analyss 8.3 Multple Regresson wth K Independent Varables 8.4 Sgnfcance tests of Parameters Populaton Regresson Model For K ndependent varables, the populaton regresson and

More information

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM

FUZZY GOAL PROGRAMMING VS ORDINARY FUZZY PROGRAMMING APPROACH FOR MULTI OBJECTIVE PROGRAMMING PROBLEM Internatonal Conference on Ceramcs, Bkaner, Inda Internatonal Journal of Modern Physcs: Conference Seres Vol. 22 (2013) 757 761 World Scentfc Publshng Company DOI: 10.1142/S2010194513010982 FUZZY GOAL

More information

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

Lossy Compression. Compromise accuracy of reconstruction for increased compression. Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost

More information

Now we relax this assumption and allow that the error variance depends on the independent variables, i.e., heteroskedasticity

Now we relax this assumption and allow that the error variance depends on the independent variables, i.e., heteroskedasticity ECON 48 / WH Hong Heteroskedastcty. Consequences of Heteroskedastcty for OLS Assumpton MLR. 5: Homoskedastcty var ( u x ) = σ Now we relax ths assumpton and allow that the error varance depends on the

More information

Exponential Type Product Estimator for Finite Population Mean with Information on Auxiliary Attribute

Exponential Type Product Estimator for Finite Population Mean with Information on Auxiliary Attribute Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 193-9466 Vol. 10, Issue 1 (June 015), pp. 106-113 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) Exponental Tpe Product Estmator

More information

Scroll Generation with Inductorless Chua s Circuit and Wien Bridge Oscillator

Scroll Generation with Inductorless Chua s Circuit and Wien Bridge Oscillator Latest Trends on Crcuts, Systems and Sgnals Scroll Generaton wth Inductorless Chua s Crcut and Wen Brdge Oscllator Watcharn Jantanate, Peter A. Chayasena, and Sarawut Sutorn * Abstract An nductorless Chua

More information

Confidence Intervals for the Overall Effect Size in Random-Effects Meta-Analysis

Confidence Intervals for the Overall Effect Size in Random-Effects Meta-Analysis Psychologcal Methods 008, Vol. 13, No. 1, 31 48 Copyrght 008 by the Amercan Psychologcal Assocaton 108-989X/08/$1.00 DOI: 10.1037/108-989X.13.1.31 Confdence Intervals for the Overall Effect Sze n Random-Effects

More information

ECONOMICS 351* -- Stata 10 Tutorial 6. Stata 10 Tutorial 6

ECONOMICS 351* -- Stata 10 Tutorial 6. Stata 10 Tutorial 6 ECONOMICS 35* -- Stata Tutoral 6 Stata Tutoral 6 TOPICS: Functonal Form and Varable Re-scalng n Smple Lnear Regresson Models, and An Introducton to Multple Lnear Regresson Models DATA: auto.dta (a Stata-format

More information

Code_Aster. Identification of the Summarized

Code_Aster. Identification of the Summarized Verson Ttre : Identfcaton du modèle de Webull Date : 2/09/2009 Page : /8 Responsable : Aurore PARROT Clé : R70209 Révson : 609 Identfcaton of the Summarzed Webull model One tackles here the problem of

More information

STK4080/9080 Survival and event history analysis

STK4080/9080 Survival and event history analysis SK48/98 Survval and event hstory analyss Lecture 7: Regresson modellng Relatve rsk regresson Regresson models Assume that we have a sample of n ndvduals, and let N (t) count the observed occurrences of

More information

Available online at Kazuhiro Tsuboi a, * Received 31 January 2010; revised 7 March 2010; accepted 21 March 2010

Available online at  Kazuhiro Tsuboi a, * Received 31 January 2010; revised 7 March 2010; accepted 21 March 2010 Avalable onlne at www.scencedrect.com Proceda Proceda Engneerng 00 (010) (009) 000 000 305 310 Proceda Engneerng www.elsever.com/locate/proceda 8 th Conference of the Internatonal Sports Engneerng Assocaton

More information

XII.3 The EM (Expectation-Maximization) Algorithm

XII.3 The EM (Expectation-Maximization) Algorithm XII.3 The EM (Expectaton-Maxzaton) Algorth Toshnor Munaata 3/7/06 The EM algorth s a technque to deal wth varous types of ncoplete data or hdden varables. It can be appled to a wde range of learnng probles

More information

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed

Irregular vibrations in multi-mass discrete-continuous systems torsionally deformed (2) 4 48 Irregular vbratons n mult-mass dscrete-contnuous systems torsonally deformed Abstract In the paper rregular vbratons of dscrete-contnuous systems consstng of an arbtrary number rgd bodes connected

More information

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI

Study on Active Micro-vibration Isolation System with Linear Motor Actuator. Gong-yu PAN, Wen-yan GU and Dong LI 2017 2nd Internatonal Conference on Electrcal and Electroncs: echnques and Applcatons (EEA 2017) ISBN: 978-1-60595-416-5 Study on Actve Mcro-vbraton Isolaton System wth Lnear Motor Actuator Gong-yu PAN,

More information

Differentiating Gaussian Processes

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

More information

Mud-rock line estimation via robust locally weighted scattering smoothing method

Mud-rock line estimation via robust locally weighted scattering smoothing method Mud-roc lne va LOESS Mud-roc lne estmaton va robust locally weghted scatterng smoothng method A. Nassr Saeed, Laurence R. Lnes, and Gary F. Margrave ABSTRACT The robust locally weghted scatterng smoothng

More information

If two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9.

If two volatile and miscible liquids are combined to form a solution, Raoult s law is not obeyed. Use the experimental data in Table 9. 9.9 Real Solutons Exhbt Devatons from Raoult s Law If two volatle and mscble lquds are combned to form a soluton, Raoult s law s not obeyed. Use the expermental data n Table 9.3: Physcal Chemstry 00 Pearson

More information