Posterior Confidence Intervals in Linear Calibration Problems: Calibrating the Thompson Ice Core Index. (Extended Version)

Size: px
Start display at page:

Download "Posterior Confidence Intervals in Linear Calibration Problems: Calibrating the Thompson Ice Core Index. (Extended Version)"

Transcription

1 TECHNOMETRICS SUPPLEMENTARY MATERIALS Poserior Confidence Inervals in Linear Calibraion Problems: Calibraing he Thomson Ice Core Inde (Eended Version) J. Huson McCulloch Economics Dearmen Ohio Sae Universi Oc Kewords: Baesian inerval credible inerval diffuse rior aleoclimaolog ogen- 8 Medieval Warm Period Hocke Sick Absrac The auhor is indebed o Sehen McInre Jussi Collin aricians in he OSU Economerics Seminar and Saisics Colloquium and several JASA and NSF referees for helful commens and suggesions. Auhor s mcculloch.@osu.edu.

2 In calibraion roblems an eogenous sae variable and an endogenous resonse variable or ro are boh observed in a se of calibraion observaions. We wish o make inferences abou he unobserved sae variable from an addiional observaion on he resonse variable under a diffuse rior. Hoadle (970) argued ha an informaive rior is required in order o obain a roer oserior disribuion. Huner and Lambo (98) roosed a soluion bu were sharl criicized a he ime. This aer resens a new derivaion of he Huner-Lambo oserior disribuion under a diffuse rior ha mees hese objecions. A he same ime i is shown ha Hoadle s aroach was based on a suble inconsisenc in he alicaion of Baes Rule. The reinsaed Huner-Lambo oserior is alied o he roblem of calibraing he ice core inde of Thomson e al. (003) o insrumenal emeraures. I is found conrar o he famous claim of Gore (006) ha his inde is in fac uninformaive abou he quesion of wheher Medieval Warm Period was warmer or cooler han he resen. I is shown ha he classical confidence inervals roosed b Fieller (954) are a good aroimaion o he oserior confidence inervals when he calibraion sloe coefficien is highl significan relaive o he desired confidence inerval ail robabili. However when he sloe is onl marginall significan he Fieller inervals become increasingl disored and hen meaningless. The roosed oserior confidence inervals siml become wider as he sloe loses significance bu remain bounded. Eensions o mulile roies sequeniall organized daa and rior resricions on he calibraion sloe coefficien are oulined bu no imlemened.

3 I. Inroducion In he classical calibraion roblem as reviewed b Osborne (99) and Brown (993) an eogenous sae variable i is measured indirecl b an endogenous resonse variable i ha has an affine bu nois relaion o i: ε. () i i i We have a se of calibraion observaions on boh for i n. In he benchmark case considered here i is assumed ha he measuremen errors ε i are i.i.d. Gaussian wih E(ε i i ) 0 and var(ε i ). We wish o make inferences abou he unobserved sae variable for an addiional observaion or observaions ouside he calibraion se for which we have onl he observed resonse. This resonse obes he same rule ' ' ε '. () Such calibraion roblems arise in such diverse fields as hermodnamics harmacolog urolog and chemical analsis (cf. Brown 993) zoolog (du Plessis and van der Merwe 996) and aleoclimaolog (e.g. Mann e al. 999 Kaufman e al. 009). In economerics alicaions could include he reconsrucion of hisorical Gross Naional Produc from fragmenar records (e.g. Romer 989) or he hisorical sandard of living from heigh daa (e.g. Seckel 995). Solving () for ields ' ( ' ε ') / (3) so ha he naural or classical calibraion esimaor of roosed b Eisenhar (939) is ˆ ( ' ˆ) / ˆ (4) I here use and for he reconsrucion values following he nomenclaure of Brown (98).

4 where ˆ and ˆ are he Ordinar Leas Squares (OLS) esimaes of and from (). The eression in (4) is also he Maimum Likelihood esimaor of (Hoadle 970 Brown 993). The resen aer resens a new derivaion of he Huner and Lambo (HL 98) oserior disribuion for under an uninformaive diffuse rior. The resuling confidence inervals (CIs) or credible inervals as he are ofen called in he Baesian cone are alwas bounded and coniguous and ma readil be consruced wih equal ail robabiliies. The radiional or classical mehod of comuing confidence regions advocaed b mos auhors is based on a heorem of Fieller (954). However his aroach leads o unsaisfacor ses ha ma be unbounded a one or even boh ends and ma even be disconiguous. Even when he are bounded inervals he classical CIs are in mos cases ecessivel skewed awa from he mean of he calibraion values of he sae variable so ha he wo ail robabiliies are unequal. Hoadle (970) aems o consruc a Baesian alernaive bu is unable o do his wihou an unnecessaril informaive rior. The reinsaed Huner-Lambo mehod is alied o he roblem of calibraing he Thomson e al. (003) ice core isooe raio inde o insrumenal emeraure. This inde was made famous in Al Gore s An Inconvenien Truh (006) where i was said o rovide definiive confirmaion ha emeraures during he Medieval Warm Period (MWP) were lower han hose a resen. I is found ha he inde does iml a oin esimae of emeraure ha is uniforml below he insrumenal average from he beginning of he series in AD unil he 930 s. However he 95% oserior CI for he emeraure anomal icall eends from a leas -.8 C o

5 3. C so ha he reconsrucion is in fac uninformaive. Alhough his inde does no oin o a MWP i b no means rules one ou conrar o he claim of Gore. The muliro emeraure reconsrucion of Loehle and McCulloch (008) insead shows a saisicall significan MWP relaive o he bimillenial average during mos of AD lus AD as well as a significanl cool Lile Ice Age (LIA) during mos of AD. Secion II of his aer resens he new derivaion of he HL oserior disribuion for he unobserved sae variable conrass i wih he aroach of Hoadle (970) and answers he criicisms ha were raised agains he HL aroach when i was firs ublished. Secion III consrucs a close aroimaion o he Thomson e al. (003) ice core isooe inde from source daa. Secion IV calibraes his inde o global emeraure and comues oserior 50% and 95% CIs. The concluding secion V comares hese inervals o he classical Fieller regions and discusses eensions of he mehod of Secion II. II. Poserior Confidence Inervals for he Classical Calibraion Esimaor The join rior disribuion for ' and he regression arameers ( ) ha governs he enire eercise ma be facored as follows: ( ) ( ) ( ) ( ). (5) Baes Rule imlies ha he disribuion of condiional on and he rue arameers is ( ' ' ) ( ' ' ) ( ' ) (6)

6 4 so ha in he benchmark case in which he rior densi comonen ( ) ' is diffuse and herefore uninformaive (cf. Zellner 97: 4-53) ( ) ( ) ( ) ( ) ( ) ( ) ( ) / / ' N ~ / / ' ' e ' e ' ' ' ' (7) whence ' N ~ '. (8) The disribuion of he OLS esimaorsˆ and condiional on he rue arameers is ˆ S S N ~ ˆ ˆ where i i i n s s s s S. Then under he sandard diffuse rior for ( ) (Zellner 97 ch. 3) (9) S S ˆ ˆ N ~ ˆ ˆ whence. (0) s s s s ˆ ˆ N ~ ˆ ˆ S

7 5 Since ' ( ' ) ' ields ε is indeenden of he regression daa combining (8) wih (0) ˆ s ' ˆ ˆ S ~ N ˆ s s s. () I follows from ' ( ') / ha where RN ( μ Σ) ( ) ˆ s s ' ˆ ˆ ~ RN ˆ s s S () indicaes he disribuion of he raio of wo normall disribued random variables wih mean vecor μ and covariance mari Σ (Fieller 93 Hinckle 969). This disribuion has as secial cases he Cauch disribuion when he numeraor and denominaor boh have mean 0 and he Gaussian disribuion when he denominaor has zero variance. In general i is heav-ailed wih ails inermediae beween hese wo cases and is skewed unless he numeraor or denominaor has mean 0. I ma even be bimodal if he numeraor is large relaive o is sandard deviaion and he denominaor small relaive o is sandard deviaion. See Marsaglia (965) for grahs of he densi in he secial case of zero correlaion. The Cumulaive Disribuion Funcion (CDF) P ( ' ' ˆ ˆ S ) as required for confidence inervals on ma easil be comued in erms of he bivariae normal CDF as shown in Aendi I of he resen aer equaion (6). I ma easil be seen ha condiional on he rue arameers he classical oin esimae ˆ in (4) has he same RN disribuion as () onl wih he OLS esimaes of and relaced b heir rue values:

8 6 s s s ( ) ˆ S ~ RN. (3) The RN oserior disribuion for in () is herefore aricularl naural and inuiive. Williams (969) noes in he cone of (3) ha he RN disribuion has undefined mean and variance. In fac i lies he domain of aracion of he Cauch disribuion so ha even he Generalized (sable domain of aracion) Law of Large Numbers canno be relied uon o make he disribuion of he average more comac han ha of he conribuions. Care should herefore be eercised in merging mulile classical esimaes b simle averaging as in Loehle and McCulloch (008) or even b GLS as in Brown (98). In racice is unknown and mus be esimaed from he variance of he residuals s s n i e i /( n ) where e i i ˆ ˆ. Under a diffuse and herefore uninformaive rior for log( ) so i (Zellner 97: 4-53) ha ( ) Then ( s ) ~ ( n ) s / χ. ˆ ˆ S s ( n) ˆ ~ BVT ˆ s S n (4) where BVT(μ V ν) reresens he elliical bivariae Suden disribuion wih mean μ covariance mari esimae V and degrees of freedom ν (cf. Zellner 97 ch. 3 aendi B.; Genz 004). I follows from () and (4) ha

9 7 where RT( Vν ) ( ) ˆ ' s s ˆ ˆ s ~ RT ˆ s n s s S (5) μ reresens he disribuion of he raio of wo elliical bivariae Suden random variables wih mean vecor μ covariance mari esimae V and ν degrees of freedom. Is CDF P( ' ' ˆ ˆ S s ) ma easil be comued in erms of he elliical bivariae Suden CDF as also shown in Aendi I of he resen aer. 3 The resen aer firs uses Baes Rule o derive ( ' ' ) ( ' ' ) under a diffuse rior for ( ' ) from. Then he regression arameers are esimaed using he calibraion daa and he rue values are inegraed ou o obain ( ' ' ˆ ˆ S s ). Hoadle (970: ) likewise considers he use of a diffuse rior for. However he insead firs esimaes he regression arameers o obain he classical Suden disribuion wih n- degrees of freedom for funcion of : as a ( ) ( ) ( ˆ ˆ ) ' ' ˆ ˆ v / s ( n ) / S (6) ( n ) v( ) where ( ) s ( s s s ) v (7) is he esimaed variance of as a funcion of. Hoadle hen noes ha holding consan and considering his eression as a funcion of i has ails ha are asmoicall roorional o / and herefore inegrae o infini when mulilied 3 Huner and Lambo (98) reor he eac formula for he densi of he raio of wo Suden random variables bu hen one mus sill inegrae his densi numericall in order o comue CIs so ha i is easier jus o use he BVT CDF o comue he RT CDF direcl.

10 8 b a uniform rior. He concludes ha an imroer uniform rior on [] leads o a nonsensical Baes esimaor and herefore rejecs he uniform rior. Brown (98; ) and Osborne (99) concur wih Hoadle ha a roer rior is required in he case considered here of a single resonse variable. Hoadle consequenl inroduces an informaive rior ha is based on he assumion ha is drawn from he marginal disribuion of he calibraion values and shows ha his leads o he inverse esimaor favored b Kruchkoff (967) based on a calibraion regression in which he indeenden variable is regressed on he deenden variable. Alhough here are alicaions in which such a rior would be jusified he resen aer is rimaril concerned wih he imoran case of an uninformaive diffuse rior. 4 However here is a suble fallac in Hoadle s reasoning: Because he Suden densi (6) is condiioned on he regression daa Hoadle s rior for mus be as well Baes Rule requires no ( ' ' ˆ ˆ S s ) ( ' ' ˆ ˆ S s ) ( ' ˆ ˆ S s ) ( ' ' ˆ ˆ S s ) ( ' ' ˆ ˆ S s ) ( ' ). I is no immediael obvious wh condiioning he rior for ' on he regression daa would make an difference. However consider generalizing he join rior (5) ha governs he enire eercise o be condiioned on an informaion se Ω ha ma or ma no be em: 4 In he case of a sequenial sae variable ha follows a sochasic rocess wih saionar incremens an informaive rior ha deends on he neares of he calibraion sae variables is aroriae as discussed below. Even hen however his rior should be alied before he regression coefficiens are inegraed ou and no aferwards.

11 9 ( Ω) ( Ω) ( Ω) ( Ω). (8) Since Hoadle s rior for imlicil deends on he regression daa his Ω mus also conain he regression daa in each of he oher erms in (8) including ( Ω). However condiioning he rior for he regression arameers on he regression resuls is inconsisen wih using he same regression resuls o inform he oserior for he arameers. Hoadle s sequence of oeraions is herefore invalid. I has been shown here ha here are no real roblems wih he Baesian aroach under a diffuse rior rovided ha a diffuse rior for he condiional densi ( ' ) ha arises from he facorizaion (5) is alied before he regression arameers are inegraed ou raher han a diffuse rior for he condiional densi ( ' ˆ ˆ S s ) afer he rue regression arameers have been eliminaed. Huner and Lambo (HL 98) derived in densi form eacl he same RT oserior disribuion for as we do above bu using a diffuse rior for η insead of for iself as in he above derivaion. In he same issue of Technomerics Hill (98) Lawless (98) and Orban (98) sharl criicized HL under he assumion ha he were reaing he rior for heir η as being indeenden of ha for he regression arameers. Hill noed ha if he rior for η is uniform he imlici rior for mus be roorional o so ha he riors for η and canno boh be indeenden of he regression arameers and argues ha he laer is he more naural assumion. However indeendence of he riors as insised uon b he HL criics is in fac nowhere required. The full join rior for he eercise ma eiher be eressed in erms of as in (5) or in erms of he HL η as

12 0 ( ) ( η ) ( ) ( ) η. (9) The assumion ha our ( ) HL assumion ha heir imlici ( ) is consan wih resec o is equivalen o he η is consan wih resec o η. I is immaerial wheher eiher of hese is roorional o or is inverse since an such consan will siml dro ou in he normalizing inegraion when Baes Rule is alied as in (7). 5 III. The Thomson Ice Core Inde The reinsaed Huner-Lambo CI mehod is illusraed b calibraing he millenial ice core isooe raio inde of Lonnie Thomson e al. (003) o recen global insrumenal emeraures. This inde was famousl described b Al Gore in his Nobel Peace Prize-winning An Inconvenien Truh (006: 60-65) as follows: Scienis Lonnie Thomson akes his eam o he os of glaciers all over he world. The dig core drills down ino he ice eracing long clinders filled wih ice ha was formed ear b ear over man cenuries. Lonnie and his eam of eers... can... measure he eac emeraure of he amoshere each ear b calculaing he raio of differen isooes of ogen (ogen-6 and ogen-8) which rovides an ingenious and highl accurae hermomeer.... The hermomeer o he righ measures emeraures in he Norhern Hemishere over he as 000 ears.... [T]he so-called global-warming skeics ofen sa ha global warming is reall an illusion reflecing naure s cclical flucuaions. To suor heir view he frequenl refer o he Medieval Warm Period. Bu as Dr. Thomson s hermomeer shows he vauned Medieval Warm Period (he hird lile red bli from he lef below) was in comared o 5 HL (98: 37) erroneousl acce ha he difference beween heir aroach and ha of Hoadle (970) is ha he assume a uniform rior on heir η whereas Hoadle considers a uniform rior for. In fac he difference is ha HL condiion heir rior for η or equivalenl on he rue arameers whereas Hoadle inaroriael condiions his rior for on he arameer esimaes as noed above.

13 he enormous increases in emeraure of he las half cenur (he red eaks a he far righ of he char). These global warming skeics... launched a fierce aack agains anoher measuremen of he 000-ear correlaion beween CO and emeraure known as he hocke sick a grahic image reresening he research of Michael Mann and his colleagues. Bu in fac scieniss have confirmed he same basic conclusions in mulile was wih Thomson s Ice core record as one of he mos definiive. (Gore 006: 60-64). McInre and McKirick ( ) had argued ha he hocke sick reconsrucion of Mann e al. (999) is based on invalid sribark brislecone reering daa as well as imroer use of Princial Comonens Analsis. As i haens he grah ha Gore (006) resens as Dr. Thomson s Thermomeer acuall was he disued hocke sick emeraure reconsrucion for which i was suosed o rovide definiive indeenden confirmaion sliced ogeher wih an insrumenal emeraure record as if he were a single series and had nohing o do wih Thomson s ice core daa. 6 Thomson who served as a member of he Science Advisor Board for An Inconvenien Truh had indeed ublished a similar-looking grah based on decadal averages of ice core ogen isooe raios in Figure 7 of Thomson e al. (003) which is reroduced in Figure below. 6 In fac wo of he series emloed b Mann e al. (999) for he crucial AD eriod were Quelccaa δ 8 O and reciiaion series so ha he hocke sick is no comleel indeenden of even he rue Thomson ice core inde.

14 Figure. Thomson e al. (003) Fig. 7 wih original caion. Panel d) of Fig. is he disued Mann e al. hocke sick reconsrucion iself overlain wih a Norhern Hemishere emeraure inde. Gore was suosed o have used anel c) bu misakenl used anel d) of he same figure insead merging is wo lines ino a single grah. 7 The slighl osiive values around 30 AD in anel d) 7 This subsiuion was confirmed b Lonnie Thomson in resonse o a quesion b he auhor a a seminar a Ohio Sae Universi Jan. 008.

15 3 became he hird lile red bli from he lef reresening he vauned Medieval Warm Period according o Gore. The acual Low Laiude comosie ice core inde in anel c) does urn u sharl in he 0 h cenur much like he insrumenall augmened hocke sick. However he aes a he o of anels a)-c) do no indicae emeraure bu merel comosie Z-scores comued from δ 8 O ice measuremens. Thomson e al. in fac did no calibrae hese Z- scores o emeraure le alone rovide CI s for such a calibraion so ha Dr. Thomson s Thermomeer was no e a hermomeer a all. The resen aer aems o fill his deficienc b calibraing his inde o insrumenal emeraure and comuing confidence inervals. Thomson e al. rovide decadall averaged daa for 5 of he 6 cores used in heir aricle in a sreadshee online a <brc.osu.edu/icecore/climaic-change-003-fig5- able.xls>. Unforunael however he Himalaan comosie Z-score (henceforh HCZ) series is missing along wih he daa for Quelccaa so ha even hough he las column is idenified as 6 core comosie i is in fac based onl on he wo included Andean cores and is definiel no he series loed in anel c) of heir Fig. 7. On March I ed Lonnie Thomson and several of his co-auhors requesing he missing series in his sreadshee bu received no rel. Neverheless i is ossible o a leas aroimael reconsruc he 6-core comosie from he five decadal δ 8 O ice series in he sreadshee lus archived annual values for he Quelccaa Summi core wih he hel of formulas embedded in he sreadshee.

16 4 Figure shows he decadall averaged δ 8 O ice values from he sreadshee for he hree 3 Himalaan cores Gulia Dunde and Dasuou. 8 Each decadal value is idenified b he dae a he o of is inerval. These daes are even mulilies of 0 so ha he decades used for hese hree series are ec. Noe ha onl Dasuou has daa for he final decade of he sud (hrough 997). Gulia ends wih and Dunde wih All hree of hese Himalaan series visuall mach hose in Figure 5 of Thomson e al. (003). Figure. Decadal averages of δ 8 O ice for he hree Himalaan cores: Gulia Dunde and Dasuou. 8 The decadall averaged daa from his sreadshee for he 3 Himalaan cores has been archived a f://f.ncdc.noaa.gov/ub/daa/aleo/icecore/ro/dasuou/dasuou-d8o. <... /dunde/dunde-d8o.> and <... /gulia/gulia-d8o.>. However he decadal Andean daa from he sread shee is no archived here ece inadverenl a he lower righ side of he Gulia file. IPCC AR4 WGI eer reviewer Sehen McInre ( ) noes ha several inconsisen versions of Gulia and Dunde have been used in he climae lieraure. The resen aer siml uses he version in he sreadshee consruced for Thomson e al. (003).

17 5 Onl wo of he hree Andean cores used Sajama and Huascarán Core are abulaed in he online sreadshee. Each of hese decadal values is associaed wih wo daes idenified as boom and o decade bu boh are even muliles of 0 saring wih (000 00) (00 00) ec. Evidenl hese Andean cores were averaged over decades beginning wih an even mulile of 0 such as ec. so ha he header o decade refers o he firs ear of he following decade. This inerreaion is reinforced b an annoaion ha he mos recen value for Huascarán Core covers he eriod 990 o 993. The hird Andean core is idenified b Thomson e al. (003) onl as Quelccaa even hough annual δ 8 O ice daa is archived a f.ncdc.noaa.gov for wo differen Quelccaa cores Core and he Summi Core. I was found ha averages of he Summi Core for decades beginning wih even muliles of 0 gave a erfec visual fi o he grah of he Quelccaa series in Fig. 6 of Thomson e al. (003) while decades ending wih even muliles of 0 had obvious differences. Core had obvious differences under eiher inerreaion and clearl was no he source of heir Fig. 6. The average of he wo cores also has obvious differences. The archived annual daa for he Quelccaa Summi Core (Thomson e al. 005) onl eend o 984 e Thomson e al. (003) use daa for Quelccaa hrough 000 based on newer shallow cores (see caion o heir Fig. 7 in Fig. above) ha are no available numericall. Accordingl hese wo values were read visuall off Fig. 6 of Thomson e al. (003) as er mil for and er mil for

18 6 Decadal averages of δ 8 O ice for he hree Andean cores Quelccaa Summi Sajama and Huascarán are shown in Figure 3 below. The Sajama series is a good visual mach o ha in Figure 6 of Thomson e al. (003). The relaive sizes of he local eaks in he Huascarán series since 790 are no quie he same as in Figure 6 of Thomson e al. (003) alhough he general shae of he wo series are ver similar. As noed above Quelccaa is a erfec mach bu onl rovided he Summi Core is used wih decades beginning wih 0. Figure 3 Decadal averages of δ 8 O ice for he hree Andean cores Quelccaa Summi Sajama and Huascarán. The ice-core indices shown in Fig. are comosie Z-scores are suosed o be derived somehow from he decadal δ 8 O ice daa shown in Figs. and 3 above bu

19 7 Thomson e al. (003) do no rovide deails of how he were comued. Were he δ 8 O ice values averaged and hen Z-scores comued or were he averages comued from Z-scores? Were he means and variances for he Z-scores comued before or afer decadal averaging? Were he comued using each enire series or jus he orion afer 000 AD used for he sud? And were he means and variances comued o he end of each series or onl using he common orions? Alhough he rovided daa sreadshee is incomlee i does conain embedded formulas ha reveal how he Z-scores would have been comued if he daa were resen: Firs each decadall averaged series was convered o Z-scores using is mean and sandard deviaion for all is available decades bu resriced o he eriod of he sud i.e. beginning wih for he Himalaan cores and wih for he Andean cores. These Z-scores are shown in Figures 4 and 5 below:

20 8 Fig. 4 Z-scores for he hree Himalaan cores.

21 9 Fig. 5 Z-scores for he hree Andean cores. Then regional comosie Z-scores were comued b averaging each region s Z- scores as available. The resuling Andean Comosie Z-score series (ACZ) shown in Fig 6 below is hus siml he average of he hree comonen Z-scores for all 00 decades. Unforunael however his series is no a erfec visual mach o Fig. a of Thomson e al. (003) (Fig. a above). Alhough mos of he high frequenc wiggles mach he local high in he 70s and he record low in he 80s are greal aenuaed in he Thomson e al. version. Furhermore alhough he 0h cenur is generall he highes in he Thomson version i is no quie as high as he h-4h cenuries in our emulaion. Perhas Thomson e al. in fac used Quelccaa Core or he average of he

22 0 wo Quelccaa cores in he consrucion of he inde even hough he Summi core was clearl wha was used in heir Fig. 6. Fig. 6 Andean Comosie Z-score series ACZ. The emulaed Himalaan Comosie Z-score (HCZ) series shown in Fig. 7 below is based on all 3 Himalaan cores for he firs 99 decades (shown in blue) and hen abrul is based onl on Dasuou for he final decade (shown in red) afer he oher wo dro ou. The line segmen connecing he oin reresening he 980s (3 cores) and ha reresening he 990s (Dasuou onl) is also in red. The resuling series is a erfec visual mach o Panel (b) of Fig..

23 Fig. 7. Himalaan Comosie Z-score series (HCZ) using all hree Himalaan cores (in blue) and Dasuou onl (in red). The line segmen connecing (3 cores) and (Dasuou onl) is also in red. Finall he sreadshee formulas show ha wo regional comosies were averaged wih equal weighs o obain a single Low Laiude Comosie Z-Score (LCZ) as shown in Figure 8 below. The final decade as well as he line segmen connecing he 980s (6 cores) o he final decade (4 cores) are loed in red. This series is a ver good visual mach o he 6-core comosie series in Panel (c) of Figure above including in aricular he unrecedenedl high value for he final decade of he 990s. Our relicaion of he series in anel c) in Figure 7 of Thomson e al. (003) is herefore reasonabl successful desie is absence from he daa sreadshee and desie some inconsisencies in he Andean sub-inde.

24 Fig. 8. Low-Laiude Comosie Z-Scores based on 6 cores for (in blue) and on 4 cores for (in red). The line segmen connecing (6 cores) and (4 cores) is also in red. I is no clear wh one would wan o average Z-scores in his manner since i makes inefficien use of he daa. However his aears o be a common racice in aleoclimaolog (see e.g. Kaufman e al. 009). The rimar goal of he resen aer is merel o relicae and calibrae he Thomson e al. (003) comosie ice core series and no o imrove uon i. Secion V below suggess how his daa migh be used more efficienl. I is singular ha alhough he LCZ series in anel (c) of Figure ends on an unrecedenedl high value for he 990s and anels (a) and (b) of Figure above shows

25 3 ha his is being driven b is HCZ comonen raher han he ACZ comonen none of he individual Himalaan Z-scores shown in Fig. 4 ends on a record high value. This onl comes abou as a saisical arifac however because wo of he Himalaan series dro ou in he las decade leaving onl Dasuou which was running higher or warmer han he oher series o reresen he region b iself. The effec is accenuaed because when HCZ and ACZ are averaged ogeher o obain LCZ he ne weigh on Dasuou suddenl changes from /6 in he firs 99 decades o / in he final decade. Fig. 9 Himalaan Comosie Z-score (HCZ) comued from all 3 Himalaan cores (in blue) and from Dasuou onl (in red). Figs. 9 and 0 show HCZ and LCZ for he full eriod AD comued boh from all 6 cores (LCZ6) as well as from onl he 4 cores ha make i o he final

26 4 decade of he sud eriod (LCZ4). Clearl he 990s are no an warmer han he 940s when eiher HCZ or LCZ is comued on a consisen basis desie he imression given in anels b) and c) of Fig.. Fig. 0 Low-Laiude Comosie Z-score (LCZ) comued from all 6 cores (in blue) and from onl he 4 cores available in he las decade (in red). The changing comosiion of LCZ would no aler is eeced resonse o emeraure if δ 8 O ice were a universall valid and linear indicaor of annual average emeraure aside from a locaion shif o comensae for differences in laiude and/or aliude. However some of he cores such as Quelccaa Summi and Dasuou are srongl correlaed wih annual average global emeraures while ohers such as Sajama and Dunde are no so ha he resonse is no universal. Some of he sies ma

27 5 be reflecing hisorical variaions in he seasonali or aliude of snow formaion or in reciiaion aerns as he air is ransored o he sie raher han average annual emeraures so ha he correlaion if an wih local or global annual emeraure can onl be deermined emiricall on a sie-b-sie basis. The makeu of he comosie series herefore can make a subsanial difference for is correlaion wih emeraure so ha in fac he final decade of LCZ (based on 4 cores wih / weigh on Dasuou) is a quie differen series han he receding 99 decades (based on 6 cores wih equal weighs on all) and canno be used in is calibraion. IV. Calibraion o Global Temeraure Thomson e al. (003) sae ha Comarison of his [LCZ] ice core comosie wih he Norhern Hemishere ro record ( A.D.) reconsruced b Mann e al. (999) and measured emeraures ( ) reored b Jones e al. (999) suggess he ice cores have caured he decadal scale variabili in he global emeraure rends. Of course i would be circular o calibrae LCZ o he hocke sick and hen o u i forward as indeenden confirmaion of he hocke sick so onl he comarison o insrumenal emeraures is relevan for our uroses. Alhough Thomson e al. (003) comare heir LCZ favorabl o a Norhern Hemishere insrumenal emeraure inde a global emeraure inde would more aroriae for univariae calibraion given ha half of he si cores are from he Souhern Hemishere. Figure below shows annual averages for of he CRUTEM3vGL global land air inde (variance adjused) of Brohan e al. (006) a

28 6 global sequel o he Jones e al. (999) series o which Thomson e al. comare heir series. 9 Fig.. CRUTEM3vGL global land air emeraure inde annual averages 850-Oc The widel used CRU series are ariall based on confidenial weaher daa ha has no been made ublicl available (see Me Office 009). Since he are no scienificall relicable he should onl be used wih cauion. However he alernaive GISSem series roduced b NASA/GISS is onl comued back o 880 and hence has hree fewer decades han he CRU series. Preliminar calculaions show ha he 9 Source <h:// The annual average loed for 009 is incomlee running onl hrough Ocober. However he calibraion onl uses he orion hrough 000.

29 7 correlaions wih GISSem are much weaker han wih CRUTEM3vGL bu his is more likel due o he smaller samle size han o an difference in behavior since 880. Since half of he Thomson e al. decadal ice core series are based on decades ending in a ear divisible b 0 while he oher half use decades beginning in a ear divisible b 0 i is a maer of indifference which convenion is used for he emeraure series. This aer arbiraril ados he former convenion i.e ec. and reas LCZ as if i were comued consisenl on his basis. Since he 4-core LCZ series LCZ4 does no necessaril have he same relaionshi o emeraure as he 6-core series LCZ6 he wo mus be calibraed searael o emeraure. Table below shows he resuls of regressing he 4 decadal values of LCZ6 for on he corresonding decadal averages of he emeraure series. For uroses of comarison he 4-core inde LCZ4 using onl Dasuou for he Himalaan region wih a / weigh as in was comued for all 5 os-850 decades and also regressed on he emeraure series as shown in he las wo columns of he same able. Table Decadal Regression of 6- and 4- core LCZ on CRUTEM3vGL and res. LCZ6 LCZ4 ˆ ˆ ˆ ˆ coef s.e sa (-sa) n 4 5 R s s s

30 8 r DW I ma be seen ha here is a osiive and weakl significan correlaion beween LCZ6 and insrumenal emeraure (-ailed 0.079) and a osiive and highl significan correlaion beween LCZ4 and emeraure (-ailed 0.008). The firs order serial correlaion coefficiens of he residuals r are esseniall 0 and he Durbin- Wason saisics DW are ver close o so ha osiive serial correlaion is no an issue wih hese regressions a leas no a his decadal frequenc. Fig. CRUTEM3vGL global emeraure inde wih reconsrucions using he 6-core LCZ and he 4-core LCZ along wih oserior median esimaes.

31 9 Fig. shows he decadal averages of CRUTEM3vGL used for he calibraion along wih he 6- and 4-core reconsruced emeraures using he classical esimaor (4). Also shown are he medians of he oserior disribuion for he 6- and 4- core reconsrucion. In mos cases he classical reconsrucion and oserior median are almos indisinguishable. Since he classical esimaor (4) is much easier o comue and elain i is naural o focus on i as he oin esimae of reconsruced emeraure. I is conjecured ha he classical esimaor coincides wih he oserior mode. Even hough LCZ4 eceeded all revious values of LCZ6 during he decade of he 990s he emeraures reconsruced from i are somewha lower han hose of he 940s using LCZ6. This comes abou because he esimaes of boh he inerce and sloe are higher for LCZ4 han for LCZ6 in Table. Boh hese facors work o reduce he reconsruced emeraure as a funcion of he inde. Figure 3 shows he emeraure reconsrucion for using he classical oin esimae (4) raher han he oserior median. The reconsrucion is based on LCZ6 for (in blue) and on LCZ4 for (in red). The line segmen connecing (6 cores) and (4 cores) is also in red. Alhough he final decade has a decidedl warm oin esimae i is no as warm as he 940s conrar o he imression given in Figs. and 8.

32 30 Fig. 3. Reconsruced emeraure oin esimaes using (4) wih LCZ6 for (in blue) and LCZ4 for (in red). The line segmen connecing (6 cores) and (4 cores) is also in red. Of course since we have relaivel recise insrumenal emeraures for 850 o he resen he reconsrucion in Fig. 3 is no our bes esimae of global emeraures for hese ears. The full reconsrucion is rovided here and wih is confidence inervals below merel o show how he full Dr. Thomson s hermomeer calibraes o emeraure. The oserior cumulaive disribuion funcions (CDFs) for he LCZ6 reconsrucion for he 99 decades o were comued in ses of.0

33 3 from o 5.00 C as shown in Fig. 4 below. 0 Because he sloe coefficien is onl weakl significan for his reconsrucion (.079) he ails of he oserior disribuions are eremel heav and slow o converge o 0 and for low and high emeraures. Fig. 4. Poserior CDFs for 99 LCZ6 emeraure reconsrucions o Fig. 5 below shows he oserior CDF for he LCZ4 reconsrucion comued for he sake of comarison for he 5 calibraion decades o Since he sloe coefficien in his regression is highl significan ( 0.008) he disribuion is 0 CRU (<h:// gives he recision of CRUTEM3vGL ( sandard errors) as aroimael 0.05 C since 95 and abou 0. C in 85 wih gradual imrovemen from 860 o 950 ece for warimes so ha ses of 0.0 C are overkill. Neverheless using Malab s mvcdf as described in Aendi I he rogram for all 00 reconsrucion daes runs in jus a few minues on a desko PC.

34 3 nearl Suden wih 3 degrees of freedom and he ails converge much more quickl o 0 and for low and high emeraures han wih LCZ6. Fig. 5. Poserior CDF for 5 LCZ4 reconsrucions o The oserior CDFs in Figs. 4 and 5 were invered b linear inerolaion a robabiliies and The resuling 50% and 95% LCZ6 Confidence Inervals (CIs) are loed in blue for in Fig. 6 below along wih he classical oin esimaes from Fig. 3. Corresonding values for LCZ4 are loed in red for The line segmens connecing (6 cores) and (4 cores) are also shown in red.

35 33 Fig. 6. LCZ6 emeraure reconsrucion for wih 50% and 95% Confidence Inervals in blue. LCZ4 emeraure reconsrucion for wih 50% and 95% Confidence Inervals in red. Line segmens connecing (6 cores) and (4 cores) are also in red. The arial decade was he warmes in he insrumenal record used a 0.64 C relaive o I ma be seen from Fig. 6 ha emeraures hroughou he eriod could have been a leas. C warmer han or a leas.7 C colder. The Medieval Warm Period (MWP) or even he Lile Ice Age (LIA) for ha maer herefore could well have been even warmer han he mos recen decade so far as his reconsrucion goes. (Since was no used in he calibraion i is no shown in Fig..)

36 34 An even beer Thomson δ 8 O ice series o calibrae should have been he 7-core inde of Thomson e al. (006). This inde adds a sevenh Himalaan core Puruogangri and eends back 000 ears where hree of he cores sill have daa. B all righs i should suersede he 6-core 000 ear inde of Thomson e al. (003). Unforunael however he suoring daa ses ha accoman ha aer onl abulae he seven cores back 400 ears as five-ear averages so ha he differen subses of he seven cores ha are acive in differen eriods canno be searael calibraed as he mus be. Furhermore he decadall averaged Z-score indices in Daa Se 3 do no mach he averages of he illusraive five-ear average Z-score indices for eiher region ha are abulaed in Daa Se so ha here is no linear relaionshi beween he archived 000- ear comosie Z-score inde and he illusraive 400-ear δ 8 O ice daa on which i is suosed o be based as documened b McCulloch (009a 009b). Three s o Lonnie Thomson and mos of his co-auhors asking for he comlee daa and clarificaion of his inconsisenc received no rel. The resen sud herefore focuses insead on he 000-ear 6-core δ 8 O ice inde of Thomson e al. (003) which is a leas aroimael relicable in erms of is full abulaed decadal comonen series. V. Furher Calibraion Issues This concluding secion discusses a number of furher calibraion issues: i) he classical join samling confidence inervals discussed b Hoadle (970) and Brown (993) ii) efficien muliro calibraion iii) a more owerful sequenial rior aroach; and iv) efficien use of rior informaion abou he sign of he calibraion sloe

37 35 coefficien. No aem is made a resen o imlemen he roosed soluions o issues ii) iv). i) The radiional or classical alernaive aroach o calibraion CI consrucion is based on Fieller s (954) confidence region for he raio of wo Suden random variables. This aroach firs considers a 00(-γ) CI for as a funcion of using he sandard Suden disribuion (6) for ( ' ' ˆ ˆ S s ) b a air of herbolic funcions of : ( ) / ˆ c ˆ v( ) / ˆ ˆ c v. These CIs bound where c ( γ / ) is he Suden -ailed criical value for es size γ Tν() is he T n sandard Suden CDF wih ν degrees of freedom and v ( ) is as defined in (7). Le ( T ( ˆ ( s ) n / s be he wo-ailed -value of he es saisic for he hohesis ha ˆ 0. If < γ so ha he sloe is significanl differen from 0 a level γ solving he boundaries for as a funcion of he observed ro value ields a quadraic equaion wih wo real roos ξ < ξ where inersecs he wo herbolas and he confidence region is he bounded inerval (ξ ξ ). Bu if γ his inerval becomes eiher (ξ ) (- ξ ) or (- ) deending on he signs of ˆ and where is he mean of he calibraion i -values. And if > γ so ha he sloe is insignificanl differen from 0 a level γ eiher wice inersecs one of he herbolas and he confidence region is he disconiguous se ( ξ ) (ξ ) or else here are no real roos in which case he confidence region is he enire real line. Hoadle (970) drl remarks ha his se ossesses inheren difficulies.

38 36 Brown (98) rooses ha (ξ ξ ) is a resecable inerval rovided he -es of he hohesis 0 is rejeced. Indeed if sa ˆ > 0 and i were known wih erfec cerain ha has he same sign hen for an < ξ he robabili ha could have been as high as is observed value would be less han or equal o γ/ while for an > ξ he robabili ha could have been as low as is observed value would also be less han or equal o γ/. The hich however is ha here in fac is alwas robabili / ha and ˆ have oosie signs. The Baesian aroach on he oher hand akes full accoun of he ossibili ha and ˆ migh have oosie signs and never gives unbounded CIs le alone disconiguous confidence regions. When > γ as in he LCZ6 eamle above he Baesian CI is siml somewha wider han i oherwise would be as is onl naural. Fig. 7 comares he oserior and Fieller 50% CIs for our LCZ calibraion. Since for LCZ6.079 <<.50 he Fieller 50% region is an inerval and he wo are ver similar. For he final decade based on LCZ so ha he similari is even closer. In general he classical CI is a good aroimaion o he Baesian CI when << γ so ha here is onl negligible chance relaive o γ ha and ˆ can have oosie signs. Indeed Huner and Lambo (98) hemselves acuall advocae i as an aroimaion.

39 37 Fig. 7. Comarison of oserior (lines) and Fieller (shaded) 50% confidence inervals. However when he sloe is onl marginall significan or even insignifican relaive o he desired γ as is ofen he case in aleoclimae cones he eac Baesian inerval is o be referred. Fig. 8 comares he oserior and Fieller 95% confidence regions for our LCZ calibraion. Since for LCZ6 we now have.079 > γ.05 he Fieller region is unreasonabl eiher he enire real line or else a air of semiinfinie inervals. For he final decade however.0008 <<.05 so ha he Fieller 95% region becomes an inerval and is a good aroimaion o he oserior inerval.

40 38 Fig. 8. Comarison of oserior (lines) and Fieller (shaded) 95% confidence regions. ii) When several roies are available and here is an a riori eecaion ha he will have he same quaniaive relaionshi o he sae variable ne of a locaion shif heir average ma siml be univariael calibraed o he sae variable. However when as in he case of he Thomson e al. (003) ice core daa he individual roies resond unequall o he sae variable and someimes no a all a muliro aroach should be more efficien han arbiraril averaging eiher he raw roies hemselves or heir Z-scores as in Thomson e al. ( ) and Kaufman e al. (009). Of course one ma no siml disregard he insignifican roies or hose ha give he wrong

41 39 sign sloe for hen here will be selecion bias ha will invalidae convenional significance ess. Suose here are q roies wih ij being he i-h calibraion-eriod observaion on he j-h ro and j being he reconsrucion-eriod observaion on he j-h ro. Assume firs ha ij ε (0) j j i ij where ε ij ~ N (0 i ) and are indeenden across roies. Le ( ) T j j K nj and T ( n K ) so ha ˆ j ˆ j and s are imlici in j and and ( ) j ˆ j ma be comued under a diffuse rior as in Secion of he e. Then under a diffuse rior for q ( ) ( ). () q q j j If he regression errors are no indeenden across roies as reliminar calculaions show o be he case wih he Thomson e al. ice core daa hen () is invalid. However if we have q n- roies we ma obain errors ha are uncorrelaed across roies and herefore indeenden under our Gaussian assumion b insead regressing each ro j on he sae variable lus roies... j- during he calibraion eriod i... n as follows: j j j ij j j j i h γ jh ih ε. ij The calibraion OLS esimaes of hese arameers are imlici in he arial calibraion daa... j and. Using he reconsrucion values of he roies define he orhogonal innovaion of ro j and is esimae b

42 40 j j j h γ jh h for j and j h ˆ ˆ γ j j jh h ˆ. Then ( ˆ... ) ma be comued as in () j j siml b adding he aroriae erms o (6) (0). Equaion () above hen generalizes o q ( ) (... q... q ˆ j... j ). () j Equaion () or () ma be evaluaed numericall eiher b calculaing he RT densiies direcl as in Huner and Lambo (98) or b finie differences from he CDF and hen numericall inegraing he roduc of oserior densiies o normalize and obain CIs. A ransformaion such as Suden wih q- degrees of freedom ma be useful in order o caure he enire ails. When he roies are significanl correlaed wih he sae variable even afer condiioning on he receding roies and if he ell a consisen sor abou he reconsrucion sae variable he muliro oserior disribuion () will end o be much igher abou is mode han he single-ro reconsrucions would be. Bu when he roies are significan e ell an inconsisen sor abou he muliro oserior densi will end o be more sread ou or even mulimodal. CIs comued from his disribuion will hen end o be wide enough o include all or a leas mos of he individual modes. As long as all he roies have he same coverage across observaions he order in which he are arranged should make no difference for he muliro oserior

43 4 disribuion. However if he have differen coverage as is he case for he calibraion eriod wih wo of he Thomson e al. (003) roies and as is he case during he 000-ear reconsrucion eriod for several of he roies in Thomson e al. (006) i is imoran o arrange hem in decreasing order of coverage so ha all he required condiioning roies are available a each se. If a ro has more coverage in he reconsrucion eriod han he ohers bu less in he calibraion eriod or vice versa i ma unforunael be necessar o give u some daa oins. In he muliro case he oin esimae of ma be aken as he oserior median since if he roies are inconsisen here ma be mulile oserior modes. The GLS esimaor suggesed b Brown (98: eq..6) which generalizes he single-ro classical esimaor in equaion (4) is far simler o comue. However i ma be roblemaic because i is esseniall a weighed average of single ro esimaes each of which has Cauch-like heav ails wih undefined mean. The Law of Large Numbers herefore does no ensure ha he disribuion of he average will be more comac han he disribuions of he conribuing values. When he roies do no have indeenden simle regression errors and q n- as is for eamle he case wih he Kaufman e al. (009) daa () becomes invalid since evenuall zero and even negaive degrees of freedom will be encounered. In his case i ma be aroriae o reduce he dimensionali of he ro daa se b careful use of Princial Comonens Analsis (Preisendorfer 988 Mann e al. 999 McInre and McKirick 005). Such an aroach ma also be referred even when q < n- if n/q is no large.

44 4 iii). Secion II above reas he roblem of inferring a single reconsrucion value for an observaion ha has no sequenial relaionshi o he calibraion daa. In a ime series or oher sequenial cone however here ma be valuable informaion from wha is known or has been inferred abou he adjacen observaions. In aricular le us suose ha he calibraion sae observaions now wih a subscri... n o indicae ime follow a random walk so ha η (3) where η ~ iid N(0 τ ) and are indeenden of he regression errors. The random walk signal variance ma be esimaed from he calibraion daa b n ˆ τ ( ) /( n ) ~ τ χ n /( n ). Suose we have a ime series of calibraion ro values... n as well as an adjoining ime series of reconsrucion ro values... T wih T ec. as is ordinaril he case in aleoclimae emeraure reconsrucions. We wish o infer he T sae values... T which also follow he random walk (3) wih T ec. Le and reresen he corresonding vecors of reconsrucion values. The random walk imlies ha he disribuion of he sae variable k eriods before or k eriods afer n condiional on he calibraion daa has a normal disribuion wih mean or n and variance kτ and ha he uncondiional disribuion has infinie variance. I herefore moivaes he imroer diffuse rior assumion of Secion II for an isolaed reconsrucion dae ha is so far from he calibraion daa ha

45 43 he calibraion values become uninformaive. However when he reconsrucion dae is eiher close o he calibraion daa or surrounded b oher reconsrucion daes condiioning on his informaion can greal imrove he recision of he reconsrucion. Condiional on he rue arameer values he reconsrucion ro values and he calibraion sae vecor he oserior disribuions of he are Gaussian wih means and variances deermined b a simle Kalman smooher (e.g. Harve 989 McCulloch 005). This meshes he informaion in a forward Kalman filer ha sars a wih a diffuse rior as in Secion II wih he informaion in a reverse Kalman filer ha sars a T wih a rior ha is governed b and works backwards. Then he rue arameers ma be inegraed ou o find he disribuions of he condiional on he arameer esimaes. To simlif he noaion he following develomen of he Kalman filer and smooher is imlicil condiional on he rue arameer values and τ unil oherwise indicaed. The forward Kalman filer is iniialized wih a diffuse rior as in secion II so ha as in (7) ( ) ~ N ( ) ( / / ) N( c d v ) wih c / d / v. / Now suose ha for an... T he mean of filer densi for - is likewise affine in : (... ) ~ N( c d v ).

46 44 Then (3) imlies ha he redicive densi for ime is ( ) ( ) ~... τ v d c N. Baes Rule hen imlies ha he filer densi for iself is also affine in : ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ( ) ~ / / ~ v d c N v d c N N τ ) wih ( ) ( ) ( ) ( ) ( ) ( ) ( ). / / / / d g g d c g g c v g v v τ The reverse Kalman filer works b backward inducion beginning wih ime T using a rior based on he observed value of T and (3): ( ) ( ) ( ) ~ τ N T T Baes Rule hen imlies ha he reverse filer for T is again affine in : ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ). / / / / / ; ~ / / ~ τ τ ( ) T T T T T T T T T T T T T T T T T T T T T g d g g c v g v v d c N N N Reasoning analogous o ha for he forward filer hen imlies ha he reverse filer is likewise affine in for all :

47 45 ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ). / / / / ; ~... T T d g g d c g g c v g v v v d c N τ Finall merging ogeher he indeenden informaion in ( )... and ( ) ( ) ~... τ T v d c N ields he Kalman smooher ( ) ( ) ( ) ( ) ( ) ( ). / / / ; ~ d G G d D c G G c C v V G v v V V D C N τ τ (4) Smooher oin esimaes ˆ analogous o he classical oinwise esimaor (4) ma be found b evaluaing he mean D C using he arameer esimaes in lace of he arameers hemselves. In order obain valid confidence inervals i remains o inegrae ou he unknown arameers. Unforunael his is no longer deermined b he RN disribuion even condiional on he wo variances since now eners in a comlicaed wa ino C D and V. However he deendence on he inerce is sill simle since i onl aears in he affine eression for he mean. Equaion (9) imlies ( ) ( )( ) ( ) ( ) ˆ / ˆ N ~ ˆ ˆ ρ τ s s s S where. Therefore ( ) / / ρ s s s Noe ha he informaion se { } is equivalen o { } ˆ ˆ since once we know boh and ˆ and add no informaion. Therefore he informaion se ˆ { } ma be aken o imlicil include { } ˆ ˆ.

48 46 ( ˆ ˆ S τ ) ~ N ( ˆ ( / )( ˆ C D s s ) V D ( ρ ) ( s ). (5) The remaining hree arameers and τ mus be inegraed ou numericall o obain ( ˆ ˆ S s ˆ τ ) using ( τ ˆ s s ˆ τ ) ( ˆ s ) ( s ) ( τ ˆ τ ) and a uniform rior on log(τ ). These inegrals ma be erformed using ernar inegraion as described in Aendi II. Alhough his is a rile inegral i is robabl adequae o use ernar inegraion wih onl a modes number of oins over each arameer e.g. m 0 for and m 4 for he wo variances making a oal of 60 combinaions o evaluae. (In fac he sequenial srucure of he sae variable means ha he ro daa conain some indirec informaion abou and τ ha is no conained direcl in s or τˆ. However he gain from eloiing his indirec informaion is likel o be small.) The resuling sequenial reconsrucion will b consrucion equal he observed sae daa during he calibraion eriod. In order o make his rivial relaionshi clear he calibraion orion should be loed wih a differen color or smbol han he nonrivial orion. A he end of he reconsrucion eriod he reconsrucion will be cenered nearl on wih a ver igh disribuion. Moving back ino he reconsrucion eriod however he oin esimaes will begin o look like a smoohed version of he oinwise reconsrucion wih smooh CIs ha graduall widen. However because he informaion in adjacen ro values is being aken ino accoun he CIs will never become as wide as hose for he oinwise reconsrucion. Siml smoohing he oinwise CIs is herefore no equivalen o comuing a CI for he smoohed

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at

BOX-JENKINS MODEL NOTATION. The Box-Jenkins ARMA(p,q) model is denoted by the equation. pwhile the moving average (MA) part of the model is θ1at BOX-JENKINS MODEL NOAION he Box-Jenkins ARMA(,q) model is denoed b he equaion + + L+ + a θ a L θ a 0 q q. () he auoregressive (AR) ar of he model is + L+ while he moving average (MA) ar of he model is

More information

Chickens vs. Eggs: Replicating Thurman and Fisher (1988) by Arianto A. Patunru Department of Economics, University of Indonesia 2004

Chickens vs. Eggs: Replicating Thurman and Fisher (1988) by Arianto A. Patunru Department of Economics, University of Indonesia 2004 Chicens vs. Eggs: Relicaing Thurman and Fisher (988) by Ariano A. Paunru Dearmen of Economics, Universiy of Indonesia 2004. Inroducion This exercise lays ou he rocedure for esing Granger Causaliy as discussed

More information

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag.

- The whole joint distribution is independent of the date at which it is measured and depends only on the lag. Saionary Processes Sricly saionary - The whole join disribuion is indeenden of he dae a which i is measured and deends only on he lag. - E y ) is a finie consan. ( - V y ) is a finie consan. ( ( y, y s

More information

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2

Check in: 1 If m = 2(x + 1) and n = find y when. b y = 2m n 2 7 Parameric equaions This chaer will show ou how o skech curves using heir arameric equaions conver arameric equaions o Caresian equaions find oins of inersecion of curves and lines using arameric equaions

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

Vector autoregression VAR. Case 1

Vector autoregression VAR. Case 1 Vecor auoregression VAR So far we have focused mosl on models where deends onl on as. More generall we migh wan o consider oin models ha involve more han one variable. There are wo reasons: Firs, we migh

More information

Computer Vision. Motion Extraction

Computer Vision. Motion Extraction Comuer Moion Eracion Comuer Alicaions of moion eracion Change / sho cu deecion Surveillance / raffic monioring Moion caure / gesure analsis HC image sabilisaion Moion comensaion e.g. medical roboics Feaure

More information

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles

Diebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance

More information

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé

Bias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070

More information

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES

FORECASTS GENERATING FOR ARCH-GARCH PROCESSES USING THE MATLAB PROCEDURES FORECASS GENERAING FOR ARCH-GARCH PROCESSES USING HE MALAB PROCEDURES Dušan Marček, Insiue of Comuer Science, Faculy of Philosohy and Science, he Silesian Universiy Oava he Faculy of Managemen Science

More information

Lecture 6 - Testing Restrictions on the Disturbance Process (References Sections 2.7 and 2.10, Hayashi)

Lecture 6 - Testing Restrictions on the Disturbance Process (References Sections 2.7 and 2.10, Hayashi) Lecure 6 - esing Resricions on he Disurbance Process (References Secions 2.7 an 2.0, Hayashi) We have eveloe sufficien coniions for he consisency an asymoic normaliy of he OLS esimaor ha allow for coniionally

More information

Description of the MS-Regress R package (Rmetrics)

Description of the MS-Regress R package (Rmetrics) Descriion of he MS-Regress R ackage (Rmerics) Auhor: Marcelo Perlin PhD Suden / ICMA Reading Universiy Email: marceloerlin@gmail.com / m.erlin@icmacenre.ac.uk The urose of his documen is o show he general

More information

Vector autoregression VAR

Vector autoregression VAR Vecor auoregression VAR So far we have focused mosly on models where y deends on as y. More generally we migh wan o consider models for more han on variable. If we only care abou forecasing one series

More information

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H.

ACE 564 Spring Lecture 7. Extensions of The Multiple Regression Model: Dummy Independent Variables. by Professor Scott H. ACE 564 Spring 2006 Lecure 7 Exensions of The Muliple Regression Model: Dumm Independen Variables b Professor Sco H. Irwin Readings: Griffihs, Hill and Judge. "Dumm Variables and Varing Coefficien Models

More information

20. Applications of the Genetic-Drift Model

20. Applications of the Genetic-Drift Model 0. Applicaions of he Geneic-Drif Model 1) Deermining he probabiliy of forming any paricular combinaion of genoypes in he nex generaion: Example: If he parenal allele frequencies are p 0 = 0.35 and q 0

More information

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

KEY. Math 334 Midterm I Fall 2008 sections 001 and 003 Instructor: Scott Glasgow 1 KEY Mah 4 Miderm I Fall 8 secions 1 and Insrucor: Sco Glasgow Please do NOT wrie on his eam. No credi will be given for such work. Raher wrie in a blue book, or on our own paper, preferabl engineering

More information

Chapter 3 Common Families of Distributions

Chapter 3 Common Families of Distributions Chaer 3 Common Families of Disribuions Secion 31 - Inroducion Purose of his Chaer: Caalog many of common saisical disribuions (families of disribuions ha are indeed by one or more arameers) Wha we should

More information

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8) I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression

More information

Unit Root Time Series. Univariate random walk

Unit Root Time Series. Univariate random walk Uni Roo ime Series Univariae random walk Consider he regression y y where ~ iid N 0, he leas squares esimae of is: ˆ yy y y yy Now wha if = If y y hen le y 0 =0 so ha y j j If ~ iid N 0, hen y ~ N 0, he

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation

The fundamental mass balance equation is ( 1 ) where: I = inputs P = production O = outputs L = losses A = accumulation Hea (iffusion) Equaion erivaion of iffusion Equaion The fundamenal mass balance equaion is I P O L A ( 1 ) where: I inpus P producion O oupus L losses A accumulaion Assume ha no chemical is produced or

More information

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT

ON DETERMINATION OF SOME CHARACTERISTICS OF SEMI-MARKOV PROCESS FOR DIFFERENT DISTRIBUTIONS OF TRANSIENT PROBABILITIES ABSTRACT Zajac, Budny ON DETERMINATION O SOME CHARACTERISTICS O SEMI MARKOV PROCESS OR DIERENT DISTRIBUTIONS O R&RATA # 2(3 ar 2 (Vol. 2 29, June ON DETERMINATION O SOME CHARACTERISTICS O SEMI-MARKOV PROCESS OR

More information

ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING

ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING Inernaional Journal of Social Science and Economic Research Volume:02 Issue:0 ESTIMATION OF DYNAMIC PANEL DATA MODELS WHEN REGRESSION COEFFICIENTS AND INDIVIDUAL EFFECTS ARE TIME-VARYING Chung-ki Min Professor

More information

Lecture Notes 2. The Hilbert Space Approach to Time Series

Lecture Notes 2. The Hilbert Space Approach to Time Series Time Series Seven N. Durlauf Universiy of Wisconsin. Basic ideas Lecure Noes. The Hilber Space Approach o Time Series The Hilber space framework provides a very powerful language for discussing he relaionship

More information

Testing for a Single Factor Model in the Multivariate State Space Framework

Testing for a Single Factor Model in the Multivariate State Space Framework esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics

More information

Distribution of Estimates

Distribution of Estimates Disribuion of Esimaes From Economerics (40) Linear Regression Model Assume (y,x ) is iid and E(x e )0 Esimaion Consisency y α + βx + he esimaes approach he rue values as he sample size increases Esimaion

More information

Study and Evaluation of Performances of the Digital Multimeter

Study and Evaluation of Performances of the Digital Multimeter Sudy and Evaluaion of Performances of he Digial Mulimeer. Inroducion Mulimeers allow measuring he mos common elecrical quaniies: volage, curren, resisance. They evolved from simle ones based on elecromechanical

More information

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests

Outline. lse-logo. Outline. Outline. 1 Wald Test. 2 The Likelihood Ratio Test. 3 Lagrange Multiplier Tests Ouline Ouline Hypohesis Tes wihin he Maximum Likelihood Framework There are hree main frequenis approaches o inference wihin he Maximum Likelihood framework: he Wald es, he Likelihood Raio es and he Lagrange

More information

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution

Physics 127b: Statistical Mechanics. Fokker-Planck Equation. Time Evolution Physics 7b: Saisical Mechanics Fokker-Planck Equaion The Langevin equaion approach o he evoluion of he velociy disribuion for he Brownian paricle migh leave you uncomforable. A more formal reamen of his

More information

Generalized Least Squares

Generalized Least Squares Generalized Leas Squares Augus 006 1 Modified Model Original assumpions: 1 Specificaion: y = Xβ + ε (1) Eε =0 3 EX 0 ε =0 4 Eεε 0 = σ I In his secion, we consider relaxing assumpion (4) Insead, assume

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS

Financial Econometrics Jeffrey R. Russell Midterm Winter 2009 SOLUTIONS Name SOLUTIONS Financial Economerics Jeffrey R. Russell Miderm Winer 009 SOLUTIONS You have 80 minues o complee he exam. Use can use a calculaor and noes. Try o fi all your work in he space provided. If

More information

Second Order Linear Differential Equations

Second Order Linear Differential Equations Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous

More information

Data Fusion using Kalman Filter. Ioannis Rekleitis

Data Fusion using Kalman Filter. Ioannis Rekleitis Daa Fusion using Kalman Filer Ioannis Rekleiis Eample of a arameerized Baesian Filer: Kalman Filer Kalman filers (KF represen poserior belief b a Gaussian (normal disribuion A -d Gaussian disribuion is

More information

GMM - Generalized Method of Moments

GMM - Generalized Method of Moments GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................

More information

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1

Vectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1 Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies

More information

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models

A Specification Test for Linear Dynamic Stochastic General Equilibrium Models Journal of Saisical and Economeric Mehods, vol.1, no.2, 2012, 65-70 ISSN: 2241-0384 (prin), 2241-0376 (online) Scienpress Ld, 2012 A Specificaion Tes for Linear Dynamic Sochasic General Equilibrium Models

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

Math 10B: Mock Mid II. April 13, 2016

Math 10B: Mock Mid II. April 13, 2016 Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.

More information

Measurement Error 1: Consequences Page 1. Definitions. For two variables, X and Y, the following hold: Expectation, or Mean, of X.

Measurement Error 1: Consequences Page 1. Definitions. For two variables, X and Y, the following hold: Expectation, or Mean, of X. Measuremen Error 1: Consequences of Measuremen Error Richard Williams, Universiy of Nore Dame, hps://www3.nd.edu/~rwilliam/ Las revised January 1, 015 Definiions. For wo variables, X and Y, he following

More information

Testing the Random Walk Model. i.i.d. ( ) r

Testing the Random Walk Model. i.i.d. ( ) r he random walk heory saes: esing he Random Walk Model µ ε () np = + np + Momen Condiions where where ε ~ i.i.d he idea here is o es direcly he resricions imposed by momen condiions. lnp lnp µ ( lnp lnp

More information

Modeling and Forecasting Volatility Autoregressive Conditional Heteroskedasticity Models. Economic Forecasting Anthony Tay Slide 1

Modeling and Forecasting Volatility Autoregressive Conditional Heteroskedasticity Models. Economic Forecasting Anthony Tay Slide 1 Modeling and Forecasing Volailiy Auoregressive Condiional Heeroskedasiciy Models Anhony Tay Slide 1 smpl @all line(m) sii dl_sii S TII D L _ S TII 4,000. 3,000.1.0,000 -.1 1,000 -. 0 86 88 90 9 94 96 98

More information

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,

More information

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes

23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals

More information

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum

Stopping Brownian Motion without Anticipation as Close as Possible to its Ultimate Maximum Theory Probab. Al. Vol. 45, No.,, (5-36) Research Reor No. 45, 999, De. Theore. Sais. Aarhus Soing Brownian Moion wihou Aniciaion as Close as Possible o is Ulimae Maximum S. E. GRAVERSEN 3, G. PESKIR 3,

More information

The Arcsine Distribution

The Arcsine Distribution The Arcsine Disribuion Chris H. Rycrof Ocober 6, 006 A common heme of he class has been ha he saisics of single walker are ofen very differen from hose of an ensemble of walkers. On he firs homework, we

More information

Econ Autocorrelation. Sanjaya DeSilva

Econ Autocorrelation. Sanjaya DeSilva Econ 39 - Auocorrelaion Sanjaya DeSilva Ocober 3, 008 1 Definiion Auocorrelaion (or serial correlaion) occurs when he error erm of one observaion is correlaed wih he error erm of any oher observaion. This

More information

Solutions to Odd Number Exercises in Chapter 6

Solutions to Odd Number Exercises in Chapter 6 1 Soluions o Odd Number Exercises in 6.1 R y eˆ 1.7151 y 6.3 From eˆ ( T K) ˆ R 1 1 SST SST SST (1 R ) 55.36(1.7911) we have, ˆ 6.414 T K ( ) 6.5 y ye ye y e 1 1 Consider he erms e and xe b b x e y e b

More information

Dirac s hole theory and the Pauli principle: clearing up the confusion.

Dirac s hole theory and the Pauli principle: clearing up the confusion. Dirac s hole heory and he Pauli rincile: clearing u he conusion. Dan Solomon Rauland-Borg Cororaion 8 W. Cenral Road Moun Prosec IL 656 USA Email: dan.solomon@rauland.com Absrac. In Dirac s hole heory

More information

Properties of Autocorrelated Processes Economics 30331

Properties of Autocorrelated Processes Economics 30331 Properies of Auocorrelaed Processes Economics 3033 Bill Evans Fall 05 Suppose we have ime series daa series labeled as where =,,3, T (he final period) Some examples are he dail closing price of he S&500,

More information

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:

Math 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities: Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial

More information

13.3 Term structure models

13.3 Term structure models 13.3 Term srucure models 13.3.1 Expecaions hypohesis model - Simples "model" a) shor rae b) expecaions o ge oher prices Resul: y () = 1 h +1 δ = φ( δ)+ε +1 f () = E (y +1) (1) =δ + φ( δ) f (3) = E (y +)

More information

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes

23.5. Half-Range Series. Introduction. Prerequisites. Learning Outcomes Half-Range Series 2.5 Inroducion In his Secion we address he following problem: Can we find a Fourier series expansion of a funcion defined over a finie inerval? Of course we recognise ha such a funcion

More information

How to Deal with Structural Breaks in Practical Cointegration Analysis

How to Deal with Structural Breaks in Practical Cointegration Analysis How o Deal wih Srucural Breaks in Pracical Coinegraion Analysis Roselyne Joyeux * School of Economic and Financial Sudies Macquarie Universiy December 00 ABSTRACT In his noe we consider he reamen of srucural

More information

Problem Set 7-7. dv V ln V = kt + C. 20. Assume that df/dt still equals = F RF. df dr = =

Problem Set 7-7. dv V ln V = kt + C. 20. Assume that df/dt still equals = F RF. df dr = = 20. Assume ha df/d sill equals = F + 0.02RF. df dr df/ d F+ 0. 02RF = = 2 dr/ d R 0. 04RF 0. 01R 10 df 11. 2 R= 70 and F = 1 = = 0. 362K dr 31 21. 0 F (70, 30) (70, 1) R 100 Noe ha he slope a (70, 1) is

More information

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006 2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)

More information

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions

Inventory Analysis and Management. Multi-Period Stochastic Models: Optimality of (s, S) Policy for K-Convex Objective Functions Muli-Period Sochasic Models: Opimali of (s, S) Polic for -Convex Objecive Funcions Consider a seing similar o he N-sage newsvendor problem excep ha now here is a fixed re-ordering cos (> 0) for each (re-)order.

More information

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H.

ACE 562 Fall Lecture 8: The Simple Linear Regression Model: R 2, Reporting the Results and Prediction. by Professor Scott H. ACE 56 Fall 5 Lecure 8: The Simple Linear Regression Model: R, Reporing he Resuls and Predicion by Professor Sco H. Irwin Required Readings: Griffihs, Hill and Judge. "Explaining Variaion in he Dependen

More information

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds

Kriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1

Chapter 5. Heterocedastic Models. Introduction to time series (2008) 1 Chaper 5 Heerocedasic Models Inroducion o ime series (2008) 1 Chaper 5. Conens. 5.1. The ARCH model. 5.2. The GARCH model. 5.3. The exponenial GARCH model. 5.4. The CHARMA model. 5.5. Random coefficien

More information

OBJECTIVES OF TIME SERIES ANALYSIS

OBJECTIVES OF TIME SERIES ANALYSIS OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging

More information

Smoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T

Smoothing. Backward smoother: At any give T, replace the observation yt by a combination of observations at & before T Smoohing Consan process Separae signal & noise Smooh he daa: Backward smooher: A an give, replace he observaion b a combinaion of observaions a & before Simple smooher : replace he curren observaion wih

More information

1. Diagnostic (Misspeci cation) Tests: Testing the Assumptions

1. Diagnostic (Misspeci cation) Tests: Testing the Assumptions Business School, Brunel Universiy MSc. EC5501/5509 Modelling Financial Decisions and Markes/Inroducion o Quaniaive Mehods Prof. Menelaos Karanasos (Room SS269, el. 01895265284) Lecure Noes 6 1. Diagnosic

More information

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito

Intermediate Macroeconomics: Mid-term exam May 30 th, 2016 Makoto Saito 1 Inermediae Macroeconomics: Mid-erm exam May 30 h, 2016 Makoo Saio Try he following hree roblems, and submi your answer in handwrien A4 aers. You are execed o dro your aers ino he mailbox assigned for

More information

Kalman filtering for maximum likelihood estimation given corrupted observations.

Kalman filtering for maximum likelihood estimation given corrupted observations. alman filering maimum likelihood esimaion given corruped observaions... Holmes Naional Marine isheries Service Inroducion he alman filer is used o eend likelihood esimaion o cases wih hidden saes such

More information

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A

Licenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In

More information

Computing with diode model

Computing with diode model ECE 570 Session 5 C 752E Comuer Aided Engineering for negraed Circuis Comuing wih diode model Objecie: nroduce conces in numerical circui analsis Ouline: 1. Model of an examle circui wih a diode 2. Ouline

More information

Technical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.

Technical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models. Technical Repor Doc ID: TR--203 06-March-203 (Las revision: 23-Februar-206) On formulaing quadraic funcions in opimizaion models. Auhor: Erling D. Andersen Convex quadraic consrains quie frequenl appear

More information

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015

Explaining Total Factor Productivity. Ulrich Kohli University of Geneva December 2015 Explaining Toal Facor Produciviy Ulrich Kohli Universiy of Geneva December 2015 Needed: A Theory of Toal Facor Produciviy Edward C. Presco (1998) 2 1. Inroducion Toal Facor Produciviy (TFP) has become

More information

Regression with Time Series Data

Regression with Time Series Data Regression wih Time Series Daa y = β 0 + β 1 x 1 +...+ β k x k + u Serial Correlaion and Heeroskedasiciy Time Series - Serial Correlaion and Heeroskedasiciy 1 Serially Correlaed Errors: Consequences Wih

More information

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.

Nature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time. Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike

More information

Time series Decomposition method

Time series Decomposition method Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

References are appeared in the last slide. Last update: (1393/08/19)

References are appeared in the last slide. Last update: (1393/08/19) SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be

More information

Testing H 0 : ρ = 0: Comparing A Single Correlation to Zero

Testing H 0 : ρ = 0: Comparing A Single Correlation to Zero Tesing H 0 : ρ = 0: Comaring A Single Correlaion o Zero February 5, 2018 Conens Tesing H 0 : ρ = 0 Examle 1: The correlaion beween heighs of your arens Examle 2: Drinking and Facebook friends Using R o

More information

fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /10/10

fakultät für informatik informatik 12 technische universität dortmund Petri Nets Peter Marwedel TU Dortmund, Informatik /10/10 2 Peri Nes Peer Marwedel TU Dormund, Informaik 2 2008/0/0 Grahics: Alexandra Nole, Gesine Marwedel, 2003 Generalizaion of daa flow: Comuaional grahs Examle: Peri nes Inroduced in 962 by Carl Adam Peri

More information

A One Line Derivation of DCC: Application of a Vector Random Coefficient Moving Average Process*

A One Line Derivation of DCC: Application of a Vector Random Coefficient Moving Average Process* A One Line Derivaion of DCC: Alicaion of a Vecor Random Coefficien Moving Average Process* Chrisian M. Hafner Insiu de saisique, biosaisique e sciences acuarielles Universié caholique de Louvain Michael

More information

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n

Module 2 F c i k c s la l w a s o s f dif di fusi s o i n Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms

More information

Errata (1 st Edition)

Errata (1 st Edition) P Sandborn, os Analysis of Elecronic Sysems, s Ediion, orld Scienific, Singapore, 03 Erraa ( s Ediion) S K 05D Page 8 Equaion (7) should be, E 05D E Nu e S K he L appearing in he equaion in he book does

More information

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance

Lecture 5. Time series: ECM. Bernardina Algieri Department Economics, Statistics and Finance Lecure 5 Time series: ECM Bernardina Algieri Deparmen Economics, Saisics and Finance Conens Time Series Modelling Coinegraion Error Correcion Model Two Seps, Engle-Granger procedure Error Correcion Model

More information

F x 3-4x 2 >25-6x + 6 = 0, whose only real root is x a. relative minimum of at x 1.404; c. 3p2 - p3

F x 3-4x 2 >25-6x + 6 = 0, whose only real root is x a. relative minimum of at x 1.404; c. 3p2 - p3 or Use Onl in Pilo Program F Seleced Answers A-7. 7. Abs min (, ) Local and abs ma (, ) Inflecion oin (, ) Inflecion oin (, ) Abs min (, ) Local and abs ma ((),.) Inflecion oin ((), 9.9) -inerce (, ) -inerce

More information

EXAMPLES OF EVALUATION OF RATE FORMS FROM KINETIC DATA IN BATCH SYSTEMS

EXAMPLES OF EVALUATION OF RATE FORMS FROM KINETIC DATA IN BATCH SYSTEMS HE 47 LETURE 5 EXMPLES OF EVLUTION OF RTE FORMS FROM INETI DT IN BTH SYSTEMS EXMPLE : Deermine he reacion order and he rae consan or a single reacion o he ye roducs based on he ollowing exerimenal inormaion

More information

And the solution to the PDE problem must be of the form Π 1

And the solution to the PDE problem must be of the form Π 1 5. Self-Similar Soluions b Dimensional Analsis Consider he diffusion problem from las secion, wih poinwise release (Ref: Bluman & Cole, 2.3): c = D 2 c x + Q 0δ(x)δ() 2 c(x,0) = 0, c(±,) = 0 Iniial release

More information

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel

15. Bicycle Wheel. Graph of height y (cm) above the axle against time t (s) over a 6-second interval. 15 bike wheel 15. Biccle Wheel The graph We moun a biccle wheel so ha i is free o roae in a verical plane. In fac, wha works easil is o pu an exension on one of he axles, and ge a suden o sand on one side and hold he

More information

Forecasting optimally

Forecasting optimally I) ile: Forecas Evaluaion II) Conens: Evaluaing forecass, properies of opimal forecass, esing properies of opimal forecass, saisical comparison of forecas accuracy III) Documenaion: - Diebold, Francis

More information

Some Basic Information about M-S-D Systems

Some Basic Information about M-S-D Systems Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,

More information

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon 3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of

More information

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates

Biol. 356 Lab 8. Mortality, Recruitment, and Migration Rates Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese

More information

PMT. Version. General Certificate of Education (A-level) January 2013 MPC4. Mathematics. (Specification 6360) Pure Core 4. Final.

PMT. Version. General Certificate of Education (A-level) January 2013 MPC4. Mathematics. (Specification 6360) Pure Core 4. Final. Version General Cerificae of Educaion (A-level) January 0 Mahemaics MPC (Secificaion 660) Pure Core Final Mark Scheme Mark schemes are reared by he Princial Examiner and considered, ogeher wih he relevan

More information

ln 2 1 ln y x c y C x

ln 2 1 ln y x c y C x Lecure 14 Appendi B: Some sample problems from Boas Here are some soluions o he sample problems assigned for Chaper 8 8: 6 Soluion: We wan o find he soluion o he following firs order equaion using separaion

More information

) were both constant and we brought them from under the integral.

) were both constant and we brought them from under the integral. YIELD-PER-RECRUIT (coninued The yield-per-recrui model applies o a cohor, bu we saw in he Age Disribuions lecure ha he properies of a cohor do no apply in general o a collecion of cohors, which is wha

More information

(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum.

(10) (a) Derive and plot the spectrum of y. Discuss how the seasonality in the process is evident in spectrum. January 01 Final Exam Quesions: Mark W. Wason (Poins/Minues are given in Parenheses) (15) 1. Suppose ha y follows he saionary AR(1) process y = y 1 +, where = 0.5 and ~ iid(0,1). Le x = (y + y 1 )/. (11)

More information

An EM algorithm for maximum likelihood estimation given corrupted observations. E. E. Holmes, National Marine Fisheries Service

An EM algorithm for maximum likelihood estimation given corrupted observations. E. E. Holmes, National Marine Fisheries Service An M algorihm maimum likelihood esimaion given corruped observaions... Holmes Naional Marine Fisheries Service Inroducion M algorihms e likelihood esimaion o cases wih hidden saes such as when observaions

More information

KINEMATICS IN ONE DIMENSION

KINEMATICS IN ONE DIMENSION KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec

More information

STA 114: Statistics. Notes 2. Statistical Models and the Likelihood Function

STA 114: Statistics. Notes 2. Statistical Models and the Likelihood Function STA 114: Saisics Noes 2. Saisical Models and he Likelihood Funcion Describing Daa & Saisical Models A physicis has a heory ha makes a precise predicion of wha s o be observed in daa. If he daa doesn mach

More information

Lecture 33: November 29

Lecture 33: November 29 36-705: Inermediae Saisics Fall 2017 Lecurer: Siva Balakrishnan Lecure 33: November 29 Today we will coninue discussing he boosrap, and hen ry o undersand why i works in a simple case. In he las lecure

More information

Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size.

Methodology. -ratios are biased and that the appropriate critical values have to be increased by an amount. that depends on the sample size. Mehodology. Uni Roo Tess A ime series is inegraed when i has a mean revering propery and a finie variance. I is only emporarily ou of equilibrium and is called saionary in I(0). However a ime series ha

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information