Choice of Thresholds. for Wavelet Shrinkage Estimate of the Spectrum. Hong-ye Gao. MathSoft, Inc Westlake Ave. N., Suite 500.
|
|
- Rosalind Edwards
- 5 years ago
- Views:
Transcription
1 Choice of Thresholds for Wavelet Shrinkage Estimate of the Spectrum Hong-ye Gao Statistical Sciences Division MathSoft, Inc Westlake Ave. N., Suite 500 Seattle, WA Abstract We study the problem of estimating the log spectrum of a stationary Gaussian time series by thresholding the empirical wavelet coecients. We propose the use of thresholds tj;n depending on sample size n, wavelet basis and resolution level j. At ne resolution levels (j = 1; 2; :::), we propose tj;n = j log n; where fjg are level-dependent constants and at coarse levels (j 1) tj;n = p 3 p log n: The purpose of this thresholding level is to make the reconstructed log-spectrum as nearly noise-free as possible. In addition to being pleasant from a visual point of view, the noise-free character leads to attractive theoretical properties over a wide range of smoothness assumptions. Previous proposals set much smaller thresholds and did not enjoy these properties. Key Words and Phrases: Log Spectrum Estimation; Orthogonal Wavelet Transformation; Shrinkage Estimator. Acknowledgements: This work was supported by NASA contract NCA2-488, by NSF DMS and by NSF DMS This is part of the author's Ph.D dissertation, under the supervision of Professor David L. Donoho, whose generous guidance and encouragement are gratefully acknowledged and appreciated.
2 1 Introduction Suppose we want to study a time series from an observed segment of the series X 1 ; :::; X N. There are two aspects to the study of time series analysis and modelling. The aim of analysis is to understand the nature of the process. The main reason for modelling a time series is to enable forecasts of future values of the process. No forecasting can be made before we understand the salient features of the process being characterized. The analysis of a time series can be done either in the time domain or in the frequency domain. In the time domain, attention is focused on the relationship between observations at dierent points in time, e.g. ARMA models; while in the frequency domain it is cyclical movements which are studied and this can be done by studying the spectral density h. The two forms of analysis are complementary rather than competitive. They give dierent insights into the nature of the process. Our goal in this paper is to estimate h from the data X 1 ; :::; X N. Here h(!) = 1 2 1X s=01 (s) cos(s!): The techniques currently used in spectral analysis are: Window methods (also known as kernel methods), x6.2.3 of Priestley (1981). A spectral window (also called weight function) is used to smooth the periodogram, or equivalently a lag window (also known as weight sequences, Priestley (1981), pp. 437, covariance window, Bentkus and Susinskas (1982)) W (1) is applied to the sample auto-covariances. The resulting estimate is: ^h(!) = 1 2 X jsj<n W (s)^(s) cos(s!) Typical choices of windows include the Fejer window: W (s) = I [jsjm ] for some 1 M N, for example M = p N; Bartlett window: W (s) = (1 0 jsj=m) + ; Daniell window, W (s) = sin(s=m)=(s=m); see x6.2.3 of Priestley (1981) for more choices of windows. Bentkus and Susinskas (1982) has studied optimal L 2 convergence rate for the window method over certain classes of smooth spectra. Autoregressive spectral estimation (AR approximation), Parzen (1974). A high order autoregressive model y t = 1 y t p y t0p + " t is t and the corresponding spectrum, a rational function, ^h(!) = ^ 2 2j1 0 ^ 1 z ^ p z p j 2 is taken as an approximation to f. 2
3 Mixtures of the above, e.g. the prewhitening technique, x7.4.1 of Priestley (1981). Maximum entropy (ME) methods, introduced by Burg (1967, 1972). This method is essentially equivalent to AR method. ARMA approximation: similar idea as AR approximation. They perform well in many cases and have been used frequently in practice. They have both advantages and disadvantages. For example, window methods have computational advantages but they may perform poorly in cases of high dynamic range; AR approximation is suitable only for relatively smooth spectra, Tukey (1978); ARMA procedures are almost always eective and furnishes attractive interpretations for the derived model; but it is computationally unrealistic to t high-order ARMA models, and rigorous theory of asymptotic properties is lacking. In this paper, we will develop a technique, based on wavelet decomposition of the periodogram and reconstruction of the spectrum; our new approach avoids shortcomings of earlier methods. Our new approach is computationally ecient; it can estimate spectra which are nonsmooth at a near-optimal rate; and it can be rigorously analyzed. It is well known that the asymptotic variances of the empirical wavelet coecients are functions of the unknown spectrum, and that non-linear Wavelet Shrinkage procedure below heavily depends on the asymptotic variances. Since the asymptotic variances are \proportional" to the coecients, i.e., the bigger the coecients at location (j; k), the higher the noise level at this location. This makes estimation more dicult in practice. In order to overcome the heteroscedasticity, we therefore consider estimating the log-spectrum. Suppose we observe a segment X 0 ; X 1 ; :::; X 2n01 of a Gaussian time series with mean zero and spectral density h, where for simplicity we assume that n = 2 m is dyadic. In this paper we discuss the following 4-step wavelet shrinkage procedure for estimation of log spectrum, g = log h: [1] Calculate the log-periodogram z l = log I(! l ); l = 0; : : : ; n 0 1; where! l = l=2n and X I(!) = 1 2n01 4n j X t e 0it! j 2 (1) t=0 [2] Take a standard periodic wavelet transform of (z l ) to get the empirical wavelet coef- cients fy j;k g j=1;2;:::m01;k=0;:::;n=2j 01. [3] Apply the soft threshold t (x) = sgn(x)(jxj 0 t) + (2) 3
4 to the empirical wavelet coecients fy j;k g, with level-dependent thresholds t = t j;n according to formula t j;n = p 3 p log n _ j log n (3) where ( 1 ; :::; 10 ) = (1:29; 1:09; 0:92; 0:77; 0:65; 0:54; 0:46; 0:39; 0:32; 0:27) for the commonly used compactly supported orthogonal wavelet bases (coiets, daublets and symmlets, Chapters 6 and 8 of Daubechies, 1992) and some selected n. [4] Invert the wavelet transform, producing an estimate (^g l 3 ) of the log-spectrum at the Fourier frequencies! l. This procedure falls in the general category of wavelet shrinkage estimates for noisy data. Donoho and Johnstone (1992a,b,c), have discussed the application of such thresholded wavelet transforms for recovering curves from noisy data, where the noise is assumed to have a Gaussian distribution, and they have proposed the level-independent threshold t n = p 2 log(n), where is the noise level. Such methods have also been developed in density estimation: Johnstone, Kerkyacharian and Picard (1992) propose the use of level-dependent thresholds p j. In general such wavelet shrinkage methods have a number of theoretical advantages, including near-optimal mean-squared error and near-ideal spatial adaptation. Recently, the author and others have studied the possibility of using wavelet shrinkage to de-noise periodogram data. Application of wavelet shrinkage to the log-periodogram is particularly attractive, Moulin (1993a, b), since the logarithm is the variance stabilizing transformation for the periodogram, Wahba (1980). It is easy to apply wavelet shrinkage software designed for Gaussian noise to the log periodogram, and in some cases acceptable reconstructions have been obtained. However, this will not always be the case. The noise in the wavelet coecients of the log-periodogram has a non-gaussian character; it reaches high levels somewhat more frequently than Gaussian theory would predict. Consequently, thresholds set based on Gaussian theory will not be high enough to completely suppress noise in the coecients. The thresholding we propose is based on a careful analysis of the non-gaussian character, and is somewhat larger than Gaussian theory would predict. As an example of the dierence between our method and other proposals, we present in Figure 1 a display for an AR(24) time series: (a) its log spectrum; (b) log periodogram; and wavelet reconstructions based on (c) our proposal and (d) based on Gaussian theory. It is visually evident that the spectral estimate based on Gaussian theory has spurious noise spikes, which are unrelated to the true underlying spectrum. More plots for this example, along with other examples, for dierent sample sizes and dierent wavelet bases, are given at the end of the paper. The pattern visually evident in this example is conrmed by two theorems which we prove in this paper. 4
5 Theorem 1 Under the Wahba approximation for log periodograms (see section 2.2 below), as n! 1, P f [ j [sup jy j;k 0 Ey j;k j > t j;n ]g! 0: (4) k In words, the thresholds are set high enough so that the noise does not surpass them. Theorem 2 Under the Wahba approximation for log periodograms and a compactly supported wavelet basis, at the nest level j =1, as n! 1, Similar results hold at levels j =2; 3; :::. p P fsup jy 1;k 0 Ey 1;k j > 2 log ng! 1: (5) k In words: noise spikes are almost surely to be absent for our proposal and almost surely to be present in a proposal based on Gaussian noise. The signicance of these noise spikes and the need to de-noise is more than cosmetic. In a recent article, Donoho, Johnstone, Kerkyacharian and Picard (1992) study the Gaussian white noise model, and show that by exploiting a noise-free property like (4) one can show that the estimator attains near-optimal reconstruction in a wide variety of norms simultaneously over a wide variety of smoothness assumptions. Their arguments are general and abstract and transfer to the present setting after one has established the de-noising property (4) If noise-free property like (4) is not maintained, then an inspection of their arguments will show that the resulting estimate does not achieve near optimal reconstruction in norms which measure smoothness of the reconstruction, although it may well be true that the error in `2 norm, which does not measure smoothness, is still acceptable. 2 The Choice of Thresholds Our choice of the thresholds (t j;n ) depends on the tail behavior of the empirical coecients fy j;k g. Let p j (t) = max P fjy j;k 0 Ey j;k j > tg (6) k be the tail probability of y j;k, where wavelet coecient y j;k is a standardized linear combination of the log periodogram ordinates. Then we aim to establish that n p j (t j;n )! 0; n! 1; (7) uniformly for j = 1; 2; :::; m 0 1 as this implies (4) in Theorem 1 (see (15). In this section, we would like to give two heuristic arguments, upon which our proposed thresholds (3) are based. 5
6 2.1 Normal Approximation For normally distributed wavelet coecients where " j;k N(0; 2 ) are iid, the threshold w j;k = j;k + " j;k t j;n t n = p 2 log n enjoys a variety of nice theoretical properties, Donoho and Johnstone (1992a, b, c). Under certain regularity conditions, for j close to m = [log 2 n] (technically j! 1 as m! 1), as n! 1, individual empirical wavelet coecients will be asymptotically normal: w j;k 0 Ew j;k N(0; 2 =6) this follows by applying results of Taniguchi (1979, 1980). Moulin (1993a, b) makes similar argument. So when the sample size n is large, for coarse levels (j close to m), y j;0 ; :::; y j;2j 01 are approximately normally distributed with the same variance 2 =6. One may argue that p j ( qlog(n)=3) 1 + o(1) n p ; m > j! 1: (8) 2 log n Therefore (7) holds and this leads to the rst part of our thresholds. Here \regularity" refers to the length of memory and the asymptotic normality for the empirical coecients not only depends on n (large) and j (large), but also the regularity of the time series. In general, the longer the memory, the weaker the asymptotic normality. 2.2 Wahba Approximation For small j, the normal approximation deteriorates. Under a compactly supported wavelet basis (chapters 6 and 8 of Daubechies, 1992), for xed j, the coecients of the nest level y j;k 's are linear combinations of xed number of flog I(! l )g. For example, with the Haar wavelet, y 1;k is just the dierence of two adjacent log I(! l )'s. Consider the following non-gaussian additive noise model: z l = g l + " l (9) l = 1; 2; :::; n 0 1, where (g l ) is the object to be estimated (e.g. log spectrum at frequency! l ) and " i = log( i =2) + (10) where f i ; i = 1; 2; :::; n 0 1g are independently 2 2 = exp(1=2) distributed and = 0:57721 is the Euler-Mascheroni constant. It can be shown that E" i = 0 and var(" i ) = 2 =6. Wahba (1980) proposed this as a model for the log periodogram, i.e., z l log I(! l ) and 6
7 g l log f(! l ) 0. Note that I(! 0 ) 2 1. Since the inuence of this term is negligible for n large, we will ignore this term in our discussion. For circulant time series, see x4.3 of Harvey (1989), log I(! l ) follows model (9){(10) exactly. For general stationary processes, the above model is only asymptotically true (see Theorem of Brillinger, 1981). The exact distribution of log periodogram can be found in Wittwer (1986). Now we will show that under the model (9){(10), our advertised thresholds are indeed adequate for the ner levels. Note that under model (9){(10), each y j;k is a standardized linear combination of z i 's. Let us use the very nest level, j = 1, as an example. Similar results hold for levels j = 2; 3; :::. Suppose a compactly supported wavelet basis with L coecients in the dilation equation. For example, for the Haar wavelet, L=2; D4 wavelet, L=4, etc., more examples can be found in Chapters 6 and 8 of Daubechies (1992). Then y 1;k = LX l=1 a l z l+2k0s where P s is some shift parameter for computational purpose (e.g. s = L=2) and P a l = 0; a 2 l = 1. Inequalities below leads to the threshold (3). Let a = max 1lL ja l j, from the proofs of the Theorems (see (16) and (22) of Appendix), we can show that for t large, Hence, for n large, there is a constant A so that 0:25e 02t=a p 1 (t) (et=al) L e 0t=a : (11) 0:25 n 2=a p A(log n)l 1( log n) n =a : (12) Hence thresholds at the nest level must have the form: log n, for some constant > a, and then np 1 ( log(2n))! 0; n! 1: (13) Note that the wavelet lters obey P L l=1 a 2 l = 1, which implies a < 1. Hence, for the nest level, taking = 1 always guarantees (13) and (7) follows for our proposal t j;n. In general, there are L j non-zero constants a i;j so that y j;k = L X j i=1 a i;j z i8k2 j (14) where 8 in subscript is interpreted modulo n. Let a j = max i ja i;j j, then we have L j 2 j (L 0 1) 0 L + 2 and a j 2 0(j01)=2 ; where L is the length of the lter and = sup j1 2 (j01)=2 a j < 1:534. In Tables 1 and 2 of Appendix, we list some L j 's and 2 (j01)=2 a j 's for the commonly used wavelets (apparently they converge to sup x j (x)j). This leads to the threshold (3). 7
8 3 Numerical Examples In this section, we study four time series and compare our proposed method with the one based on Gaussian theory. Figure 1 contains an AR(24) signal, with coecients: Figure 2 contains a white noise siganl, the true log spectral density in this case is g(!) = log( 2 =2), a constant. Figure 3 contains an MA(15001) time series with coecients a 1 = =4; AR method is compared in this example. a n+1 = sin(n=2) ; n = 1; 2; :::; 15000: n Figure 4 contains the Sunspots signal, WaveShrink estimate with comparisons with an AR estimate from SPLUS function spec.ar. An AR(2) model, suggested by Priestley (1981), pp. 882, is also plotted. 8
9 (a) (b) (c) (d) Figure 1: (a). log spectrum of an AR(24) time series. (b). log periodogram, (c). WaveShrink estimate based the new thresholding scheme, and (d). WaveShrink estimate based on Gaussian theory. 9
10 (a) (b) (c) (d) Figure 2: (a). log spectrum of a white noise signal, (b). log periodogram, (c). WaveShrink estimate based the new thresholding scheme, and (d). WaveShrink estimate based on Gaussian theory. 10
11 (a) (b) (c) (d) Figure 3: (a). log spectrum of a long MA series, (b). log periodogram, (c). WaveShrink estimate based the new thresholding scheme, and (d). AR estimate from SPLUS function spec.ar. 11
12 (a) (b) (c) (d) Figure 4: (a). sunspots data, (b). log periodogram, (c). WaveShrink estimate (solid line) and AR(2) estimate (dotted line), and (d). AR (of order 27) estimate. 12
13 4 Appendix 4.1 Tables L1 H1 L2 H2 L3 H3 L4 H4 L5 H5 L6 H6 C C C C C D D D D D D D D D D S S S S S S S Table 1: Maximum Coecients, Lj low-pass cascade lter at level j, Hj high-pass cascade lter at level j. C for Coiets; D for Daublets, Daubechies' original wavelets; S for Symmlet, Daubechies' near-symmetric wavelets. 13
14 L=2 L=4 L=6 L=8 L=10 L=12 L=14 L=16 L=18 L=20 j= j= j= j= j= j= Table 2: Number of Non-zero Coecients, L is the length of the lter and j = level. 4.2 Proofs of the Theorems Proof of Theorem 1: First of all, P f [ j [sup jy j;k 0Ey j;k j>t j;n ]g k m01 X j=1 m01 X j=1 P fsup jy j;k 0 Ey j;k j > t j;n g k 2 j p j (t j;n ) n max p j(t j;n ) (15) Under model (9) and (10), from (14), coecients in level j can be written as y j;k 0 Ey j;k = L X j i=1 a i;j " i8k2 j with P i a i;j = 0 and P i a 2 i;j = 1. Let a j = max i ja i;j j. The moment generating function for " log( 2 2=2) is where 0 is Euler's Gamma function. Let 3 3 j(t) = inf s f0st + Lj M(u) Ee u" = 0(u + 1) X i=1 (log 0(1 + a i;j s) + log 0(1 0 a i;j s))g It can be shown that for any random variable X with EX = 0, for all a > 0 and t. Then and therefore, P (X > a) e 0at Ee tx p j (t) e 33 j (t) max p j(t) sup e 33(t) j = exp( sup 3 3 j (t)) (16) 14
15 Let `n = (log n) 7=6 and consider t 3 j = 8 < : log n L 1=4 j p log n L j < `n L j `n : (17) For the thresholds we have t j;n t 3 j. Therefore, t j;n = log n L 1=4 j _ p log n and sup p j (t j;n ) sup p j (t 3 j) exp( sup 3 3 j (t 3 j)) _ sup 3 3 j (t3 j ) sup 3 3 j (t3 j ) sup 3 3 j (t3 j ) 1L j <`n `nl j n = sup 1L j <`n When L j < `n = (log n) 7=6 and n large, we have log n > L 1=4 j 3 3 j( log n _ p ) sup 3 3 j( log n) (18) L 1=4 `nl j j n q L j a j L j Consider the following non-positive function H, dened on R + by Lemma 1 For t a j L j, and H(t) = 0t log t 3 3 t j(t) L j H( ) a j L j L H( log n ) # in L for large n and 1 L `n: al3=4 The proof will be given later. Combining (16) and Lemma 1 yields the second part of (11). From the Lemma, sup 1L j <`n 3 3 j( log n ) sup L 1=4 1L j j <`n L j H( log n ) al 3=4 j H( a log n) = 0 log n (1 + o(1)) (19) a and When L j `n, we have p log n L 3=7 j 15
16 Lemma 2 For j 2 (0; 1) and 0 t 2 j =6a j, t j(t) 0 : 2 =3 + 6 j In particular, for n large and j 6a j L 3=7 j = 2 L 01=14 j (log n) 01=12, 3 3 j The proof will be given later. From the Lemma, (t) 03t2 (1 + o(1)) t 2 L3=7 j : p sup 3 3 j( log n) 0 32 log n (1 + o(1)): (20) `nl j n 2 Combining (18), (19) and (20), we have sup 3 3 j(t 3 j ) 0( a Therefore, for any > a and > = p 3, n sup and this completes the proof of Theorem 1. ^ 3 ) log n (1 + o(1)): 2 p j (t j;n ) n exp( sup 3 3 j (t j;n )) = o(1): Proof of Lemma 1: First of all, it can be shown that for 0 x < 1, 0(1 + x) 0(1 0 x) and 0(x) 1=x. (In fact, 0(x) = Q x 01 1 (1+1=k) x k=1 1+x=k and (1 + 1=k) x 1 + x=k when 0 x < 1.) Therefore, 3 3 j (t) inf jsj<1=a j f0st + L X j i=1 log 0(1 + a i;j s)g L j inf f0 st 0 log(1 0 s)g 0s<1 a j L j and the rst part of the Lemma can be easily solved from the last expression. For n large, such that C n a 01 log n > 11 and `n = (log n) 7=6 (C n =11) 4=3. Let h(x) = xh(c n x 03=4 ), then for 1 x < `n, h 0 (x) is increasing in x and h 0 (`n) < 0. Therefore, h 0 (x) < 0 for 1 x < `n and this means h(x) is decreasing in x. The proof is completed. Proof of Lemma 2: Let K(x) = log 0(1 + x), then for x > 01, by Taylor expansion, K(x) = K(0) + K 0 (0)x + x 2 K 00 (x)=2 = 0x + x 2 K 00 (x)=2; (21) 16
17 where = 0:5772:::: is the Euler constant and = (x) 2 (0; 1). Note that K 0 (x 0 1) is the famous Psi function. By the properties of Psi function (Davis (1935), pp. 11), K 00 (x) = P 1 k=1 (k + x) 02 and K 00 (0) = 2 =6. Hence, for any 2 (0; 1), we have K 00 (0) K 00 (0) + 3 K and for x 0, Recall that P i a i;j = 0 and P i a 2 i;j = 1, K(x) 0x + x 2 K =2: X 3 3 j (t) = inf f0st + s i K(a i;j s)g and Lemma 2 follows immediately. X inf f0st + (0a i;j s + a 2 i;j jsj< j =a s2 K =2)g j i = inf f0st + s 2 K =2g jsj< j =a j Proof of Theorem 2: To prove Theorem 2, we need a lower bound for p j (t). It is easy to show that if X; Y are independent, then for any s; t > 0, P fjx + Y j tg P fjxj < sgp fjy j t + sg: Suppose ja 1 j = max ja l j = a. There exists t 0 > 0 such that for t t 0, P fj LX l=2 a l log( l =2)j < tg 1 2 : Also notice that for 0 x 1, 1 0 e 0x x=2. Then p 1 (t) = P fj LX l=1 a l log( l =2)j tg P fja 1 log( 1 =2)j 2tgP fj 1 2 P fj log( 1=2)j 2t a g 1 2 P flog( 1=2) 0 2t a g and this completes the proof of Lemma 1. LX l=2 a l log( l =2)j < tg = 1 2 f1 0 exp(0e02t=a )g 1 4 e02t=a : (22) From (12) we can see that for any constant C < a=2, np 1 (C log n)! 1; n! 1: 17
18 In particular, p s 2 np 1 ( 2 log n) = np 1 ( log n)! 1; n! 1: log n There exists K > 0 such that fy 1;iK g are independent (linear combinations over disjoint intervals). Therefore, P fsup k p p jy 1;k 0 Ey 1;k j > 2 log ng P fsup jy 1;iK 0 Ey 1;iK j > 2 log ng i Y and this implies the Theorem 2. n=2k = 1 0 i=0 (1 0 p 1 ( p 2 log n))! 1; References [1] Bentkus, R. Yu. and Susinskas, Ju. V. (1982) On Optimal Statistical Estimators of the Spectral Density. Soviet Math. Dokl., Vol. 25, No. 2, [2] Brillinger, D. R. (1981). Time Series: Data Analysis and Theory, Holden-Day, San Francisco. [3] Burg, J. P. (1967). Maximum Entropy Spectral Analysis. Paper Presented at the 37th Annual international S.E.G. Meeting, Oklahoma City, Oklahoma. [4] Burg, J. P. (1972). The Relationship between Maximum Entropy Spectra and Maximum Likelihood Spectra. Geophysice, 37, [5] Daubechies, I. (1992). Ten Lectures on Wavelets. SIAM. Philadelphia. [6] Davis, H. T. (1935). Tables of The Mathematical Functions, Vol II, Trinity University. [7] Donoho, D. (1992). De-Noising By Soft-Thresholding. Technical Report, Department of Statistics, Stanford University. [8] Donoho, D. and Johnstone, I (1992a). Minimax Estimation via Wavelet Shrinkage. Technical Report, Department of Statistics, Stanford University. [9] Donoho, D. and Johnstone, I (1992b). Adapting to Unknown Smoothness via Wavelet Shrinkage. Technical Report, Department of Statistics, Stanford University. [10] Donoho, D. and Johnstone, I (1992c). Ideal Spatial Adaptation by Wavelet Shrinkage. Technical Report, Department of Statistics, Stanford University. [11] Donoho, D., Johnstone, I., Kerkyacharian, G. and Picard, D. (1992). Wavelet Shrinkage: Asymptopia? Technical Report, Department of Statistics, Stanford University. 18
19 [12] Gao, Hong-Ye (1992). Spectral Density Estimation by Using Wavelets. Technical Report, Department of Statistics, Stanford University. [13] Harvey, A. C. (1981). Time Series Models, Philip Allan. [14] Harvey, A. C. (1989). Forecasting, Structural Time Series Models and the Kalman Filter. Cambridge University Press. [15] Johnstone, I., Kerkyacharian, G. and Picard, D. (1992). Estimation d'une densite de probabilite par methode d'ondelettes. Comptes Rendus Acad. Sciences Paris (A), [16] Moulin, P. (1993a) A Wavelet Regularization Method for Diuse Radar-Target Imaging and Speckle-Noise Reduction. Journal of Math. Imaging and Vision, Special Issue on Wavelet, Vol. 3, [17] Moulin, P. (1993b) Wavelet Thresholding Techniques for Power Spectrum Estimation. Submitted to IEEE Trans. on Signal Processing. [18] Parzen, E. (1974). Some Recent Advances in Time Series Modelling. IEEE Trans. Automatic Control, AC-19, [19] Priestley, M. B. (1981). Spectral Analysis and Time Series, Academic Press, London. [20] Taniguchi, M. (1979). On Estimation of Parameters of Gaussian Stationary Processes. J. Appl. Probab., 16, [21] Taniguchi, M. (1980). On Estimation of the Integrals of Certain Functions of Spectral Density. J. Appl. Probab., 17, [22] Tukey, J. W. (1978). When Should Which Spectrum Approach Be Used? Paper Presented at the Institute of Statisticians Conference, King's College, Cambridge, July, [23] Wahba, G. (1980). Automatic Smoothing of the Log Periodogram. Journal of the American Statistical Association, 75, [24] Wittwer, G. (1986) On the Distribution of the Periodogram for Stationary Random Sequences (with Discussion). Math. Opera. und Stat., 17,
Introduction Wavelet shrinage methods have been very successful in nonparametric regression. But so far most of the wavelet regression methods have be
Wavelet Estimation For Samples With Random Uniform Design T. Tony Cai Department of Statistics, Purdue University Lawrence D. Brown Department of Statistics, University of Pennsylvania Abstract We show
More information(a). Bumps (b). Wavelet Coefficients
Incorporating Information on Neighboring Coecients into Wavelet Estimation T. Tony Cai Bernard W. Silverman Department of Statistics Department of Mathematics Purdue University University of Bristol West
More informationThe Lifting Wavelet Transform for Periodogram Smoothing
ISSN : 976-8491 (Online) ISSN : 2229-4333 (Print) IJCST Vo l. 3, Is s u e 1, Ja n. - Ma r c h 212 The Lifting for Periodogram Smoothing 1 M.Venakatanarayana, 2 Dr. T.Jayachandra Prasad 1 Dept. of ECE,
More informationWavelet Thresholding Techniques for Power Spectrum Estimation
Wavelet Thresholding Techniques for Power Spectrum Estimation Pierre Moulin Bell Communications Research 445 South Street Morristown, NJ 07960 email: moulin@bellcore.com tel (201) 829-5021 fax (201) 829-4391
More informationWavelet Shrinkage for Nonequispaced Samples
University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 1998 Wavelet Shrinkage for Nonequispaced Samples T. Tony Cai University of Pennsylvania Lawrence D. Brown University
More informationDiscussion of Regularization of Wavelets Approximations by A. Antoniadis and J. Fan
Discussion of Regularization of Wavelets Approximations by A. Antoniadis and J. Fan T. Tony Cai Department of Statistics The Wharton School University of Pennsylvania Professors Antoniadis and Fan are
More informationComments on \Wavelets in Statistics: A Review" by. A. Antoniadis. Jianqing Fan. University of North Carolina, Chapel Hill
Comments on \Wavelets in Statistics: A Review" by A. Antoniadis Jianqing Fan University of North Carolina, Chapel Hill and University of California, Los Angeles I would like to congratulate Professor Antoniadis
More informationAdaptive Bayesian Wavelet Shrinkage. Hugh A. Chipman, Eric D. Kolaczyk, Robert E. McCulloch
Adaptive Bayesian Wavelet Shrinkage Hugh A. Chipman, Eric D. Kolaczyk, Robert E. McCulloch To appear in the Journal of the American Statistical Association. Hugh A. Chipman and Robert E. McCulloch are
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationIntroduction to Spectral and Time-Spectral Analysis with some Applications
Introduction to Spectral and Time-Spectral Analysis with some Applications A.INTRODUCTION Consider the following process X t = U cos[(=3)t] + V sin[(=3)t] Where *E[U] = E[V ] = 0 *E[UV ] = 0 *V ar(u) =
More informationAutomatic Local Smoothing for Spectral Density. Abstract. This article uses local polynomial techniques to t Whittle's likelihood for spectral density
Automatic Local Smoothing for Spectral Density Estimation Jianqing Fan Department of Statistics University of North Carolina Chapel Hill, N.C. 27599-3260 Eva Kreutzberger Department of Mathematics University
More informationTime series denoising with wavelet transform
Paper Time series denoising with wavelet transform Bartosz Kozłowski Abstract This paper concerns the possibilities of applying wavelet analysis to discovering and reducing distortions occurring in time
More informationSize properties of wavelet packets generated using finite filters
Rev. Mat. Iberoamericana, 18 (2002, 249 265 Size properties of wavelet packets generated using finite filters Morten Nielsen Abstract We show that asymptotic estimates for the growth in L p (R- norm of
More informationTime Series: Theory and Methods
Peter J. Brockwell Richard A. Davis Time Series: Theory and Methods Second Edition With 124 Illustrations Springer Contents Preface to the Second Edition Preface to the First Edition vn ix CHAPTER 1 Stationary
More informationREAL-TIME TIME-FREQUENCY BASED BLIND SOURCE SEPARATION. Scott Rickard, Radu Balan, Justinian Rosca. Siemens Corporate Research Princeton, NJ 08540
REAL-TIME TIME-FREQUENCY BASED BLIND SOURCE SEPARATION Scott Rickard, Radu Balan, Justinian Rosca Siemens Corporate Research Princeton, NJ 84 fscott.rickard,radu.balan,justinian.roscag@scr.siemens.com
More informationResidual Bootstrap for estimation in autoregressive processes
Chapter 7 Residual Bootstrap for estimation in autoregressive processes In Chapter 6 we consider the asymptotic sampling properties of the several estimators including the least squares estimator of the
More informationGaussian distributions and processes have long been accepted as useful tools for stochastic
Chapter 3 Alpha-Stable Random Variables and Processes Gaussian distributions and processes have long been accepted as useful tools for stochastic modeling. In this section, we introduce a statistical model
More information[km]
SEGMENTATION INDEPENDENT ESTIMATES OF TURBULENCE PARAMETERS G. C. Papanicolaou a, K. Solna b and D. Washburn c adepartment of Mathematics, Stanford University, Stanford, CA 9435 bdepartment of Mathematics,
More informationcovariance function, 174 probability structure of; Yule-Walker equations, 174 Moving average process, fluctuations, 5-6, 175 probability structure of
Index* The Statistical Analysis of Time Series by T. W. Anderson Copyright 1971 John Wiley & Sons, Inc. Aliasing, 387-388 Autoregressive {continued) Amplitude, 4, 94 case of first-order, 174 Associated
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY & Contents PREFACE xiii 1 1.1. 1.2. Difference Equations First-Order Difference Equations 1 /?th-order Difference
More informationSmooth functions and local extreme values
Smooth functions and local extreme values A. Kovac 1 Department of Mathematics University of Bristol Abstract Given a sample of n observations y 1,..., y n at time points t 1,..., t n we consider the problem
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7
More informationAn algorithm for robust fitting of autoregressive models Dimitris N. Politis
An algorithm for robust fitting of autoregressive models Dimitris N. Politis Abstract: An algorithm for robust fitting of AR models is given, based on a linear regression idea. The new method appears to
More informationFull terms and conditions of use:
This article was downloaded by:[smu Cul Sci] [Smu Cul Sci] On: 28 March 2007 Access Details: [subscription number 768506175] Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered
More informationAggregation of Spectral Density Estimators
Aggregation of Spectral Density Estimators Christopher Chang Department of Mathematics, University of California, San Diego, La Jolla, CA 92093-0112, USA; chrchang@alumni.caltech.edu Dimitris Politis Department
More information524 Jan-Olov Stromberg practice, this imposes strong restrictions both on N and on d. The current state of approximation theory is essentially useless
Doc. Math. J. DMV 523 Computation with Wavelets in Higher Dimensions Jan-Olov Stromberg 1 Abstract. In dimension d, a lattice grid of size N has N d points. The representation of a function by, for instance,
More informationAdaptive Wavelet Estimation: A Block Thresholding and Oracle Inequality Approach
University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 1999 Adaptive Wavelet Estimation: A Block Thresholding and Oracle Inequality Approach T. Tony Cai University of Pennsylvania
More informationOn the usefulness of wavelet-based simulation of fractional Brownian motion
On the usefulness of wavelet-based simulation of fractional Brownian motion Vladas Pipiras University of North Carolina at Chapel Hill September 16, 2004 Abstract We clarify some ways in which wavelet-based
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationMathematical Institute, University of Utrecht. The problem of estimating the mean of an observed Gaussian innite-dimensional vector
On Minimax Filtering over Ellipsoids Eduard N. Belitser and Boris Y. Levit Mathematical Institute, University of Utrecht Budapestlaan 6, 3584 CD Utrecht, The Netherlands The problem of estimating the mean
More informationWavelet Shrinkage for Unequally Spaced Data Sylvain Sardy, Donald B Percival 3, Andrew G Bruce, Hong-Ye Gao 4 and Werner Stuetzle MathSoft, Inc, 700 W
MathSoft Research Report No 4 Wavelet Shrinkage for Unequally Spaced Data Sylvain Sardy, Donald B Percival, Andrew G Bruce, Hong-Ye Gao and Werner Stuetzle Last revised April 8, 998 Acknowledgements: This
More informationNonparametric Estimation of Functional-Coefficient Autoregressive Models
Nonparametric Estimation of Functional-Coefficient Autoregressive Models PEDRO A. MORETTIN and CHANG CHIANN Department of Statistics, University of São Paulo Introduction Nonlinear Models: - Exponential
More informationMultivariate Bayes Wavelet Shrinkage and Applications
Journal of Applied Statistics Vol. 32, No. 5, 529 542, July 2005 Multivariate Bayes Wavelet Shrinkage and Applications GABRIEL HUERTA Department of Mathematics and Statistics, University of New Mexico
More informationSYSTEM RECONSTRUCTION FROM SELECTED HOS REGIONS. Haralambos Pozidis and Athina P. Petropulu. Drexel University, Philadelphia, PA 19104
SYSTEM RECOSTRUCTIO FROM SELECTED HOS REGIOS Haralambos Pozidis and Athina P. Petropulu Electrical and Computer Engineering Department Drexel University, Philadelphia, PA 94 Tel. (25) 895-2358 Fax. (25)
More informationFast Angular Synchronization for Phase Retrieval via Incomplete Information
Fast Angular Synchronization for Phase Retrieval via Incomplete Information Aditya Viswanathan a and Mark Iwen b a Department of Mathematics, Michigan State University; b Department of Mathematics & Department
More informationGMM-based inference in the AR(1) panel data model for parameter values where local identi cation fails
GMM-based inference in the AR() panel data model for parameter values where local identi cation fails Edith Madsen entre for Applied Microeconometrics (AM) Department of Economics, University of openhagen,
More informationRobust Nonlinear Wavelet Transform based on Median-Interpolation
Robust Nonlinear Wavelet Transform based on Median-Interpolation David L. Donoho & Thomas P.Y. Yu Department of Statistics Stanford University Stanford, CA 943, USA donoho@stat.stanford.edu pyu@stat.stanford.edu
More informationA Data-Driven Block Thresholding Approach To Wavelet Estimation
A Data-Driven Block Thresholding Approach To Wavelet Estimation T. Tony Cai 1 and Harrison H. Zhou University of Pennsylvania and Yale University Abstract A data-driven block thresholding procedure for
More informationDesign of Image Adaptive Wavelets for Denoising Applications
Design of Image Adaptive Wavelets for Denoising Applications Sanjeev Pragada and Jayanthi Sivaswamy Center for Visual Information Technology International Institute of Information Technology - Hyderabad,
More informationMultiresolution Models of Time Series
Multiresolution Models of Time Series Andrea Tamoni (Bocconi University ) 2011 Tamoni Multiresolution Models of Time Series 1/ 16 General Framework Time-scale decomposition General Framework Begin with
More informationROYAL INSTITUTE OF TECHNOLOGY KUNGL TEKNISKA HÖGSKOLAN. Department of Signals, Sensors & Systems
The Evil of Supereciency P. Stoica B. Ottersten To appear as a Fast Communication in Signal Processing IR-S3-SB-9633 ROYAL INSTITUTE OF TECHNOLOGY Department of Signals, Sensors & Systems Signal Processing
More informationLinear Regression and Its Applications
Linear Regression and Its Applications Predrag Radivojac October 13, 2014 Given a data set D = {(x i, y i )} n the objective is to learn the relationship between features and the target. We usually start
More informationdata lam=36.9 lam=6.69 lam=4.18 lam=2.92 lam=2.21 time max wavelength modulus of max wavelength cycle
AUTOREGRESSIVE LINEAR MODELS AR(1) MODELS The zero-mean AR(1) model x t = x t,1 + t is a linear regression of the current value of the time series on the previous value. For > 0 it generates positively
More information1. Introduction Over the last three decades a number of model selection criteria have been proposed, including AIC (Akaike, 1973), AICC (Hurvich & Tsa
On the Use of Marginal Likelihood in Model Selection Peide Shi Department of Probability and Statistics Peking University, Beijing 100871 P. R. China Chih-Ling Tsai Graduate School of Management University
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationWavelet Neural Networks for Nonlinear Time Series Analysis
Applied Mathematical Sciences, Vol. 4, 2010, no. 50, 2485-2495 Wavelet Neural Networks for Nonlinear Time Series Analysis K. K. Minu, M. C. Lineesh and C. Jessy John Department of Mathematics National
More informationSea Surface. Bottom OBS
ANALYSIS OF HIGH DIMENSIONAL TIME SERIES: OCEAN BOTTOM SEISMOGRAPH DATA Genshiro Kitagawa () and Tetsuo Takanami (2) () The Institute of Statistical Mathematics, 4-6-7 Minami-Azabu, Minato-ku, Tokyo 06-8569
More informationA New Method for Varying Adaptive Bandwidth Selection
IEEE TRASACTIOS O SIGAL PROCESSIG, VOL. 47, O. 9, SEPTEMBER 1999 2567 TABLE I SQUARE ROOT MEA SQUARED ERRORS (SRMSE) OF ESTIMATIO USIG THE LPA AD VARIOUS WAVELET METHODS A ew Method for Varying Adaptive
More informationNonparametric Regression
Adaptive Variance Function Estimation in Heteroscedastic Nonparametric Regression T. Tony Cai and Lie Wang Abstract We consider a wavelet thresholding approach to adaptive variance function estimation
More informationApplied Time. Series Analysis. Wayne A. Woodward. Henry L. Gray. Alan C. Elliott. Dallas, Texas, USA
Applied Time Series Analysis Wayne A. Woodward Southern Methodist University Dallas, Texas, USA Henry L. Gray Southern Methodist University Dallas, Texas, USA Alan C. Elliott University of Texas Southwestern
More information1 Introduction This work follows a paper by P. Shields [1] concerned with a problem of a relation between the entropy rate of a nite-valued stationary
Prexes and the Entropy Rate for Long-Range Sources Ioannis Kontoyiannis Information Systems Laboratory, Electrical Engineering, Stanford University. Yurii M. Suhov Statistical Laboratory, Pure Math. &
More informationA Wavelet-based Method for Overcoming the Gibbs Phenomenon
AMERICAN CONFERENCE ON APPLIED MATHEMATICS (MATH '08), Harvard, Massachusetts, USA, March 4-6, 008 A Wavelet-based Method for Overcoming the Gibbs Phenomenon NATANIEL GREENE Department of Mathematics Computer
More informationWhich wavelet bases are the best for image denoising?
Which wavelet bases are the best for image denoising? Florian Luisier a, Thierry Blu a, Brigitte Forster b and Michael Unser a a Biomedical Imaging Group (BIG), Ecole Polytechnique Fédérale de Lausanne
More informationModel Selection and Geometry
Model Selection and Geometry Pascal Massart Université Paris-Sud, Orsay Leipzig, February Purpose of the talk! Concentration of measure plays a fundamental role in the theory of model selection! Model
More informationEstimation Error Bounds for Frame Denoising
Estimation Error Bounds for Frame Denoising Alyson K. Fletcher and Kannan Ramchandran {alyson,kannanr}@eecs.berkeley.edu Berkeley Audio-Visual Signal Processing and Communication Systems group Department
More informationPOINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES
November 1, 1 POINT VALUES AND NORMALIZATION OF TWO-DIRECTION MULTIWAVELETS AND THEIR DERIVATIVES FRITZ KEINERT AND SOON-GEOL KWON,1 Abstract Two-direction multiscaling functions φ and two-direction multiwavelets
More informationWavelet Thresholding for Non Necessarily Gaussian Noise: Functionality
Wavelet Thresholding for Non Necessarily Gaussian Noise: Functionality By R. Averkamp 1 and C. Houdré 2 Freiburg University and Université Paris XII and Georgia Institute of Technology For signals belonging
More informationGoodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach
Goodness-of-Fit Tests for Time Series Models: A Score-Marked Empirical Process Approach By Shiqing Ling Department of Mathematics Hong Kong University of Science and Technology Let {y t : t = 0, ±1, ±2,
More informationTime Since Burst Trigger (Seconds) TTS Data, Energy Channel 3, kev
SLAC{PUB{8364 February 2 Properties of Gamma-Ray Burst Time Proles Using Pulse Decomposition Analysis Andrew Lee and Elliott D. Bloom Stanford Linear Accelerator Center, Stanford University, Stanford,
More informationA general theory of discrete ltering. for LES in complex geometry. By Oleg V. Vasilyev AND Thomas S. Lund
Center for Turbulence Research Annual Research Briefs 997 67 A general theory of discrete ltering for ES in complex geometry By Oleg V. Vasilyev AND Thomas S. und. Motivation and objectives In large eddy
More informationTwo Denoising Methods by Wavelet Transform
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 47, NO. 12, DECEMBER 1999 3401 Two Denoising Methods by Wavelet Transform Quan Pan, Lei Zhang, Guanzhong Dai, and Hongcai Zhang Abstract Two wavelet-based noise
More informationWavelet Regression Estimation in Longitudinal Data Analysis
Wavelet Regression Estimation in Longitudinal Data Analysis ALWELL J. OYET and BRAJENDRA SUTRADHAR Department of Mathematics and Statistics, Memorial University of Newfoundland St. John s, NF Canada, A1C
More informationKey Words and Phrases. Segmented Multi-Resolution analysis. Edge- Preserving Image processing methods. Edge detection. Sub-pixel resolution. De-Noisin
On Minimum Entropy Segmentation David L. Donoho Department of Statistics Stanford University September 1993 Revised March 1994 Abstract We describe segmented multiresolution analyses of [0; 1]. Such multiresolution
More informationStatistics of Stochastic Processes
Prof. Dr. J. Franke All of Statistics 4.1 Statistics of Stochastic Processes discrete time: sequence of r.v...., X 1, X 0, X 1, X 2,... X t R d in general. Here: d = 1. continuous time: random function
More informationLikelihood Ratio Tests and Intersection-Union Tests. Roger L. Berger. Department of Statistics, North Carolina State University
Likelihood Ratio Tests and Intersection-Union Tests by Roger L. Berger Department of Statistics, North Carolina State University Raleigh, NC 27695-8203 Institute of Statistics Mimeo Series Number 2288
More informationNon-Gaussian Maximum Entropy Processes
Non-Gaussian Maximum Entropy Processes Georgi N. Boshnakov & Bisher Iqelan First version: 3 April 2007 Research Report No. 3, 2007, Probability and Statistics Group School of Mathematics, The University
More informationBayesian Regularization
Bayesian Regularization Aad van der Vaart Vrije Universiteit Amsterdam International Congress of Mathematicians Hyderabad, August 2010 Contents Introduction Abstract result Gaussian process priors Co-authors
More information1/sqrt(B) convergence 1/B convergence B
The Error Coding Method and PICTs Gareth James and Trevor Hastie Department of Statistics, Stanford University March 29, 1998 Abstract A new family of plug-in classication techniques has recently been
More informationResearch Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance
Advances in Decision Sciences Volume, Article ID 893497, 6 pages doi:.55//893497 Research Article Least Squares Estimators for Unit Root Processes with Locally Stationary Disturbance Junichi Hirukawa and
More information446 SCIENCE IN CHINA (Series F) Vol. 46 introduced in refs. [6, ]. Based on this inequality, we add normalization condition, symmetric conditions and
Vol. 46 No. 6 SCIENCE IN CHINA (Series F) December 003 Construction for a class of smooth wavelet tight frames PENG Lizhong (Λ Π) & WANG Haihui (Ξ ) LMAM, School of Mathematical Sciences, Peking University,
More informationAdditive functionals of infinite-variance moving averages. Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535
Additive functionals of infinite-variance moving averages Wei Biao Wu The University of Chicago TECHNICAL REPORT NO. 535 Departments of Statistics The University of Chicago Chicago, Illinois 60637 June
More informationWavelet Footprints: Theory, Algorithms, and Applications
1306 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 5, MAY 2003 Wavelet Footprints: Theory, Algorithms, and Applications Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli, Fellow, IEEE Abstract
More informationFourier Analysis of Stationary and Non-Stationary Time Series
Fourier Analysis of Stationary and Non-Stationary Time Series September 6, 2012 A time series is a stochastic process indexed at discrete points in time i.e X t for t = 0, 1, 2, 3,... The mean is defined
More informationTime Series I Time Domain Methods
Astrostatistics Summer School Penn State University University Park, PA 16802 May 21, 2007 Overview Filtering and the Likelihood Function Time series is the study of data consisting of a sequence of DEPENDENT
More informationPCA with random noise. Van Ha Vu. Department of Mathematics Yale University
PCA with random noise Van Ha Vu Department of Mathematics Yale University An important problem that appears in various areas of applied mathematics (in particular statistics, computer science and numerical
More informationThe Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-Zero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In
The Best Circulant Preconditioners for Hermitian Toeplitz Systems II: The Multiple-ero Case Raymond H. Chan Michael K. Ng y Andy M. Yip z Abstract In [0, 4], circulant-type preconditioners have been proposed
More information9 Brownian Motion: Construction
9 Brownian Motion: Construction 9.1 Definition and Heuristics The central limit theorem states that the standard Gaussian distribution arises as the weak limit of the rescaled partial sums S n / p n of
More informationEstimating the Number of Common Factors in Serially Dependent Approximate Factor Models
Estimating the Number of Common Factors in Serially Dependent Approximate Factor Models Ryan Greenaway-McGrevy y Bureau of Economic Analysis Chirok Han Korea University February 7, 202 Donggyu Sul University
More informationTIME SERIES ANALYSIS. Forecasting and Control. Wiley. Fifth Edition GWILYM M. JENKINS GEORGE E. P. BOX GREGORY C. REINSEL GRETA M.
TIME SERIES ANALYSIS Forecasting and Control Fifth Edition GEORGE E. P. BOX GWILYM M. JENKINS GREGORY C. REINSEL GRETA M. LJUNG Wiley CONTENTS PREFACE TO THE FIFTH EDITION PREFACE TO THE FOURTH EDITION
More informationParametric Method Based PSD Estimation using Gaussian Window
International Journal of Engineering Trends and Technology (IJETT) Volume 29 Number 1 - November 215 Parametric Method Based PSD Estimation using Gaussian Window Pragati Sheel 1, Dr. Rajesh Mehra 2, Preeti
More informationA Lower Bound Theorem. Lin Hu.
American J. of Mathematics and Sciences Vol. 3, No -1,(January 014) Copyright Mind Reader Publications ISSN No: 50-310 A Lower Bound Theorem Department of Applied Mathematics, Beijing University of Technology,
More informationOle Christensen 3. October 20, Abstract. We point out some connections between the existing theories for
Frames and pseudo-inverses. Ole Christensen 3 October 20, 1994 Abstract We point out some connections between the existing theories for frames and pseudo-inverses. In particular, using the pseudo-inverse
More informationConstruction of Multivariate Compactly Supported Orthonormal Wavelets
Construction of Multivariate Compactly Supported Orthonormal Wavelets Ming-Jun Lai Department of Mathematics The University of Georgia Athens, GA 30602 April 30, 2004 Dedicated to Professor Charles A.
More informationAN INEQUALITY FOR THE NORM OF A POLYNOMIAL FACTOR IGOR E. PRITSKER. (Communicated by Albert Baernstein II)
PROCDINGS OF TH AMRICAN MATHMATICAL SOCITY Volume 9, Number 8, Pages 83{9 S -9939()588-4 Article electronically published on November 3, AN INQUALITY FOR TH NORM OF A POLYNOMIAL FACTOR IGOR. PRITSKR (Communicated
More informationSpatial autoregression model:strong consistency
Statistics & Probability Letters 65 (2003 71 77 Spatial autoregression model:strong consistency B.B. Bhattacharyya a, J.-J. Ren b, G.D. Richardson b;, J. Zhang b a Department of Statistics, North Carolina
More informationStatistics 349(02) Review Questions
Statistics 349(0) Review Questions I. Suppose that for N = 80 observations on the time series { : t T} the following statistics were calculated: _ x = 10.54 C(0) = 4.99 In addition the sample autocorrelation
More informationSuper-resolution via Convex Programming
Super-resolution via Convex Programming Carlos Fernandez-Granda (Joint work with Emmanuel Candès) Structure and Randomness in System Identication and Learning, IPAM 1/17/2013 1/17/2013 1 / 44 Index 1 Motivation
More informationTime Series Analysis -- An Introduction -- AMS 586
Time Series Analysis -- An Introduction -- AMS 586 1 Objectives of time series analysis Data description Data interpretation Modeling Control Prediction & Forecasting 2 Time-Series Data Numerical data
More informationoutput dimension input dimension Gaussian evidence Gaussian Gaussian evidence evidence from t +1 inputs and outputs at time t x t+2 x t-1 x t+1
To appear in M. S. Kearns, S. A. Solla, D. A. Cohn, (eds.) Advances in Neural Information Processing Systems. Cambridge, MA: MIT Press, 999. Learning Nonlinear Dynamical Systems using an EM Algorithm Zoubin
More informationParametric Inference on Strong Dependence
Parametric Inference on Strong Dependence Peter M. Robinson London School of Economics Based on joint work with Javier Hualde: Javier Hualde and Peter M. Robinson: Gaussian Pseudo-Maximum Likelihood Estimation
More informationHigher order correlation detection in nonlinear aerodynamic systems using wavelet transforms
Higher order correlation detection in nonlinear aerodynamic systems using wavelet transforms K. Gurley Department of Civil and Coastal Engineering, University of Florida, USA T. Kijewski & A. Kareem Department
More informationA Modified Confidence Set for the S Title Date in Linear Regression Models.
A Modified Confidence Set for the S Title Date in Linear Regression Models Author(s) Yamamoto, Yohei Citation Issue 24-5-7 Date Type Technical Report Text Version publisher URL http://hdl.handle.net/86/26678
More informationDigital Filter Design Using Non-Uniform Sampling
Proceedings of the 6th WSEAS International Conference on Applications of Electrical Engineering, Istanbul, Turkey, May 27-29, 2007 120 Digital Filter Design Using Non-Uniform Sampling K. BUSAWON, WAI-PANG
More informationSubmitted to the Brazilian Journal of Probability and Statistics
Submitted to the Brazilian Journal of Probability and Statistics Multivariate normal approximation of the maximum likelihood estimator via the delta method Andreas Anastasiou a and Robert E. Gaunt b a
More informationObjective Functions for Tomographic Reconstruction from. Randoms-Precorrected PET Scans. gram separately, this process doubles the storage space for
Objective Functions for Tomographic Reconstruction from Randoms-Precorrected PET Scans Mehmet Yavuz and Jerey A. Fessler Dept. of EECS, University of Michigan Abstract In PET, usually the data are precorrected
More informationH. Krim and J.-C. Pesquet Following a section of preliminaries and denitions of notations, we derive in Section 3 a procedure for an estimation enhanc
On the Statistics of Best Bases Criteria? H. Krim and J.-C. Pesquet Stochastic Systems Group, LIDS, Massachusetts Institute of Technology, Cambridge MA 039, USA Laboratoire des Signaux et Systemes, CNRS/UPS
More informationRelative Irradiance. Wavelength (nm)
Characterization of Scanner Sensitivity Gaurav Sharma H. J. Trussell Electrical & Computer Engineering Dept. North Carolina State University, Raleigh, NC 7695-79 Abstract Color scanners are becoming quite
More informationDiagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations
Diagnostic Test for GARCH Models Based on Absolute Residual Autocorrelations Farhat Iqbal Department of Statistics, University of Balochistan Quetta-Pakistan farhatiqb@gmail.com Abstract In this paper
More informationA Subspace Approach to Estimation of. Measurements 1. Carlos E. Davila. Electrical Engineering Department, Southern Methodist University
EDICS category SP 1 A Subspace Approach to Estimation of Autoregressive Parameters From Noisy Measurements 1 Carlos E Davila Electrical Engineering Department, Southern Methodist University Dallas, Texas
More informationStatistical inference on Lévy processes
Alberto Coca Cabrero University of Cambridge - CCA Supervisors: Dr. Richard Nickl and Professor L.C.G.Rogers Funded by Fundación Mutua Madrileña and EPSRC MASDOC/CCA student workshop 2013 26th March Outline
More information