Discussion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli
|
|
- Aubrey Hill
- 5 years ago
- Views:
Transcription
1 1 Introdction Discssion of The Forward Search: Theory and Data Analysis by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen Department of Economics, University of Copenhagen and CREATES, University of Aarhs Bent Nielsen, Department of Economics, University of Oxford and N eld College 15 Janary 010 The paper by Atkinson, Riani and Ceroli, henceforth ARC, is concerned with detection of otliers and nsspected strctres which is rather important in practice. This is done throgh a Forward Search Algorithm. The statistical analysis of sch algorithms poses many challenging problems, and we wold like to contribte to the theory of the algorithm in this discssion. We establish some reslts in a simple case for a single iteration of the algorithm sing empirical process theory. This wold then have to be extended to more models, and developed frther to nderstand the properties of the fll algorithm. The established reslts sggest that the heristic reslts from ARC cold be correct if the parameters were known, bt not when the parameters are estimated. A general reference for the empirical process theory is the monograph by Kol (00) which analyses weighted empirical processes. Some frther developments are made in Johansen and Nielsen (009) which we exploit here. For simplicity we only consider a simple location-scale problem bt the reslts wold generalize to regressions and time-series regressions. The discssion is organized so that the algorithm is described in. Some potential reslts are described in 3. The analysis of a single step of the algorithm is then provided for the known parameter case in 4 and for the nknown parameter case in 5. We conclde in 6 and leave some proofs to an Appendix. The algorithm We consider the regression problem y i = +" i ; i = 1; : : : ; n; where the i.i.d. " i follow a known distribtion fnction F and symmetric density f, with mean zero and variance 1. The distribtion is assmed to be continos and to satisfy some reglarity conditions for smoothness which are met by the standard normal distribtion, see Johansen and Nielsen (009, Assmption A). The absolte standardized error j" i j has distribtion 1
2 fnction G satisfying G() = P(j" i j ) = F() F( ); and = G 1 ( ) is the -qantile of G and its density is g() = f():.1 Description of the algorithm We start with some initial robst location estimator ^ and an initial observation set of size m 0 : This set is constrcted by calclating the absolte residals ^r i = jy i ^j; nding the order statistics ^r (i), and de ning the initial observation index set of size m 0 as the m 0 observations closest to ^; that is The algorithm then proceeds in the steps S (m 0) = fi : jy i ^j ^r (m0 )g: 1. Given an index set S (m) calclate estimators ^ (m) = m 1P y is (m) i ; (^ (m) ) = m 1P is m (y i ^ (m) ) ; residals ^r (m) i = jy i (m) j; and their order statistics ^r (m) (i) for i = 1; : : : ; n:. Test whether the residal nearest to the observations in S (m) to an otlier. In ARC the test is based pon bt we sggest to se (1) test = min ^r (m) i6s (m) i =^ (m) ; () test = ^r (m) (m+1) =^(m) : 3. In the next step we either contine with step 1 or stop the algorithm. does not correspond (a) If the test based on the nearest residal does not reject, then de ne and retrn to 1. S (m+1) = fi : ^r (m) i ^r (m) (m+1) g (b) If the test rejects, then set ^m = m and de ne the terminal estimators ^ n = ^ ( ^m) ; ^ n = ^ ( ^m) and the observation set ^S n = S ( ^m) : An important problem is to determine the distribtion of the test statistic and the estimators in case of a sample withot otliers. ARC sggest that the initial estimator ^ cold be chosen as the least trimmed sqares estimator, see Rosseew (1984). This is constrcted by choosing some m > n= and nding ^ (LT S;n;m) = arg min P n i=1 (r() i ) 1 () (r i r () ); (m) r() i = jy i j: (.1)
3 . Comments on the choice of test statistic Comment.1 The motivation for the test statistic ^ (m) () test = ^r (m) (m+1) is that it is the largest of the residals with index in S (m+1) : The rank of the observation ^r (m) (m+1) may, however, not enter S (m) : Comment. The index sets S (m+1) test statistic, (1) test or () test: is constrcted independently of the choice of Comment.3 In general the test statistics (1) test and () test will be di erent. Three reslts follow. 1. the statistic ^ (m) (1) test = min (m) i6s ^r (m) i sggested by ARC is not an order statistic of the residals r (m) i, becase S (m) (m 1) is based on the previos set of residals ^r i :. it holds (1) test () test: Indeed, if S (m) statistics are eqal. If S (m) mst inclde one of the ranks of r (m) (1) ; : : : ; r(m) (m) : is the ranks of r (m) (1) ; : : : ; r(m) (m) then these does not have this form then the complement of S (m) 3. The di erence between the two test statistics relates to the distance j^ (m 1) ^ (m) j and is therefore likely to be nimportant. As a nmerical example consider the data set ( 13; 8; 5; 5; 6; 7; 8): The initial estimator is chosen as the sample average ^ = y = 0: The absolte residals are then (^r i ) n i=1 = (13; 8; 5; 5; 6; 7; 8): An initial index set of size 3 is then the ranks S (3) = (3; 4; 5) pointing at the observations ( 5; 5; 6): Now, in the rst step of the algorithm compte the estimator and absolte residals Then (1) test is based on ^ (3) = ; (^r (3) i ) n i=1 = (15; 10; 7; 3; 4; 5; 6): min ^r (3) i6s (3) i = min(15; 10; 5; 6) = 5; whereas () test is based on the order statistic ^r (3) (4) = 6: Regardless of the choice of test statistic the pdated index set is S (4) = (4; 5; 6; 7) pointing at the observations (5; 6; 7; 8): Note that so the sets S (m) (3; 4; 5) = S (3) 6 S (4) = (4; 5; 6; 7); are not in general increasing. 3
4 3 Some potential reslts In order that the Forward Search Algorithm can be applied with con dence it is important to derive the distribtions of the test statistics and the estimators. It wold also be of interest to see if the sets S (m) are monotone in m: 3.1 Consistency One wold like to have a reslt as the following Potential Theorem 3.1 If the initial estimator (^; ^ ) are consistent for n! 1, and if m = n + O(1) for 0 < < 1; or if m = ^m; then it holds that ^ (m) and ^ (m) have a probability limit. Comment 3.1 We do not have a general proof of sch a reslt, bt some examples are given in 5. One can note here that the di erence between order statistics from i.i.d. observations are of the order of O P (n 1 ); so averaging order statistics which are that close to the initial (consistent estimator) will at most give a deviation from this of the order of O P (n 1 ); and hence shold not distrb the consistency. If we cold nd the probability limit of ^ (m) and ^ (m) ; we cold correct the estimators to give consistent estimators, see ARC 3.3, and Theorem Monotonicity The next reslt is allded to in a nmber of places in ARC, even thogh it is realized that it does not hold withot some conditions. Potential Theorem 3. The sets S (m) are monotone S (m) S (m+1) ; m = ; : : : ; ^m Comment 3. As it stands it is nfortnately false, as seen in the example in Comment.3. If in general we have fond S (m) for some m; and want to add one point to ; call it y m+1 ; then the new average becomes S (m) y (m+1) = y (m) m (y m+1 y (m) ): Ths, if for instance y m+1 > y (m) ; the average is moved p by (y m+1 y (m) )=(m + 1) which is of the order of m 1 n 1 : The distance between order statistics is of the order of n 1 ; so if one of the observations in S (m) is close to the lower bondary of the band de ning S (m) ; then it cold easily happen that it falls otside when the band is moved p by (y m+1 y (m) )=(m + 1). Ths a point can leave the band when another is added even for large m: This event may have small probability, however, so the following reslt cold hold, bt we have no proof. 4
5 Potential Theorem 3.3 The sets S (m) are monotone with probability tending to one, that is for m = n + O(1) as n! 1 then P(S (m) S (m+1) )! 1; or, perhaps, n 1 P n P m=m m 0 1 (m) (S! 1: 0 S (m+1) ) 3.3 The test statistics In ARC it is sggested to nd the distribtion of the test statistics by simlation, and it is arged that it cold be costly measred by compter time. An approximation is sggested in 3.3 of ARC, sing the theory of order statistics, and in a previos paper (Atkinson and Riani 006) another approximation based on order statistics is sggested. One wold like to show that the test statistics are asymptotically normal. Potential Theorem 3.4 If the initial estimators (^; ^ ) are consistent for n! 1, and m = n + O(1); or m = ^m; and if (^ (m) ; ^ (m) ) is consistent, then the test statistic n 1= f(^ (m) ) is asymptotically normally distribted. 1^r (m) (m+1) g Comment We can prove this reslts nder sitable conditions and we have collected these contribtions in the 5, where we otline a general strategy for nding these limit distribtions. 4 The case of known location and scale When the parameters are known, the residals are r i = jy i as r (1) r (m). Then j = j" i j, which we order S (m) = fi : j i j r (m) g; m = ; 3; : : : ; ^m: In this case we clearly have S (m) S (m+1) so Theorem 3. is correct. The empirical distribtion fnction of j" i j = jy i j= is denoted G n () = n 1 nx i=1 1 (j"i j): The order statistics r (m) have the well known relations f 1 r (m) g = f m n G n()g; G n ( 1 r (m) ) = m n ; (4.1) 5
6 which transforms expressions in order statistics into expression involving the empirical distribtion fnction. Moreover, (1) test = () test = 1 r (m) and the distribtion is given by the expression (5) in ARC by applying (4.1) as P( 1 r (m) ) = Pf m n G n()g = P n n j=m G() j f1 G()g n j : (4.) j Ths, the distribtion sggested as an approximation in ARC is in fact the exact distribtion if the parameters were known. In Atkinson and Riani (006) another approximation to the distribtion of the test statistic is sggested. It is arged that if m is proportional to n; then we are essentially working in a trncated distribtion where 100 % have been inclded in the sample. Ths, it is sggested to approximate the distribtion of the test statistic (j) test = 1 r (m) by the distribtion of the largest of m observations, z i say, where z i are drawn from the trncated distribtion G()= for 0. By the same methodology as in 3.3 of ARC we then get the approximate reslt P( 1 r (m) ) P(max im z i ) = 1 G() m = h 1 G() i n ; (4.3) which is sggested as an approximation, which is clearly di erent from (4.). Yet another way of arging is that when we have already sed r (1) ; : : : ; r (m) to constrct S (m) and the estimator y (m) ; then the signi cance of r (m+1) cold be evalated in the conditional distribtion given r (1) ; : : : ; r (m) ; and that is given by P( 1 1 G() r (m+1) jr (1) ; : : : ; r (m) ) = f 1 G( 1 r (m) ) gn m ; 1 r (m) : This distribtion can be interpreted as the distribtion of the smallest observation among n m observations r i = jy i j with density 1 g( 1 r)=f1 G( 1 r (m) )g for r r (m) and cold be sed for a conditional test of the next residal given what has already been sed. 5 The case of nknown location and scale We assme we have n 1= -consistent estimators (^; ^ ) = (; ) + O P (n 1= ) and de ne ^r i = jy i ^j; with order statistics ^r (i) : We de ne ~r = ^r (m) and nx ~ = m 1 y i 1 (jyi ^j~r); (5.1) i=1 nx ~ = m 1 (y i ~) 1 (jyi ^j~r): (5.) i=1 6
7 We nd stochastic expansions and limit distribtions of ~r and the one-step estimators (~; ~ ), which are based pon the m observations closest to the initial estimator ^. When applied to the rst iteration of the algorithm then ^; ^ ; ^r (m) represent the initial estimator ^ and residal ^r (m) along with some sitable variance estimator so ~ = ^ (m0) and ~ = ^ (m0) : When applied to a later iteration of the algorithm then ^; ^; ^r (m) represent ^ (m 1) ; ^ (m 1) (m 1) ; ^r (m) so ~ = ^ (m) and ~ = ^ (m). 5.1 The asymptotic expansions We rst nd an expansion of the test statistic () test = ^ 1 ~r which can be applied to nd the asymptotic distribtion of the test for di erent choices of estimators (^; ^ ), and then give expansions of the one-estimators ~ and ~ ; which show how the estimator is changed from the initial estimators (^; ^ ) to one-step estimators (~; ~ ). The proofs are given in the Appendix. Theorem 5.1 If m = n + O(1) and if (^ ; ^ ) O P (n 1= ) then n 1= (^ 1 ~r ) = 1 f( ) n 1=P n i=1 f1 (j" i j ) g n1= ( ^ 1) + o P (1): Theorem 5. If m = n+o(1) and if (^ ; ^ ) O P (n 1= ) then the estimator ~ de ned in (5.1) satis es n 1= (~ ) = n 1=P n i=1 " i1 (j"i j< ) + f( ) n 1= (^ ) + o P (1): Theorem 5.3 If m = n+o(1) and if (^ ; ^ ) O P (n 1= ) then the estimator ~ de ned in (5.) satis es ~! P 1 ; where = R f()d: We therefore de ne ~ corrected = ~ = : It holds n 1= f ~ corrected 1g = ( ) 1 n 1=P n i=1 (" i 1 )1 (j"i j< ) ( ) 1 ( 4 1 )n 1=P n i=1 f1 (j" i j ) g + o P (1): Comment 5.1 Note that these reslt are all derived for symmetric distribtions. If this assmption is dropped, a bias term will appear in some asymptotic distribtions and terms inclding n 1= (^ ) will appear in Theorems 5.1, Examples We illstrate these reslts by nding the asymptotic distribtion of the test statistic for di erent choices of initial estimators (^; ^ ). 7
8 First, for comparison we give the reslt for the test statistic for known parameters, where we re-discover a classical reslt on the asymptotic distribtion of order statistics, see for instance David (1981, Theorem 9., p. 55). Corollary 5.1 If m = n + O(1) and if ^ = ; ^ = then n 1= ( 1 ~r ) D! N[0; (1 ) ff( )g ]: Proof of Corollary 5.1. In the expansion of Theorem 5.1 the second term drops ot since ^ = : Then apply the Central Limit Theorem to the rst term. Alternatively, the initial estimators cold be chosen as the fll sample estimators ^ = n P 1 n i=1 y i and ^ = n P 1 n i=1 (y i y). This changes the limit distribtion. Corollary 5. If m = n + O(1) and if (^; ^ ) are the fll sample estimators then n 1= (^ 1 ~r ) D! N[0; (1 ) ff( )g + f1 + f( ) g]: Proof of Corollary 5.. Apply Theorem 5.1 by inserting the expansion n 1= ( ^ 1) = n 1=P n i=1 (" i 1) + o P (1) to get that n 1= (^ 1 ~r ) eqals 1 f( ) n 1=P n i=1 f1 (j" i j ) g n 1=P n i=1 (" i 1) + o P (1): This is asymptotically normal with a variance as indicated. Comment 5. In general we get a di erent limit distribtion for the test statistic when the variance is estimated. Note the crios reslt that for the standard normal distribtion we nd = Z " f(")d" = Z f(")d" ["f(")] "= "= = f( ): so the asymptotic variance in Corollary 5. becomes ff( )g (1 ) =; which is less than the variance we get for known parameters. Finally, we shall see what happens when we choose ^ as the least trimmed sqares estimator ^ (LT S;n;m) de ned in (.1). A stochastic expansion of ^ (LT S;n;m) is given by Víšek (006, Theorem 1, p. 15) as n 1= (^ (LT S;n;m) ) = f( ) n 1=P n i=1 " i1 (j"i j ) + o P (1): (5.3) 8
9 Corollary 5.3 If m = n + O(1) and ^ = ^ (LT S;n;m) we nd the expansion n 1= (~ ) = f( ) n 1=P n i=1 " i1 (j"i j ) + o P (1); so that the limit distribtion is n 1= (~ ) D! N[0; f f( )g ]: If F is standard normal then = f( ) so the variance is = : Proof of Corollary 5.3. Insert Víšek s expansion (5.3) in Theorem 5. to get n 1= (~ ) = f1 + f( ) f( ) gn 1=P n i=1 " i1 (j"i j ) + o P (1) = f( ) n 1=P n i=1 " i1 (j"i j ) + o P (1): This converges to a normal distribtion with the variance as indicated. Comment 5.3 Ths, if the initial estimator ^ has the expansion of the least trimmed sqares estimator, then so does the estimator ~: In this sense the least trimmed sqares estimator is a " xed point" in the mapping from the initial estimator to the one-step estimator. ARC do indeed sggest, possibly for other reasons, to start with the least trimmed sqares estimator. A similar reslt holds if the initial estimator is the Hber skip, which has the same expansion and limit distribtion as least trimmed sqares estimator, see Johansen and Nielsen (009). 6 Some nal comments 6.1 The simlation method The above theoretical reslts indicate that the idea of jdging signi cance sing the exact theory of order statistics, seems to be ne if the parameters are known. Bt if the parameters are estimated the (asymptotic) distribtions change, depending on the choice of initial estimator. Ths it wold be very helpfl with some simlations of the algorithm, as it is sed, to check if the asymptotic distribtions can describe the variation of estimators and tests. We have seen that di erent initial estimators give di erent (limit) distribtions. There may also be a problem for very large ; where a di erent asymptotic theory may be needed. 9
10 6. Generalizations The reslts of Johansen and Nielsen (009) cover models with general xed or random regressors, bt also time series regressions, both stationary and non-stationary. The tools to stdy the general theory of empirical processes for residals of sch models is otlined in Engle and Nielsen (009). So it is possible to extend the theory of the Forward Search Algorithm mch beyond what we have indicated in this discssion. 6.3 Algorithms for time series In 7 of ARC it is sggested However, many things remain to be developed in the application of the Forward Search in the time series context, sch as... the constrction of an algorithm which can atomatically distingish among the di erent types of otliers and level shifts. While it cold prove very sefl to develop extensions of the Forward Search to time series one shold bear in mind that the algorithm Atometrics by Doornik (009) bilding on the work of Hoover and Perez (1999) and the PcGets algorithm of Hendry and Krolzig (005) has been developed for this prpose. In any case, it is a bit srprising to see in ARC that the reslts sggested for regression analysis be applied to ozone data where the residals are likely to be atocorrelated. A Proofs of main theorems The proofs exploit Theorem 1.17 of Johansen and Nielsen (009); see also their eqation For (^a; ^b) = O P (n 1= ) and ` = 0; 1; ; it was shown nder reglarity conditions that n 1=P n i=1 "ìf1 (j"i ^bja+^a) 1 (j"i ja)g = a` 1 n 1= (^a `+1 + ^b ` ) + o P (1); (A.1) where, for a symmetric density j = 0 and j+1 = j+1 f( ) for j = 0; 1; : : : Proof of Theorem 5.1. In order to nd the asymptotic distribtion of the test statistic ^ 1 ~r we expand as n 1= ( ~r^ ) = n 1= ^ f( ~r ) ( ^ 1)g: Since ^ = O P (n 1= ) then 1^ 1 = ( ^ 1)= + o P (n 1= ) and ^ 1 = 1 + o P (1): It follows that n 1= ( ~r^ ) = n 1= f( ~r ) ( ^ 1)gf1 + o P (1)g + o P (1): (A.) 10
11 The rst term in (A.), is now shown to be asymptotically normal n 1= ( ~r ) D! Nf0; (1 ) f( ) g: (A.3) The qantile ~r satis es b G n (~r=) = m=n where b G n (r) = n 1 P n i=1 1 (j" i 1 (^ )j<r) is an empirical distribtion fnction; see (4.1). Ths it holds P n def = Pfn 1= ( ~r ) zg = Pf ~r + n 1= zg where G n = n 1= f b G n ( + n 1= z) This can be expanded as = Pf m n b G n ( + n 1= z)g = P(0 G n ); m n g = n 1=P n i=1 [1 fj" i 1 (^ )j< +n 1= zg m n ]: G n = n 1=P n i=1 f1 (j" i j< ) g + n 1=P n i=1 [1 fj" i 1 (^ )j< +n 1= zg 1 (j"i j< )] + n 1= ( m n ): (A.4) Here, the rst term is asymptotically Nf0; (1 )g: Applying (A.1) with ` = 0; a = ; ^a = n 1= z; ^b = 1 (^ ) the second term is f( )z + o P (1): The third term vanishes by assmption. Ths, G n is asymptotically Nff( )z; (1 )g: In particlar, it follows that P n def = Pfn 1= ( ~r ) zg! f f( )z p (1 ) g: In trn, ~r is asymptotically normal as stated in (A.3). Next, an expansion for ~r is derived. Since b G n (~r=) = m=n then n 1= f b G n ( ~r ) g = n1= ( m n ) = o P (1); by the assmption to m: Expanding the left hand term as in (A.4) then gives n 1=P n i=1 f1 (j" i j< ) g + n 1=P n i=1 [1 fj" i 1 (^ )j< 1^r (m) g 1 (j"i j< )] = o P (1): Applying (A.1) with ` = 0; a = ; ^a = 1 ~r ; ^b = 1 (^ ) then shows n 1=P n i=1 f1 (j" i j< ) g + f( )n 1= ( 1 ~r ) = o P (1): Solving this eqation for n 1= ( 1 ~r into (A.). ) gives the desired expansion when inserted 11
12 Proof of Theorem 5.. The estimator ~ satis es ~ = m 1P n i=1 (y i )1 (jyi ^j~r) = n m n 1P n i=1 (y i )1 (j"i 1 (^ )j< 1 ~r): Adding and sbtracting 1 (j"i j< ) and sing y i = + " i it holds m n 1 n 1= (~ ) = n 1=P n i=1 " i1 (j"i j< ) + n 1=P n i=1 " if1 (j"i 1 (^ )j< 1 ~r) 1 (j"i j< )g: Here, m 1 n! by assmption. The rst term on the right converges in distribtion by the central limit theorem. Applying (A.1) with ` = 1; a = ; ^a = 1 ~r ; ^b = 1 (^ ) the second term is f( ) 1 n 1= (^ ) + o P (1): Inserting these reslts gives the desired expansion. Proof of Theorem 5.3. The estimator ~ satis es ~ 1 = m 1P n i=1 f (y i ^) 1 g1 (jyi ^j~r): Using y i ^ = " i (^ ) then (y i ^) 1 = (" i 1 ) + (^ ) 1 (^ )" i ; while 1 (jyi ^j~r) = 1 fj"i 1 (^ )j< 1 ~rg as before. We de ne the fnctions h () = 1 ; h 1 () = "; and h 0 () = 1; and, for ` = 0; 1; ; S` = n 1=P n i=1 h`(" i )1 fj"i 1 (^ )j< 1 ~rg = n 1=P n i=1 h`(" i )1 fj"i j< g + R`; R` = n 1=P n i=1 h`(" i )[1 fj"i 1 (^ )j< 1 ~rg 1 fj"i j< g]; so that m n n1= f ~ 1 g = S + (^ ) S 0 1 (^ )S 1 : (A.5) For S the rst term converges in distribtion by the central limit theorem. Applying (A.1) with ` = and ` = 0, a = ; ^a = 1 ~r ; ^b = 1 (^ ); the second term is R = f( )( 4 1 )n 1= ( 1 ~r ) + o P (1): It follows that S = n 1=P n i=1 (" i S = n 1=P n i=1 (" i )1 fj"i j< g + f( )( 4 )n 1= ( 1 ~r ) + o P (1): Inserting the expression for ~r from Theorem 5.1 with ^ = then shows )1 fj"i j< g ( 4 )n 1=P n i=1 f1 (j" i j ) g + o P (1): The terms S 0 and S 1 are O P (1) by a similar argment. These terms, are however, pre-mltiplied by vanishing terms in that = O P (n 1= ) by assmption. Inserting these reslts in (A.5) noting that m 1 n! by assmption shows n 1= f ~ 1 g = S + o P (1): With ~ corrected = ~ = the left hand side becomes n 1= f ~ corrected 1g and the desired reslt follows. 1
13 References Atkinson, A.C. and Riani, M. (006) Distribtion theory and simlations for tests of otliers in regression. Jornal of Comptational and Graphical Statistics 15, Atkinson, A.C., Riani, M. and Ceroli, A. (009) The forward search: Theory and data analysis. Discssion paper Jornal of Korean Statistical Society. David, H.A. (1981) Order Statistics. nd ed. New York: Wiley. Doornik, J. (009) Atometrics. In Castle, J.L. and Shephard, N. (eds.) The Methodology and Practice of Econometrics: A Festschrift in Honor of David F. Hendry, pp Oxford: Oxford University Press. Engler, E. and Nielsen, B. (009) The empirical process of atoregressive residals. Econometrics Jornal 1, Hendry, D.F. and Krolzig, H.-M. (005) The properties of atomatic Gets modelling. Economic Jornal 115, C3 C61. Hoover, K.D. and Perez, S. J. (1999) Data mining reconsidered: Encompassing and the general to speci c approach to speci cation search. Econometric Jornal, Johansen, S. and Nielsen, B. (009) An analysis of the indicator satration estimator as a robst regression estimator. In Castle, J.L. and Shephard, N. (eds.) The Methodology and Practice of Econometrics: A Festschrift in Honor of David F. Hendry, pp Oxford: Oxford University Press. Kol, H.L. (00) Weighted Empirical Processes in Dynamic Nonlinear Models, nd edition. New York: Springer. Rosseew, P.J. (1984) Least median of sqares regression. Jornal of the American Statistical Association 79, Víšek, J.Á. (006) The least trimmed sqares. Part III: Asymptotic normality. Kybernetika 4,
Discussion Papers Department of Economics University of Copenhagen
Discssion Papers Department of Economics University of Copenhagen No. 10-06 Discssion of The Forward Search: Theory and Data Analysis, by Anthony C. Atkinson, Marco Riani, and Andrea Ceroli Søren Johansen,
More informationSources of Non Stationarity in the Semivariogram
Sorces of Non Stationarity in the Semivariogram Migel A. Cba and Oy Leangthong Traditional ncertainty characterization techniqes sch as Simple Kriging or Seqential Gassian Simlation rely on stationary
More informationi=1 y i 1fd i = dg= P N i=1 1fd i = dg.
ECOOMETRICS II (ECO 240S) University of Toronto. Department of Economics. Winter 208 Instrctor: Victor Agirregabiria SOLUTIO TO FIAL EXAM Tesday, April 0, 208. From 9:00am-2:00pm (3 hors) ISTRUCTIOS: -
More informationDiscontinuous Fluctuation Distribution for Time-Dependent Problems
Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation
More informationThe Linear Quadratic Regulator
10 The Linear Qadratic Reglator 10.1 Problem formlation This chapter concerns optimal control of dynamical systems. Most of this development concerns linear models with a particlarly simple notion of optimality.
More informationSection 7.4: Integration of Rational Functions by Partial Fractions
Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to
More informationWe automate the bivariate change-of-variables technique for bivariate continuous random variables with
INFORMS Jornal on Compting Vol. 4, No., Winter 0, pp. 9 ISSN 09-9856 (print) ISSN 56-558 (online) http://dx.doi.org/0.87/ijoc.046 0 INFORMS Atomating Biariate Transformations Jeff X. Yang, John H. Drew,
More informationPulses on a Struck String
8.03 at ESG Spplemental Notes Plses on a Strck String These notes investigate specific eamples of transverse motion on a stretched string in cases where the string is at some time ndisplaced, bt with a
More informationThe Dual of the Maximum Likelihood Method
Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All
More informationSubcritical bifurcation to innitely many rotating waves. Arnd Scheel. Freie Universitat Berlin. Arnimallee Berlin, Germany
Sbcritical bifrcation to innitely many rotating waves Arnd Scheel Institt fr Mathematik I Freie Universitat Berlin Arnimallee 2-6 14195 Berlin, Germany 1 Abstract We consider the eqation 00 + 1 r 0 k2
More informationA New Approach to Direct Sequential Simulation that Accounts for the Proportional Effect: Direct Lognormal Simulation
A ew Approach to Direct eqential imlation that Acconts for the Proportional ffect: Direct ognormal imlation John Manchk, Oy eangthong and Clayton Detsch Department of Civil & nvironmental ngineering University
More informationBLOOM S TAXONOMY. Following Bloom s Taxonomy to Assess Students
BLOOM S TAXONOMY Topic Following Bloom s Taonomy to Assess Stdents Smmary A handot for stdents to eplain Bloom s taonomy that is sed for item writing and test constrction to test stdents to see if they
More informationLecture Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2
BIJU PATNAIK UNIVERSITY OF TECHNOLOGY, ODISHA Lectre Notes On THEORY OF COMPUTATION MODULE - 2 UNIT - 2 Prepared by, Dr. Sbhend Kmar Rath, BPUT, Odisha. Tring Machine- Miscellany UNIT 2 TURING MACHINE
More informationAn Investigation into Estimating Type B Degrees of Freedom
An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information
More informationSetting The K Value And Polarization Mode Of The Delta Undulator
LCLS-TN-4- Setting The Vale And Polarization Mode Of The Delta Undlator Zachary Wolf, Heinz-Dieter Nhn SLAC September 4, 04 Abstract This note provides the details for setting the longitdinal positions
More informationA Single Species in One Spatial Dimension
Lectre 6 A Single Species in One Spatial Dimension Reading: Material similar to that in this section of the corse appears in Sections 1. and 13.5 of James D. Mrray (), Mathematical Biology I: An introction,
More informationA Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time
A Theory of Markovian Time Inconsistent Stochastic Control in Discrete Time Tomas Björk Department of Finance, Stockholm School of Economics tomas.bjork@hhs.se Agatha Mrgoci Department of Economics Aarhs
More informationA Model-Free Adaptive Control of Pulsed GTAW
A Model-Free Adaptive Control of Plsed GTAW F.L. Lv 1, S.B. Chen 1, and S.W. Dai 1 Institte of Welding Technology, Shanghai Jiao Tong University, Shanghai 00030, P.R. China Department of Atomatic Control,
More information3.1 The Basic Two-Level Model - The Formulas
CHAPTER 3 3 THE BASIC MULTILEVEL MODEL AND EXTENSIONS In the previos Chapter we introdced a nmber of models and we cleared ot the advantages of Mltilevel Models in the analysis of hierarchically nested
More informationTechnical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty
Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte
More informationA Macroscopic Traffic Data Assimilation Framework Based on Fourier-Galerkin Method and Minimax Estimation
A Macroscopic Traffic Data Assimilation Framework Based on Forier-Galerkin Method and Minima Estimation Tigran T. Tchrakian and Sergiy Zhk Abstract In this paper, we propose a new framework for macroscopic
More informationUniversal Scheme for Optimal Search and Stop
Universal Scheme for Optimal Search and Stop Sirin Nitinawarat Qalcomm Technologies, Inc. 5775 Morehose Drive San Diego, CA 92121, USA Email: sirin.nitinawarat@gmail.com Vengopal V. Veeravalli Coordinated
More informationSensitivity Analysis in Bayesian Networks: From Single to Multiple Parameters
Sensitivity Analysis in Bayesian Networks: From Single to Mltiple Parameters Hei Chan and Adnan Darwiche Compter Science Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwiche}@cs.cla.ed
More informationStability of Model Predictive Control using Markov Chain Monte Carlo Optimisation
Stability of Model Predictive Control sing Markov Chain Monte Carlo Optimisation Elilini Siva, Pal Golart, Jan Maciejowski and Nikolas Kantas Abstract We apply stochastic Lyapnov theory to perform stability
More information4 Exact laminar boundary layer solutions
4 Eact laminar bondary layer soltions 4.1 Bondary layer on a flat plate (Blasis 1908 In Sec. 3, we derived the bondary layer eqations for 2D incompressible flow of constant viscosity past a weakly crved
More informationCHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK. Wassim Jouini and Christophe Moy
CHANNEL SELECTION WITH RAYLEIGH FADING: A MULTI-ARMED BANDIT FRAMEWORK Wassim Joini and Christophe Moy SUPELEC, IETR, SCEE, Avene de la Bolaie, CS 47601, 5576 Cesson Sévigné, France. INSERM U96 - IFR140-
More informationE ect Of Quadrant Bow On Delta Undulator Phase Errors
LCLS-TN-15-1 E ect Of Qadrant Bow On Delta Undlator Phase Errors Zachary Wolf SLAC Febrary 18, 015 Abstract The Delta ndlator qadrants are tned individally and are then assembled to make the tned ndlator.
More informationFRÉCHET KERNELS AND THE ADJOINT METHOD
PART II FRÉCHET KERNES AND THE ADJOINT METHOD 1. Setp of the tomographic problem: Why gradients? 2. The adjoint method 3. Practical 4. Special topics (sorce imaging and time reversal) Setp of the tomographic
More informationGeneralized Jinc functions and their application to focusing and diffraction of circular apertures
Qing Cao Vol. 20, No. 4/April 2003/J. Opt. Soc. Am. A 66 Generalized Jinc fnctions and their application to focsing and diffraction of circlar apertres Qing Cao Optische Nachrichtentechnik, FernUniversität
More informationDecoder Error Probability of MRD Codes
Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehigh.ed Zhiyan Yan Department of Electrical
More informationModelling by Differential Equations from Properties of Phenomenon to its Investigation
Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University
More informationImproved Jive Estimators for Overidenti ed Linear Models with and without Heteroskedasticity
Improved Jive Estimators for Overidenti ed Linear Models with and withot Heteroskedasticity Daniel A. Ackerberg and Pal J. Deverex` Febrary 11, 2003 Abstract This paper examines and extends work on estimators
More informationDecoder Error Probability of MRD Codes
Decoder Error Probability of MRD Codes Maximilien Gadolea Department of Electrical and Compter Engineering Lehigh University Bethlehem, PA 18015 USA E-mail: magc@lehighed Zhiyan Yan Department of Electrical
More information1 The space of linear transformations from R n to R m :
Math 540 Spring 20 Notes #4 Higher deriaties, Taylor s theorem The space of linear transformations from R n to R m We hae discssed linear transformations mapping R n to R m We can add sch linear transformations
More informationDecision Oriented Bayesian Design of Experiments
Decision Oriented Bayesian Design of Experiments Farminder S. Anand*, Jay H. Lee**, Matthew J. Realff*** *School of Chemical & Biomoleclar Engineering Georgia Institte of echnology, Atlanta, GA 3332 USA
More informationQUANTILE ESTIMATION IN SUCCESSIVE SAMPLING
Jornal of the Korean Statistical Society 2007, 36: 4, pp 543 556 QUANTILE ESTIMATION IN SUCCESSIVE SAMPLING Hosila P. Singh 1, Ritesh Tailor 2, Sarjinder Singh 3 and Jong-Min Kim 4 Abstract In sccessive
More informationPart II. Martingale measres and their constrctions 1. The \First" and the \Second" fndamental theorems show clearly how \mar tingale measres" are impo
Albert N. Shiryaev (Stelov Mathematical Institte and Moscow State University) ESSENTIALS of the ARBITRAGE THEORY Part I. Basic notions and theorems of the \Arbitrage Theory" Part II. Martingale measres
More informationStudy on the impulsive pressure of tank oscillating by force towards multiple degrees of freedom
EPJ Web of Conferences 80, 0034 (08) EFM 07 Stdy on the implsive pressre of tank oscillating by force towards mltiple degrees of freedom Shigeyki Hibi,* The ational Defense Academy, Department of Mechanical
More informationON THE SHAPES OF BILATERAL GAMMA DENSITIES
ON THE SHAPES OF BILATERAL GAMMA DENSITIES UWE KÜCHLER, STEFAN TAPPE Abstract. We investigate the for parameter family of bilateral Gamma distribtions. The goal of this paper is to provide a thorogh treatment
More informationarxiv: v1 [physics.flu-dyn] 11 Mar 2011
arxiv:1103.45v1 [physics.fl-dyn 11 Mar 011 Interaction of a magnetic dipole with a slowly moving electrically condcting plate Evgeny V. Votyakov Comptational Science Laboratory UCY-CompSci, Department
More informationAnalytical Value-at-Risk and Expected Shortfall under Regime Switching *
Working Paper 9-4 Departamento de Economía Economic Series 14) Universidad Carlos III de Madrid March 9 Calle Madrid, 16 893 Getafe Spain) Fax 34) 91649875 Analytical Vale-at-Risk and Expected Shortfall
More informationThe Cryptanalysis of a New Public-Key Cryptosystem based on Modular Knapsacks
The Cryptanalysis of a New Pblic-Key Cryptosystem based on Modlar Knapsacks Yeow Meng Chee Antoine Jox National Compter Systems DMI-GRECC Center for Information Technology 45 re d Ulm 73 Science Park Drive,
More informationNonlinear parametric optimization using cylindrical algebraic decomposition
Proceedings of the 44th IEEE Conference on Decision and Control, and the Eropean Control Conference 2005 Seville, Spain, December 12-15, 2005 TC08.5 Nonlinear parametric optimization sing cylindrical algebraic
More informationA Characterization of Interventional Distributions in Semi-Markovian Causal Models
A Characterization of Interventional Distribtions in Semi-Markovian Casal Models Jin Tian and Changsng Kang Department of Compter Science Iowa State University Ames, IA 50011 {jtian, cskang}@cs.iastate.ed
More informationLecture: Corporate Income Tax
Lectre: Corporate Income Tax Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak atoregressive cash
More informationFigure 1 Probability density function of Wedge copula for c = (best fit to Nominal skew of DRAM case study).
Wedge Copla This docment explains the constrction and properties o a particlar geometrical copla sed to it dependency data rom the edram case stdy done at Portland State University. The probability density
More informationUNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL
8th International DAAAM Baltic Conference "INDUSTRIAL ENGINEERING - 19-1 April 01, Tallinn, Estonia UNCERTAINTY FOCUSED STRENGTH ANALYSIS MODEL Põdra, P. & Laaneots, R. Abstract: Strength analysis is a
More informationConvergence analysis of ant colony learning
Delft University of Technology Delft Center for Systems and Control Technical report 11-012 Convergence analysis of ant colony learning J van Ast R Babška and B De Schtter If yo want to cite this report
More informationStep-Size Bounds Analysis of the Generalized Multidelay Adaptive Filter
WCE 007 Jly - 4 007 London UK Step-Size onds Analysis of the Generalized Mltidelay Adaptive Filter Jnghsi Lee and Hs Chang Hang Abstract In this paper we analyze the bonds of the fixed common step-size
More informationFEA Solution Procedure
EA Soltion Procedre (demonstrated with a -D bar element problem) EA Procedre for Static Analysis. Prepare the E model a. discretize (mesh) the strctre b. prescribe loads c. prescribe spports. Perform calclations
More informationEssentials of optimal control theory in ECON 4140
Essentials of optimal control theory in ECON 4140 Things yo need to know (and a detail yo need not care abot). A few words abot dynamic optimization in general. Dynamic optimization can be thoght of as
More informationClassify by number of ports and examine the possible structures that result. Using only one-port elements, no more than two elements can be assembled.
Jnction elements in network models. Classify by nmber of ports and examine the possible strctres that reslt. Using only one-port elements, no more than two elements can be assembled. Combining two two-ports
More informationFinite Difference Method of Modelling Groundwater Flow
Jornal of Water Resorce and Protection, 20, 3, 92-98 doi:0.4236/warp.20.33025 Pblished Online March 20 (http://www.scirp.org/ornal/warp) Finite Difference Method of Modelling Grondwater Flow Abstract Magns.
More informationTed Pedersen. Southern Methodist University. large sample assumptions implicit in traditional goodness
Appears in the Proceedings of the Soth-Central SAS Users Grop Conference (SCSUG-96), Astin, TX, Oct 27-29, 1996 Fishing for Exactness Ted Pedersen Department of Compter Science & Engineering Sothern Methodist
More informationLarge Panel Test of Factor Pricing Models
Large Panel est of Factor Pricing Models Jianqing Fan, Yan Liao and Jiawei Yao Department of Operations Research and Financial Engineering, Princeton University Bendheim Center for Finance, Princeton University
More informationThe spreading residue harmonic balance method for nonlinear vibration of an electrostatically actuated microbeam
J.L. Pan W.Y. Zh Nonlinear Sci. Lett. Vol.8 No. pp.- September The spreading reside harmonic balance method for nonlinear vibration of an electrostatically actated microbeam J. L. Pan W. Y. Zh * College
More informationFOUNTAIN codes [3], [4] provide an efficient solution
Inactivation Decoding of LT and Raptor Codes: Analysis and Code Design Francisco Lázaro, Stdent Member, IEEE, Gianligi Liva, Senior Member, IEEE, Gerhard Bach, Fellow, IEEE arxiv:176.5814v1 [cs.it 19 Jn
More informationLecture: Corporate Income Tax - Unlevered firms
Lectre: Corporate Income Tax - Unlevered firms Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak
More informationComputational Geosciences 2 (1998) 1, 23-36
A STUDY OF THE MODELLING ERROR IN TWO OPERATOR SPLITTING ALGORITHMS FOR POROUS MEDIA FLOW K. BRUSDAL, H. K. DAHLE, K. HVISTENDAHL KARLSEN, T. MANNSETH Comptational Geosciences 2 (998), 23-36 Abstract.
More informationLecture Notes: Finite Element Analysis, J.E. Akin, Rice University
9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)
More informationOptimal Control of a Heterogeneous Two Server System with Consideration for Power and Performance
Optimal Control of a Heterogeneos Two Server System with Consideration for Power and Performance by Jiazheng Li A thesis presented to the University of Waterloo in flfilment of the thesis reqirement for
More informationAdmissibility under the LINEX loss function in non-regular case. Hidekazu Tanaka. Received November 5, 2009; revised September 2, 2010
Scientiae Mathematicae Japonicae Online, e-2012, 427 434 427 Admissibility nder the LINEX loss fnction in non-reglar case Hidekaz Tanaka Received November 5, 2009; revised September 2, 2010 Abstract. In
More informationMove Blocking Strategies in Receding Horizon Control
Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it
More informationPropagation of measurement uncertainty in spatial characterisation of recreational fishing catch rates using logistic transform indicator kriging
st International Congress on Modelling and Simlation, Gold Coast, Astralia, 9 Nov to 4 Dec 05 www.mssan.org.a/modsim05 Propagation of measrement ncertainty in spatial characterisation of recreational fishing
More informationJ.A. BURNS AND B.B. KING redced order controllers sensors/actators. The kernels of these integral representations are called fnctional gains. In [4],
Jornal of Mathematical Systems, Estimation, Control Vol. 8, No. 2, 1998, pp. 1{12 c 1998 Birkhaser-Boston A Note on the Mathematical Modelling of Damped Second Order Systems John A. Brns y Belinda B. King
More informationOptimal search: a practical interpretation of information-driven sensor management
Optimal search: a practical interpretation of information-driven sensor management Fotios Katsilieris, Yvo Boers and Hans Driessen Thales Nederland B.V. Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationUsing average control chart with subgroups to monitor turnaround time
Jornal of Scientific KAO: USING & Indstrial AVERAGE Research CONTROL CHART WITH SUBGROUPS TO MONITOR TURNAROUND TIME Vol. 71, Janary 1, pp. 7-5 7 Using average control chart with sbgrops to monitor trnarond
More informationCharacterizations of probability distributions via bivariate regression of record values
Metrika (2008) 68:51 64 DOI 10.1007/s00184-007-0142-7 Characterizations of probability distribtions via bivariate regression of record vales George P. Yanev M. Ahsanllah M. I. Beg Received: 4 October 2006
More informationOptimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications
Optimization via the Hamilton-Jacobi-Bellman Method: Theory and Applications Navin Khaneja lectre notes taken by Christiane Koch Jne 24, 29 1 Variation yields a classical Hamiltonian system Sppose that
More informationApplying Fuzzy Set Approach into Achieving Quality Improvement for Qualitative Quality Response
Proceedings of the 007 WSES International Conference on Compter Engineering and pplications, Gold Coast, stralia, Janary 17-19, 007 5 pplying Fzzy Set pproach into chieving Qality Improvement for Qalitative
More informationOptimal Control, Statistics and Path Planning
PERGAMON Mathematical and Compter Modelling 33 (21) 237 253 www.elsevier.nl/locate/mcm Optimal Control, Statistics and Path Planning C. F. Martin and Shan Sn Department of Mathematics and Statistics Texas
More informationSTEP Support Programme. STEP III Hyperbolic Functions: Solutions
STEP Spport Programme STEP III Hyperbolic Fnctions: Soltions Start by sing the sbstittion t cosh x. This gives: sinh x cosh a cosh x cosh a sinh x t sinh x dt t dt t + ln t ln t + ln cosh a ln ln cosh
More informationarxiv:quant-ph/ v4 14 May 2003
Phase-transition-like Behavior of Qantm Games arxiv:qant-ph/0111138v4 14 May 2003 Jiangfeng D Department of Modern Physics, University of Science and Technology of China, Hefei, 230027, People s Repblic
More informationOn the scaling ranges of detrended fluctuation analysis for long-memory correlated short series of data
On the scaling ranges of detrended flctation analysis for long-memory correlated short series of data arxiv:126.17v1 [physics.data-an] 5 Jn 212 Darisz Grech (1) and Zygmnt Mazr (2) (1) Institte of Theoretical
More informationA Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model
entropy Article A Characterization of the Domain of Beta-Divergence and Its Connection to Bregman Variational Model Hyenkyn Woo School of Liberal Arts, Korea University of Technology and Edcation, Cheonan
More informationOPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIELD OF A POLYHEDRAL BODY WITH LINEARLY INCREASING DENSITY 1
OPTIMUM EXPRESSION FOR COMPUTATION OF THE GRAVITY FIEL OF A POLYHERAL BOY WITH LINEARLY INCREASING ENSITY 1 V. POHÁNKA2 Abstract The formla for the comptation of the gravity field of a polyhedral body
More informationSareban: Evaluation of Three Common Algorithms for Structure Active Control
Engineering, Technology & Applied Science Research Vol. 7, No. 3, 2017, 1638-1646 1638 Evalation of Three Common Algorithms for Strctre Active Control Mohammad Sareban Department of Civil Engineering Shahrood
More informationA Note on Irreducible Polynomials and Identity Testing
A Note on Irrecible Polynomials an Ientity Testing Chanan Saha Department of Compter Science an Engineering Inian Institte of Technology Kanpr Abstract We show that, given a finite fiel F q an an integer
More informationBridging the Gap Between Multigrid, Hierarchical, and Receding-Horizon Control
Bridging the Gap Between Mltigrid, Hierarchical, and Receding-Horizon Control Victor M. Zavala Mathematics and Compter Science Division Argonne National Laboratory Argonne, IL 60439 USA (e-mail: vzavala@mcs.anl.gov).
More informationLecture 17 Errors in Matlab s Turbulence PSD and Shaping Filter Expressions
Lectre 7 Errors in Matlab s Trblence PSD and Shaping Filter Expressions b Peter J Sherman /7/7 [prepared for AERE 355 class] In this brief note we will show that the trblence power spectral densities (psds)
More informationRemarks on strongly convex stochastic processes
Aeqat. Math. 86 (01), 91 98 c The Athor(s) 01. This article is pblished with open access at Springerlink.com 0001-9054/1/010091-8 pblished online November 7, 01 DOI 10.1007/s00010-01-016-9 Aeqationes Mathematicae
More informationThe Replenishment Policy for an Inventory System with a Fixed Ordering Cost and a Proportional Penalty Cost under Poisson Arrival Demands
Scientiae Mathematicae Japonicae Online, e-211, 161 167 161 The Replenishment Policy for an Inventory System with a Fixed Ordering Cost and a Proportional Penalty Cost nder Poisson Arrival Demands Hitoshi
More informationJoint Transfer of Energy and Information in a Two-hop Relay Channel
Joint Transfer of Energy and Information in a Two-hop Relay Channel Ali H. Abdollahi Bafghi, Mahtab Mirmohseni, and Mohammad Reza Aref Information Systems and Secrity Lab (ISSL Department of Electrical
More informationDEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS (H/U < 20) AND CONSEQUENCES ON CRITICALITY SAFETY
DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS ( < 20) AND CONSEQUENCES ON CRITICALITY SAFETY N. Leclaire, S. Evo, I.R.S.N., France Introdction In criticality stdies, the blk density
More informationSecond-Order Wave Equation
Second-Order Wave Eqation A. Salih Department of Aerospace Engineering Indian Institte of Space Science and Technology, Thirvananthapram 3 December 016 1 Introdction The classical wave eqation is a second-order
More informationActive Flux Schemes for Advection Diffusion
AIAA Aviation - Jne, Dallas, TX nd AIAA Comptational Flid Dynamics Conference AIAA - Active Fl Schemes for Advection Diffsion Hiroaki Nishikawa National Institte of Aerospace, Hampton, VA 3, USA Downloaded
More informationSimplified Identification Scheme for Structures on a Flexible Base
Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles
More informationImprecise Continuous-Time Markov Chains
Imprecise Continos-Time Markov Chains Thomas Krak *, Jasper De Bock, and Arno Siebes t.e.krak@.nl, a.p.j.m.siebes@.nl Utrecht University, Department of Information and Compting Sciences, Princetonplein
More informationRELIABILITY ASPECTS OF PROPORTIONAL MEAN RESIDUAL LIFE MODEL USING QUANTILE FUNC- TIONS
RELIABILITY ASPECTS OF PROPORTIONAL MEAN RESIDUAL LIFE MODEL USING QUANTILE FUNC- TIONS Athors: N.UNNIKRISHNAN NAIR Department of Statistics, Cochin University of Science Technology, Cochin, Kerala, INDIA
More informationOn relative errors of floating-point operations: optimal bounds and applications
On relative errors of floating-point operations: optimal bonds and applications Clade-Pierre Jeannerod, Siegfried M. Rmp To cite this version: Clade-Pierre Jeannerod, Siegfried M. Rmp. On relative errors
More informationElectron Phase Slip in an Undulator with Dipole Field and BPM Errors
CS-T--14 October 3, Electron Phase Slip in an Undlator with Dipole Field and BPM Errors Pal Emma SAC ABSTRACT A statistical analysis of a corrected electron trajectory throgh a planar ndlator is sed to
More informationVIBRATION MEASUREMENT UNCERTAINTY AND RELIABILITY DIAGNOSTICS RESULTS IN ROTATING SYSTEMS
VIBRATIO MEASUREMET UCERTAITY AD RELIABILITY DIAGOSTICS RESULTS I ROTATIG SYSTEMS. Introdction M. Eidkevicite, V. Volkovas anas University of Technology, Lithania The rotating machinery technical state
More informationA Regulator for Continuous Sedimentation in Ideal Clarifier-Thickener Units
A Reglator for Continos Sedimentation in Ideal Clarifier-Thickener Units STEFAN DIEHL Centre for Mathematical Sciences, Lnd University, P.O. Box, SE- Lnd, Sweden e-mail: diehl@maths.lth.se) Abstract. The
More informationComputational Fluid Dynamics Simulation and Wind Tunnel Testing on Microlight Model
Comptational Flid Dynamics Simlation and Wind Tnnel Testing on Microlight Model Iskandar Shah Bin Ishak Department of Aeronatics and Atomotive, Universiti Teknologi Malaysia T.M. Kit Universiti Teknologi
More informationAPPLICATION OF MICROTREMOR MEASUREMENTS TO EARTHQUAKE ENGINEERING
170 APPLICATION OF MICROTREMOR MEASUREMENTS TO EARTHQUAKE ENGINEERING I. M. Parton* and P. W. Taylor* INTRODUCTION Two previos papers pblished in this Blletin (Refs 1, 2) described the methods developed
More informationNon-Linear Models Topic Outlines
Non-Linear Models Mohammad Ehsanl Karim Institte of Statistical Research and training; University of Dhaka, Dhaka, Bangladesh Topic Otlines Intrinsically Linear Regression Models Intrinsically
More informationEXCITATION RATE COEFFICIENTS OF MOLYBDENUM ATOM AND IONS IN ASTROPHYSICAL PLASMA AS A FUNCTION OF ELECTRON TEMPERATURE
EXCITATION RATE COEFFICIENTS OF MOLYBDENUM ATOM AND IONS IN ASTROPHYSICAL PLASMA AS A FUNCTION OF ELECTRON TEMPERATURE A.N. Jadhav Department of Electronics, Yeshwant Mahavidyalaya, Ned. Affiliated to
More informationQueueing analysis of service deferrals for load management in power systems
Qeeing analysis of service deferrals for load management in power systems Andrés Ferragt and Fernando Paganini Universidad ORT Urgay Abstract With the advent of renewable sorces and Smart- Grid deployments,
More informationCalculations involving a single random variable (SRV)
Calclations involving a single random variable (SRV) Example of Bearing Capacity q φ = 0 µ σ c c = 100kN/m = 50kN/m ndrained shear strength parameters What is the relationship between the Factor of Safety
More informationInternational Journal of Physical and Mathematical Sciences journal homepage:
64 International Jornal of Physical and Mathematical Sciences Vol 2, No 1 (2011) ISSN: 2010-1791 International Jornal of Physical and Mathematical Sciences jornal homepage: http://icoci.org/ijpms PRELIMINARY
More information