A STATISTICAL SURVEY OF EARTHQUAKES IN THE MAIN SEISMIC REGION OF NEW ZEALAND

Size: px
Start display at page:

Download "A STATISTICAL SURVEY OF EARTHQUAKES IN THE MAIN SEISMIC REGION OF NEW ZEALAND"

Transcription

1 Reprinted frm the NEW ZEALAND JOURNAL OF GEOLOGY AND GEOPHYSICS, Vl. 9, N.3 A STATISTICAL SURVEY OF EARTHQUAKES IN THE MAIN SEISMIC REGION OF NEW ZEALAND PART 2-TIME SERIES ANALYSES D. VERE-JONES* and R. B. DAVIES Applied Mathematics Divisin, Department f Scientific and Industrial Research, Wellingtn (Receit1ed fr publicatin 24 March 1965) ABSTRACT Time series analyses are carried ut n earthquake data frm the main seismic regin f New Zealand fr the years Origin times nly are cnsidered, the energies and exact psitins f the shcks being largely ignred. The relevant statis tical thery fr the first and secnd rder prperties f the prcess is described, and simple prbability mdels fr earthquake ccurrence are put frward. On the basis f these results, the. data are examined fr peridic and gruping effects. N significant peridic effects are fund, either amng the shallw shcks (depths up t 100 km) r amng the deep shcks (depths 100 km r greater). Bth cmpnents shw strng evidence f gruping, and several alternative mdels t describe this effect are put frward and cmpared. l-introduction Time series methds have been used fr a lng time in the analysis f earthquake frequencies, but many f the papers make little use f the statistical techniques which are nw available. Our aim in this paper is t apply sme f these techniques, especially spectral analysis, t earthquake recrds frm New Zealand. Any statistical analysis relies ultimately n a prbability mdel f the prcess under cnsideratin, and in the curse f the paper we shall develp sme simple mdels fr earthquake ccurrence. Like the chemical prcesses emplyed by the Curies, these methds can be used t extract a small quantity f active infrmatin frm a large bdy f inert data. We hpe that in the present cntext they may prve useful t the seismlgist. In particular, it is pssible that the study f suitable stchastic mdels may ultimately thrw sme further light n earthquake mechanism. ;'The recrds frm a small, self-cntained seismic regin (f which the lli-ln seismic regin f New Zealand appears t furnish a reasnable cxlunple) reveal time series f extremely cmplex structure. Large fiuctua 'Jins in the numbers f shcks per year, cmplicated sequences f related shcks, dependence f statistical parameters n depth and perhaps magnitude, and fluctuatins f activity n a larger time-scale, all appear t be haracteristic features f such recrds. An earlier paper (Vere-Jnes, *Nw with Department f Statistics, Australian Natinal University, Canberra. N.z. J. Gel. GeaphyJ. 9:

2 J 252 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL. 9 Turnvsky, and Eiby, 1964) attempted t ascertain the slw fluctuatins in activity {in the frm f. linear trends) frm data cmpiled by the New Zealand Seismlgical Ob$ervatry ver the perid In this paper we shall use the same data as the basis f a time-series study. In general, we shall ignr: the, lcatins and magnihldes f the shcks, except in s far as they determmed the fbur classes f shcks described in the earlier paper: (a) The ttal shcks. (b) The shallw shcks (depths less than 100 km).* (c) The deep shcks (depths greater than r equal t 100 km). (d) The mdified list f shallw shcks, btained after the large swarms and aftershck sequences have been remved. We refer the reade'r t the earlier paper fr a mre detailed accunt f the data, and a discussin f its limitatins. It is restricted t the main seismic regin f New Zealand, and t shcks which have an instrumental magnitude f 4'5 r greater. After remval f linear trends, the data frm each class were assumed t represent data frm a statinary time series. There is, f curse, n reasn t suppse that all the nn-statinarity is cntained in the linear trends, but we culd regard ur prcedure as the first step in a series f apprximatins. Thus we ca:n use the data after remval f linear trends t frm a rugh estimate f the cvariance structure; knwing this, we can then estimate a mre cmplex trend, and s n. The techniques we shall use are standard methds f time-series analysis, mdified slightly t take int accunt the discrete nature f earthquake ccurrence. Sme f the ideas may be new in the present cntext, s that we have discussed them in mre detail than wuld usually be necessary. Sectin 2 cntains backgrund material, while sectins are cncerned with particular mdels fr earthquake ccurrence. The technique f spectral analysis is described in sectin 2.3, and the fllwing three sectins cntain sme applicatins. In sectin 2.7 we cmment briefly n energies. Sectins 3' cntain evaluatins f the results fr each f the fur classes f earthquakes in turn. 2.0-STOCHASTIC POINT PROCESSES The stchastic prcesses relevant t a descriptin f earthquake ccurrence are knwn as "stchastic pint prcesses". In these prcesses events ccur singly, r perhaps in small grups, at (t a first apprximatin) instants f time. They are t be distinguished frm the stchastic prcesses which develp cntinuusly in time. The prbability structure f a pint prcess is described by the distributins f the number f events in particular intervals r cllectins f intervals. Fr practical purpses, the first arid secnd rder prperties f these distributins (the means, variances, and cvariances), are especially imprtant. These quantities can, in general, be written dwn in terms f tw basic functins, the instantaneus rate met), *We have included the 100 km shcks with the deep shcks instead f the shallw shcks as in the earlier paper.

3 7 N.3 VERE-JONES & DAVIES - EARTIlQUAKE ANALYSES 253 and cvariance.dsity C(u, t). The instantaneus rate at a given time t is efined as the limit f the mean rate f ccurrence f events in a small time mterval abut t, as the length f the interval tends t zer, m(l) = lim E[N(t - a) t +,8]/(a +,8) (1) a,/ (t1,.t2) her rresents the number f events (a randm variable) ccurnng ill the tllne mtervals (f 1, f 2 ). The prcess is assumed t be sufficiently regular fr the limit t exist and t be finite (and fr ur purpses bunded) in the whle f the time interval under cnsideratin. The cvariance density is defined as the limit f the cvariance f the numbers f events in tw small time intervals, divided by the lengths f the time intervals, as the tw lengths tend t zer, C(u, t) [N(t- a, t + {3), N(u -y, u + 8)] lim Cv a,/3,'y, (a+{3)(y+8) (2) Again it is assed that this limit exists and is finite fr all relevant pairs f values f u ana f.. The mean and variance f the numbers f events in any interval (t1' t2) are given by (3) I r t 2 C(u, t) du dt (4) J t1 and the cvariance f the numbers f events in tw disjint intervals (t1' t2) and (t8,/4) by (5) When the prcess is statinary-when its prbablity struture is unchanging with time-these equatins assume a partlrularly Simple frm. In this case m(t) must be a cnstant (say m) and C(u, t) mst reduce.t? a functin C(v) f the difference v = (t - u) alne. In particular (wntmg T = 11 - t 2 ), Equati()ns (3) and (4) take the frm (6) B[N(T)] = mt Var [.N(T)] = mt + 2 J: (T - u) C(u) du (7) Inset-l

4 254 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL THE POISSON MODEL The natural starting pint fr a discussin'f stchastic mdels fr earthquke ccurrenc is the Pissn' prcess. The fundamental prperty f the POIssn prcess IS that the numbers f events in any tw disjint intervals are inqependent randm variables. This prperty and the supplementary cnditins that the prbability shuld be asympttically prprtinal t the length f the interval, and that the prbability f mre than ne event in the small interval shuld be asympttically negligible ((dt) ), are sufficient t specify the whle prbability structure f the prcess. The prbability f exactly n events ccurring in the interval (t1' t 2 ) given by the frmula. P.. = A(t1' t2 ).. exp { - A(t1' t2)}in! (8) where A(t1,t2) = ft2 ;\(t)dt and ;\(t) is the instantaneus rate defined in 11 equatin (1). Frm equatins (3) and (4) we have E[N(t1' t2)} = Var [N(t1' t2)} = A(t!) t2) (9) If the prcess is statinary we can write ;\(t) = ;\ (a cnstant) and A(tv t 2 ) = (t2-11);\, Then E[N(t1' t2)} = Var [N(t1' t2)} = (t2 - t1);\ (10) Many prcesses apprximate quite well t' the Pissn prcess in spite f its restrictive definitin. One reasn fr this is that a prcess built up f a number f small, independent cmpnents will resemble a Pissn prcess mre and mre clsely as the number f. cmpnents increases, and the cntributin frm each decreases (Khinchin, 1955). It might be thught, therefre, that the prcess cnsidered here, representing the sum f cntributins frm many different areas and depths, wuld be f this type. We shall see that this is nt the case. A cnsiderable prprtin f the shcks in any interval appear t cme frm nly a small number f active surces. Sme reasns fr suppsing that the Pissn mdel is nt satisfactry were put frward in the earlier paper. This cnclusin is strengthened by a direct examinatin f the cvariance structure and the distributin f interval lengths. T determine the cvariance structure, an autcrrelatin analysis was carried ut n the numbers f shcks ccurring in successive 0'1 year perids. A separate analysis was carried ut fr the fur classes f shcks listed in sectin 1. In each case a linear trend was remved frm the uncrrected data and the autcvariances cmputed frm the frmula 1 N-s N- r YrYr+B is (11)

5 N.3 VERE-JONES & DAVIES - EARTltQUAKE ANALYSES 255 where N is the ttal number f intervals (here 200), the Yr are the transfrmed variables (after the trend and mean has been remved), and the "v" n the symbl dentes an estimated value. Values fr the first 10 autcvariances are given in Table 1a. Fully efficient tests fr dependence are hard t establish, but with 200 bservatins, the distributin f the cjf; shuld be apprximately that f the rdinary crrelatin cefficient, and n this assumptin, the first f the ratis cj/c can be used t test whether there is sme degree f crrelatin present. A ne-sided test is apprpriate, and the 5%, 2%, and 1% values fr each grup (scaled up by the factr C in each case) are listed in Table lb. In all cases the value f c 1 is significant at the 5% lel, and fr the shallw shacks and fr all shcks at the 1 % level. This evidence is f mt interest fr the deep shcks, fr which it has usually been assumed that there are few r n grups f related shcks. Hwever, as pinted ut in the earlier paper, the printed lists shw numerus grups f deep shcks, ccurring verlerids f a year r s, within very narrwly defined limits f latitude an depth. In this paper als, we find repeated evidence f dependence amng the deep shcks. The high value f c 1 in the furth rw f Table 1a pints in a similar manner t a residual dependence amng the shallw shcks, even after the bvius grups f swarms and after shcks have been remved. A mre sensitive indicatr f deviatins frm the Pissn mdel is the rati f the variance t the mean (the Pissn index f dispersin). This has the expect value-f a unity fr a Pissn prcess (Equatin 9), and fr a large number f bservatins its distributin is asympttically prp.rtinal t a X2 distributin with (N - 1) degrees f freedm. The values f this index are shwn in Table 1c. Since the 0'1 % pint is 1'33, fr all grups, the Pissn hypthesis is rejected at a very high level. The large fluctuatins are even mre striking if we take a larger perid. v Fig. 1 shws the bserved variance (actually V (T)IT) fr the shallw shcks fr a perid f length T fr different values f T. The variances have been calculated frm the frmula (Murthy, 1961)., 1 JT - 7' T2 VeT) = N2(t, t + T)dt- _N2(0, T)(12) T - T 0 T2 where T is the ttal time f the sequence cnsidered. The lwer hrizntal line n the same figure represents the variance fr a Pissn prcess with the same rate (V(T) = AT). Final evidence f gruping cmes frm the distributin f the interval lengths between shcks. In Table 2 we have set ut the distributin f interval lengths between shcks fr tw series f 101 shcks, ne series fr a quiet perid ( ), and ne fr an active pid' ( ). In bth cases the numbers in each cell were cmpared Wlth the expected numbers frm an expnential distributin having th same mean. The X 2 value fr the quiet perid is. smewhat reaer than lts expect:d vale (8), bu nt significantly s. ThlS seems t mdlcate that fr a qulet pend the POlssn

6 IV VI ;. 0'> [ ' IOO , ;;-:;-r---- 0, // 0// Pissn,(,) I- -7 Expnentil mdel, / / st pwer law mdel I // 2nd pwer law mdel I / / z / /// N 50 '0' 1/ '/ ---L--125 Iwk -05yr -Iyr (- z l'" > fii ii / > Z ;;: 100 (/P-(f =-= =---.:: ==... t C) t'1 't:i I"'l en I I I I FIG. 1--Plt. f C.(r) against perid l' fr different mdels (shallw shcks). Years <: 0 r '0

7 N.3 VERE-JONES & DAVIES - EARTHQUAKE ANALYSES 257 TABLE 1a-Autcvariances fr Numbers f Shcks in 0 1 Year Perids Shck Types C C, C. C, C. C. C. C, C. C Ct. All shcks Shallw shcks Deep shcks Mdified shallw shcks TABLE 1b--Significance Levels fr Ct TABLE Ie-Mean and Pissn Index Shck Types 5% 2% 1% Mean Pissn Index All shcks Shallw shcks ; Deep shcks Mdified shallw shcks TABLE 2-Distributin f Time Intervals between Shcks SERIES 1-( ) Intervals (days) Observed N Expected N. S l S S Mean length = 8 0 days. x' (8 degrees f freedm) = 11. SERIES 11-( ) Intervals (days) S Observed N S :; 7 Expected N ; Mean length = 4 5 days. X (8 degrees f freedm) = The 5% and 1% vlues f x with 8 degrees f freedm are 15 5 and 20 1, respectively.

8 258 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL. 9 apprximatin may nt be unreasnable. Fr the active perid, hwever, the X 2 value is significant at the 1% level indicating that the Pissn apprximatin is then quite inadequate. 1n each the largest deviatin is at zerthere are mre clsely spaced shcks than there shuld be n the Pissn mdel. The psibi1i/ f describing earthquake ccurrenc by a Pissn press has been InvestIgated by a number f authrs, especially by V. N. GalSky, wh has studied the recrds f Central Asian earthquakes. These investigatins have chiefly cncerned the ccurrence f large shcks, fr which n appreciabte divergence frm the Pissn law is fund (Gaisky, 1961; Gaisky and Katk, 1960; Gaisky and Birman, 1962). These results d nt frm a real cntradictin t thse f the present paper since we are nt examining the distributin f large shcks, but the crrelatins amng shcks f small r mderate size. Indeed such results give an added plausibility t the mdels discussed in the next chapter, where it is suppsed' that grups f small shcks are initiated by majr events beying the Pissn Law. An earlier paper by Ljapunv and Fandjusina (1950) shwed that with Central Asian earthquakes, als, discrepancies frm the simple Pissn law may ccur when earthquakes f smaller size are cnsidered., 2.2-TRIGGER MODELS The results s far established pint t a significant degree f dependence amng earthquakes. There is a psitive crrelatin amng the numbers f earthquakes in successive time intervals, which prduces (In accrdance with Equatin 4) a high variance fr the numbers f shcks in an interval f fixed length. Such prperties place the earthquake prcess amng the general class f "cntagius prcesses" which have been widely studied in applicatin t eclgy and frestry (Matern, 1960), t csmlgy (Neyman and Sctt, 1958), and in ther cntexts. Fr earthquakes it is natural t suppse that the crrelatin effects arise frm the ccurrence f grups f statistically related shcks. The mdels that we shall cnsider in this sectin are all based n this idea. We shall nt, hwever, immediately identify these grups with earthquake swarms, after-shck sequences, and ther grups which are knwn t llave ccurred. In view f the practical difficulties f identifying all members f a given swarm r after-shck sequence we shall at first make n assumptins as t the nature f the grups, but attempt t ascertain indirectly, by estimating the larameters f certain mdels, the extent and character f the clustering e ect. The simplest mdel t allw fr gruping is the cmpund Pissn prcess, in which grups f events ccur tgether at the instants f a simple Pissn prcess. This mdel is nt satisfactry as it stands, since in practice the grup will spread ver a finite and cnceivably large segment f the time I axis. A mre suitable class ()If mdel can be built up in the fllwing way. Suppse that events initiating grups f shcks ("trigger events"-nt necessarily themselves shcks)* ccur at the instants f a *The case in which the trigger event is itself a shck can be brught within the scpe f the analysis f allwing A},,(t) t have a discrete cmpnent (atm) at t = O.

9 N.3 VERE-JONES & DAVIES-EARTIIQUAKE ANALYSES 259 si';dple Pissn trcess with cnstant rate J-t. Let '\(t) be a decay functin with the prperhes '\(t) 0 (t < 0) '\(t) > 0 (t >0) (13) f '\(t)dt = 1 and suppse that the cnditinal prbability that a shck will ccur in the small time interval (t + x, t +.'JC + dx), given that a trigger event ccurred at time t, is independent f t and equal t A'\(x). We shall assume als that A itself, the quantity characterising the size f the grups, is a randm variable with distributin F(a), finite mean a and variance v a= E(A) = f adf(a) (14) v = Var (A) = f (a - a)2 df(a) Finally we shall assume that separate trigger events generate independent sequences f shcks, and that the randm variables A, assciated with distinct trigger events, are als independent. It is nt difficult t shw that these assumptins define a statinary pint prcess, fr which the instantaneus rate is given by m=p.di (15) and cvariance density by C(u) =J-t(a 2 + v) f '\(t)'\(1 + u)dt (16) Hence (frm Equatins 6 and 7) the mean and variance f the numbers f events in an interval f length T are given by B[N(T)] pflt (17) _ C(r) = pat + 2.J-t(tP + v) rt (r - u)c(ti)du (18) where c(u) - J '\(t) '\(t + ti)dt. I J 0

10 260 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL. 9 Similarly the rth cvariance (the cvariance between the numbers f shcks in tw intervals rt apart, each f length T) is given by t er(t) = (a2 + v) (7 - u)[e(rt + u) + c(rt - u)]du (19a) Fr T large, (20) The right-hand side f this expressin is the variance f a cmpund Pissn prcess fr which grups f events ccur at rate ft, and the size f the grup has the distributin (21) This is als the distributin f the numbers f shcks in a grup initiated by a single trigger event. Thus Equatin (20) expresses the fact that when the interval becmes very large, we may ignre effects due t the shcks within a single grup spreading ut alng the time axis. It is imprtant in practice t knw hw quickly the asympttic frm (Equatin 20) is reached. The faster the functin A(t) appraches zer, the faster the asympttic frm is reached. The bserved variances (Fig. 1) suggest that the asympttic frm is reached nly slwly, and hence that the decay functin itself is nt a rapid ne.. Such mdels have a fair degree f bth intuitive appeal and f generality. Their secnd-rder prperties depend n the tw parameters p.a and p.(a2 + v) and n the decay functin A(t). Different chices f these quantities give rise. t particular mdels that can be tested fr fit against the bserved prcess. We shall cnsider in detail tw chices f A(t), ne with an expnential decay A(t) = pe-pf (p > 0) (22) and the ther with an inverse pwer-law decay AU) = pepi(e + t)p + 1 (e > 0) (23) 2.3-THE SPECTRUM Estimates fr the parameters f these mdels and the tests f gdness f fit will nt be based directly n the cvariances described in the previus sectin, but n their Furier transfrm, the spectrum,

11 N.3 VERE-JONES & DAVIES - EARTIiQUAKE ANALYSES 261 t (<0) = 00 C,,(T) cs nw (24) 'T -00 The bserved quantity crrespnding t the theretical spectrum is the peridgram I (w) T (25) (25a) where the C,,(T) = C-,,(T) are the bserved cvariances defined in Equatin (11), Xt is the number f shcks in the tth interval, N is the ttal number f intervals cnsidered, and xt is the average f the x t. Accrding t Equatins (24) and (25) the spectrum may be interpreted as an analysis f the ttal variance int cmpnents crrespnding t different frequencies. In the case f a cntinuus variable--say a vltage varying in time--the spectrum may be regarded as giving the average pwer transferred at different frequencies. In the case f a stchastic pint-prcess this interpretatin n lnger hlds, but the analgy may give sme idea f the meaning f the spectrum t thse nt familiar with the thery. In particular, a cyclic effect in the Xt with perid T (measured in units f the basic interval T) willgiverisetapeakini (w) atw = 27r/T. T Under wide cnditins, the rdinates Ij = 1(00;) (Wj = 27rj/N, j = 0, 1,..., N - 1) are apprximately uncrrelated. When the Xt are independent, the I j are als apprximately identically distributed, and in particular, the expected value f If is then apprximately equal t the cnstant value u = 2 vat (Xt) fr all values f j. Fr such a prcess the theretical spectrum is just a hrizntalline at height u2 These results are the basis f mst f the simple tests fr cyclic effects. In the case f a prcess with spectrum few), the rati 1/ = Ij/f f has prperties similar t thse f the peridgram arising frm independent Xt (Hannan, 1960). We shall use Ii' fr tests f gdness f fit f ur prpsed mdels. The peridgram is nt a cnsistent estimatr f the spectrum, and fr mst purpses it is necessary t smth it. Several methds have been suggested; we shall use Bartlett's frmula Inset-2 M-l l(w} 2 (1 -) c. cs (26) -M+ 1

12 " '--''''' i6i N.Z. JOURNAL OF GEOLOGY AND GEOPttYStCS VOL. 9 where M < N is chsen with regard t the degree f smthing required. This frm may be regarded as (apprximately) the result f averaging the values f I (w) in the vicinity f w, r alternatively (again apprximately) f dividing the ttal recrd f data in t N / M sectins f equal length, and averaging the resulting peridgrams. The use f spectral methds in the analysis f stchastic pint prcesses has been discussed recently by Bartlett (1%3). Bartlett uses the Furier transfrm f the cvariance density <I>(w) = J: C(t)e iwt dt +- m (27) as the basis f his discussin, where C(t) and m are defined as in sectin 2'0, and we may take C(O) =. f (w) is related t <I>(w) by T T sin 2(1w ) (28) where T -00 is the fixed interval length assciated with f (w). Fr the trigger mdel f the previus sectin, 1 foo 12 <I>(w) = f1j1 + p,(a + 2 v) A(t)eiwt dt I (29) I 0, The idea f a dubly stchastic prcess mtined in Bartlett's paper is als relevant. Such a prcess may be described as a Pissn prcess, whse parameter is itself a randm functin, say t(t).* The trigge,r mdels are f this type, where (t) is a "saw-tth" functin, whse peaks are lcated at the trigger events. The decay frm each peak is determined by the functin '\(t). Bartlett estimates <I> ( w) directly by the statistic T (30) where t. is the time f the sth event, and T is the ttal time cnsidered. Fr ur purpses f (w) has the advantage f the finite length, and is mre T simply estimated, but suffers frm "aliasing" (e.g., a peak at w in f (w) culd be due t peaks at any f x =,/T, x = (2n7T + W)/T, x = (2n7T - w)/r in <I>(x), althugh nrmally nly the first fur terms r s need be cnsidered). In practice the difference between Equatins (25) and (30) amunts t a difference in the chice f r. Perhaps the mst relevant factr here is the time-scale f the effects we wish t study. If we are interested in the *It is interesting t nte that the pectrum f HI) is just 4>«(0) - m (d. 27). T

13 N.3 VERE-JONES & DAVIES - EARTHQUAKE ANALYSES 263 fine structure f the prcess-that is t say, in phenmena having a timescale f the same rder as the typical interval between events-then we shuld chse a value f 7 which is small with respect t this interval, and such a: chice crrespnds t the statistic (Equatin 30). On the ther hand, little infrmatin will be lst by taking a lnger value f T if we are interested in phenmena having a larger time-scale. The aliasing effect remains as a disadvantage f chsing intervals f equal length, but we shall nt enter int a discussin f ther pssibilities. Tw analyses will be described belw. Fr the first f these we have taken 7 = 0 1 year, and fr the secnd 7 = 1 day. The first analysis was made with the primary aim f investigating a pssible cyclic effect with a perid f tw t three years. It als yields infrmatin abut the cvariance structure. f the prces ver re.iatively lng perids. The secnd analysis yields cnstderably mre mfrmattn abut the fine structure f the prcess. In this case f «(0) is a reasnable apprximatin t <p(w). T. Te smthed peridgrams resulting frm the first analysis are shw m Ftgs. 2-5, and thse frm the secnd analysis in Figs Fr the first set N = 200 and M = 40 (in the ntatin f Equatin 26), and fr the secnd set N = 7,304 and M = 200. We have als marked in the values f the first few (unsmthed) periqgram rdinates; the value CO(7), which gives the mean height f the peridgram; theretical spectra frm several prpsed mdels; and apprximate significance levels against the hypthesis f independent x t with variance (12 = C( T) These levels refer t an verall maximum fr the peridgram, and nt t a peak at a particular value chsen in advance. The marked value f 20% is a cnservative estimate, btained by regarding the smthed peridgram as the average f N!M independent cmpnents. This appears t give a rather crude apprximatin, particularly when N!M is large, and a value f 5% (based n the frmulae in Hannan, 1960, sectin 3.3) might be mre apprpriate fr the ne-day spectra. In any case, the significance levels fr a prcess with such cmplex structure are likely t be very apprximate, and shuld be taken nly as a guide. Fr the interpretatin f theseperidgrams it is helpful t bear in mind the frms f the spectra fr the simple and cmpund Pissn prcesses mentined previusly. Fr the simple Pissn prcess with mean rate,\= p.a, f (w) =.p.th (31) 'T whereas fr the cmpund prcess t which ur trigger mdel apprximates fr lng intervals, f (w) =p.(a + a2 + v)-r (32) 'T Hence fr (0 crrespnding t perids lng cmpared with the lengths f grups f shcks (w 0) we might expect Equatin (32) t hld, and fr shrt perids «(O 71") Equatin (31) t hld. Is is apparent that this gives th general frm f the peridgrams-,..a regin f high values near the

14 264 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSiCS VOL. 9 rigin, subsequently.decreasing twards a value crrespnding t tl:le mean rte. We shall cnsier the peridgrs in mre detail in the ensuing sectins, wh.er they will be cmpared with the spectra f sme sme trigger mdels; It IS dear at nce, hwever, that n prnunced cyclic effects are present. In cnclusin it shuld be nted that the variance curve, the c<>variances, a!ld the spectrum are merely different ways f expressing the same infrmatin, the spectrum having the advantages f simpler interpretatin, and f transfrming simply under simple transfrmatins f the basic prcess TRIGGER MODEL WITH EXPONENTIAL DECAY If the decay functin has the expnential frm (Equatin 22), the cvariance density is als expnential, and frmulae (Equatin 18) and (Equatin 19) becme CO(T) -ppt + p.(az + V)T[l - (1 - e-pt)/pt} (33) C,,(T) - p.(a 2 + v)t e- npt (csh pt - l)/pt (34) The spectrum takes the frm cs w - f «1) = C(T) + 2C1 (T) (35) T 1 - e-2pt - 2e-PT cs (1) e-pt 2p. sin2 (ia» sinh (pt) p.(at + a 2 + V)T - - (a2 + v) (35a) P csb(p1') - cs w In rder t estimate the parameters f the mdel we shall first suppse that the mean rate (!JIlT) is determined directly as the ttal number f shcks divided by the ttal time. The remaining parameters will be estimated by using the methd f least uares t fit a functin f the frm (Equatin 35) t the peridgram rdlllates I j subject t the cnditin These methds are a cmprmise between the demands f technique and cnvenience: The peridgram rdinates, althugh apprximately independent, are expnentially rather than nrmally distributed, while least-squares methds have ptimal prperties fr nrmal variates. A general thery f maximum-likelihd estimatin f. the spectrum has been develped by Whittle and thers (e.g., Whittle, 1952), but their frmulae are nt easy t apply. With as many as 100 rdinates it seems likely that any tw reasnable methds f estimatin will lead t similarestimates--althugh the leastsquares methd will tend t give mre weight t pints fr which few) is large, than wuld Whittle's methd.

15 N. 3 VERE-JONES & DAVIES - EARnIQUAKE ANALYSES 265 The least squares equatins lead, after sme reductin, t the frmulae and ----lgy A(y). where 'I = e-fj'f, A(y) = lg Y -(1 - 'I)' B(y) = 1/y- 2 lgy - 'I. N-l the dash dentes differentiatin, and PN('I) is the plynmial an 'In-I 1 with cefficients an = (1 - n/n)cn + (n/n)c N -... Since 'I is f the rder f 0'2 these equatins are readily slved by numerical methds, the plynmials n the L.H.S. being largely determined by their first tw r three terms. Althugh the spectrum deeends n bth a and 11, as well as n 'I' there are nly tw parameters in these equatins, 'I and the factrp.(a2 + 11)7. The parameters a and 11 cannt be separated Wlless sme particular functinal frm is assumed fr the distributin fwlcti F(a). This pint will be discussed in mre detail in a later section. We have shwn in Table 3 the three parameters piit'l.j-(tp + li)t and 'I = rfj'f, tgether with the half-life (descrlbed belw), and the rati R = (a + tp + 1I)/a. Thes are shwn fr the ttal, shallw, and deep shcks. The crrespnding spectra are shown by the dtted lines in Figs. 2-4 and 6-8. T assess the "gdness. f fit" f the mdel, we first transfrm the peridgram by frming Ii'..:... Ii/Ii where Ii is the fitted- spectrum at III = 27rJ/N. We then examine the fwlctin This functin fr the shallw shcks is shwn in Fig. 10. The maximum deviatin f this functin frm the line y = x will,. the hypthesis that theel is crrect, fllw the Klmgrv-Smirnv law (e.g., Hannail, 1960). On this basis 20% and 5% significance levels are shwn n the diagram. This test des nt indicate disagreement with the mdel. Hwever in n case is the agreement very gd near the rigin (Figs. 2-5), r with the ne-day spectra.. The behaviur f the spectrum near the

16 266 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL. 9 rigin has a reciprcal relatin with the behaviur f the cvariance functin fr large values f t; a sharp peak at the rigin crrespnds t a slw rate f decay f the cvariance density and vice versa. The expnential spectrum in Figs. 2-5 is t flat near the rigin, crrespnding t the fact that the bserved wvariances fall ff much mre slwly than expnentially. Indeed, the values f C" calculated n the basis f Equatin (34) are negligable fr n greater than abut seven, while -several f the bserved values, fr example C 10, are quite large. It is hard t assess the significance f these later values, but it seems likely that they culd apprach the 50/0 level, if the expnential mdel was taken as the null hypthesis. Similarly, it can be seen frm Fig. 1 that the expnential mdel underestimates V (T) fr large values f T. The disagreement with the ne-day spectra is f the ppsite type. Fr the shrt perid, the expnential functin is t fiat. This effect is partly due t the fact that we have nrmalised ur mdel.s that CO(T) = CO(T), when T = 0'1 year. Hwever, these results nly enhance the general cnclusin that the expnential mdel cannt match the bserved prcess ver an extended range. Fr shrt perids the expnential functin rises at the rigin t slwly; fr lnger perids itdecays t fast. T summarise, it appears that the expnential mdel des nt affrd a gd descriptin f the data, althugh the disagreement is nt sufficiently prnunced fr the mdel t be rejected by the Klmgrv-Smirnv test applied t the 0'1 year spectra. These results have been btained n the basis f the data fr the shallw shcks, but they are supprted als fr the data fr the deep shcks. Indeed, the peak at the rigin in Fig. 3 is even sharper than that in Fig. 2, indicating that fr large T the cvariance density decays even mre slwly fr the deep shcks than it des fr the shallw shcks. This remark pints t an imprtant qualitative difference between the deep and shallw shcks, which is cnfirmed by the results set ut in Table 3. A "half life" f a grup f shcks may be defined as the time after which, n the average, half the shcks f the grup have ccurred. With the expnential mdel this equals (lge 2) / p. The higher value f the half-life fr the deep shcks suggests that fr these the time scale is lnger than fr the shallw shcks. TABLE 3-Parameters fr Trigger Mdel with Expnential Decay IJ4T p.(a' + V)T Y = r=pi R i-life (days) All shcks 5' '14 3'2 12') Shallw shcks 3'2 10'5 0'11 4'3 11'5 Deep shcks 2'5 1' '8 16'

17 N.3 VERE-JONES & DAVIES - EARTHQUAKE ANALYSES WESTPORT EARTHQUAKE SEQUENCE-1962 S far we have nt islated any particular earthquake sequence fr examinatin, but analysed the list as a whle. In this sectin we cnsider a particular earthquake sequence that ccurred near Westprt, New Zealand, in 1962, in rder t btain sme justificatin for ur next chice f decay functin. Our data is frm Adams and Le Frt (1963). The cumuliative graph,f the number f shcks f magnitude greater than r equal t three, against time, is shwn in Fig. 11, tgether with a graph f the functin N=A[l - ()P ] with A = 100, p = 0'25, and c = 0'25 days. This suggests a decay functin f the frm.\(t) = pcp/(c + 1)1 +p (23) might give abetter fit than the expnential functin. It shuld be nted that we have included here earthquakes f a lwer magnitude than were were used in the verall study, and that we have been cncerned with a relatively shrt perid. The functin (Equatin 23) decays extremely slwly; even at the end f the perid shwn, the cumulative graph has reached little mre than 75% f its asympttic value. The use f such a functin in relatin t after-shck sequences is by n means new,and in particular it has been the subject f an extensive study by Utsu (1961) f after-shck sequences in Japan. Whereas Utsu was largely cncerned with the analysis f individual sequences, in the sectin belw we shall examine hw far the same hypthesis is capable f explaining verall features f the data, as reflected in the peridgrams f Figs As a final remark, it is wrth pinting ut that in Fig. 11 there is a tendency fr the shcks t ccur in grups within the aftershck sequence itself. 2'6-TRIGGER MODEL WITH INVERSE POWER-LAW DECAY where k is numerical cnstant, apprximately equal t 0'9 when p = 0'25. The cvariance's 'T ) 1-p -(Hp) Cn(T),..., p.(a 2 + v)cp (. -c- n (37)

18 268 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL. 9 f.r nand () large. There des nt apear tg be a simple fgrmula fr the spectrum, but the ngn-aliased fgrm <I>(x) = p.a + p.(a2 + v) exp (- ip7r/2) r(1 -,p) (ex)p,..., pix + p.(a! + v)p2/(ex)2 fgr ex large. 00 (ex)i +p-- (38) (j _ p) j! The main feature.f <I>(x) and hence f (x) is the very sharp peak at the T.rigin-see Fig. 3, This cannt be reprduced n a smthed perid.gram f.r tw reasns, Firstly, frm EquatiGn (25a), the effect.f centring the.bservatins at the.bserved mean Xt is tg reduce the value.f the peridgram at the rigin tg zer, Remval.f the linear trend will have a similar effect. Secndly, the smthing itself will tend t.. remgve any sharp peaks. The unsmthed I; shuld ngt be greatly affected by this, and S. we have sh.wn IclIO in Figs These values have the cmpensating disadvantage that. their standard deviatin is f the same.rder as their expected value, SG that little significance can be placed.n their psiti.ns as individuals. Estimating the parameters.f the mdel is again difficult. The previgus article suggests taking p = 0'25. In this case the variance can be expressed in terms.f elliptic integrals, VeT) = CO(T) =!JIlT + p.(a 2 + V)T (1 - I(q) ] (39) 4 where I(q) = - 2( 2a )1 [2P(ep,1T/4) - B(ep, 7r/4) } 3 1 a 4 4 (1 - a)a 4 + and a =.COS ep = {e/(e + t)}i, and E (ep, 1T/4), P(ep,7I'/4) are elliptic integrals.f the first and secgnd kinds, respectively. If we then equatep.at X C = C C1 = c1 we find fgr the shallgw shcks, p.(a 2 + v)t = 31'1, e = 2'3 days, and f.r the deep shcks, p.(# + V)T = 6'2, e = 3'4 days. The c.rresp.nding "half-lives" are 35 days and 51 days respectively, althugh fgr such pwerlaw mdels the "half-life" gives a pogr indicati.n.f the time range, and the dependence remains appreciable fr several years.

19 N.3 VERE-JONES & DAVIES - EARTHQUAKE ANALYSES 269 These values f c are cnsiderably larger than the value fr the Westp()ft sequence, and fr this reasn a secnd pwer-lllw mdel fr the shallw shcks was taken, with c = 0'37 days. Fr this mdel we have set pfr = x and (in rder t secure a gd fit with the ne-day spectrum ) /L(d? + v)t = 18'3. The "half-life" crrespnding t his value f c is 5.2 days.. The spectra fr the first mdels are shwn by the dashed lines in Figs. 3 and 6 (shallw shcks), and Figs. 4 and 7 (deep shcks). The secnd pwer-law mdel fr the shallw shcks is shwn by the half-dashed line in Figs. 3a.nd 6. The crrespnding graphs f V('T)/'T are shwn fr the shallw shcks in Fig. 1. Taking the shallw shcks first, we see that bth the pwer-law spectra fit the 0'1 year peridgrams reasnably well. As cmpared with the expnential mdel, there is an imprvement f fit near the rigin, but the verall shape f the graph is nt greatly altered. The Klmgrv-Smirnv test again indicates n significant disagreement with the mdel. The results shwn in Fig. 1 and the ne-day peridgram are mre interesting. Frm Fig. 1 it can be seen that the pwer law with the larger value f c gives a better picture f the grwth f the cvariance fr mderate. values f 'T. On the ther hand, the behaviur f the prcess fr small values f 'T appears t be better described by the pwer law with the smaller value f c. This is shwn up particularly well by the ne"day spectra (Fig. 6). Again it appears that n smgle mdel will give a satisfactry descriptin f the prcess ver the- whle range f values f 'T. The ne-day spectra shw further discrepancies with the bserved peridgram near W = 7/'. This effect shws up mst strngly in the scaled peridgram I/(w) ::::: J(w)/!(IJl), which is shwn in Fig. 12. The peak at 7T indicates a marked deviatin frm the behaviur t be expected frm any f the simple trigger mdels. It is significant at abut the 21% level.* The interpretatin f this peak will be taken up in Sectin 3. Figs. 4 and 7 shw that fr the deep shcks the pwer-law mdel gives a reasnable fit bth with the 0'1 year and the ne-day peridgrams. A further mdel (with a "blck decay" functin) will be described in sectin 3. The deep shcks shw n sign f the peaks near 7T n the ne-day peridgram. 2.7-ENERGY SPECTRA It might be thughtpsible t use a peridgram calculated frm the energy released by earthquakes in successive perids. In fact it is unsatisfac *Care has t be taken in assessing the significance f this peak, since the peridgram rdinates at 0, 7T are (tughly speaking) based n distributins with half the usual number f degrees f fr«<lm. A further difficulty is that the theretical spectrum f(",) has nt been scaled s that it enclses the same area as the bserved peridgram (i.e., the variances have nt been exactly matched). Hwever, fr any simple mdel, the spectrum shuld apprach the value pa as '" appraches 7T, and. rescaling (by altering the value fp.(ti' + V)T) will mdify the spectrum near", = 0, nt near '" = 7T. Thus the significance f the peak at '" = 7T will nt he greatly affected by such changes.

20 270 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL. 9 try t d this directly because f the very large fluctuatins in the energy released. The peridgram cnveys very little infrmatin since it is almst entirely determined by the few earthquakes f large magnitude; thus n gruping effects are apparent in such a study. The lgarithm f the energy in each interval (plus a cnstant t keep the lgarithm finite when n earthquakes are recrded) is better behaved, and its peridgram shws evidence f sme cvariance structure. Other pssibilities (fractinal pwers f the energy released) might als be cnsidered, but their. physical interpretatin is nt clear. Sme recent studies (Keilis-BrOk and Malinvskaya, 1964) have suggested that this might prve a fruitful field fr further investigatin. We shall nt be able t cnsider it further in the present paper. EVALUATION OF THE RESULTS 3.1-SHALLOW SHOCKS N significant peridic effects are apparent, at least within the ranges cvered by the 0'1 year and ne-day spectra (say three mnths t three years and 3-15 days), The gruping effect is certainly significant. This is shwn bth bya direct examinatin f the cvariances and by the spectra. Thus the simple Pissn mdel cannt be accepted as an adequate descriptin f earthquake ccurrence amng the shallw shcks. Nne f the prpsed mdels seems t give a fully satisfactry accunt f the gruping effect. There are disagreements bth fr very lng perids and fr very shrt perids. The cntinued grwth f the variance fr large values f T tells heavily aga:inst the expnential mdel, and suggests that the dependence amng shcks may last fr many mnths. The first pwerlaw mdel (c = 2'6) gives a reasnable fit fr mderate values f T, but even this mdel appears t underestimate the variance when.,. is very large. This fact pints t the pssibility that sme f the basic assumptins under lying the mdel may be at fault. The fluctuatins in activity. frm year t year may be t large t be represented in terms f any simple statinary mdel, even after the remval f linear trends. Alternatively we may be faced with a prcess whse variance diverges. The validity f any secnd. rder analysis wuld have t be examined carefully in this case. The disagreements fr shrt perids are als very interesting. A feature f the cvariances e" fr the ne-day analysis is that they drp away quite slwly frm a very high initial value C. Fr this reasn the expnential and first pwer.law mdels, which give a reasnable fit fr the 0'1 year spectra, are here quite inapprpriate. A much better fit is given by the secnd pwerlaw, which has a smaller time scale (the parameter c being reduced by a factr f 6-7). On. the ther hand this mdel agrees less well with the behaviur f the prcess fr larger values f T, particularly as regards the variance curve shwn in Fig. 1. A secnd remarkable feature f the ne-day analysis is the peak in the peridgram at <t) = 7r. This peak may be partly a genuine peridic effect, but it seems rather t mark a failure f the peridgram t decay t the mean rate as quickly as the theretical spectra. This is mst likely an indica

21 N.3 VERE-JONES & DAVIES - EARTIIQUAKE ANALYSES 271 tin f a secndary gruping effect within the main prcess. Such an eff wuld' have a characteristic time scale f the rder f ne t five days, as cmpared with characteristic time scale f ne t six mnths fr the main effect The recrds f after-shck sequences such as the Westprt sequence lend sme direct supprt t this hypthesis. All f these ideas need further investigatin. As was pinted ut in a previus sectin, the quantities p" a, and v cannt be separated by the present analysis, withut makmg sme further hyptheses abut the distributin f the number f shcks in a grup. Hwever, It wuld seem unwise t try and give a particular frm t this distributin, because f the extremely skew frm suggested by a direct examinatin f the grups (see Table 7 in the earlier paper). In fact the appeaqnce f this distributin raises again the unpleasant pssibility f infinite variance. Despite these difficulties it may' be f sme interest t cmpare different estimates f the rati R = p,(a + a2 + v)/p,a which can be interpreted as the rati f the mean-square grup-size t the mean grup-size. Fr the first pwer-law mdel R = 11, and fr the secnd R = 6. An estimate can als be btained directly frm the spectrum by using the first few peridgram rdinates t estimate p, (a + a2 + v) (Equatin 32). This methd suggests R = 8. These values are larger than wuld be btained frm the distributin f grup sizes given in the earlier paper, r frm the expnential mdel DEEP SHOCKS Again there is n evidence f significant peridic effects. Bth the 0'1 year and ne-day peridgrams shw marked differences frm the crrespnding peridgrams fr the shallw shcks. The 0'1 year spectrum has a high narrw peak at the rigin. This may be partly due t secndrder r higher trends which have nt been remved, but the cvariances alne leave little dubt that a significant gruping effect is present. Thus the Pissn mdel is rejected fr this class als. The high, narrw peak suggests a small dependence extending ver very lng perids. The. expnential and pwer-law mdels suggest that the- time scale fr the deep shcks may be anything frm ne and a half t three times lnger than fr the shallw shcks. As an alternative methd f btaining a rugh estimate f this time scale, we have fitted t the data a trigger mdel having a' decay functin f "blck" type, In this case A(t) lit (0 < t < T) therwise. <I> (x) pil + 4p,(a2 + v) sin2 (jxt)lx2ta

22 272 N.Z. JOURNAL OF GEOLOGY AND GEOPHYSICS VOL. 9 The parameters were estimated frm the height f the first tw peri.dgram rdinates, and the pint at which the spectrum drps dwn t Its lwer value ptj. The results are T - 4 years p.(a2 + V)T 12 R 6. The value f T may seem rather lng cmpared with thse given by the ther mdels, but it shuld perhaps be cmpared with the "75% life" rather than the "half-life", the frmer having the values 34 days and 2'4 years fr the expnential and pwer-law mdels, respectively. : The ne-day peridgram shws little evidence f the fine structure.which is present with the shallw earthquakes. The tw peaks crrespndmg t perids near three and fur days apprach significance, and may be due t a phenmenn similar t that which causes the peak at a perid f tw days fr the shallw shcks. Such an interpretatin wuld cnfrm t the idea f a dilated time scale at greater depths. 3.3-TTAL SHOCKS We will nt prpse a mdel here because f the different frms Of clustering, which apparently predminate at different depths. It is f interest t cmpare the sum f the shallw- and deep-shck peridgrams with that f the ttal shcks. If the shallw and deep shcks are independent the sum f their peridgrams shuld be equal t that f the tcal shcks. N significant difference shws up, althugh there is sme difference near the rigin. This may pint twards a crrelatin with a cnsiderable time delay, but a crss-crrelatin analysis wuld be needed befre any mre definite statement culd be made. 3.4-MDIFIED LIST OF SHALLOW SHOCKS As pinted ut earlier and in the previus paper, the Pissn mdel is still unsatisfactry here. This indicates that sme gruping is still present. The ne-day peridgram shws that the clse gruping has been largely remved, -althugh there are still sme signs f a peak near 7r crrespnding t a peri"d f tw days. Mre extensive and varied types f "purging" might prvide anther methd leading t infrmatin abut the gruping effect. 4-DISCUSSION OF METHOD AND CONCLUSIONS (a) Te. thery f stchastic pint prcesses prvides a suitable setting fr a statistical study f earthquake ccurrence. Useful infrmatin can be O?tl1ed. frm te cvariances, peridgrams, variance-time curves, and the distnbutins f mterval lengths between shcks. Fr a preliminary survey

23 N.3 VERE-JONES & DAVIES - EAR11IQUAKE ANALYSES 273 perhaps the mst useful tl is a peridgram analysis based n the numbers f shcks ccurring in successive intervals f equal length. The chice f interval-length shuld be related t the time scale f the effect being studied. (b) The methds have the limitatins f utilising nly the time crdinates f the shcks, and nly first- and secnd-rder mments. Infrmatin is lst by neglecting the ther crdinates and the higher-rder mments. In additin, the techniques used require that the first-, secnd, and third-rder. prperties f the prcess are finite, but the results s far btained suggest that the pssibility f divergent mments must be brne in mind. (c) The "cntagius prcesses" prvide a suitable classf mdels which can be applied t earthquake ccurrence. In particular, earthquake ccurrence canbe described in a general way by trigger mdels, where the (cnditinal) prbability f a shck ccurring at a time t after a trigger event is prprtinal t a decay functin.\.(t). The pwer-law frm fr.\.(t) gives a better general fit than the expnential frm, but the tests used did nt shw up any significant difference between them. Nne f the mdels discussed can be cnsidered fully satisfactry, as they fail t predict the cntinued grwth f the variance fr lng intervals, and the details f the shrt-term structure, (d) The New Zealand data shw n evidence f any genuine peridic effect with a perid f mderate length (say ne mnth t tw years). On the ther hand, each class f data shws strng evidence f clustering, and fr each class the Pissn mdel must be rejected in favur f sme mre cmplex prcess. The time-perid cvered by the analysis is barely lng enugh t give an accurate idea f the gruping effect, especially with the deep shcks, and certainly n effects invlving times cmparable with r greater than 20 years culd have been detected. (e) There are imprtant qualitative differences between the gruping effects fr the deep ana shallw shcks. Estimates f the quantity R (the rati f the mean-square grup-size t the mean grup-size) vary frm 4'3 t 11 fr the shallw shcks, and frm 1'8 t 6 fr the deep shcks, Mrever the, time scale f the gruping effect appears t be lnger fr the deep shcks than fr the shallw shcks, being f the rder ne t six mnths fr the shallw shcks and up t fur years fr the deep shcks. A secndary gruping effect with a perid f tw days appears in the shallw shcks, but is nt present in the deep shcks, althugh these shw sme signs f secndry effects with perids f three t fur days. ACKNOWLEDGMENTS We are again grateful t the members f the Applied Mathematics Divisin staff wh have helped with the calculatins, and t Dr F. F. Evisn and the staff f the Sesmlgical Observatry, Wellingtn, f their cmments nd cntinued interest in this wrk.. We shyld hk t thank especially, Mr G. E. Dlckmsn fr handling the prgrammlog asciated With the spectral analyses, and Dr H. Thmpsn fr his help 10 he preparation f tl;1e fil drafts f the manuscript. We are indebted t Prfessr Galsky fr cmmuoicahng his wrk n the statistics f Central Asian earthquakes.

, which yields. where z1. and z2

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

More information

AP Statistics Notes Unit Two: The Normal Distributions

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

More information

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

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

More information

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

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

More information

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

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

More information

MATHEMATICS SYLLABUS SECONDARY 5th YEAR

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

More information

Lead/Lag Compensator Frequency Domain Properties and Design Methods

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

More information

CHAPTER 4 DIAGNOSTICS FOR INFLUENTIAL OBSERVATIONS

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

More information

LHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers

LHS Mathematics Department Honors Pre-Calculus Final Exam 2002 Answers LHS Mathematics Department Hnrs Pre-alculus Final Eam nswers Part Shrt Prblems The table at the right gives the ppulatin f Massachusetts ver the past several decades Using an epnential mdel, predict the

More information

Physics 2010 Motion with Constant Acceleration Experiment 1

Physics 2010 Motion with Constant Acceleration Experiment 1 . Physics 00 Mtin with Cnstant Acceleratin Experiment In this lab, we will study the mtin f a glider as it accelerates dwnhill n a tilted air track. The glider is supprted ver the air track by a cushin

More information

Differentiation Applications 1: Related Rates

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

More information

ENSC Discrete Time Systems. Project Outline. Semester

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

More information

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

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

More information

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

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

More information

Determining the Accuracy of Modal Parameter Estimation Methods

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

More information

B. Definition of an exponential

B. Definition of an exponential Expnents and Lgarithms Chapter IV - Expnents and Lgarithms A. Intrductin Starting with additin and defining the ntatins fr subtractin, multiplicatin and divisin, we discvered negative numbers and fractins.

More information

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

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

More information

Preparation work for A2 Mathematics [2017]

Preparation work for A2 Mathematics [2017] Preparatin wrk fr A2 Mathematics [2017] The wrk studied in Y12 after the return frm study leave is frm the Cre 3 mdule f the A2 Mathematics curse. This wrk will nly be reviewed during Year 13, it will

More information

Hypothesis Tests for One Population Mean

Hypothesis Tests for One Population Mean Hypthesis Tests fr One Ppulatin Mean Chapter 9 Ala Abdelbaki Objective Objective: T estimate the value f ne ppulatin mean Inferential statistics using statistics in rder t estimate parameters We will be

More information

Computational modeling techniques

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

More information

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic.

Sections 15.1 to 15.12, 16.1 and 16.2 of the textbook (Robbins-Miller) cover the materials required for this topic. Tpic : AC Fundamentals, Sinusidal Wavefrm, and Phasrs Sectins 5. t 5., 6. and 6. f the textbk (Rbbins-Miller) cver the materials required fr this tpic.. Wavefrms in electrical systems are current r vltage

More information

Module 4: General Formulation of Electric Circuit Theory

Module 4: General Formulation of Electric Circuit Theory Mdule 4: General Frmulatin f Electric Circuit Thery 4. General Frmulatin f Electric Circuit Thery All electrmagnetic phenmena are described at a fundamental level by Maxwell's equatins and the assciated

More information

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

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

More information

Lab 1 The Scientific Method

Lab 1 The Scientific Method INTRODUCTION The fllwing labratry exercise is designed t give yu, the student, an pprtunity t explre unknwn systems, r universes, and hypthesize pssible rules which may gvern the behavir within them. Scientific

More information

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

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

More information

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

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

More information

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b

MODULE 1. e x + c. [You can t separate a demominator, but you can divide a single denominator into each numerator term] a + b a(a + b)+1 = a + b . REVIEW OF SOME BASIC ALGEBRA MODULE () Slving Equatins Yu shuld be able t slve fr x: a + b = c a d + e x + c and get x = e(ba +) b(c a) d(ba +) c Cmmn mistakes and strategies:. a b + c a b + a c, but

More information

Inference in the Multiple-Regression

Inference in the Multiple-Regression Sectin 5 Mdel Inference in the Multiple-Regressin Kinds f hypthesis tests in a multiple regressin There are several distinct kinds f hypthesis tests we can run in a multiple regressin. Suppse that amng

More information

Study Group Report: Plate-fin Heat Exchangers: AEA Technology

Study Group Report: Plate-fin Heat Exchangers: AEA Technology Study Grup Reprt: Plate-fin Heat Exchangers: AEA Technlgy The prblem under study cncerned the apparent discrepancy between a series f experiments using a plate fin heat exchanger and the classical thery

More information

NUMBERS, MATHEMATICS AND EQUATIONS

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

More information

Homology groups of disks with holes

Homology groups of disks with holes Hmlgy grups f disks with hles THEOREM. Let p 1,, p k } be a sequence f distinct pints in the interir unit disk D n where n 2, and suppse that fr all j the sets E j Int D n are clsed, pairwise disjint subdisks.

More information

Distributions, spatial statistics and a Bayesian perspective

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

More information

Dead-beat controller design

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

More information

Biplots in Practice MICHAEL GREENACRE. Professor of Statistics at the Pompeu Fabra University. Chapter 13 Offprint

Biplots in Practice MICHAEL GREENACRE. Professor of Statistics at the Pompeu Fabra University. Chapter 13 Offprint Biplts in Practice MICHAEL GREENACRE Prfessr f Statistics at the Pmpeu Fabra University Chapter 13 Offprint CASE STUDY BIOMEDICINE Cmparing Cancer Types Accrding t Gene Epressin Arrays First published:

More information

Math Foundations 20 Work Plan

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

More information

1 The limitations of Hartree Fock approximation

1 The limitations of Hartree Fock approximation Chapter: Pst-Hartree Fck Methds - I The limitatins f Hartree Fck apprximatin The n electrn single determinant Hartree Fck wave functin is the variatinal best amng all pssible n electrn single determinants

More information

ABSORPTION OF GAMMA RAYS

ABSORPTION OF GAMMA RAYS 6 Sep 11 Gamma.1 ABSORPTIO OF GAMMA RAYS Gamma rays is the name given t high energy electrmagnetic radiatin riginating frm nuclear energy level transitins. (Typical wavelength, frequency, and energy ranges

More information

A Matrix Representation of Panel Data

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

More information

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

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

More information

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

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

More information

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

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

More information

Thermodynamics Partial Outline of Topics

Thermodynamics Partial Outline of Topics Thermdynamics Partial Outline f Tpics I. The secnd law f thermdynamics addresses the issue f spntaneity and invlves a functin called entrpy (S): If a prcess is spntaneus, then Suniverse > 0 (2 nd Law!)

More information

BASD HIGH SCHOOL FORMAL LAB REPORT

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

More information

Interference is when two (or more) sets of waves meet and combine to produce a new pattern.

Interference is when two (or more) sets of waves meet and combine to produce a new pattern. Interference Interference is when tw (r mre) sets f waves meet and cmbine t prduce a new pattern. This pattern can vary depending n the riginal wave directin, wavelength, amplitude, etc. The tw mst extreme

More information

Chemistry 20 Lesson 11 Electronegativity, Polarity and Shapes

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

More information

OTHER USES OF THE ICRH COUPL ING CO IL. November 1975

OTHER USES OF THE ICRH COUPL ING CO IL. November 1975 OTHER USES OF THE ICRH COUPL ING CO IL J. C. Sprtt Nvember 1975 -I,," PLP 663 Plasma Studies University f Wiscnsin These PLP Reprts are infrmal and preliminary and as such may cntain errrs nt yet eliminated.

More information

Resampling Methods. Chapter 5. Chapter 5 1 / 52

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

More information

Department of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets

Department of Economics, University of California, Davis Ecn 200C Micro Theory Professor Giacomo Bonanno. Insurance Markets Department f Ecnmics, University f alifrnia, Davis Ecn 200 Micr Thery Prfessr Giacm Bnann Insurance Markets nsider an individual wh has an initial wealth f. ith sme prbability p he faces a lss f x (0

More information

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

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

More information

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

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

More information

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

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

More information

ENGI 4430 Parametric Vector Functions Page 2-01

ENGI 4430 Parametric Vector Functions Page 2-01 ENGI 4430 Parametric Vectr Functins Page -01. Parametric Vectr Functins (cntinued) Any nn-zer vectr r can be decmpsed int its magnitude r and its directin: r rrˆ, where r r 0 Tangent Vectr: dx dy dz dr

More information

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

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

More information

Computational modeling techniques

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

More information

How do scientists measure trees? What is DBH?

How do scientists measure trees? What is DBH? Hw d scientists measure trees? What is DBH? Purpse Students develp an understanding f tree size and hw scientists measure trees. Students bserve and measure tree ckies and explre the relatinship between

More information

Preparation work for A2 Mathematics [2018]

Preparation work for A2 Mathematics [2018] Preparatin wrk fr A Mathematics [018] The wrk studied in Y1 will frm the fundatins n which will build upn in Year 13. It will nly be reviewed during Year 13, it will nt be retaught. This is t allw time

More information

CHAPTER 2 Algebraic Expressions and Fundamental Operations

CHAPTER 2 Algebraic Expressions and Fundamental Operations CHAPTER Algebraic Expressins and Fundamental Operatins OBJECTIVES: 1. Algebraic Expressins. Terms. Degree. Gruping 5. Additin 6. Subtractin 7. Multiplicatin 8. Divisin Algebraic Expressin An algebraic

More information

Eric Klein and Ning Sa

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

More information

CHM112 Lab Graphing with Excel Grading Rubric

CHM112 Lab Graphing with Excel Grading Rubric Name CHM112 Lab Graphing with Excel Grading Rubric Criteria Pints pssible Pints earned Graphs crrectly pltted and adhere t all guidelines (including descriptive title, prperly frmatted axes, trendline

More information

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

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

More information

Experiment #3. Graphing with Excel

Experiment #3. Graphing with Excel Experiment #3. Graphing with Excel Study the "Graphing with Excel" instructins that have been prvided. Additinal help with learning t use Excel can be fund n several web sites, including http://www.ncsu.edu/labwrite/res/gt/gt-

More information

Admissibility Conditions and Asymptotic Behavior of Strongly Regular Graphs

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

More information

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

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

More information

A study on GPS PDOP and its impact on position error

A study on GPS PDOP and its impact on position error IndianJurnalfRadi& SpacePhysics V1.26,April1997,pp. 107-111 A study n GPS and its impact n psitin errr P Banerjee,AnindyaBse& B SMathur TimeandFrequencySectin,NatinalPhysicalLabratry,NewDelhi110012 Received19June

More information

Supporting information

Supporting information Electrnic Supplementary Material (ESI) fr Physical Chemistry Chemical Physics This jurnal is The wner Scieties 01 ydrgen perxide electrchemistry n platinum: twards understanding the xygen reductin reactin

More information

ChE 471: LECTURE 4 Fall 2003

ChE 471: LECTURE 4 Fall 2003 ChE 47: LECTURE 4 Fall 003 IDEL RECTORS One f the key gals f chemical reactin engineering is t quantify the relatinship between prductin rate, reactr size, reactin kinetics and selected perating cnditins.

More information

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

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

More information

Module 3: Gaussian Process Parameter Estimation, Prediction Uncertainty, and Diagnostics

Module 3: Gaussian Process Parameter Estimation, Prediction Uncertainty, and Diagnostics Mdule 3: Gaussian Prcess Parameter Estimatin, Predictin Uncertainty, and Diagnstics Jerme Sacks and William J Welch Natinal Institute f Statistical Sciences and University f British Clumbia Adapted frm

More information

Function notation & composite functions Factoring Dividing polynomials Remainder theorem & factor property

Function notation & composite functions Factoring Dividing polynomials Remainder theorem & factor property Functin ntatin & cmpsite functins Factring Dividing plynmials Remainder therem & factr prperty Can d s by gruping r by: Always lk fr a cmmn factr first 2 numbers that ADD t give yu middle term and MULTIPLY

More information

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

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

More information

Kinetic Model Completeness

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

More information

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

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

More information

A solution of certain Diophantine problems

A solution of certain Diophantine problems A slutin f certain Diphantine prblems Authr L. Euler* E7 Nvi Cmmentarii academiae scientiarum Petrplitanae 0, 1776, pp. 8-58 Opera Omnia: Series 1, Vlume 3, pp. 05-17 Reprinted in Cmmentat. arithm. 1,

More information

Name: Block: Date: Science 10: The Great Geyser Experiment A controlled experiment

Name: Block: Date: Science 10: The Great Geyser Experiment A controlled experiment Science 10: The Great Geyser Experiment A cntrlled experiment Yu will prduce a GEYSER by drpping Ments int a bttle f diet pp Sme questins t think abut are: What are yu ging t test? What are yu ging t measure?

More information

ECE 5318/6352 Antenna Engineering. Spring 2006 Dr. Stuart Long. Chapter 6. Part 7 Schelkunoff s Polynomial

ECE 5318/6352 Antenna Engineering. Spring 2006 Dr. Stuart Long. Chapter 6. Part 7 Schelkunoff s Polynomial ECE 538/635 Antenna Engineering Spring 006 Dr. Stuart Lng Chapter 6 Part 7 Schelkunff s Plynmial 7 Schelkunff s Plynmial Representatin (fr discrete arrays) AF( ψ ) N n 0 A n e jnψ N number f elements in

More information

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System

Flipping Physics Lecture Notes: Simple Harmonic Motion Introduction via a Horizontal Mass-Spring System Flipping Physics Lecture Ntes: Simple Harmnic Mtin Intrductin via a Hrizntal Mass-Spring System A Hrizntal Mass-Spring System is where a mass is attached t a spring, riented hrizntally, and then placed

More information

Chapter Summary. Mathematical Induction Strong Induction Recursive Definitions Structural Induction Recursive Algorithms

Chapter Summary. Mathematical Induction Strong Induction Recursive Definitions Structural Induction Recursive Algorithms Chapter 5 1 Chapter Summary Mathematical Inductin Strng Inductin Recursive Definitins Structural Inductin Recursive Algrithms Sectin 5.1 3 Sectin Summary Mathematical Inductin Examples f Prf by Mathematical

More information

Hubble s Law PHYS 1301

Hubble s Law PHYS 1301 1 PHYS 1301 Hubble s Law Why: The lab will verify Hubble s law fr the expansin f the universe which is ne f the imprtant cnsequences f general relativity. What: Frm measurements f the angular size and

More information

OF SIMPLY SUPPORTED PLYWOOD PLATES UNDER COMBINED EDGEWISE BENDING AND COMPRESSION

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

More information

Writing Guidelines. (Updated: November 25, 2009) Forwards

Writing Guidelines. (Updated: November 25, 2009) Forwards Writing Guidelines (Updated: Nvember 25, 2009) Frwards I have fund in my review f the manuscripts frm ur students and research assciates, as well as thse submitted t varius jurnals by thers that the majr

More information

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018

Physics 2B Chapter 23 Notes - Faraday s Law & Inductors Spring 2018 Michael Faraday lived in the Lndn area frm 1791 t 1867. He was 29 years ld when Hand Oersted, in 1820, accidentally discvered that electric current creates magnetic field. Thrugh empirical bservatin and

More information

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture

A Few Basic Facts About Isothermal Mass Transfer in a Binary Mixture Few asic Facts but Isthermal Mass Transfer in a inary Miture David Keffer Department f Chemical Engineering University f Tennessee first begun: pril 22, 2004 last updated: January 13, 2006 dkeffer@utk.edu

More information

PSU GISPOPSCI June 2011 Ordinary Least Squares & Spatial Linear Regression in GeoDa

PSU GISPOPSCI June 2011 Ordinary Least Squares & Spatial Linear Regression in GeoDa There are tw parts t this lab. The first is intended t demnstrate hw t request and interpret the spatial diagnstics f a standard OLS regressin mdel using GeDa. The diagnstics prvide infrmatin abut the

More information

ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS

ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS ENGINEERING COUNCIL CERTIFICATE LEVEL THERMODYNAMIC, FLUID AND PROCESS ENGINEERING C106 TUTORIAL 5 THE VISCOUS NATURE OF FLUIDS On cmpletin f this tutrial yu shuld be able t d the fllwing. Define viscsity

More information

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

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

More information

Dispersion Ref Feynman Vol-I, Ch-31

Dispersion Ref Feynman Vol-I, Ch-31 Dispersin Ref Feynman Vl-I, Ch-31 n () = 1 + q N q /m 2 2 2 0 i ( b/m) We have learned that the index f refractin is nt just a simple number, but a quantity that varies with the frequency f the light.

More information

This section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving.

This section is primarily focused on tools to aid us in finding roots/zeros/ -intercepts of polynomials. Essentially, our focus turns to solving. Sectin 3.2: Many f yu WILL need t watch the crrespnding vides fr this sectin n MyOpenMath! This sectin is primarily fcused n tls t aid us in finding rts/zers/ -intercepts f plynmials. Essentially, ur fcus

More information

General Chemistry II, Unit I: Study Guide (part I)

General Chemistry II, Unit I: Study Guide (part I) 1 General Chemistry II, Unit I: Study Guide (part I) CDS Chapter 14: Physical Prperties f Gases Observatin 1: Pressure- Vlume Measurements n Gases The spring f air is measured as pressure, defined as the

More information

AN INTERMITTENTLY USED SYSTEM WITH PREVENTIVE MAINTENANCE

AN INTERMITTENTLY USED SYSTEM WITH PREVENTIVE MAINTENANCE J. Operatins Research Sc. f Japan V!. 15, N. 2, June 1972. 1972 The Operatins Research Sciety f Japan AN INTERMITTENTLY USED SYSTEM WITH PREVENTIVE MAINTENANCE SHUNJI OSAKI University f Suthern Califrnia

More information

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

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

More information

Tree Structured Classifier

Tree Structured Classifier Tree Structured Classifier Reference: Classificatin and Regressin Trees by L. Breiman, J. H. Friedman, R. A. Olshen, and C. J. Stne, Chapman & Hall, 98. A Medical Eample (CART): Predict high risk patients

More information

FIELD QUALITY IN ACCELERATOR MAGNETS

FIELD QUALITY IN ACCELERATOR MAGNETS FIELD QUALITY IN ACCELERATOR MAGNETS S. Russenschuck CERN, 1211 Geneva 23, Switzerland Abstract The field quality in the supercnducting magnets is expressed in terms f the cefficients f the Furier series

More information

Performance Bounds for Detect and Avoid Signal Sensing

Performance Bounds for Detect and Avoid Signal Sensing Perfrmance unds fr Detect and Avid Signal Sensing Sam Reisenfeld Real-ime Infrmatin etwrks, University f echnlgy, Sydney, radway, SW 007, Australia samr@uts.edu.au Abstract Detect and Avid (DAA) is a Cgnitive

More information

arxiv:hep-ph/ v1 2 Jun 1995

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

More information

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function

Modeling the Nonlinear Rheological Behavior of Materials with a Hyper-Exponential Type Function www.ccsenet.rg/mer Mechanical Engineering Research Vl. 1, N. 1; December 011 Mdeling the Nnlinear Rhelgical Behavir f Materials with a Hyper-Expnential Type Functin Marc Delphin Mnsia Département de Physique,

More information

NOTE ON APPELL POLYNOMIALS

NOTE ON APPELL POLYNOMIALS NOTE ON APPELL POLYNOMIALS I. M. SHEFFER An interesting characterizatin f Appell plynmials by means f a Stieltjes integral has recently been given by Thrne. 1 We prpse t give a secnd such representatin,

More information

Lyapunov Stability Stability of Equilibrium Points

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

More information

Lecture 24: Flory-Huggins Theory

Lecture 24: Flory-Huggins Theory Lecture 24: 12.07.05 Flry-Huggins Thery Tday: LAST TIME...2 Lattice Mdels f Slutins...2 ENTROPY OF MIXING IN THE FLORY-HUGGINS MODEL...3 CONFIGURATIONS OF A SINGLE CHAIN...3 COUNTING CONFIGURATIONS FOR

More information

Application of ILIUM to the estimation of the T eff [Fe/H] pair from BP/RP

Application of ILIUM to the estimation of the T eff [Fe/H] pair from BP/RP Applicatin f ILIUM t the estimatin f the T eff [Fe/H] pair frm BP/RP prepared by: apprved by: reference: issue: 1 revisin: 1 date: 2009-02-10 status: Issued Cryn A.L. Bailer-Jnes Max Planck Institute fr

More information

7 TH GRADE MATH STANDARDS

7 TH GRADE MATH STANDARDS ALGEBRA STANDARDS Gal 1: Students will use the language f algebra t explre, describe, represent, and analyze number expressins and relatins 7 TH GRADE MATH STANDARDS 7.M.1.1: (Cmprehensin) Select, use,

More information