Decomposing compositional data: minimum chi-squared reduced-rank approximations on the simplex
|
|
- MargaretMargaret Long
- 5 years ago
- Views:
Transcription
1 Decomposin compositional data: minimum chi-squared reduced-rank approximations on the simplex Gert Jan Welte Department of Applied Earth Sciences Delft University of Technoloy PO Box 508 NL-600 GA Delft The Netherlands Introduction The loratio transformation (Aitchison, 1981, 1986 opened the way to statistically riorous analysis of compositional data. Most statistical problems in compositional data analysis can be formulated in terms of loratios, and solved accordinly. However, some problems are more easily described and solved in terms of raw compositional variables. Raw compositions with k components are restricted to a k-1 dimensional subspace of the real space, a so-called simplex, defined by the followin two constraints: k p i i1 1 and p i 0. Compositional variation within suites of rocks and rock-formin minerals is often believed to reflect linear mixin. The correspondin mathematical model, expressed in terms of raw compositional variables may be represented in matrix notation as: P MB E where P is a (n k matrix of observed compositions, believed to represent a set of n mixtures, M is a (n q matrix of mixin coefficients, and B a (q k matrix of end members. The number of end members in the mixture is defined as q. E is a (n k matrix of error terms. Note that M and B are also compositions, so that the rows of E are forced to sum to zero. In many cases of linear mixin, B is assumed to be known exactly. For instance, in a procedure known as normative partitionin, chemical compositions of rocks are expressed as mixtures of perfect model minerals, whose chemical compositions are predicted by theory. However, typical mixin problems in sedimentary eoloy are of an explicit nature, characterised by a near total lack of information about the end members. The nonlinear inverse mixin problem in sedimentary eoloy relies on the existence of an unobservable matrix of perfect mixtures, which is defined as an exact convex linear combination of matrices M and B: X q MB
2 where X q is a (n k matrix of perfect mixtures of exact rank q, implyin that the end members are linearly independent. It thus follows that q min( n, k and: P X q E In other words, X q is a reduced-rank approximation of the compositional data P. The main purpose of this contribution is to present a method for estimatin X q which permits statistically riorous testin of the oodness-of-fit. Current estimation methods The estimation methods developed in the earth sciences (Renner, 1993; Welte, 1997 combine factorisation of P by means of a sinular value decomposition with constrained least squares approximation to ensure that all elements of X q are nonneative. This numerical approach is computationally efficient, because there is as yet no need to estimate M and B (Welte, Instead, use can be made of the unidentifiability (nonuniqueness of model parameters to postulate the existence of a (q q transformation matrix T: 1 Xq MB AT TF AF The estimation problem is now reduced to findin X q AF, where F is a (q k matrix of the first q transposed eienvectors of the data matrix P scaled to unit sum, and A is a (n q matrix of coordinates onto the eienvectors. This problem must be solved for each row of P separately: Minimise T T F a p T Subect to a 1 1 Exact linear equality constraints (unit sum Subect to F T a T h Exact linear inequality constraints (nonneativity/positivity Where h is a k-vector of lower bounds on X values: h 0. It follows that the elements of Xq k are restricted to the interval h ; 1 h. 1 This method of reduced-rank approximation may be used to estimate the number of linearly independent end members q needed to model P independently of the actual end-member compositions B. The weakness of this approach is the lack of a formal oodness-of-fit test for X q. Althouh tests with synthetic data provided support for the use of certain empirical measures (Welte, 1995, 1997; Prins & Welte, 1999a, 1999b; Van der Ark, 1999 this problem was not fully solved. If the number of linearly independent end members has been estimated, a physically realistic mixin model of the data can be produced by solvin the followin set of equations: Minimise MT A Subect to 1T 1 Exact linear equality constraints (unit sum 1 Subect to T B F Exact linear equality constraints (unit sum: nonlinear in T Subect to B 0 Exact linear inequality constraints (nonneativity
3 Subect to M 0 Exact linear inequality constraints (nonneativity A detailed account of the procedure for estimatin M and B is beyond the scope of this contribution (for details see Welte, However, it is important to note that it may not be possible to approximate P as a product of two entirely nonneative matrices M and B. Hence, a set of q linearly independent end members may not exist, even if the fit of X q to P is ood. The implicit estimation of X q thus allows the detection of possible row deeneracy amon the true end members B. An estimation method developed for latent budet analysis (De Leeuw & Van der Heden, 1988, an end-member model used in the social sciences, is based on the EM alorithm for aumentation of incomplete data (Dempster et al., This method employs a maximum likelihood criterion to simultaneously estimate X q, M and B, based on the assumption that P follows a product-multinomial distribution. Accordin to this model, each composition (row of P is considered a random sample from an unknown multinomial distribution. Althouh this approach potentially allows a statistically riorous assessment of the discrepancies between X q and P, it is notoriously slow to convere, and it does not allow testin of the rank of X q without estimatin M and B. In the followin sections a new estimation method will be presented that combines the advantaes of the above approaches: Detection of row deeneracy amon end members by implicit estimation of X q ; Minimisation of a statistically riorous measure of discrepancy between X q and P. Distributional assumptions A statistically riorous measure of discrepancy between P and X q may be formulated by reardin compositions as enumerative data, i.e., counts of discrete units. This is strictly true for many types of eoloical data, such as point counts of mineral rains or map reions. However, other types of compositional data, with units of measurement not correspondin to discrete entities, may also be treated as enumerative. The "count lenth" (sample size N of compositional data can be estimated from the spread of a series of replicate analyses, because samples from a population with fixed composition follow a multinomial distribution if all variation is attributable to random error. Goodness-of-fit statistics for the multinomial distribution have been studied extensively by Cressie & Read (1984 and Read & Cressie (1988, who showed that all commonly used statistics can be rearded as members of a one-parameter family of power-diverence statistics. For lare values of N, as commonly encountered in eoloical data, the behaviour of all power-diverence statistics tends towards that of the well-known Pearson's chi-squared statistic (PCS. Hence, for all practical purposes, minimisation of PCS may be rearded as equivalent to maximum likelihood estimation. Estimation of sample size N from a series of m replicate compositions with k constituents proceeds as follows: m k ( p p m ˆ 1 PCS N, where pˆ p, the vector of arithmetic means. i1 1 pˆ m i1 The fit of the replicates to the model is ood by definition:
4 PCS E (Pr0.5, with (m-1(k-1 derees of freedom. By combinin the above expressions we obtain: (Pr0.5 N m k ( p pˆ i1 1 pˆ Strictly speakin, every matrix of compositional data follows a product-multinomial rather than a multinomial distribution, because of the constant row-sum constraint. However, the oodness-of-fit tests for both types of distributions are alebraically equivalent (Read & Cressie, 1988, provided that the number of derees of freedom is reduced to account for this constraint. Every compositional data matrix thus may be treated as a series of samples from one or more fixed populations. In theory, each composition represents a sample that could have been drawn from a different population, but this need not be the case. The main purpose of end-member modellin is to quantify the relations between these populations. For example, if sediments have been sampled across an entire sedimentary basin fed by multiple sediment sources, each sample can be thouht of as representin the local sediment population (Welte, 001. Each of these local populations in turn represents a mixture of a limited number of end-member populations, whose mixin proportions vary continuously in space. There are no riid assumptions with respect to the distribution of mixin proportions on the simplex spanned by the end members. The Minimum Chi-Squared Decomposition alorithm The basic equation for columnwise scalin of the matrix of compositional data is: 1 W PG where G is a diaonal (k k matrix of nonneative weihts. The obective function we seek to minimise is Pearson's chi-squared statistic (PCS: k n ( w wˆ OPCS N 1 i1 wˆ k n ( p x N x 1 i1 We prefer least-squares estimation to maximum likelihood methods, because the latter are not as computationally efficient (Van der Ark, The method of constrained weihted least squares (CWLS minimises the obective function: O WLS k 1 i1 k n 1 i1 n ( w wˆ ( p x subect to nonneativity and unit-sum constraints on the estimated composition vectors x. By equatin these two obective functions we obtain the followin expression for the minimum chi-squared (MCS scalin factor of the -th column: N n i1 n ( p i1 ( p x x x
5 The above expression shows that the reduced-rank estimate X q is required to calculate the MCS scalin factors, implyin that they must be estimated by an iterative process. An initial uess of the MCS scalin factors is based on the fact that the case q = 1 has a simple solution known from chi-squared testin theory (as outlined above for the analysis of replicates. The null hypothesis implies compositional homoeneity of P: H : q 1 0 H1 : q 1 The trivial mixin model with one end member that best fits the data is simply the vector of arithmetic column means. Substitutin the one end-member solution into the eneral expression ives: n p i1 Nn The obective functions are now identical and approximately distributed as χ with (n-1(k-1 derees of freedom. If the null hypothesis is not accepted, more advanced mixin models are needed to account for the variation in P, implyin that q > 1. For the purpose of estimatin X in case q > 1, we define the vector of MCS scalin factors as ~, the vector of diaonal elements of G normalised to unit sum. In other words, we treat ~ as a composition because the oodness-of-fit of the model is solely determined by the relative values of its elements: ~ k 1 The iterative estimation starts by scalin the data usin the current estimate of MCS scalin factors: p w ~ Then the eienvectors V (riht sinular vectors of W are extracted by a sinular value decomposition: T W USV A set of reference vectors F is calculated from the first q columns of V in three steps: First the column sums of vector coefficients are constrained to be nonneative: k 1 v m 0 v m v m 1,,..., k, m 1,,..., q Then a fully nonneative matrix is formed by the linear transformation T qq The mm are calculated from the coefficients of each column of V: VT, where:
6 I v m v 1 mm max v m Where I represents the indicator function: I v0 0 I v0 1 The final step combines scalin and transposition to form F: m fm k ~ 1 m The ~ reference vectors F now obey the followin constraints: F G 1 Linear equality (unit sum F 0 Linear inequality (nonneativity ~ An important result of the above calculations is that X q AFG with unknown coefficients A 1 all of the same order of manitude. It follows that X G ~ AF, implyin that the problem of reduced-rank approximation may be formulated in unweihted form as: Minimise T T F a w T Subect to a 1 1 Linear equality (unit sum Subect to F T a T G ~ 1 h Linear inequality (nonneativity/positivity This is a so-called LSIE problem (least-squares with inequality and equality constraints, which may be solved by applyin the followin transformations: Transformation of a LSIE to a LSI problem by removal of linear equality constraints (LeMaitre, 1979; Transformation of a LSI to a LDP (least-distance prorammin problem of eneral form: Minimise x subect to Ax b (Lawson & Hanson, 1974; Menke, 1989 Transformation of a LDP to a NNLS (nonneative least-squares problem of eneral form: Minimise Cy d subect to y 0 (Lawson & Hanson, 1974; Menke, 1989 The solution to the LSIE problem is obtained by application of the inverse transformations to the NNLS solution. This completes the reduced-rank estimation for the current ~. The vector of MCS scalin coefficients is estimated by Powell's Direction Set Method, a multidimensional minimisation method that does not require calculation of derivatives (Press et al., Multiple random starts may be used to investiate the behaviour of PCS as a function of the composition of ~. Limited experimentation has shown the presence of several local minima, which could be a matter of concern in applications to lare data sets. In the data analysed thus far, multiple minima have comparable PCS values, which merely suests that "valley floors" are not quite flat. The resultin minimum chi-squared estimate is approximately distributed as χ with (n-q(k-q derees of freedom, which allows for riorous oodness-of-fit testin of the reduced rank approximation. The probability of the model X under the null hypothesis as estimated from q
7 PCS is known to deviate from the theoretical χ in cases with low expected frequencies (Read & Cressie, 1988; Greenwood & Nikulin, 1996, which is not uncommon for eoloical data. In such cases an exact chi-squared test should be used (Romesbur et al., References Aitchison. J. (1981 A new approach to null correlations of proportions. Journal of Mathematical Geoloy, 13, Aitchison, J. (1986 The statistical analysis of compositional data. Chapman and Hall, London. Cressie, N. & Read, T.R.C. (1984 Multinomial oodness-of-fit tests. Journal of the Royal Statistical Society, Series B, 46, De Leeuw, J., and Van der Heden, P.G.M. (1988 The analysis of time budets with a latent time-budet model. In: E. Diday et al., Eds., Data analysis and informatics 5. North Holland Publishers, Amsterdam, p Dempster, A.P., Laird, N.M. & Rubin, D.B. (1977 Maximum Likelihood from incomplete data via the EM alorithm. Journal of the Royal Statistical Society, Series B, 39, Greenwood, P.E. & Nikulin, M.S. (1996 A uide to Chi-Squared testin. Wiley & Sons, New York. Lawson, C.L., and Hanson, R.J. (1974 Solvin least squares problems. Prentice-Hall, Enlewood Cliffs. LeMaitre, R.W. (1979 A new eneralised petroloical mixin model. Contributions to Mineraloy and Petroloy, 71, Menke, W. (1989 Geophysical data analysis: discrete inverse theory, revised edition. Academic Press, San Dieo. Press, W.H., Teukolsky, S.A., Vetterlin, W.T. & Flannery, B.P. (1994 Numerical recipes in Fortran: The art of scientific computin, second edition. University Press, Cambride. Prins, M.A. & Welte, G.J. (1999a End-member modelin of siliciclastic rain-size distributions: The late Quaternary record of eolian and fluvial sediment supply to the Arabian Sea and its paleoclimatic sinificance. In: Numerical Experiments in Stratiraphy: Recent Advances in Stratiraphic / Sedimentoloic Computer Simulations (Ed. by J.W. Harbauh, L. Watney, et al., SEPM Special Publication 6, p Prins, M.A. & Welte, G.J. (1999b End-member modellin of rain-size distributions of sediment mixtures. In: Prins, M.A. (1999 Pelaic, hemipelaic and turbidite deposition in the Arabian Sea durin the late Quaternary [PhD Thesis Utrecht University]. Geoloica Ultraiectina 168, p Read, T.R.C. & Cressie, N. (1988 Goodness-of-fit statistics for discrete multivariate data. Spriner Verla, New York. Renner, R.M. (1993 The resolution of a compositional data set into mixtures of fixed source compositions. Applied Statistics, Journal of the Royal Statistical Society, Series C, 4, Romesbur, H.C., Marshall, K. & Mauk, T.P. (1981 FITEST: a computer proram for "exact chi-square" oodness-of-fit sinificance tests. Computers and Geosciences, 7, Van der Ark, L.A. (1999 Contributions to latent budet analysis: a tool for the analysis of compositional data [Ph.D. Thesis Utrecht University], DSWO Press, Leiden. Welte, G.J. (1995 Unravellin mixed provenance of coastal sands: The Po delta and adacent beaches of the northern Adriatic sea as a test case. In: Geoloy of deltas (M.A. Oti & G. Postma, Eds., Balkema, Rotterdam, p
8 Welte, G.J. (1997 End-member modellin of compositional data: Numerical-statistical alorithms for solvin the explicit mixin problem. Journal of Mathematical Geoloy, 9, Welte, G.J. (001 Quantitative analysis of detrital modes: statistically riorous confidence reions in ternary diarams and their use in sedimentary petroloy. Earth-Science Reviews, in press.
Improved transformation of ϕ-divergence goodness-of-fit test statistics based on minimum ϕ -divergence estimator for GLIM of binary data
SUT Journal of Mathematics Vol. 52, No. 2 (216), 193 214 Improved transformation of ϕ-diverence oodness-of-fit test statistics based on minimum ϕ -diverence estimator for GLIM of binary data Nobuhiro Taneichi,
More informationMODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION
= = = MODEL SELECTION CRITERIA FOR ACOUSTIC SEGMENTATION Mauro Cettolo and Marcello Federico ITC-irst - Centro per la Ricerca Scientifica e Tecnoloica I-385 Povo, Trento, Italy ABSTRACT Robust acoustic
More informationCorrespondence Analysis of Longitudinal Data
Correspondence Analysis of Longitudinal Data Mark de Rooij* LEIDEN UNIVERSITY, LEIDEN, NETHERLANDS Peter van der G. M. Heijden UTRECHT UNIVERSITY, UTRECHT, NETHERLANDS *Corresponding author (rooijm@fsw.leidenuniv.nl)
More informationGeneralized Least-Squares Regressions V: Multiple Variables
City University of New York (CUNY) CUNY Academic Works Publications Research Kinsborouh Community Collee -05 Generalized Least-Squares Reressions V: Multiple Variables Nataniel Greene CUNY Kinsborouh Community
More informationGenetic Algorithm Approach to Nonlinear Blind Source Separation
Genetic Alorithm Approach to Nonlinear Blind Source Separation F. Roas, C.G. Puntonet, I.Roas, J. Ortea, A. Prieto Dept. of Computer Architecture and Technoloy, University of Granada. fernando@atc.ur.es,
More informationLinearized optimal power flow
Linearized optimal power flow. Some introductory comments The advantae of the economic dispatch formulation to obtain minimum cost allocation of demand to the eneration units is that it is computationally
More informationStatistics Hotelling s T Gary W. Oehlert School of Statistics 313B Ford Hall
\ C f C E A tatistics 5041 11 Hotellin s T Gary W ehlert chool of tatistics 313B Ford Hall 612-625-155 ary@statumnedu Let s think about the univariate -test First recall that there are one-sample tests,
More informationAn EM Algorithm for the Student-t Cluster-Weighted Modeling
An EM Alorithm for the Student-t luster-weihted Modelin Salvatore Inrassia, Simona. Minotti, and Giuseppe Incarbone Abstract luster-weihted Modelin is a flexible statistical framework for modelin local
More informationStat260: Bayesian Modeling and Inference Lecture Date: March 10, 2010
Stat60: Bayesian Modelin and Inference Lecture Date: March 10, 010 Bayes Factors, -priors, and Model Selection for Reression Lecturer: Michael I. Jordan Scribe: Tamara Broderick The readin for this lecture
More informationAsymptotic Behavior of a t Test Robust to Cluster Heterogeneity
Asymptotic Behavior of a t est Robust to Cluster Heteroeneity Andrew V. Carter Department of Statistics University of California, Santa Barbara Kevin. Schnepel and Doulas G. Steierwald Department of Economics
More informationA Multigrid-like Technique for Power Grid Analysis
A Multirid-like Technique for Power Grid Analysis Joseph N. Kozhaya, Sani R. Nassif, and Farid N. Najm 1 Abstract Modern sub-micron VLSI desins include hue power rids that are required to distribute lare
More informationAdjustment of Sampling Locations in Rail-Geometry Datasets: Using Dynamic Programming and Nonlinear Filtering
Systems and Computers in Japan, Vol. 37, No. 1, 2006 Translated from Denshi Joho Tsushin Gakkai Ronbunshi, Vol. J87-D-II, No. 6, June 2004, pp. 1199 1207 Adjustment of Samplin Locations in Rail-Geometry
More informationNonlinear Model Reduction of Differential Algebraic Equation (DAE) Systems
Nonlinear Model Reduction of Differential Alebraic Equation DAE Systems Chuili Sun and Jueren Hahn Department of Chemical Enineerin eas A&M University Collee Station X 77843-3 hahn@tamu.edu repared for
More information1 Introduction. 2 A regression model
Regression Analysis of Compositional Data When Both the Dependent Variable and Independent Variable Are Components LA van der Ark 1 1 Tilburg University, The Netherlands; avdark@uvtnl Abstract It is well
More informationLecture 5 Processing microarray data
Lecture 5 Processin microarray data (1)Transform the data into a scale suitable for analysis ()Remove the effects of systematic and obfuscatin sources of variation (3)Identify discrepant observations Preprocessin
More informationMatrix multiplication: a group-theoretic approach
CSG399: Gems of Theoretical Computer Science. Lec. 21-23. Mar. 27-Apr. 3, 2009. Instructor: Emanuele Viola Scribe: Ravi Sundaram Matrix multiplication: a roup-theoretic approach Given two n n matrices
More informationKeywords: Solid-Fluid Equilibrium, Solid-Fluid-Fluid Equilibrium, Continuation Method, binary n-alkane systems.
APPLICATION OF CONTINUATION METHOD TO THE CALCULATION OF OLID -FLUID -FLUID AND RELATED OLID -FLUID EQUILIBRIA OF BINARY AYMMETRIC MIXTURE IN WIDE RANGE OF CONDITION.B. Rodriuez-Reartes (1), M. Cismondi
More informationA Mathematical Model for the Fire-extinguishing Rocket Flight in a Turbulent Atmosphere
A Mathematical Model for the Fire-extinuishin Rocket Fliht in a Turbulent Atmosphere CRISTINA MIHAILESCU Electromecanica Ploiesti SA Soseaua Ploiesti-Tiroviste, Km 8 ROMANIA crismihailescu@yahoo.com http://www.elmec.ro
More informationNumerical and Experimental Investigations of Lateral Cantilever Shaft Vibration of Passive-Pitch Vertical-Axis Ocean Current
R. Hantoro, et al. / International Enery Journal 1 (011) 191-00 191 Numerical and Experimental Investiations of Lateral Cantilever Shaft Vibration of Passive-Pitch Vertical-Axis Ocean Current R. Hantoro
More informationTwo new spectral conjugate gradient algorithms based on Hestenes Stiefel
Research Article Two new spectral conjuate radient alorithms based on Hestenes Stiefel Journal of Alorithms & Computational Technoloy 207, Vol. (4) 345 352! The Author(s) 207 Reprints and permissions:
More informationEfficient crystal size distribution estimation approach for growth dominated crystallisation processes
Louhborouh University Institutional Repository Efficient crystal size distribution estimation approach for rowth dominated crystallisation processes This item was submitted to Louhborouh University's Institutional
More informationSimple Optimization (SOPT) for Nonlinear Constrained Optimization Problem
(ISSN 4-6) Journal of Science & Enineerin Education (ISSN 4-6) Vol.,, Pae-3-39, Year-7 Simple Optimization (SOPT) for Nonlinear Constrained Optimization Vivek Kumar Chouhan *, Joji Thomas **, S. S. Mahapatra
More informationarxiv: v2 [math.oc] 11 Jun 2018
: Tiht Automated Converence Guarantees Adrien Taylor * Bryan Van Scoy * 2 Laurent Lessard * 2 3 arxiv:83.673v2 [math.oc] Jun 28 Abstract We present a novel way of eneratin Lyapunov functions for provin
More informationAMERICAN INSTITUTES FOR RESEARCH
AMERICAN INSTITUTES FOR RESEARCH LINKING RASCH SCALES ACROSS GRADES IN CLUSTERED SAMPLES Jon Cohen, Mary Seburn, Tamas Antal, and Matthew Gushta American Institutes for Research May 23, 2005 000 THOMAS
More information10log(1/MSE) log(1/MSE)
IROVED MATRI PENCIL METHODS Biao Lu, Don Wei 2, Brian L Evans, and Alan C Bovik Dept of Electrical and Computer Enineerin The University of Texas at Austin, Austin, T 7872-84 fblu,bevans,bovik@eceutexasedu
More informationA Performance Comparison Study with Information Criteria for MaxEnt Distributions
A Performance Comparison Study with nformation Criteria for MaxEnt Distributions Ozer OZDEMR and Aslı KAYA Abstract n statistical modelin, the beinnin problem that has to be solved is the parameter estimation
More informationNumber of cases (objects) Number of variables Number of dimensions. n-vector with categorical observations
PRINCALS Notation The PRINCALS algorithm was first described in Van Rickevorsel and De Leeuw (1979) and De Leeuw and Van Rickevorsel (1980); also see Gifi (1981, 1985). Characteristic features of PRINCALS
More informationInvestigation of ternary systems
Investiation of ternary systems Introduction The three component or ternary systems raise not only interestin theoretical issues, but also have reat practical sinificance, such as metallury, plastic industry
More informationRobust Semiparametric Optimal Testing Procedure for Multiple Normal Means
Veterinary Dianostic and Production Animal Medicine Publications Veterinary Dianostic and Production Animal Medicine 01 Robust Semiparametric Optimal Testin Procedure for Multiple Normal Means Pen Liu
More informationAG DANK/BCS Meeting 2013 in London University College London, 8/9 November 2013
AG DANK/S Meetin 3 in London University ollee London 8/9 November 3 MDELS FR SIMULTANEUS LASSIFIATIN AND REDUTIN F THREE-WAY DATA Roberto Rocci University Tor erata Rome A eneral classification model:
More informationStability Analysis of Nonlinear Systems using Dynamic-Routh's Stability Criterion: A New Approach
Stability Analysis of Nonlinear Systems usin Dynamic-Routh's Stability Criterion: A New Approach Bast Kumar Sahu, St. M., IEEE Centre for Industrial Electronics d Robotics National Institute of Technoloy
More informationTermination criteria in the Moore-Skelboe Algorithm for Global Optimization by Interval Arithmetic
To appear in series Nonconvex Global Optimization and Its Applications Kluwer Academic Publishers Termination criteria in the Moore-Skelboe Alorithm for Global Optimization by Interval Arithmetic M.H.
More informationGeophysical Data Analysis: Discrete Inverse Theory
Geophysical Data Analysis: Discrete Inverse Theory MATLAB Edition William Menke Lamont-Doherty Earth Observatory and Department of Earth and Environmental Sciences Columbia University. ' - Palisades, New
More informationLevenberg-Marquardt-based OBS Algorithm using Adaptive Pruning Interval for System Identification with Dynamic Neural Networks
Proceedins of the 29 IEEE International Conference on Systems, Man, and Cybernetics San Antonio, X, USA - October 29 evenber-marquardt-based OBS Alorithm usin Adaptive Prunin Interval for System Identification
More informationPHY 133 Lab 1 - The Pendulum
3/20/2017 PHY 133 Lab 1 The Pendulum [Stony Brook Physics Laboratory Manuals] Stony Brook Physics Laboratory Manuals PHY 133 Lab 1 - The Pendulum The purpose of this lab is to measure the period of a simple
More informationNormed Distances and Their Applications in Optimal Circuit Design
Optimization and Enineerin, 4, 197 213, 2003 c 2003 Kluwer Academic Publishers. Manufactured in he Netherlands Normed Distances and heir Applications in Optimal Circuit Desin ABDEL-KARIM S.O. HASSAN Department
More informationRound-off Error Free Fixed-Point Design of Polynomial FIR Predictors
Round-off Error Free Fixed-Point Desin of Polynomial FIR Predictors Jarno. A. Tanskanen and Vassil S. Dimitrov Institute of Intellient Power Electronics Department of Electrical and Communications Enineerin
More informationSolution to the take home exam for ECON 3150/4150
Solution to the tae home exam for ECO 350/450 Jia Zhiyan and Jo Thori Lind April 2004 General comments Most of the copies we ot were quite ood, and it seems most of you have done a real effort on the problem
More informationChapter 5 Confidence Intervals
Chapter 5 Confidence Intervals Confidence Intervals about a Population Mean, σ, Known Abbas Motamedi Tennessee Tech University A point estimate: a single number, calculated from a set of data, that is
More informationDigital Electronics Paper-EE-204-F SECTION-A
B.Tech 4 th Semester (AEIE) F Scheme, May 24 Diital Electronics Paper-EE-24-F Note : Attempt five questions. Question is compulsory and one question from each of the four sections.. What is a loic ate?
More informationLasso-based linear regression for interval-valued data
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session CPS065) p.5576 Lasso-based linear regression for interval-valued data Giordani, Paolo Sapienza University of Rome, Department
More informationover the parameters θ. In both cases, consequently, we select the minimizing
MONOTONE REGRESSION JAN DE LEEUW Abstract. This is an entry for The Encyclopedia of Statistics in Behavioral Science, to be published by Wiley in 200. In linear regression we fit a linear function y =
More informationThe Computational Complexity Analysis of a MINLP-Based Chemical Process Control Design
J.Oto.Ktrl.Inst (J.Auto.Ctrl.Inst Vol (, 00 ISSN : 085-57 Abstract The Computational Complexity Analysis of a MINLP-Based Chemical Process Control Desin E. Ekawati, S. Yuliar Center for Instrumentation
More informationRATE OPTIMIZATION FOR MASSIVE MIMO RELAY NETWORKS: A MINORIZATION-MAXIMIZATION APPROACH
RATE OPTIMIZATION FOR MASSIVE MIMO RELAY NETWORKS: A MINORIZATION-MAXIMIZATION APPROACH Mohammad Mahdi Nahsh, Mojtaba Soltanalian, Petre Stoica +, Maryam Masjedi, and Björn Ottersten Department of Electrical
More informationExtended Target Poisson Multi-Bernoulli Filter
1 Extended Taret Poisson Multi-Bernoulli Filter Yuxuan Xia, Karl Granström, Member, IEEE, Lennart Svensson, Senior Member, IEEE, and Maryam Fatemi arxiv:1801.01353v1 [eess.sp] 4 Jan 2018 Abstract In this
More informationTwo-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix
Two-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix René Vidal Stefano Soatto Shankar Sastry Department of EECS, UC Berkeley Department of Computer Sciences, UCLA 30 Cory Hall,
More informationThis paper is not to be removed from the Examination Halls
~~MT1173 ZA d0 This paper is not to be removed from the Examination Halls UNIVERSITY OF LONDON MT1173 ZA BSc degrees and Diplomas for Graduates in Economics, Management, Finance and the Social Sciences,
More informationPhase Diagrams: construction and comparative statics
1 / 11 Phase Diarams: construction and comparative statics November 13, 215 Alecos Papadopoulos PhD Candidate Department of Economics, Athens University of Economics and Business papadopalex@aueb.r, https://alecospapadopoulos.wordpress.com
More informationThe Expectation Maximization Algorithm
The Expectation Maximization Algorithm Frank Dellaert College of Computing, Georgia Institute of Technology Technical Report number GIT-GVU-- February Abstract This note represents my attempt at explaining
More informationNumerical Analysis: Solving Systems of Linear Equations
Numerical Analysis: Solving Systems of Linear Equations Mirko Navara http://cmpfelkcvutcz/ navara/ Center for Machine Perception, Department of Cybernetics, FEE, CTU Karlovo náměstí, building G, office
More informationLame equation, quantum top and elliptic Bernoulli polynomials.
Louhborouh University Institutional Repository Lame equation, quantum top and elliptic Bernoulli polynomials. This item was submitted to Louhborouh University's Institutional Repository by the/an author.
More informationA Social Interactions Model with Endogenous Friendship Formation and Selectivity*
A Social Interactions Model with Endoenous Friendship Formation and Selectivity* Chih-Shen Hsieh and Lun-Fei Lee Department of Economics Ohio State University Columbus, Ohio, USA Date: Auust, 2011 Abstract
More informationEquivalent rocking systems: Fundamental rocking parameters
Equivalent rockin systems: Fundamental rockin parameters M.J. DeJon University of Cambride, United Kindom E.G. Dimitrakopoulos The Hon Kon University of Science and Technoloy SUMMARY Early analytical investiations
More informationModel Estimation Example
Ronald H. Heck 1 EDEP 606: Multivariate Methods (S2013) April 7, 2013 Model Estimation Example As we have moved through the course this semester, we have encountered the concept of model estimation. Discussions
More informationUSE OF FILTERED SMITH PREDICTOR IN DMC
Proceedins of the th Mediterranean Conference on Control and Automation - MED22 Lisbon, Portual, July 9-2, 22. USE OF FILTERED SMITH PREDICTOR IN DMC C. Ramos, M. Martínez, X. Blasco, J.M. Herrero Predictive
More informationFLAWS IN BLACK HOLE THEORY AND GENERAL RELATIVITY
FLAWS IN BLACK HOLE THEORY AND GENERAL RELATIVITY Stephen J. Crothers Queensland, Australia http://vixra.or/author/stephen_j_crothers 1 /15 Black Hole and Bi Ban Mutually Exclusive BLACK HOLE UNIVERSE
More informationAN ALTERNATING MINIMIZATION ALGORITHM FOR NON-NEGATIVE MATRIX APPROXIMATION
AN ALTERNATING MINIMIZATION ALGORITHM FOR NON-NEGATIVE MATRIX APPROXIMATION JOEL A. TROPP Abstract. Matrix approximation problems with non-negativity constraints arise during the analysis of high-dimensional
More informationREPORT DOCUMENTATION PAGE
REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-088 Public reportin burden for this collection of information is estimated to averae hour per response, includin the time for reviewin instructions,
More informationModeling for control of a three degrees-of-freedom Magnetic. Levitation System
Modelin for control of a three derees-of-freedom Manetic evitation System Rafael Becerril-Arreola Dept. of Electrical and Computer En. University of Toronto Manfredi Maiore Dept. of Electrical and Computer
More informationMODULE - 2 LECTURE NOTES 3 LAGRANGE MULTIPLIERS AND KUHN-TUCKER CONDITIONS
Water Resources Systems Plannin and Manaement: Introduction to Optimization: arane Multipliers MODUE - ECTURE NOTES 3 AGRANGE MUTIPIERS AND KUHN-TUCKER CONDITIONS INTRODUCTION In the previous lecture the
More informationFirst- and Second Order Phase Transitions in the Holstein- Hubbard Model
Europhysics Letters PREPRINT First- and Second Order Phase Transitions in the Holstein- Hubbard Model W. Koller 1, D. Meyer 1,Y.Ōno 2 and A. C. Hewson 1 1 Department of Mathematics, Imperial Collee, London
More informationNonlinear Support Vector Machines through Iterative Majorization and I-Splines
Nonlinear Support Vector Machines through Iterative Majorization and I-Splines P.J.F. Groenen G. Nalbantov J.C. Bioch July 9, 26 Econometric Institute Report EI 26-25 Abstract To minimize the primal support
More informationANALYZE In all three cases (a) (c), the reading on the scale is. w = mg = (11.0 kg) (9.8 m/s 2 ) = 108 N.
Chapter 5 1. We are only concerned with horizontal forces in this problem (ravity plays no direct role). We take East as the +x direction and North as +y. This calculation is efficiently implemented on
More information1 Introduction. Korea Atomic Energy Research Institute, South Korea.
Tenth International Conference on Computational Fluid Dynamics (ICCFD10), Barcelona, Spain, July 9-13, 2018 ICCFD10-0206 Multi-physics approach for nuclear reactor analysis usin thermal-hydraulics and
More informationSmarandache Curves According to Sabban Frame on
Smarandache Curves Accordin to Sabban Frame on S Kemal Taşköprü, Murat Tosun Faculty of Arts and Sciences, Department of Mathematics Sakarya University, Sakarya 5487 TURKEY Abstract: In this paper, we
More informationStatistical methods and data analysis
Statistical methods and data analysis Teacher Stefano Siboni Aim The aim of the course is to illustrate the basic mathematical tools for the analysis and modelling of experimental data, particularly concerning
More informationStrong Interference and Spectrum Warfare
Stron Interference and Spectrum Warfare Otilia opescu and Christopher Rose WILAB Ruters University 73 Brett Rd., iscataway, J 8854-86 Email: {otilia,crose}@winlab.ruters.edu Dimitrie C. opescu Department
More informationEmbedded ADMM-based QP Solver for MPC with polytopic constraints
25 European Control Conference (ECC) July 5-7, 25. Linz, Austria Embedded ADMM-based QP Solver for MPC with polytopic constraints Thuy V. Dan, K.V. Lin and J.M. Maciejowsi Abstract We propose an alorithm
More informationSigmaplot di Systat Software
Sigmaplot di Systat Software SigmaPlot Has Extensive Statistical Analysis Features SigmaPlot is now bundled with SigmaStat as an easy-to-use package for complete graphing and data analysis. The statistical
More informationMODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM. Apurba Sakti EGEE 520, Mathematical Modeling of EGEE systems Spring 2007
MODELING A METAL HYDRIDE HYDROGEN STORAGE SYSTEM Apurba Sakti EGEE 520, Mathematical Modelin of EGEE systems Sprin 2007 Table of Contents Abstract Introduction Governin Equations Enery balance Momentum
More informationExperimental and Computational Studies of Gas Mixing in Conical Spouted Beds
Refereed Proceedins The 1th International Conference on Fluidization - New Horizons in Fluidization Enineerin Enineerin Conferences International Year 007 Experimental and Computational Studies of Gas
More informationPERFORMANCE OF THE EM ALGORITHM ON THE IDENTIFICATION OF A MIXTURE OF WAT- SON DISTRIBUTIONS DEFINED ON THE HYPER- SPHERE
REVSTAT Statistical Journal Volume 4, Number 2, June 2006, 111 130 PERFORMANCE OF THE EM ALGORITHM ON THE IDENTIFICATION OF A MIXTURE OF WAT- SON DISTRIBUTIONS DEFINED ON THE HYPER- SPHERE Authors: Adelaide
More informationECEN 5005 Crystals, Nanocrystals and Device Applications Class 18 Group Theory For Crystals
ECEN 5005 Crystals, Nanocrystals and Device Applications Class 18 Group Theory For Crystals Multi-Electron Crystal Field Theory Weak Field Scheme Stron Field Scheme Tanabe-Suano Diaram 1 Notation Convention
More informationDynamic sequential analysis of careers
Dynamic sequential analysis of careers Fulvia Pennoni Department of Statistics and Quantitative Methods University of Milano-Bicocca http://www.statistica.unimib.it/utenti/pennoni/ Email: fulvia.pennoni@unimib.it
More informationMULTIDIMENSIONAL COUPLED PHOTON-ELECTRON TRANSPORT SIMULATIONS USING NEUTRAL PARTICLE S N CODES
Computational Medical Physics Workin Group Workshop II ep 30 Oct 3 007 University of Florida (UF) Gainesville Florida UA on CD-ROM American Nuclear ociety LaGrane Park IL (007) MULTIDIMNIONAL COUPLD PHOTON-LCTRON
More informationA New Nonlinear Conjugate Gradient Coefficient. for Unconstrained Optimization
Applie Mathematical Sciences, Vol. 9, 05, no. 37, 83-8 HIKARI Lt, www.m-hiari.com http://x.oi.or/0.988/ams.05.4994 A New Nonlinear Conjuate Graient Coefficient for Unconstraine Optimization * Mohame Hamoa,
More informationarxiv:quant-ph/ v3 20 May 2005
For Distinuishin Conjuate Hidden Subroups, the Pretty Good Measurement is as Good as it Gets arxiv:quant-ph/0501177v3 20 May 2005 Cristopher Moore moore@cs.unm.edu Department of Computer Science University
More informationDynamical Diffraction
Dynamical versus Kinematical Di raction kinematical theory valid only for very thin crystals Dynamical Diffraction excitation error s 0 primary beam intensity I 0 intensity of other beams consider di racted
More informationModeling Multiscale Differential Pixel Statistics
Modeling Multiscale Differential Pixel Statistics David Odom a and Peyman Milanfar a a Electrical Engineering Department, University of California, Santa Cruz CA. 95064 USA ABSTRACT The statistics of natural
More informationA turbulence closure based on the maximum entropy method
Advances in Fluid Mechanics IX 547 A turbulence closure based on the maximum entropy method R. W. Derksen Department of Mechanical and Manufacturing Engineering University of Manitoba Winnipeg Canada Abstract
More informationStatistics. Lent Term 2015 Prof. Mark Thomson. 2: The Gaussian Limit
Statistics Lent Term 2015 Prof. Mark Thomson Lecture 2 : The Gaussian Limit Prof. M.A. Thomson Lent Term 2015 29 Lecture Lecture Lecture Lecture 1: Back to basics Introduction, Probability distribution
More information2.1.3 The Testing Problem and Neave s Step Method
we can guarantee (1) that the (unknown) true parameter vector θ t Θ is an interior point of Θ, and (2) that ρ θt (R) > 0 for any R 2 Q. These are two of Birch s regularity conditions that were critical
More informationOptimization in Predictive Control Algorithm
Latest rends in Circits, Systems, Sinal Processin and Atomatic Control Optimization in Predictive Control Alorithm JAN ANOŠ, MAREK KUBALČÍK omas Bata University in Zlín, Faclty of Applied Informatics Nám..
More informationBidirectional Clustering of Weights for Finding Succinct Multivariate Polynomials
IJCSNS International Journal of Computer Science and Network Security, VOL.8 No.5, May 28 85 Bidirectional Clusterin of Weihts for Findin Succinct Multivariate Polynomials Yusuke Tanahashi and Ryohei Nakano
More informationSYSTEMS ANALYSIS HOPF BIFURCATION IN COASTAL ECOGEOSYSTEMS. I. T. Selezov, 1 Iu. G. Kryvonos, 2 and V. M. Moskovkin 3 UDC
DOI.7/s559-6-9856-8 Cybernetics and Systems Analysis Vol. 5 No. 4 July 6 SYSTEMS ANALYSIS HOPF BIFURCATION IN COASTAL ECOGEOSYSTEMS I. T. Selezov Iu. G. Kryvonos and V. M. Moskovkin 3 UDC 53.59 Abstract.
More informationBuckling analysis of thin-walled members via semi-analytical finite strip transfer matrix method
Special Issue Article Bucklin analysis of thin-walled members via semi-analytical finite strip transfer matrix method Advances in Mechanical Enineerin 16, ol. 8() 1 11 Ó The Author(s) 16 DOI: 1.11/168811661
More informationRecovery of Piecewise Smooth Images from Few Fourier Samples
Recovery of Piecewise Smooth Imaes from Few Fourier Samples Gre Onie Department of Mathematics University of Iowa Iowa City, Iowa 54 Email: reory-onie@uiowa.edu Mathews Jacob Department of Electric and
More informationExpanded Knowledge on Orifice Meter Response to Wet Gas Flows
32 nd International North Sea Flow Measurement Workshop 21-24 October 2014 Expanded Knowlede on Orifice Meter Response to Wet Gas Flows Richard Steven, Colorado Enineerin Experiment Station Inc Josh Kinney,
More informationLinear Systems. Carlo Tomasi. June 12, r = rank(a) b range(a) n r solutions
Linear Systems Carlo Tomasi June, 08 Section characterizes the existence and multiplicity of the solutions of a linear system in terms of the four fundamental spaces associated with the system s matrix
More informationIEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING 1
IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING 1 Intervention in Gene Reulatory Networks via a Stationary Mean-First-Passae-Time Control Policy Golnaz Vahedi, Student Member, IEEE, Babak Faryabi, Student
More informationAn Optimization Problem With a Surprisingly Simple Solution
An Optimization Problem With a Surprisinly Simple Solution D. Drinen, K. G. Kennedy, and W. M. Priestley 1 INTRODUCTION. Suppose you and n of your friends play the followin ame. A random number from the
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 4 for Applied Multivariate Analysis Outline 1 Eigen values and eigen vectors Characteristic equation Some properties of eigendecompositions
More informationNing Wu Institute for Traffic Engineering Ruhr University Bochum, Germany Tel: ; Fax: ;
MODELLING THE IMPACT OF SIDE-STREET TRAFFIC VOLUME ON MAJOR- STREET GREEN TIME AT ISOLATED SEMI-ACTUATED INTERSECTIONS FOR SIGNAL COORDINATION DECISIONS Donmei Lin, Correspondin Author Center for Advanced
More informationarxiv: v1 [hep-th] 16 Feb 2018
Fast Summation of Diverent Series and Resurent Transseries in Quantum Field Theories from Meijer-G Approximants arxiv:1802.06034v1 [hep-th] 16 Feb 2018 Héctor Mera, 1, 2, 3, Thomas G. Pedersen, 2, 3 and
More informationGENOME WIDE ASSOCIATION ANALYSES BASED ON BROADLY DIFFERENT SPECIFICATIONS FOR PRIOR DISTRIBUTIONS, GENOMIC WINDOWS, AND ESTIMATION METHODS
Genetics: Early Online, published on June 1, 017 as 10.1534/enetics.117.059 1 3 4 GENOME WIDE ASSOCIATION ANALYSES BASED ON BROADLY DIFFERENT SPECIFICATIONS FOR PRIOR DISTRIBUTIONS, GENOMIC WINDOWS, AND
More informationUniversity of Luxembourg. Master in Mathematics. Student project. Compressed sensing. Supervisor: Prof. I. Nourdin. Author: Lucien May
University of Luxembourg Master in Mathematics Student project Compressed sensing Author: Lucien May Supervisor: Prof. I. Nourdin Winter semester 2014 1 Introduction Let us consider an s-sparse vector
More informationOn the Unified Approach of the Least Squares Method
On the Unified Approach of the Least Squares Method Dimitrios A. Rossikopoulos Department of Geodesy and Surveyin, School of Rural and Surveyin Enineerin, he Aristotle University of hessaloniki. E-mail:
More informationAdvanced Methods Development for Equilibrium Cycle Calculations of the RBWR. Andrew Hall 11/7/2013
Advanced Methods Development for Equilibrium Cycle Calculations of the RBWR Andrew Hall 11/7/2013 Outline RBWR Motivation and Desin Why use Serpent Cross Sections? Modelin the RBWR Axial Discontinuity
More informationLeast squares multidimensional scaling with transformed distances
Least squares multidimensional scaling with transformed distances Patrick J.F. Groenen 1, Jan de Leeuw 2 and Rudolf Mathar 3 1 Department of Data Theory, University of Leiden P.O. Box 9555, 2300 RB Leiden,
More informationPackage lsei. October 23, 2017
Package lsei October 23, 2017 Title Solving Least Squares or Quadratic Programming Problems under Equality/Inequality Constraints Version 1.2-0 Date 2017-10-22 Description It contains functions that solve
More information