Robust estimation of the external drift and the variogram of spatial data
|
|
- Shanon Crawford
- 5 years ago
- Views:
Transcription
1 Researh Colletion Conferene Paper Robust estimation of the external drift and the variogram of spatial data Author(s): Künsh, Hansruedi; Papritz, Andreas Jürg; Shwierz, Cornelia; Stahel, Werner A. Publiation Date: 2013 Permanent Link: Rights / Liense: In Copyright - Non-Commerial Use Permitted This page was generated automatially upon download from the ETH Zurih Researh Colletion. For more information please onsult the Terms of use. ETH Library
2 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August Robust estimation of the external drift and the variogram of spatial data Künsh, Hans R. ETH Zurih, Seminar for Statistis Rämistrasse Zürih, Switzerland hans-ruedi.kuensh@stat.math.ethz.h Papritz, Andreas ETH Zurih, Institute of Terrestrial Eosystems Universitätstrasse Zürih, Switzerland papritz@env.ethz.h Shwierz, Cornelia ETH Zurih, Seminar for Statistis & Institute of Terrestrial Eosystems Rämistrasse Zürih, Switzerland ornelia.shwierz@stat.math.ethz.h Stahel, Werner A. ETH Zurih, Seminar for Statistis Rämistrasse Zürih, Switzerland stahel@stat.math.ethz.h Introdution Most of the geostatistial software tools, available either ommerially in GIS and in statistial software pakages or as open soure software (e.g. R, R Development Core Team, 2011) rely on non-robust algorithms. This is unfortunate, beause outlying observations are rather the rule than the exeption, in partiular in environmental data sets. Outlying observations may results from errors (e.g. in data transription) or from loal perturbations in the proesses that are responsible for a given pattern of spatial variation. As an example, the spatial distribution of some trae metal in the soils of a region may be distorted by emissions of loal anthropogeni soures. Outliers affet the modelling of the large-sale spatial variation, the so-alled external drift or trend, the estimation of the spatial dependene of the residual variation and the preditions by kriging. Identifying outliers manually is umbersome and requires expertise. A better approah is to use robust algorithms that prevent automatially that outlying observations have undue influene. Former studies on robust geostatistis foused on robust estimation of the sample variogram (f. Lark, 2000, for a review), kriging without external drift (Hawkins and Cressie, 1984) and on robust estimation of parameters of a linear model for the external drift (Militino and Ugarte, 1997). Furthermore, Rihardson and Welsh (1995) proposed a robustified version of (restrited) maximum likelihood ([REML) estimation for the variane omponents of a linear mixed model, whih was later used by Marhant and Lark (2007) for robust REML estimation of the variogram. We propose here a novel method for robust REML estimation of the variogram of a Gaussian random field that is possibly ontaminated by independent errors from a long-tailed distribution. Besides robust estimates of the parameters of the external drift and of the variogram, the method also provides robustified kriging
3 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August preditions. Theory We use the model Y (s) = x(s) T β+z(s)+ε(s) where x(s) T β is the external drift, Z(s) is a stationary Gaussian field with mean zero and ovariane R(h; σ0 2, α), and ε(s) is an independently distributed error with a probability density funtion proportional to 1/σ exp( ρ(ε/σ)) so that ρ(x) = x 2 /2 gives the ase where ε(s) has a Gaussian distribution. σ 2 is the squared sale parameter of ε and is ommonly alled the nugget, and θ T = (σ0 2, α) are the sill and range parameters of the ovariane funtion (or variogram). For some variogram models θ ontains further parameters. Using the observations, y T = (y(s 1 ), y(s 2 ),..., y(s n )) we want to estimate β, σ 2 and θ and to predit Z(s) for some s where s an be equal to one of the sampled loations s i or different. The starting point for developing our new proedure is the Gaussian log-likelihood (1) l(β, σ 2, θ y) = 1 2 log ( det(σ 2 I + V θ )) 1 2 (y Xβ) T (σ 2 I + V θ ) 1 (y Xβ) that arises when ε(s) is Gaussian. V θ = σ 2 0 V α is the ovariane and V α the orrelation matrix of z T = ( z(s 1 ), z(s 2 ),..., z(s n ) ). We define the pseudo log-likelihood that depends in addition to the parameters on the latent variable z (2) l (β, σ 2, θ, z y) = 1 2 log ( det(σ 2 I + V θ ) ) 1 2 n ( y(si ) x(s i ) T ) 2 β z(s i ) 1 σ 2 zt V 1 θ z. i=1 The Gaussian log-likelihood may be onsidered as a profile log-likelihood of l (β, σ 2, θ, z y) that has been maximized with respet to z. Hene, l(β, σ 2, θ y) = l (β, σ 2, θ, z y) where ( z = argmax l (β, σ 2, θ, z y) ) = V θ (σ 2 I + V θ ) 1 (y Xβ). z Maximizing (1) with respet to β, σ 2 and θ is thus equivalent to solving the expanded system of ML estimating equations (3) (4) (5) (6) (7) = 1ˆσ y X ˆβ ẑ V 1 ẑ = 0, z ˆσ ˆθ β = y X ˆβ ẑ XT = 0, ˆσ σ 2 = trae ( (I + 1ˆσ 2 Vˆθ) 1) σ 2 0 ( n i=1 = n trae ( (I + 1ˆσ 2 Vˆθ) 1) ẑ T V 1 ẑ = 0, ˆθ α = trae( ( ˆσ2 ˆσ 2 0 I + Vˆα ) 1 V α α α=ˆα ) 1 ˆσ 2 0 ) 2 y(s i ) x(s i ) T ˆβ ẑ(si ) = 0, ˆσ ẑ T V 1 V α ˆα V 1 α ˆα ẑ = 0. α=ˆα Closed-form expressions exist in general only for the generalized least squares (GLS) estimates of β ˆβ = (X T (ˆσ 2 I + Vˆθ) 1 X) 1 X T (ˆσ 2 I + Vˆθ) 1 y, and the plug-in mean square preditions of z ẑ = Vˆθ(ˆσ 2 I + Vˆθ) 1 (y X ˆβ), whih are alled universal (UK) or external drift kriging preditions in geostatistis.
4 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August Equations (3) (7) show that the Gaussian ML estimates and the UK preditions depend on the standardized residuals ˆr/ˆσ = (y X ˆβ ẑ)/ˆσ, and it is obvious that both estimates and preditions are sensitive to outliers in y. Welsh and Rihardson (1997) proposed several approahes to make ML estimates in linear mixed models robust to outlying observations. A first approah onsists of robustifying the Gaussian log-likelihood. One has then to replae (y(s i ) x(s i ) T β z(s i )) 2 /σ 2 in (2) by a funtion ρ ( ) of the standardized errors that grows more slowly than the quadrati ( denotes here a tuning parameter that ontrols the amount of robustness of the estimates). Initially, we tried suh a strategy and approximated the restrited likelihood funtion (Harville, 1974) of a model with long-tailed ε(s) by a Laplae approximation and maximized the approximation to get robust REML estimates of the drift and variogram parameters (Shwierz et al., 2010). However, to bound the influene of outliers, we had to use a bounded ρ ( ) funtion, whih made the numerial maximization of the restrited log-likelihood more diffiult beause the related system of estimating equations has then multiple roots. Furthermore, we had diffiulties to find appropriate bias orretions to obtain Fisher onsistent estimates for Gaussian ε(s). Therefore, we eventually pursued another suggestion by Welsh and Rihardson (1997) and robustified (3) (7) by replaing ˆr/ˆσ by a bounded and odd funtion ψ ( ) of the standardized residuals, yielding thereby a system of robustified estimating equations (8) (9) (10) (11) (12) trae ( (I + 1 σ 2 V θ) 1) n i=1 trae ( ( σ2 σ 0 2 I + V α ) 1 V α ) 1 α σ 0 2 α= α 1 σ ψ (y X β z ) V 1 z + a σ θ 1 = 0, X T (y X ψ β z ) + a2 σ = 0, ψ 2 (y(s i ) x(s i ) T β z(si )) + a3 σ = 0, n trae ( (I + 1 σ 2 V θ) 1) z T V θ 1 z + a 4 = 0, z T V α 1 V α V α 1 α z + a 5 = 0. α= α In the above equations, a 1 to a 5 are adjustments for Fisher onsisteny at the Gaussian model. They are determined by the ondition that the expetations of (8) (12) must vanish for Gaussian ε(s) (e.g. Maronna et al., 2006, p. 67). Hene, a 1 = a 2 = 0 sine ψ ( ) is odd, and (13) a 3 = a 3 trae ( (I + 1 σ 2 V θ) 1) (14) (15) a 4 = a 4 ( n trae ( (I + 1 σ 2 V θ) 1)), a 5 = a 5 trae ( ( σ2 σ 0 2 I + V α ) 1 V α ), α α= α where (using (8) to derive (16)) n a 3 = E [ψ 2 (y(s i ) x(s i ) T β z(si )) (16) = σ 2 trae ( V 2 σ θ E [ z, zt ), (17) (18) i=1 a 4 = E [ z T V θ 1 z = trae ( V 1 a 5 = E [ 1 σ 2 0 z T V 1 α θ E [ z, zt ), V α α α= α V 1 α z = 1 σ 2 0 trae ( V α 1 V α V α 1 α E [ z, zt ). α= α To orret the biases we thus need the ovariane matrix of z, whih we approximated in the following way: The starting point was a first-order Taylor series approximation for ψ ( r/σ) at ε/σ ψ ( r/σ) ψ (ε/σ) 1 σ diag( ψ (ε/σ) )( z z + X( β β) )
5 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August where r T = y X β z and ε T = (ε(s 1 ), ε(s 2 ),..., ε(s n )). Then we replaed ψ (ε(s i )/σ) by its expetation, say b = E [ψ (ε(s i )/σ) (for Gaussian ε(s)), and substituted the modified Taylor approximation in (8) (9) for ψ ( r/σ, whih resulted in [ z β β M 1 [ I X T I X T [ b z σψ(ε/σ) where M = [ b I + σ 2 V 1 θ b X T b X b X T X. Sine z and ψ (ε/σ) are independent, and letting a = E [ψ 2 (ε(s)/σ) (for Gaussian ε(s)), we finally obtained the approximation [ Cov [ z, z T Cov [ z, (19) β T [ Cov [ β, z T Cov [ β, β T M 1 Λ ΛX X T Λ X T M 1, ΛX where Λ = b 2 V θ + a σ 2 I. Solving (8) (12) subjet to (13) (15) and (19) for the non-robust ase yields the ustomary REML estimates of σ 2 and θ and the related plug-in GLS estimate of β and the UK predition of z. Hene, apart from the Bayesian interpretation of the REML estimates (Harville, 1974) and the interpretation as ML estimates obtained from error ontrasts (Patterson and Thompson, 1971), the REML estimates an also be onsidered as ML estimates, fored to be Fisher onsistent for finite sample size. Estimation proedure We wrote an R funtion that solves the robustified estimating equations for z, β, σ 2 and θ. Instead of the popular Huber ψ -funtion (Maronna et al., 2006, p. 26), we used the ontinuously differentiable funtion (20) ψ (x) = exp( 2x/), whih orresponds to a shifted and saled logisti umulative distribution funtion. Note that lim ψ (x) = ±. We omputed the roots of the estimating equations by a ombination of iteratively re-weighted least squares (IRWLS, e.g. Maronna et al., 2006, p. 105) and a Broyden s sheme. x ± For given σ 2 and θ, z and β were determined by IRWLS. The estimates were then inserted into (10) (12), and the roots of these equations were omputed by Broyden s method as implemented in the R pakage nleqslv. The estimation proedure needs initial values for z and all the parameters. An initial guess of β was obtained by MM-estimation (Maronna et al., 2006, pp. 124). The robust regression residuals were then spatially smoothed by loess to get an initial z. The MAD of the loess residuals was taken as initial σ. An initial guess of θ was obtained by fitting the sill and range parameters of a parametri variogram model (with nugget fixed at the initial σ 2 ) by weighted non-linear least squares (Cressie, 1993, p. 96) to the sample variogram of the regression residuals that had been omputed by the MAD estimator (Dowd, 1984).
6 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August Table 1: Statistis of relative errors (%) of the approximations of a 3, a 4 and a 5 in simulations of n = Y (s) (s uniformly distributed on the unit square) for tuning parameter = 1. The data were generated with a linear external drift in the oordinates and an exponential variogram with ranges varying from 0.01 to 0.2 and nugget:sill ratios from 0.25 to 4. minimum 1st quartile median 3rd quartile maximum a a a Simulation study We explored the properties of our robust REML estimation method by simulations. First, we heked approximation (19) for Cov [ z, z T. We simulated for several external drifts and exponential or ubi variograms (with various ombinations of nugget, sill and range) data sets onsisting of Gaussian Y (s) on grids or for uniformly distributed s in one- and two-dimensional domains. For eah senario, we simulated realizations and omputed for eah realization a 3, a 4 and a 5 using (19) and the true σ2 and θ. Then we ompared the mean of these a i values with the true value that was omputed with the empirial ovarianes of the estimates z, obtained by interept slope easting slope northing sill nugget range Figure 1: Boxplots of estimates β, σ 2 and θ omputed from 200 Gaussian data sets that were simulated with a linear external drift in the oordinates and an exponential variogram with σ 2 = 0.5, σ0 2 = 2, and α = 0.05 for n = 200 uniformly distributed s in the unit square (dotted lines mark true parameter values, for = 1000 estimates orrespond to ustomary REML estimates).
7 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August interept Z ~ (s) at ontam. lo s=0.895 Z ~ (s) next to ontam. lo. (s=0.905) Z ~ (s) far from ontam. lo. (s=0.099) = 1000 = 2 = 1.67 = 1.33 = = 1000 = 2 = 1.67 = 1.33 = = 1000 = 2 = 1.67 = 1.33 = = 1000 = 2 = 1.67 = 1.33 = 1 sill nugget range = 1000 = 2 = 1.67 = 1.33 = = 1000 = 2 = 1.67 = 1.33 = = 1000 = 2 = 1.67 = 1.33 = 1 Figure 2: Standardized urves for one realization of the simulation senario 500 Gaussian Y (s) on entred regular grid in [0, 1; no trend (β = 0.5); exponential variogram with σ 2 = 0.5, σ0 2 = 2, α = The observation at s = was ontaminated by adding x ( 20, 19.75, 19.5,..., 19.75, 20). Then we estimated for eah x the model parameters and z, and plotted n( φ(x) φ(0)) vs x, where φ denotes a parameter or an element of z. Estimates and preditions orrespond for = 1000 to ustomary REML estimates and UK preditions. solving (8) (9) with the true ovariane parameters. Table 1 lists as example the relative errors of the approximations for the simulations with the spatially uniformly distributed s on the unit square. On average, the approximation underestimated the a i, but, exept for the very short ranges, the modulus of the relative errors rarely exeeded 10 %. Errors in the approximations of a 3 to a 5 lead to some bias in the parameter estimates. Figure 1 shows boxplots of the estimates β, σ 2 and θ for a simulation senario where the approximation resulted for = 1 in relatively large errors (a 3 : %, a 4 : -8.3 %, a %). The medians of the variane estimates deviate with dereasing inreasingly from the true values, but ompared to the biases that we faed when maximizing the robustified log-likelihood these systemati errors seem tolerable. Seond, we omputed for a single realization of one simulation senario urves (Maronna et al., 2006, p. 55). Figure 2 shows that the urves remained bounded for 2, both for the parameter estimates and the preditions. As expeted from theory, the dereased with dereasing. In ontrast, the non-robust estimates ( = 1000) got inreasingly orrupted by inreasing the severity of the outlier. Third, we simulated 200 realizations of a ontaminated Gaussian field where 5 % of the ε(s) had been replaed by normal variates with ten times as large standard deviation. The boxplots of the estimates are shown in Figure 3. The ontamination strongly inflated the non-robust ( = 1000)
8 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August interept sill nugget range Figure 3: Boxplots of estimates β, σ 2 and θ omputed from 200 ontaminated Gaussian data sets (500 Y (s) on entred regular grid in [0, 1; no trend (β = 0.5); exponential variogram with σ 2 = 0.5, σ0 2 = 2, α = 0.05; 5 % of ε(s) randomly replaed by normal variates with mean zero and standard deviation 7.9) (dotted lines mark true parameter values, for = 1000 estimates orrespond to ustomary REML estimates). REML estimate of the nugget, biased (negatively) the REML estimate of θ and inreased the varianes of the estimates of all three ovariane parameters. The robust estimates of θ showed no systemati errors and σ 2 was muh less biased and more effiient than the REML estimate. The systemati error in σ 2 dereased with dereasing. Summary and onlusions We developed a novel robust REML method to estimate the external drift and the variogram parameters of Gaussian spatial data that are possibly ontaminated by outliers. The methods also provides robustified universal kriging preditions. Simulations showed that our proedure bounds the influene of outliers on parameter estimates and preditions, and yields approximately Fisher onsistent estimates at the Gaussian model. For large, the estimates and preditions oinide with the ustomary REML estimates and the plug-in UK preditions. So far, we did not yet hek whether (19) provides an aurate approximation of Cov [ β, β T and Cov [( z z), ( z z) T. These quantities are required for omputing the mean square predition error at non-sampled s. Furthermore, we have to develop robust proedures for testing hypotheses about β. Some of these further developments will be presented at the onferene, along with an analysis of data on heavy metals in the soils around a metal smelter in Switzerland where removal and displaement of ontaminated soil has led to unpreditable irregularities in the spatial distribution of the metals. Aknowledgements We gratefully aknowledge the support of the Swiss Federal Offie for the Environment (FOEN), whih funded a part of this researh. REFERENCES Cressie, N. A. C. (1993). Statistis for Spatial Data. John Wiley & Sons, New York, revised edition. Dowd, P. A. (1984). The variogram and kriging: Robust and resistant estimators. In G. Verly, M. David,
9 ISI 58th World Statistis Congress of the International Statistial Institute, Dublin, August A. Journel, and A. Maréhal, editors, Geostatistis for Natural Resoures Charaterization, volume 1, pages , Dordreht. D. Reidel Publishing Company. Harville, D. A. (1974). Bayesian inferene for variane omponents using only error ontrasts. Biometrika, 61, Hawkins, D. M. and Cressie, N. (1984). Robust kriging a proposal. Mathematial Geology, 16, Lark, R. M. (2000). A omparison of some robust estimators of the variogram for use in soil survey. European Journal of Soil Siene, 51, Marhant, B. P. and Lark, R. M. (2007). Robust estimation of the variogram by residual maximum likelihood. Geoderma, 140, Maronna, R. A., Martin, R. D., and Yohai, V. J. (2006). Robust Statistis Theory and Methods. John Wiley & Sons, Chihester. Militino, A. F. and Ugarte, M. D. (1997). A GM estimation of the loation parameters in a spatial linear model. Communiations in Statistis, Theory and Methods, 26(7), Patterson, H. D. and Thompson, R. (1971). Reovery of inter-blok information when blok sizes are unequal. Biometrika, 58(3), R Development Core Team (2011). R: A Language and Environment for Statistial Computing. R Foundation for Statistial Computing, Vienna, Austria. Rihardson, A. M. and Welsh, A. H. (1995). Robust restrited maximum likelihood in mixed linear models. Biometris, 51, Shwierz, C., Papritz, A., Stahel, W. A., and Künsh, H. R. (2010). Robust REML estimation of the variogram and robust kriging. In Proeedings of the International Conferene on Robust Statistis Prague 2010, Book of Abstrats, pages 91 92, Prague. Welsh, A. H. and Rihardson, A. M. (1997). Approahes to the robust estimation of mixed models. In Robust Inferene, volume 15 of Handbook of Statistis, pages Elsevier.
Likelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationMethods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationSome recent developments in probability distributions
Proeedings 59th ISI World Statistis Congress, 25-30 August 2013, Hong Kong (Session STS084) p.2873 Some reent developments in probability distributions Felix Famoye *1, Carl Lee 1, and Ayman Alzaatreh
More informationOptimization of Statistical Decisions for Age Replacement Problems via a New Pivotal Quantity Averaging Approach
Amerian Journal of heoretial and Applied tatistis 6; 5(-): -8 Published online January 7, 6 (http://www.sienepublishinggroup.om/j/ajtas) doi:.648/j.ajtas.s.65.4 IN: 36-8999 (Print); IN: 36-96 (Online)
More informationLOGISTIC REGRESSION IN DEPRESSION CLASSIFICATION
LOGISIC REGRESSIO I DEPRESSIO CLASSIFICAIO J. Kual,. V. ran, M. Bareš KSE, FJFI, CVU v Praze PCP, CS, 3LF UK v Praze Abstrat Well nown logisti regression and the other binary response models an be used
More informationChapter 8 Hypothesis Testing
Leture 5 for BST 63: Statistial Theory II Kui Zhang, Spring Chapter 8 Hypothesis Testing Setion 8 Introdution Definition 8 A hypothesis is a statement about a population parameter Definition 8 The two
More informationFiber Optic Cable Transmission Losses with Perturbation Effects
Fiber Opti Cable Transmission Losses with Perturbation Effets Kampanat Namngam 1*, Preeha Yupapin 2 and Pakkinee Chitsakul 1 1 Department of Mathematis and Computer Siene, Faulty of Siene, King Mongkut
More informationExploring the feasibility of on-site earthquake early warning using close-in records of the 2007 Noto Hanto earthquake
Exploring the feasibility of on-site earthquake early warning using lose-in reords of the 2007 Noto Hanto earthquake Yih-Min Wu 1 and Hiroo Kanamori 2 1. Department of Geosienes, National Taiwan University,
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationCONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION. Smail Mahdi
Serdia Math. J. 30 (2004), 55 70 CONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION Smail Mahdi Communiated by N. M. Yanev Abstrat. A two-sided onditional onfidene interval
More informationPerforming Two-Way Analysis of Variance Under Variance Heterogeneity
Journal of Modern Applied Statistial Methods Volume Issue Artile 3 5--003 Performing Two-Way Analysis of Variane Under Variane Heterogeneity Sott J. Rihter University of North Carolina at Greensboro, sjriht@ung.edu
More informationControl Theory association of mathematics and engineering
Control Theory assoiation of mathematis and engineering Wojieh Mitkowski Krzysztof Oprzedkiewiz Department of Automatis AGH Univ. of Siene & Tehnology, Craow, Poland, Abstrat In this paper a methodology
More informationDevelopment of Fuzzy Extreme Value Theory. Populations
Applied Mathematial Sienes, Vol. 6, 0, no. 7, 58 5834 Development of Fuzzy Extreme Value Theory Control Charts Using α -uts for Sewed Populations Rungsarit Intaramo Department of Mathematis, Faulty of
More informationREFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS. 1. Introduction
Version of 5/2/2003 To appear in Advanes in Applied Mathematis REFINED UPPER BOUNDS FOR THE LINEAR DIOPHANTINE PROBLEM OF FROBENIUS MATTHIAS BECK AND SHELEMYAHU ZACKS Abstrat We study the Frobenius problem:
More informationSensitivity Analysis in Markov Networks
Sensitivity Analysis in Markov Networks Hei Chan and Adnan Darwihe Computer Siene Department University of California, Los Angeles Los Angeles, CA 90095 {hei,darwihe}@s.ula.edu Abstrat This paper explores
More informationAn iterative least-square method suitable for solving large sparse matrices
An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method
More informationRESEARCH ON RANDOM FOURIER WAVE-NUMBER SPECTRUM OF FLUCTUATING WIND SPEED
The Seventh Asia-Paifi Conferene on Wind Engineering, November 8-1, 9, Taipei, Taiwan RESEARCH ON RANDOM FORIER WAVE-NMBER SPECTRM OF FLCTATING WIND SPEED Qi Yan 1, Jie Li 1 Ph D. andidate, Department
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More information23.1 Tuning controllers, in the large view Quoting from Section 16.7:
Lesson 23. Tuning a real ontroller - modeling, proess identifiation, fine tuning 23.0 Context We have learned to view proesses as dynami systems, taking are to identify their input, intermediate, and output
More informationAssessing the Performance of a BCI: A Task-Oriented Approach
Assessing the Performane of a BCI: A Task-Oriented Approah B. Dal Seno, L. Mainardi 2, M. Matteui Department of Eletronis and Information, IIT-Unit, Politenio di Milano, Italy 2 Department of Bioengineering,
More informationCopyright 2018 Society of Photo-Optical Instrumentation Engineers (SPIE). One print or electronic copy may be made for personal use only.
Copyright 018 Soiety of Photo-Optial Instrumentation Engineers (SPIE) One print or eletroni opy may be made for personal use only Systemati reprodution and distribution, dupliation of any material in this
More informationAnalysis of discretization in the direct simulation Monte Carlo
PHYSICS OF FLUIDS VOLUME 1, UMBER 1 OCTOBER Analysis of disretization in the diret simulation Monte Carlo iolas G. Hadjionstantinou a) Department of Mehanial Engineering, Massahusetts Institute of Tehnology,
More informationarxiv:gr-qc/ v2 6 Feb 2004
Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this
More informationProbabilistic Graphical Models
Probabilisti Graphial Models David Sontag New York University Leture 12, April 19, 2012 Aknowledgement: Partially based on slides by Eri Xing at CMU and Andrew MCallum at UMass Amherst David Sontag (NYU)
More informationSupplementary Materials
Supplementary Materials Neural population partitioning and a onurrent brain-mahine interfae for sequential motor funtion Maryam M. Shanehi, Rollin C. Hu, Marissa Powers, Gregory W. Wornell, Emery N. Brown
More informationA Spatiotemporal Approach to Passive Sound Source Localization
A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,
More informationQCLAS Sensor for Purity Monitoring in Medical Gas Supply Lines
DOI.56/sensoren6/P3. QLAS Sensor for Purity Monitoring in Medial Gas Supply Lines Henrik Zimmermann, Mathias Wiese, Alessandro Ragnoni neoplas ontrol GmbH, Walther-Rathenau-Str. 49a, 7489 Greifswald, Germany
More informationComputer Science 786S - Statistical Methods in Natural Language Processing and Data Analysis Page 1
Computer Siene 786S - Statistial Methods in Natural Language Proessing and Data Analysis Page 1 Hypothesis Testing A statistial hypothesis is a statement about the nature of the distribution of a random
More informationModel-based mixture discriminant analysis an experimental study
Model-based mixture disriminant analysis an experimental study Zohar Halbe and Mayer Aladjem Department of Eletrial and Computer Engineering, Ben-Gurion University of the Negev P.O.Box 653, Beer-Sheva,
More informationDanielle Maddix AA238 Final Project December 9, 2016
Struture and Parameter Learning in Bayesian Networks with Appliations to Prediting Breast Caner Tumor Malignany in a Lower Dimension Feature Spae Danielle Maddix AA238 Final Projet Deember 9, 2016 Abstrat
More informationarxiv:quant-ph/ v1 18 Dec 2000
Effets of Imperfet Gate Operations in Shor s rime Fatorization Algorithm Hao Guo 1,, Gui-Lu Long 1,,3,4,5 and Yang Sun 1,,6,7 1 Department of hysis, Tsinghua University, Beijing 100084 Key Laboratory for
More informationEstimation for Unknown Parameters of the Extended Burr Type-XII Distribution Based on Type-I Hybrid Progressive Censoring Scheme
Estimation for Unknown Parameters of the Extended Burr Type-XII Distribution Based on Type-I Hybrid Progressive Censoring Sheme M.M. Amein Department of Mathematis Statistis, Faulty of Siene, Taif University,
More informationRisk Analysis in Water Quality Problems. Souza, Raimundo 1 Chagas, Patrícia 2 1,2 Departamento de Engenharia Hidráulica e Ambiental
Risk Analysis in Water Quality Problems. Downloaded from aselibrary.org by Uf - Universidade Federal Do Ceara on 1/29/14. Coyright ASCE. For ersonal use only; all rights reserved. Souza, Raimundo 1 Chagas,
More informationSensitivity analysis for linear optimization problem with fuzzy data in the objective function
Sensitivity analysis for linear optimization problem with fuzzy data in the objetive funtion Stephan Dempe, Tatiana Starostina May 5, 2004 Abstrat Linear programming problems with fuzzy oeffiients in the
More informationRobust Flight Control Design for a Turn Coordination System with Parameter Uncertainties
Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo
More informationApplying CIECAM02 for Mobile Display Viewing Conditions
Applying CIECAM2 for Mobile Display Viewing Conditions YungKyung Park*, ChangJun Li*, M.. Luo*, Youngshin Kwak**, Du-Sik Park **, and Changyeong Kim**; * University of Leeds, Colour Imaging Lab, UK*, **
More informationPhase Diffuser at the Transmitter for Lasercom Link: Effect of Partially Coherent Beam on the Bit-Error Rate.
Phase Diffuser at the Transmitter for Laserom Link: Effet of Partially Coherent Beam on the Bit-Error Rate. O. Korotkova* a, L. C. Andrews** a, R. L. Phillips*** b a Dept. of Mathematis, Univ. of Central
More informationThe Distributors of Change Points in Long 11ennoryProcesses
Summer Researh Projet 25-26 The Distributors of Change Points in Long 11ennoryProesses Guan Yu Zheng Department of Mathematis and Statistis University of Canterbury The distributions of hange points in
More informationarxiv: v1 [physics.soc-ph] 26 Nov 2014
A study of the influene of the mobility on the phase transitions of the synhronous SIR model Roberto da Silva 1, Henrique A. Fernandes 2 1 Instituto de Físia, Universidade Federal do Rio Grande do Sul,
More informationModels for the simulation of electronic circuits with hysteretic inductors
Proeedings of the 5th WSEAS Int. Conf. on Miroeletronis, Nanoeletronis, Optoeletronis, Prague, Czeh Republi, Marh 12-14, 26 (pp86-91) Models for the simulation of eletroni iruits with hystereti indutors
More informationThe gravitational phenomena without the curved spacetime
The gravitational phenomena without the urved spaetime Mirosław J. Kubiak Abstrat: In this paper was presented a desription of the gravitational phenomena in the new medium, different than the urved spaetime,
More informationCOMPARISON OF ESTIMATORS FOR STRESS-STRENGTH RELIABILITY IN THE GOMPERTZ CASE
Haettepe Journal of Mathematis and Statistis olume 383 29, 339 349 COMPARISON OF ESTIMATORS FOR STRESS-STRENGTH RELIABILITY IN THE GOMPERTZ CASE Buğra Saraçoğlu, Mehmet Fedai Kaya and A. M. Abd-Elfattah
More informationSingular Event Detection
Singular Event Detetion Rafael S. Garía Eletrial Engineering University of Puerto Rio at Mayagüez Rafael.Garia@ee.uprm.edu Faulty Mentor: S. Shankar Sastry Researh Supervisor: Jonathan Sprinkle Graduate
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationOn the Licensing of Innovations under Strategic Delegation
On the Liensing of Innovations under Strategi Delegation Judy Hsu Institute of Finanial Management Nanhua University Taiwan and X. Henry Wang Department of Eonomis University of Missouri USA Abstrat This
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationA Unified View on Multi-class Support Vector Classification Supplement
Journal of Mahine Learning Researh??) Submitted 7/15; Published?/?? A Unified View on Multi-lass Support Vetor Classifiation Supplement Ürün Doğan Mirosoft Researh Tobias Glasmahers Institut für Neuroinformatik
More informationMath 151 Introduction to Eigenvectors
Math 151 Introdution to Eigenvetors The motivating example we used to desrie matrixes was landsape hange and vegetation suession. We hose the simple example of Bare Soil (B), eing replaed y Grasses (G)
More informationA new projection estimate for multivariate location with minimax bias.
A new projetion estimate for multivariate loation with minimax bias. Jorge G. Adrover Universidad Naional de Cordoba and CONICET and Vtor J. Yohai y Universidad de Buenos Aires and CONICET TECHNICAL REPORT
More informationCavity flow with surface tension past a flat plate
Proeedings of the 7 th International Symposium on Cavitation CAV9 Paper No. ## August 7-, 9, Ann Arbor, Mihigan, USA Cavity flow with surfae tension past a flat plate Yuriy Savhenko Institute of Hydromehanis
More informationUPPER-TRUNCATED POWER LAW DISTRIBUTIONS
Fratals, Vol. 9, No. (00) 09 World Sientifi Publishing Company UPPER-TRUNCATED POWER LAW DISTRIBUTIONS STEPHEN M. BURROUGHS and SARAH F. TEBBENS College of Marine Siene, University of South Florida, St.
More informationIMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION
09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE
More informationTests of fit for symmetric variance gamma distributions
Tests of fit for symmetri variane gamma distributions Fragiadakis Kostas UADPhilEon, National and Kapodistrian University of Athens, 4 Euripidou Street, 05 59 Athens, Greee. Keywords: Variane Gamma Distribution,
More informationOn dealing with spatially correlated residuals in remote sensing and GIS
On dealing with spatially correlated residuals in remote sensing and GIS Nicholas A. S. Hamm 1, Peter M. Atkinson and Edward J. Milton 3 School of Geography University of Southampton Southampton SO17 3AT
More informationAn I-Vector Backend for Speaker Verification
An I-Vetor Bakend for Speaker Verifiation Patrik Kenny, 1 Themos Stafylakis, 1 Jahangir Alam, 1 and Marel Kokmann 2 1 CRIM, Canada, {patrik.kenny, themos.stafylakis, jahangir.alam}@rim.a 2 VoieTrust, Canada,
More informationONLINE APPENDICES for Cost-Effective Quality Assurance in Crowd Labeling
ONLINE APPENDICES for Cost-Effetive Quality Assurane in Crowd Labeling Jing Wang Shool of Business and Management Hong Kong University of Siene and Tehnology Clear Water Bay Kowloon Hong Kong jwang@usthk
More information10.5 Unsupervised Bayesian Learning
The Bayes Classifier Maximum-likelihood methods: Li Yu Hongda Mao Joan Wang parameter vetor is a fixed but unknown value Bayes methods: parameter vetor is a random variable with known prior distribution
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationDIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS
CHAPTER 4 DIGITAL DISTANCE RELAYING SCHEME FOR PARALLEL TRANSMISSION LINES DURING INTER-CIRCUIT FAULTS 4.1 INTRODUCTION Around the world, environmental and ost onsiousness are foring utilities to install
More informationImproved likelihood inference for discrete data
Improved likelihood inferene for disrete data A. C. Davison Institute of Mathematis, Eole Polytehnique Fédérale de Lausanne, Station 8, 1015 Lausanne, Switzerland. D. A. S. Fraser and N. Reid Department
More informationGeneral probability weighted moments for the three-parameter Weibull Distribution and their application in S-N curves modelling
General probability weighted moments for the three-parameter Weibull Distribution and their appliation in S-N urves modelling Paul Dario Toasa Caiza a,, Thomas Ummenhofer a a KIT Stahl- und Leihtbau, Versuhsanstalt
More information18.05 Problem Set 6, Spring 2014 Solutions
8.5 Problem Set 6, Spring 4 Solutions Problem. pts.) a) Throughout this problem we will let x be the data of 4 heads out of 5 tosses. We have 4/5 =.56. Computing the likelihoods: 5 5 px H )=.5) 5 px H
More informationChapter 4 - Fundamentals of spatial processes Lecture notes
Chapter 4 - Fundamentals of spatial processes Lecture notes Geir Storvik January 21, 2013 STK4150 - Intro 2 Spatial processes Typically correlation between nearby sites Mostly positive correlation Negative
More information36-720: Log-Linear Models: Three-Way Tables
36-720: Log-Linear Models: Three-Way Tables Brian Junker September 10, 2007 Hierarhial Speifiation of Log-Linear Models Higher-Dimensional Tables Digression: Why Poisson If We Believe Multinomial? Produt
More informationMODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP OF CONCRETE IN UNIAXIAL COMPRESSION
VIII International Conferene on Frature Mehanis of Conrete and Conrete Strutures FraMCoS-8 J.G.M. Van Mier, G. Ruiz, C. Andrade, R.C. Yu and X.X. Zhang Eds) MODELLING THE POSTPEAK STRESS DISPLACEMENT RELATIONSHIP
More informationTwo Points Hybrid Block Method for Solving First Order Fuzzy Differential Equations
Journal of Soft Computing and Appliations 2016 No.1 (2016) 43-53 Available online at www.ispas.om/jsa Volume 2016, Issue 1, Year 2016 Artile ID jsa-00083, 11 Pages doi:10.5899/2016/jsa-00083 Researh Artile
More informationthe following action R of T on T n+1 : for each θ T, R θ : T n+1 T n+1 is defined by stated, we assume that all the curves in this paper are defined
How should a snake turn on ie: A ase study of the asymptoti isoholonomi problem Jianghai Hu, Slobodan N. Simić, and Shankar Sastry Department of Eletrial Engineering and Computer Sienes University of California
More informationLightpath routing for maximum reliability in optical mesh networks
Vol. 7, No. 5 / May 2008 / JOURNAL OF OPTICAL NETWORKING 449 Lightpath routing for maximum reliability in optial mesh networks Shengli Yuan, 1, * Saket Varma, 2 and Jason P. Jue 2 1 Department of Computer
More information( ) ( ) Volumetric Properties of Pure Fluids, part 4. The generic cubic equation of state:
CE304, Spring 2004 Leture 6 Volumetri roperties of ure Fluids, part 4 The generi ubi equation of state: There are many possible equations of state (and many have been proposed) that have the same general
More informationThe Hanging Chain. John McCuan. January 19, 2006
The Hanging Chain John MCuan January 19, 2006 1 Introdution We onsider a hain of length L attahed to two points (a, u a and (b, u b in the plane. It is assumed that the hain hangs in the plane under a
More informationResearch Collection. Mismatched decoding for the relay channel. Conference Paper. ETH Library. Author(s): Hucher, Charlotte; Sadeghi, Parastoo
Researh Colletion Conferene Paper Mismathed deoding for the relay hannel Author(s): Huher, Charlotte; Sadeghi, Parastoo Publiation Date: 2010 Permanent Link: https://doi.org/10.3929/ethz-a-005997152 Rights
More informationUniversity of Bristol - Explore Bristol Research
Chen, J., Ravey, E., Hallett, S. R., Wisnom, M. R., & Grassi, M. (2009). Predition of delamination in braided omposite T-piee speimens. Composites Siene and Tehnology, 69(14), 2363-2367. 10.1016/j.ompsiteh.2009.01.027
More informationA Time-Dependent Model For Predicting The Response Of A Horizontally Loaded Pile Embedded In A Layered Transversely Isotropic Saturated Soil
IOSR Journal of Mehanial and Civil Engineering (IOSR-JMCE) e-issn: 2278-1684,p-ISSN: 232-334X, Volume 16, Issue 2 Ser. I (Mar. - Apr. 219), PP 48-53 www.iosrjournals.org A Time-Dependent Model For Prediting
More informationMean Activity Coefficients of Peroxodisulfates in Saturated Solutions of the Conversion System 2NH 4. H 2 O at 20 C and 30 C
Mean Ativity Coeffiients of Peroxodisulfates in Saturated Solutions of the Conversion System NH 4 Na S O 8 H O at 0 C and 0 C Jan Balej Heřmanova 5, 170 00 Prague 7, Czeh Republi balejan@seznam.z Abstrat:
More informationCOMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION
4 TH INTERNATIONAL CONGRESS OF THE AERONAUTICAL SCIENCES COMBINED PROBE FOR MACH NUMBER, TEMPERATURE AND INCIDENCE INDICATION Jiri Nozika*, Josef Adame*, Daniel Hanus** *Department of Fluid Dynamis and
More informationDiscrete Bessel functions and partial difference equations
Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat
More informationA Heuristic Approach for Design and Calculation of Pressure Distribution over Naca 4 Digit Airfoil
IOSR Journal of Engineering (IOSRJEN) ISSN (e): 2250-3021, ISSN (p): 2278-8719 PP 11-15 www.iosrjen.org A Heuristi Approah for Design and Calulation of Pressure Distribution over Naa 4 Digit Airfoil G.
More informationSINCE Zadeh s compositional rule of fuzzy inference
IEEE TRANSACTIONS ON FUZZY SYSTEMS, VOL. 14, NO. 6, DECEMBER 2006 709 Error Estimation of Perturbations Under CRI Guosheng Cheng Yuxi Fu Abstrat The analysis of stability robustness of fuzzy reasoning
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
More informationCalculation of Desorption Parameters for Mg/Si(111) System
e-journal of Surfae Siene and Nanotehnology 29 August 2009 e-j. Surf. Si. Nanoteh. Vol. 7 (2009) 816-820 Conferene - JSSS-8 - Calulation of Desorption Parameters for Mg/Si(111) System S. A. Dotsenko, N.
More informationSPLINE ESTIMATION OF SINGLE-INDEX MODELS
SPLINE ESIMAION OF SINGLE-INDEX MODELS Li Wang and Lijian Yang University of Georgia and Mihigan State University Supplementary Material his note ontains proofs for the main results he following two propositions
More informationWhen p = 1, the solution is indeterminate, but we get the correct answer in the limit.
The Mathematia Journal Gambler s Ruin and First Passage Time Jan Vrbik We investigate the lassial problem of a gambler repeatedly betting $1 on the flip of a potentially biased oin until he either loses
More informationCollinear Equilibrium Points in the Relativistic R3BP when the Bigger Primary is a Triaxial Rigid Body Nakone Bello 1,a and Aminu Abubakar Hussain 2,b
International Frontier Siene Letters Submitted: 6-- ISSN: 9-8, Vol., pp -6 Aepted: -- doi:.8/www.sipress.om/ifsl.. Online: --8 SiPress Ltd., Switzerland Collinear Equilibrium Points in the Relativisti
More informationTransformation to approximate independence for locally stationary Gaussian processes
ransformation to approximate independene for loally stationary Gaussian proesses Joseph Guinness, Mihael L. Stein We provide new approximations for the likelihood of a time series under the loally stationary
More informationSOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS
SOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS Prepared by S. Broverman e-mail 2brove@rogers.om website http://members.rogers.om/2brove 1. We identify the following events:. - wathed gymnastis, ) - wathed baseball,
More informationDynamic Simulation and Composition Control in A 10 L Mixing Tank
roeeding of International Conferene on Chemial and Material Engineering ISBN : 978--97-8-7 SE.4 - Dynami Simulation and Composition Control in A L Mixing Tank Yulius Deddy Hermawan Chemial Engineering
More informationThe shape of a hanging chain. a project in the calculus of variations
The shape of a hanging hain a projet in the alulus of variations April 15, 218 2 Contents 1 Introdution 5 2 Analysis 7 2.1 Model............................... 7 2.2 Extremal graphs.........................
More informationDiscrete Generalized Burr-Type XII Distribution
Journal of Modern Applied Statistial Methods Volume 13 Issue 2 Artile 13 11-2014 Disrete Generalized Burr-Type XII Distribution B. A. Para University of Kashmir, Srinagar, India, parabilal@gmail.om T.
More informationAvailable online Journal of Scientific and Engineering Research, 2017, 4(7): Research Article
Available online www.jsaer.om Journal of ientifi and Engineering Researh, 2017, 4(7):299-308 Researh Artile IN: 2394-2630 CODEN(UA): JERBR Bayesian and E-Bayesian estimations for the ompound Rayleigh distribution
More informationINFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE BEHAVIOR OF HYDRAULIC CYLINDER SUBJECTED TO JERKY MOTION
Proeedings of ICFDP 8: 8 th International Congress of Fluid Dynamis & Propulsion Deember 14-17, 006, Sharm El-Shiekh, Sinai, Egypt ICFDP8-EG-154 INFLUENCE OF OPERATING AND CONSTRUCTION PARAMETERS ON THE
More informationA Characterization of Wavelet Convergence in Sobolev Spaces
A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More informationRandom Effects Logistic Regression Model for Ranking Efficiency in Data Envelopment Analysis
Random Effets Logisti Regression Model for Ranking Effiieny in Data Envelopment Analysis So Young Sohn Department of Information & Industrial Systems Engineering Yonsei University, 34 Shinhon-dong, Seoul
More informationSupporting Information for
Eletroni Supplementary Material (ESI) for Nanosale This journal is The Royal Soiety of Chemistry 013 Supporting Information for Exitation polarization modulation in loalization mirosopy allows to resolve
More informationSensor management for PRF selection in the track-before-detect context
Sensor management for PRF seletion in the tra-before-detet ontext Fotios Katsilieris, Yvo Boers, and Hans Driessen Thales Nederland B.V. Haasbergerstraat 49, 7554 PA Hengelo, the Netherlands Email: {Fotios.Katsilieris,
More informationLikelihood-Based Methods
Likelihood-Based Methods Handbook of Spatial Statistics, Chapter 4 Susheela Singh September 22, 2016 OVERVIEW INTRODUCTION MAXIMUM LIKELIHOOD ESTIMATION (ML) RESTRICTED MAXIMUM LIKELIHOOD ESTIMATION (REML)
More informationGreen s function for the wave equation
Green s funtion for the wave equation Non-relativisti ase January 2019 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 43 and 44): 1 2 A 2 2 2 A = µ 0
More informationCalibrating a nanoindenter for very shallow depth indentation using equivalent contact radius
Calibrating a nanoindenter for very shallow depth indentation using equivalent ontat radius Damir R Tadjiev, Russell J Hand, Simon A Hayes To ite this version: Damir R Tadjiev, Russell J Hand, Simon A
More informationRemark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.
Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need
More information