PROXSCAL. Notation. W n n matrix with weights for source k. E n s matrix with raw independent variables F n p matrix with fixed coordinates
|
|
- Andra Garrett
- 6 years ago
- Views:
Transcription
1 PROXSCAL PROXSCAL perfors ultidiensional scaling of proxiity data to find a leastsquares representation of the obects in a low-diensional space. Individual differences odels can be specified for ultiple sources. A aorization algorith guarantees onotone convergence for optionally transfored, etric and nonetric data under a variety of odels and constraints. Detailled atheatical derivations concerning the algorith can be found in Coandeur and Heiser (1993). Notation he following notation is used throughout this chapter, unless stated otherwise. For the diensions of the vectors and atrices are used: n Nuber of obects Nuber of sources p Nuber of diensions s Nuber of independent variables h axiu(s, p) l Length of transforation vector r Degree of spline t Nuber of interior nots for spline he input and input-related variables are: n n atrix with raw proxiities for source W n n atrix with weights for source E n s atrix with raw independent variables F n p atrix with fixed coordinates Output and output-related variables are: 1
2 D ˆ n n atrix with transfored proxiities for source Z n p atrix with coon space coordinates A p p atrix with space weights for source X n p atrix with individual space coordinates for source Q n h atrix with transfored independent variables B h p atrix with regression weights for independent variables S l ( r+ t) atrix of coefficients for the spline basis Special atrices and functions are: J D( X ) n n V n n B( X ) n n I 11 / 1 1, centering atrix of appropriate size d, atrix with distances, with eleents { i } where d i = ( xi x )( xi x ) atrix with eleents { v i } atrix with eleents { b i } wi f ( δi ) ( X ) wi for i n, where vi = wil for i = l i, where if di > 0 and i di ( X ) bi = 0 if di ( X ) = 0 and i n bil if i = l i Introduction he following loss function is iniized by PROXSCAL, n 1 σ ˆ w d i di i X, (1.1) = 1 i< which is the weighted ean squared error between the transfored proxiities and the distances of n obect within sources. he transforation function for the
3 3 Preliinaries proxiities provides nonnegative, onotonically nondecreasing values for the transfored proxiities d ˆi. he distances d i ( X ) are siply the Euclidean distances between the obect points, with the coordinates in the rows of X. he ain algorith consists of the following aor steps: 1. find initial configurations X, and evaluate the loss function;. find an update for the configurations X ; 3. find an update for the transfored proxiities d ˆi ; 4. evaluate the loss function; if soe predefined stop criterion is satisfied, stop; otherwise, go to step. At the start of the procedure, several preliinary coputations are perfored to handle issing weights or proxiities, and initialize the raw proxiities. Missings On input, issing values ay occur for both weights and proxiities. If a weight is issing, it is set equal to zero. If a proxiity is issing, the corresponding weight is set equal to zero. Proxiities Only the upper or lower triangular part (without the diagonal) of the proxiity atrix is needed. In case both triangles are given, the weighted ean of both triangles is used. Next, the raw proxiities are transfored such that siilarities becoe dissiilarities by ultiplying with 1, taing into account the conditionality, and setting the sallest dissiilarity equal to zero. ransforations For ordinal transforations, the nonissing proxiities are replaced by their ascending ran nubers, also taing into account the conditionality. For spline transforations, the spline basis S is coputed.
4 4 Noralization he proxiities are noralized such that the weighted squared proxiities equal the su of the weights, again, taing into account the conditionality. Step 1: Initial Configuration PROXSCAL allows for several initial configurations. Before deterining the initial configuration, issings are handled, and the raw proxiities are initialized. Finally, after one of the starts described below, the coon space Z is centered on the origin and optially dilated in accordance with the noralized proxiities. Siplex Start he siplex start consists of a ran-p approxiation of the atrix VBJ. Set H, an n p colunwise orthogonal atrix, satisfying HH= I p equal to I p, where I p denotes the atrix with the first p coluns of the identity atrix. he nonzero rows are selected in such a way that the first Z = B( J) H contains the p coluns of BJ with the largest diagonal eleents. he following steps are coputed in turn, until convergence is reached: 1. For a fixed Z, H = PQ, where PQ is taen fro the singular value decoposition BJZ = PLQ ;. For a fixed H, Z = 1/ V B J H, where V is the pseudo-inverse of V. For a restricted coon space Z, the second step is adusted in order to fullfill the restictions. his procedure was introduced in Heiser (1985). orgerson Start he proxiities are aggregated over sources, squared, double centered and ultiplied with 0.5, after which an eigenvalue decoposition is used to deterine the coordinate values, thus * 0.5JD J = Q 4, where eleents of * D are defined as * ˆ di = widi wi = 1 = 1 1
5 5 1/ followed by Z = Q, where only the first p positive ordered eigenvalues ( λ 1 λ K λ n ) and eigenvectors are used. his technique, classical scaling, is due to orgerson (195, 1958) and Gower (1966) and also nown under the naes orgerson scaling or orgerson-gower scaling. (Multiple) Rando Start he coordinate values are randoly generated fro a unifor distribution using the default rando nuber generator fro the SPSS syste. User-Provided Start he coordinate values provided by the user are used. Step : Configuration Update Update for the Coon Space X he coon space Z is related to the individual spaces = 1,..., through the odel X = ZA, where A are atrices containing space weights. Assue that weight atrix A is of full ran. Only considering Z defines (1.1) as σ z = c + z Hz z t, (3.1) where ( Z) z vec, 1 H ( AA V), = 1 1 vec, = 1 t B X X A (3.) for which a solution is found as z = H t (3.3) Several special cases exist for which (3.3) can be siplified. First, the weights atrices W ay all be equal, or even all equal to one. In these cases H will siplify, as will the pseudo-inverse of H. Another siplification is concerned with the different odels, reflected in restrictions for the space weights. Equation
6 6 (3.3) provides a solution for Z, where A is of full ran. his odel is the generalized Euclidean odel, also nown as IDIOSCAL (Carroll and Chang, 197). he weighted Euclidean odel, or INDSCAL, restricts A to be diagonal, which does siplify H, but not the pseudo-inverse. he identity odel requires A I for all, and does siplify H and its pseudo-inverse, for the ronecer = product vanishes. o avoid coputing the pseudo-inverse of a large atrix, PROXSCAL uses three technical siplifications when appropriate. First, the pseudo-inverse can be replaced by a proper inverse by adding the nullspace, taing the proper inverse and then subtracting the nullspace again as 1 H = H+ N N N= 11 / 1 1. where Furtherore, a diensionwise approach (Heiser and Stoop, 1986) is used which results in a solution for diension a of Z as z = V z, a a a where 1 V V e A A e, a = a a = 1 where e a is the a-th colun of an identity atrix, and 1 za = ( ) a a, B X X A V P A A e = 1 with P a an n p atrix equal to Z, but with the a-th colun containing zeros. Still, the proper inverse of a n n atrix is required. he final siplification is concerned with a aorization function in which the largest eigenvalue of V allows for an easy update (Heiser, 1987; Groenen, Heiser, and Meulan, 1999). Instead of the largest eigenvalue itself, an upper bound is used for this scalar (Wolowicz and Styan, 1980). Update for the Space Weights An update for the space weights ( = 1,..., ) odel is given by A for the generalized Euclidean -1 ( ) A = Z V Z Z B X X. (3.4)
7 7 Suppose, PLQ is the singular value decoposition of A, for which the diagonal atrix with singular values L is in nonincreasing order. hen, for the reduced ran odel, the best r( r < p) ran approxiation of A is given by R, where R contains the first r coluns of PL, and contains the first r coluns of Q. For the weighted Euclidean odel, (3.4) reduces to a diagonal atrix -1 ( ) A = diag Z V Z diag Z B X X. he space weights for the identity odel need no update, since A = I for all. Siplifications can be obtained if all weights W are equal to one and for the reduced ran odel, which can be done in r diensions, as explained in Heiser and Stoop (1986). Restrictions Fixed coordinates If soe of the coordinates of Z are fixed by the user, then only the free coordinates of Z need to be updated. he diensionwise approach is taen one step further, which results in an update for obect i on diension a as p z = e i a i B X X A e = 1 a e A A e V = 1 z e V z % eve eve i a i ia i a a i a ia where the ath colun of Z is divided into za = z% ia +ziaei, with ei the ith colun of 1 the identity atrix, and Va = eaaeav. = 1 his update procedure will only locally iniize (3.1) and repeatedly cycling through all free coordinates until convergence is reached, will provide global optiization. After all free coordinates have been updated, Z is centered on the origin. On output, the configuration is adapted as to coincide with the initial fixed coordinates. Independent variables Independent variables Q are used to express the coordinates of the coon space Z as a weighted su of these independent variables as
8 8 Z QB q b. h = = = 1 An update for Z is found by perforing the following calculations for = 1,..., h: 1.. U h = qb 1 = = 1 = 1 1 C V U A A, where C B( X ) X A = 3. update b as 1 1 = = 1 b q V q A A q 4. optionally, copute optially transfored variables by regressing 1 1 q = b + I V q %, where is greater than or equal to the largest eigenvalue of variable V 1 = = 1 baabv and V, on the original q. Missing eleents in the original variable are replaced with the corresponding values fro Finally, set h = 1 q%. Z QB q b. = = Independent variables restrictions were introduced for the MDS odel in Bentler and Wees (1978), Bloxo (1978), de Leeuw and Heiser (1980) and Meulan and Heiser (1984). If there are ore diensions (p) than independent variables (s), p-s duy variables are created and treated copletely free in the analysis. he transforations for the independent variables fro Step 4 are identical to the transforations of the proxiities, except that the nonnegativety constraint does not apply. After transforation, the variables q are centered on the origin, noralized on n, and the reverse noralization is applied to the regression weights b.
9 9 Step 3: ransforation Update Conditionality wo types of conditionalities exist in PROXSCAL. Conditionality refers to the possible coparison of proxiities in the transforation step. For unconditional transforations, all proxiities are allowed to be copared with each other, irrespective of the source. Matrix-conditional transforations only allow for coparison of proxiities within one atrix, in PROXSCAL refered to as one source. Here, the transforation is coputed for each source seperately (thus ties). ransforation Functions All transforation functions in PROXSCAL result in nonnegative values for the transfored proxiities. After the transforation, the transfored proxiities are noralized and the coon space is optially dilated accordingly. he following transforations are available. Ratio D ˆ =. No transforation is necessary, since the scale of ˆD is adusted in the noralization step. Interval D ˆ = α + β. Both α and β are coputed using linear regression, in such a way that both paraeters are nonnegative. Ordinal D ˆ = WMON, :. Weighted onotone regression (WMON) is coputed using the up-and-down-blocs iniu violators algorith (Krusal, 1964; Barlow et al., 197). For the secondary approach to ties, ties are ept tied, the proxiities within tieblocs are first contracted and expanded afterwards.
10 10 Spline ( D ˆ ) vec = Sb. PROXSCAL uses onotone spline transforations (Rasay, 1988). In this case, the spline transforation gives a sooth nondecreasing piecewise polynoial transforation. It is coputed as a weighted regression of D on the spline basiss. Regression weights b are restricted to be nonnegative and coputed using nonnegative alternating least squares (Groenen, van Os and Meulan, 000). Noralization After transforation, the transfored proxiities are noralized such that the su-of-squares of the weighted transfored proxiities are equal to n(n-1)/ in the unconditional case and equal to n(n-1)/ in the atrix-conditional case. Step 4: erination After evaluation of the loss function, the old function value and new function values are used to decide whether iterations should continue. If the new function value is saller than or equal to the iniu Stress value MINSRESS, provided by the user, iterations are terinated. Also, if the difference in consecutive Stress values is saller than or equal to the convergence criterion DIFFSRESS, provided by the user, iterations are terinated. Finally, iterations are terinated if the current nuber of iterations, exceeds the axiu nuber of iterations MAXIER, also provided by the user. In all other cases, iterations continue. Reaining Issues Acceleration For the identity odel without further restictions, the coon space can be new update old updated with acceleration as Z = Z Z, also refered to as the relaxed update. Lowering diensionality For a restart in p-1 diensions, the p-1 ost iportant diensions need to be identified. For the identity odel, the first p-1 principal axes are used. For the
11 11 weighted Euclidean odel, the p-1 ost iportant space weights are used, and for the generalized Euclidean and reduced ran odels, the p-1 largest singular values of the space weights deterine the reaining diensions. Stress easures he following statistics are used for the coputation of the Stress easures: η η η η ρ ρ κ ( Dˆ ) n = ˆ widi = 1 i< n 4 4 ( Dˆ ) = ˆ widi = 1 i< n i i = 1 i< n 4 4 i i = 1 i< n ( X) = w d ( X ) ( X) = w d ( X ) ( X) = w dˆ d ( X ) i i i = 1 i< n i i i = 1 i< ( X) = w dˆ d ( X ) n ( X) = w d ( X ) d( X) = 1 i< i i where d ( X ) is the average distance. he loss function iniized by PROXSCAL, noralized raw Stress, is given by: σ = ( Dˆ ) + ( X) ( X) η ( Dˆ ) η η α ρ α, with ρ α = η ( X) X. Note that at a local iniu of X, α is equal to one. he other Fit and Stress easures provided by PROXSCAL are given by: η ( Dˆ ) + η ( αx) ρ( αx) η ( Dˆ ) Stress-I:, with α = η αx ρ X.
12 1 Stress-II: ( Dˆ ) + ( X) ( X) η η α ρ α κ ( αx) S-Stress: η 4 ( ˆ ) + η 4 ( α ) ρ ( α ), with ( D) η ˆ α = ρ ρ D X X, with α = η Dispersion Accounted For (DAF): ucer s coefficient of congruence: 1 σ. 1-σ. ( X) 4 X. X. Decoposition of noralized raw Stress Each part of noralized raw Stress, as described before, is assigned to obects and sources. Either su over obects or su over sources are equal to total noralized raw Stress. ransforations on output On output, whenever fixed coordinates or independent variables do not apply, the odels are not unique. In these cases transforations of the coon space and the space weights are in order. For the identity odel, the coon space Z is rotated to principal axes. For the weighted Euclidean odel, = n ( diag ) 1/ diag ( ZZ ) = n Z Z Z Z so that I, and reverse tranforations are applied to the space weights A. Further, the su over sources of the squared space weights are put in descending order as to specify the iportance of the diensions. For the generalized Euclidean odel, the Cholesy decoposition ZZ= LL specifies the coon space on output as 1 Z= nz L, so that ZZ= ni.
13 13 References Barlow, R. E., Bortholoew, D. J., Brener, J. M. and Brun, H. D. (197). Statistical inference under order restrictions. New Yor: Wiley. Bentler, P. M. and Wees, D. G. (1978). Restricted ultidiensional scaling odels. Journal of Matheatical Psychology, 17, Bloxo, B. (1978). Constrained ultidiensional scaling in n spaces. Psychoetria, 43, Carroll, J. D. and Chang, J. J. (197, March). IDIOSCAL (Individual Differences In Orientation SCALing): A generalization of INDSCAL allowing idiosyncratic reference systes as well as analytic approxiation to INDSCAL. Paper presented at the spring eeting of the Psychoetric Society. Princeton, NJ. Coandeur, J. J. F. and Heiser, W. J. (1993). Matheatical derivations in the proxiity scaling (PROXSCAL) of syetric data atrices (ech. Rep. No. RR ). Leiden, he Netherlands: Departent of Data heory, Leiden University. De Leeuw, J. (1977). Applications of convex analysis to ultidiensional scaling. In J. R. Barra, F. Brodeau, G. Roier, and B. van Cutse (Eds.), Recent developents in statistics (pp ). Asterda, he Netherlands: North- Holland. De Leeuw, J. and Heiser, W. J. (1980). Multidiensional scaling with restrictions on the configuration. In P. R. Krishnaiah (Ed.), Multivariate analysis (Vol. V, pp ). Asterda, he Netherlands: North-Holland. Gower, J. C. (1966). Soe distance properties of latent root and vector ethods used in ultivariate analysis. Bioetria, 53, Groenen, P. J. F., Heiser, W. J. and Meulan, J. J. (1999). Global optiization in least squares ultidiensional scaling by distance soothing. Journal of Classification, 16, Groenen, P. J. F., van Os, B. and Meulan, J. J. (000). Optial scaling by alternating length-constained nonnegative least squares, with application to distance-based analysis. Psychoetria, 65,
14 14 Heiser, W. J. (1985). A general MDS initialization procedure using the SMACOF algorith-odel with constraints (ech. Rep. No. RR-85-3). Leiden, he Netherlands: Departent of Data heory, Leiden University. Heiser, W. J. (1987). Joint ordination of species and sites: he unfolding technique. In P. Legendre and L. Legendre (Eds.), Developents in nuerical ecology (pp ). Berlin, Heidelberg: Springer-Verlag. Heiser, W. J. and De Leeuw, J. (1986). SMACOF-I (ech. Rep. No. UG-86-0). Leiden, he Netherlands: Departent of Data heory, Leiden University. Heiser, W. J. and Stoop, I. (1986). Explicit SMACOF algoriths for individual differences scaling (ech. Rep. No. RR-86-14). Leiden, he Netherlands: Departent of Data heory, Leiden University. Krusal, J. B. (1964). Nonetric ultidiensional scaling: A nuerical ethod. Psychoetria, 9, 8-4. Meulan, J. J. and Heiser, W. J. (1984). Constrained Multidiensional Scaling: ore Directions than Diensions. In. Havrane et al (Eds.), COMPSA 1984 (pp ). Wien: Physica Verlag. Rasay, J. O. (1988). Monotone regression splines in action. Statistical Science, 3(4), Stoop, I., Heiser, W. J. and De Leeuw, J. (1981). How to use SMACOF-IA. Leiden, he Netherlands: Departent of Data heory, Leiden University. Stoop, I. and De Leeuw, J. (198). How to use SMACOF-IB. Leiden, he Netherlands: Departent of Data heory, Leiden University. orgerson, W. S. (195). Multidiensional scaling: I. heory and ethod. Psychoetria, 17, orgerson, W. S. (1958). heory and ethods of scaling. New Yor: Wiley. Wolowicz, H. and Styan, G. P. H. (1980). Bounds for eigenvalues using traces. Linear algebra and its applications, 9,
Feature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationAN ITERATIVE MAXIMUM EIGENVALUE METHOD FOR SQUARED DISTANCE SCALING 1. INTRODUCTION
AN ITERATIVE MAXIMUM EIGENVALUE METHOD FOR SQUARED DISTANCE SCALING JAN DE LEEUW ABSTRACT. Meet the abstract. This is the abstract. 1. INTRODUCTION The proble we consider in this paper is a etric ultidiensional
More informationBlock designs and statistics
Bloc designs and statistics Notes for Math 447 May 3, 2011 The ain paraeters of a bloc design are nuber of varieties v, bloc size, nuber of blocs b. A design is built on a set of v eleents. Each eleent
More informationA note on the multiplication of sparse matrices
Cent. Eur. J. Cop. Sci. 41) 2014 1-11 DOI: 10.2478/s13537-014-0201-x Central European Journal of Coputer Science A note on the ultiplication of sparse atrices Research Article Keivan Borna 12, Sohrab Aboozarkhani
More informationON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD
PROCEEDINGS OF THE YEREVAN STATE UNIVERSITY Physical and Matheatical Sciences 04,, p. 7 5 ON THE TWO-LEVEL PRECONDITIONING IN LEAST SQUARES METHOD M a t h e a t i c s Yu. A. HAKOPIAN, R. Z. HOVHANNISYAN
More informationNon-Parametric Non-Line-of-Sight Identification 1
Non-Paraetric Non-Line-of-Sight Identification Sinan Gezici, Hisashi Kobayashi and H. Vincent Poor Departent of Electrical Engineering School of Engineering and Applied Science Princeton University, Princeton,
More informationA Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)
1 A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine (1900 words) Contact: Jerry Farlow Dept of Matheatics Univeristy of Maine Orono, ME 04469 Tel (07) 866-3540 Eail: farlow@ath.uaine.edu
More informationSupport Vector Machine Classification of Uncertain and Imbalanced data using Robust Optimization
Recent Researches in Coputer Science Support Vector Machine Classification of Uncertain and Ibalanced data using Robust Optiization RAGHAV PAT, THEODORE B. TRAFALIS, KASH BARKER School of Industrial Engineering
More informationDetermining OWA Operator Weights by Mean Absolute Deviation Minimization
Deterining OWA Operator Weights by Mean Absolute Deviation Miniization Micha l Majdan 1,2 and W lodziierz Ogryczak 1 1 Institute of Control and Coputation Engineering, Warsaw University of Technology,
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationAn Improved Particle Filter with Applications in Ballistic Target Tracking
Sensors & ransducers Vol. 72 Issue 6 June 204 pp. 96-20 Sensors & ransducers 204 by IFSA Publishing S. L. http://www.sensorsportal.co An Iproved Particle Filter with Applications in Ballistic arget racing
More informationLecture 13 Eigenvalue Problems
Lecture 13 Eigenvalue Probles MIT 18.335J / 6.337J Introduction to Nuerical Methods Per-Olof Persson October 24, 2006 1 The Eigenvalue Decoposition Eigenvalue proble for atrix A: Ax = λx with eigenvalues
More informationMultivariate Methods. Matlab Example. Principal Components Analysis -- PCA
Multivariate Methos Xiaoun Qi Principal Coponents Analysis -- PCA he PCA etho generates a new set of variables, calle principal coponents Each principal coponent is a linear cobination of the original
More informationPage 1 Lab 1 Elementary Matrix and Linear Algebra Spring 2011
Page Lab Eleentary Matri and Linear Algebra Spring 0 Nae Due /03/0 Score /5 Probles through 4 are each worth 4 points.. Go to the Linear Algebra oolkit site ransforing a atri to reduced row echelon for
More informationIntelligent Systems: Reasoning and Recognition. Perceptrons and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley osig 1 Winter Seester 2018 Lesson 6 27 February 2018 Outline Perceptrons and Support Vector achines Notation...2 Linear odels...3 Lines, Planes
More informationChapter 6 1-D Continuous Groups
Chapter 6 1-D Continuous Groups Continuous groups consist of group eleents labelled by one or ore continuous variables, say a 1, a 2,, a r, where each variable has a well- defined range. This chapter explores:
More informationGeneralized AOR Method for Solving System of Linear Equations. Davod Khojasteh Salkuyeh. Department of Mathematics, University of Mohaghegh Ardabili,
Australian Journal of Basic and Applied Sciences, 5(3): 35-358, 20 ISSN 99-878 Generalized AOR Method for Solving Syste of Linear Equations Davod Khojasteh Salkuyeh Departent of Matheatics, University
More informationApproximation in Stochastic Scheduling: The Power of LP-Based Priority Policies
Approxiation in Stochastic Scheduling: The Power of -Based Priority Policies Rolf Möhring, Andreas Schulz, Marc Uetz Setting (A P p stoch, r E( w and (B P p stoch E( w We will assue that the processing
More informationAlgorithms for parallel processor scheduling with distinct due windows and unit-time jobs
BULLETIN OF THE POLISH ACADEMY OF SCIENCES TECHNICAL SCIENCES Vol. 57, No. 3, 2009 Algoriths for parallel processor scheduling with distinct due windows and unit-tie obs A. JANIAK 1, W.A. JANIAK 2, and
More informationA method to determine relative stroke detection efficiencies from multiplicity distributions
A ethod to deterine relative stroke detection eiciencies ro ultiplicity distributions Schulz W. and Cuins K. 2. Austrian Lightning Detection and Inoration Syste (ALDIS), Kahlenberger Str.2A, 90 Vienna,
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationResearch Article Approximate Multidegree Reduction of λ-bézier Curves
Matheatical Probles in Engineering Volue 6 Article ID 87 pages http://dxdoiorg//6/87 Research Article Approxiate Multidegree Reduction of λ-bézier Curves Gang Hu Huanxin Cao and Suxia Zhang Departent of
More informationTopic 5a Introduction to Curve Fitting & Linear Regression
/7/08 Course Instructor Dr. Rayond C. Rup Oice: A 337 Phone: (95) 747 6958 E ail: rcrup@utep.edu opic 5a Introduction to Curve Fitting & Linear Regression EE 4386/530 Coputational ethods in EE Outline
More informationLecture 9 November 23, 2015
CSC244: Discrepancy Theory in Coputer Science Fall 25 Aleksandar Nikolov Lecture 9 Noveber 23, 25 Scribe: Nick Spooner Properties of γ 2 Recall that γ 2 (A) is defined for A R n as follows: γ 2 (A) = in{r(u)
More informationInteractive Markov Models of Evolutionary Algorithms
Cleveland State University EngagedScholarship@CSU Electrical Engineering & Coputer Science Faculty Publications Electrical Engineering & Coputer Science Departent 2015 Interactive Markov Models of Evolutionary
More informationBipartite subgraphs and the smallest eigenvalue
Bipartite subgraphs and the sallest eigenvalue Noga Alon Benny Sudaov Abstract Two results dealing with the relation between the sallest eigenvalue of a graph and its bipartite subgraphs are obtained.
More informationarxiv: v1 [cs.ds] 3 Feb 2014
arxiv:40.043v [cs.ds] 3 Feb 04 A Bound on the Expected Optiality of Rando Feasible Solutions to Cobinatorial Optiization Probles Evan A. Sultani The Johns Hopins University APL evan@sultani.co http://www.sultani.co/
More informationSupplementary Material for Fast and Provable Algorithms for Spectrally Sparse Signal Reconstruction via Low-Rank Hankel Matrix Completion
Suppleentary Material for Fast and Provable Algoriths for Spectrally Sparse Signal Reconstruction via Low-Ran Hanel Matrix Copletion Jian-Feng Cai Tianing Wang Ke Wei March 1, 017 Abstract We establish
More informationDistributed Subgradient Methods for Multi-agent Optimization
1 Distributed Subgradient Methods for Multi-agent Optiization Angelia Nedić and Asuan Ozdaglar October 29, 2007 Abstract We study a distributed coputation odel for optiizing a su of convex objective functions
More informationThe Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Parameters
journal of ultivariate analysis 58, 96106 (1996) article no. 0041 The Distribution of the Covariance Matrix for a Subset of Elliptical Distributions with Extension to Two Kurtosis Paraeters H. S. Steyn
More informationExperimental Design For Model Discrimination And Precise Parameter Estimation In WDS Analysis
City University of New York (CUNY) CUNY Acadeic Works International Conference on Hydroinforatics 8-1-2014 Experiental Design For Model Discriination And Precise Paraeter Estiation In WDS Analysis Giovanna
More informationASSUME a source over an alphabet size m, from which a sequence of n independent samples are drawn. The classical
IEEE TRANSACTIONS ON INFORMATION THEORY Large Alphabet Source Coding using Independent Coponent Analysis Aichai Painsky, Meber, IEEE, Saharon Rosset and Meir Feder, Fellow, IEEE arxiv:67.7v [cs.it] Jul
More informationA new type of lower bound for the largest eigenvalue of a symmetric matrix
Linear Algebra and its Applications 47 7 9 9 www.elsevier.co/locate/laa A new type of lower bound for the largest eigenvalue of a syetric atrix Piet Van Mieghe Delft University of Technology, P.O. Box
More informationA remark on a success rate model for DPA and CPA
A reark on a success rate odel for DPA and CPA A. Wieers, BSI Version 0.5 andreas.wieers@bsi.bund.de Septeber 5, 2018 Abstract The success rate is the ost coon evaluation etric for easuring the perforance
More informationLeast Squares Fitting of Data
Least Squares Fitting of Data David Eberly, Geoetric Tools, Redond WA 98052 https://www.geoetrictools.co/ This work is licensed under the Creative Coons Attribution 4.0 International License. To view a
More informationPolygonal Designs: Existence and Construction
Polygonal Designs: Existence and Construction John Hegean Departent of Matheatics, Stanford University, Stanford, CA 9405 Jeff Langford Departent of Matheatics, Drake University, Des Moines, IA 5011 G
More informationThe Methods of Solution for Constrained Nonlinear Programming
Research Inventy: International Journal Of Engineering And Science Vol.4, Issue 3(March 2014), PP 01-06 Issn (e): 2278-4721, Issn (p):2319-6483, www.researchinventy.co The Methods of Solution for Constrained
More informationKernel Methods and Support Vector Machines
Intelligent Systes: Reasoning and Recognition Jaes L. Crowley ENSIAG 2 / osig 1 Second Seester 2012/2013 Lesson 20 2 ay 2013 Kernel ethods and Support Vector achines Contents Kernel Functions...2 Quadratic
More informationSharp Time Data Tradeoffs for Linear Inverse Problems
Sharp Tie Data Tradeoffs for Linear Inverse Probles Saet Oyak Benjain Recht Mahdi Soltanolkotabi January 016 Abstract In this paper we characterize sharp tie-data tradeoffs for optiization probles used
More informationINNER CONSTRAINTS FOR A 3-D SURVEY NETWORK
eospatial Science INNER CONSRAINS FOR A 3-D SURVEY NEWORK hese notes follow closely the developent of inner constraint equations by Dr Willie an, Departent of Building, School of Design and Environent,
More informationA model reduction approach to numerical inversion for a parabolic partial differential equation
Inverse Probles Inverse Probles 30 (204) 250 (33pp) doi:0.088/0266-56/30/2/250 A odel reduction approach to nuerical inversion for a parabolic partial differential equation Liliana Borcea, Vladiir Drusin
More informationarxiv: v1 [math.na] 10 Oct 2016
GREEDY GAUSS-NEWTON ALGORITHM FOR FINDING SPARSE SOLUTIONS TO NONLINEAR UNDERDETERMINED SYSTEMS OF EQUATIONS MÅRTEN GULLIKSSON AND ANNA OLEYNIK arxiv:6.395v [ath.na] Oct 26 Abstract. We consider the proble
More informationThe Simplex Method is Strongly Polynomial for the Markov Decision Problem with a Fixed Discount Rate
The Siplex Method is Strongly Polynoial for the Markov Decision Proble with a Fixed Discount Rate Yinyu Ye April 20, 2010 Abstract In this note we prove that the classic siplex ethod with the ost-negativereduced-cost
More informationChaotic Coupled Map Lattices
Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each
More informationVariations on Backpropagation
2 Variations on Backpropagation 2 Variations Heuristic Modifications Moentu Variable Learning Rate Standard Nuerical Optiization Conjugate Gradient Newton s Method (Levenberg-Marquardt) 2 2 Perforance
More informationESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS. A Thesis. Presented to. The Faculty of the Department of Mathematics
ESTIMATING AND FORMING CONFIDENCE INTERVALS FOR EXTREMA OF RANDOM POLYNOMIALS A Thesis Presented to The Faculty of the Departent of Matheatics San Jose State University In Partial Fulfillent of the Requireents
More informationA1. Find all ordered pairs (a, b) of positive integers for which 1 a + 1 b = 3
A. Find all ordered pairs a, b) of positive integers for which a + b = 3 08. Answer. The six ordered pairs are 009, 08), 08, 009), 009 337, 674) = 35043, 674), 009 346, 673) = 3584, 673), 674, 009 337)
More information. The univariate situation. It is well-known for a long tie that denoinators of Pade approxiants can be considered as orthogonal polynoials with respe
PROPERTIES OF MULTIVARIATE HOMOGENEOUS ORTHOGONAL POLYNOMIALS Brahi Benouahane y Annie Cuyt? Keywords Abstract It is well-known that the denoinators of Pade approxiants can be considered as orthogonal
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informationMulti-Dimensional Hegselmann-Krause Dynamics
Multi-Diensional Hegselann-Krause Dynaics A. Nedić Industrial and Enterprise Systes Engineering Dept. University of Illinois Urbana, IL 680 angelia@illinois.edu B. Touri Coordinated Science Laboratory
More informationHESSIAN MATRICES OF PENALTY FUNCTIONS FOR SOLVING CONSTRAINED-OPTIMIZATION PROBLEMS
R 702 Philips Res. Repts 24, 322-330, 1969 HESSIAN MATRICES OF PENALTY FUNCTIONS FOR SOLVING CONSTRAINED-OPTIMIZATION PROBLEMS by F. A. LOOTSMA Abstract This paper deals with the Hessian atrices of penalty
More informationData-Driven Imaging in Anisotropic Media
18 th World Conference on Non destructive Testing, 16- April 1, Durban, South Africa Data-Driven Iaging in Anisotropic Media Arno VOLKER 1 and Alan HUNTER 1 TNO Stieltjesweg 1, 6 AD, Delft, The Netherlands
More informationFinding Rightmost Eigenvalues of Large Sparse. Non-symmetric Parameterized Eigenvalue Problems. Abstract. Introduction
Finding Rightost Eigenvalues of Large Sparse Non-syetric Paraeterized Eigenvalue Probles Applied Matheatics and Scientific Coputation Progra Departent of Matheatics University of Maryland, College Par,
More informationOptical Properties of Plasmas of High-Z Elements
Forschungszentru Karlsruhe Techni und Uwelt Wissenschaftlishe Berichte FZK Optical Properties of Plasas of High-Z Eleents V.Tolach 1, G.Miloshevsy 1, H.Würz Project Kernfusion 1 Heat and Mass Transfer
More informationarxiv: v1 [cs.ds] 29 Jan 2012
A parallel approxiation algorith for ixed packing covering seidefinite progras arxiv:1201.6090v1 [cs.ds] 29 Jan 2012 Rahul Jain National U. Singapore January 28, 2012 Abstract Penghui Yao National U. Singapore
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Notes for EE227C (Spring 2018): Convex Optiization and Approxiation Instructor: Moritz Hardt Eail: hardt+ee227c@berkeley.edu Graduate Instructor: Max Sichowitz Eail: sichow+ee227c@berkeley.edu October
More informatione-companion ONLY AVAILABLE IN ELECTRONIC FORM
OPERATIONS RESEARCH doi 10.1287/opre.1070.0427ec pp. ec1 ec5 e-copanion ONLY AVAILABLE IN ELECTRONIC FORM infors 07 INFORMS Electronic Copanion A Learning Approach for Interactive Marketing to a Custoer
More information13.2 Fully Polynomial Randomized Approximation Scheme for Permanent of Random 0-1 Matrices
CS71 Randoness & Coputation Spring 018 Instructor: Alistair Sinclair Lecture 13: February 7 Disclaier: These notes have not been subjected to the usual scrutiny accorded to foral publications. They ay
More informationPriorities in Thurstone Scaling and Steady-State Probabilities in Markov Stochastic Modeling
Journal of Modern Applied Statistical Methods Volue Issue Article 5--03 Priorities in Thurstone Scaling and Steady-State Probabilities in Markov Stochastic Modeling Stan Lipovetsky GfK Custo Research North
More informationEnvelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Damping Characteristics
Copyright c 2007 Tech Science Press CMES, vol.22, no.2, pp.129-149, 2007 Envelope frequency Response Function Analysis of Mechanical Structures with Uncertain Modal Daping Characteristics D. Moens 1, M.
More informationDERIVING PROPER UNIFORM PRIORS FOR REGRESSION COEFFICIENTS
DERIVING PROPER UNIFORM PRIORS FOR REGRESSION COEFFICIENTS N. van Erp and P. van Gelder Structural Hydraulic and Probabilistic Design, TU Delft Delft, The Netherlands Abstract. In probles of odel coparison
More informationANALYTICAL INVESTIGATION AND PARAMETRIC STUDY OF LATERAL IMPACT BEHAVIOR OF PRESSURIZED PIPELINES AND INFLUENCE OF INTERNAL PRESSURE
DRAFT Proceedings of the ASME 014 International Mechanical Engineering Congress & Exposition IMECE014 Noveber 14-0, 014, Montreal, Quebec, Canada IMECE014-36371 ANALYTICAL INVESTIGATION AND PARAMETRIC
More informationpaper prepared for the 1996 PTRC Conference, September 2-6, Brunel University, UK ON THE CALIBRATION OF THE GRAVITY MODEL
paper prepared for the 1996 PTRC Conference, Septeber 2-6, Brunel University, UK ON THE CALIBRATION OF THE GRAVITY MODEL Nanne J. van der Zijpp 1 Transportation and Traffic Engineering Section Delft University
More informationPattern Recognition and Machine Learning. Learning and Evaluation for Pattern Recognition
Pattern Recognition and Machine Learning Jaes L. Crowley ENSIMAG 3 - MMIS Fall Seester 2017 Lesson 1 4 October 2017 Outline Learning and Evaluation for Pattern Recognition Notation...2 1. The Pattern Recognition
More informationW-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS
W-BASED VS LATENT VARIABLES SPATIAL AUTOREGRESSIVE MODELS: EVIDENCE FROM MONTE CARLO SIMULATIONS. Introduction When it coes to applying econoetric odels to analyze georeferenced data, researchers are well
More informationSupplementary Materials: Proofs and Technical Details for Parsimonious Tensor Response Regression Lexin Li and Xin Zhang
Suppleentary Materials: Proofs and Tecnical Details for Parsionious Tensor Response Regression Lexin Li and Xin Zang A Soe preliinary results We will apply te following two results repeatedly. For a positive
More informationKonrad-Zuse-Zentrum für Informationstechnik Berlin Heilbronner Str. 10, D Berlin - Wilmersdorf
Konrad-Zuse-Zentru für Inforationstechnik Berlin Heilbronner Str. 10, D-10711 Berlin - Wilersdorf Folkar A. Borneann On the Convergence of Cascadic Iterations for Elliptic Probles SC 94-8 (Marz 1994) 1
More informationNecessity of low effective dimension
Necessity of low effective diension Art B. Owen Stanford University October 2002, Orig: July 2002 Abstract Practitioners have long noticed that quasi-monte Carlo ethods work very well on functions that
More information(t, m, s)-nets and Maximized Minimum Distance, Part II
(t,, s)-nets and Maxiized Miniu Distance, Part II Leonhard Grünschloß and Alexander Keller Abstract The quality paraeter t of (t,, s)-nets controls extensive stratification properties of the generated
More informationHomework 3 Solutions CSE 101 Summer 2017
Hoework 3 Solutions CSE 0 Suer 207. Scheduling algoriths The following n = 2 jobs with given processing ties have to be scheduled on = 3 parallel and identical processors with the objective of iniizing
More informationUsing EM To Estimate A Probablity Density With A Mixture Of Gaussians
Using EM To Estiate A Probablity Density With A Mixture Of Gaussians Aaron A. D Souza adsouza@usc.edu Introduction The proble we are trying to address in this note is siple. Given a set of data points
More informationACTIVE VIBRATION CONTROL FOR STRUCTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAKE EXCITATION
International onference on Earthquae Engineering and Disaster itigation, Jaarta, April 14-15, 8 ATIVE VIBRATION ONTROL FOR TRUTURE HAVING NON- LINEAR BEHAVIOR UNDER EARTHQUAE EXITATION Herlien D. etio
More informationTesting Properties of Collections of Distributions
Testing Properties of Collections of Distributions Reut Levi Dana Ron Ronitt Rubinfeld April 9, 0 Abstract We propose a fraework for studying property testing of collections of distributions, where the
More informationCharacter analysis on linear elementary algebra with max-plus operation
Available online at www.worldscientificnews.co WSN 100 (2018) 110-123 EISSN 2392-2192 Character analysis on linear eleentary algebra with ax-plus operation ABSTRACT Kalfin 1, Jufra 2, Nora Muhtar 2, Subiyanto
More informationA Markov Framework for the Simple Genetic Algorithm
A arkov Fraework for the Siple Genetic Algorith Thoas E. Davis*, Jose C. Principe Electrical Engineering Departent University of Florida, Gainesville, FL 326 *WL/NGS Eglin AFB, FL32542 Abstract This paper
More informationRECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS MEMBRANE
Proceedings of ICIPE rd International Conference on Inverse Probles in Engineering: Theory and Practice June -8, 999, Port Ludlow, Washington, USA : RECOVERY OF A DENSITY FROM THE EIGENVALUES OF A NONHOMOGENEOUS
More informationProbability Distributions
Probability Distributions In Chapter, we ephasized the central role played by probability theory in the solution of pattern recognition probles. We turn now to an exploration of soe particular exaples
More informationREDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION
ISSN 139 14X INFORMATION TECHNOLOGY AND CONTROL, 008, Vol.37, No.3 REDUCTION OF FINITE ELEMENT MODELS BY PARAMETER IDENTIFICATION Riantas Barauskas, Vidantas Riavičius Departent of Syste Analysis, Kaunas
More informationComplex Quadratic Optimization and Semidefinite Programming
Coplex Quadratic Optiization and Seidefinite Prograing Shuzhong Zhang Yongwei Huang August 4 Abstract In this paper we study the approxiation algoriths for a class of discrete quadratic optiization probles
More informationOPTIMIZATION OF SPECIFIC FACTORS TO PRODUCE SPECIAL ALLOYS
5 th International Conference Coputational Mechanics and Virtual Engineering COMEC 2013 24-25 October 2013, Braşov, Roania OPTIMIZATION OF SPECIFIC FACTORS TO PRODUCE SPECIAL ALLOYS I. Milosan 1 1 Transilvania
More informationi ij j ( ) sin cos x y z x x x interchangeably.)
Tensor Operators Michael Fowler,2/3/12 Introduction: Cartesian Vectors and Tensors Physics is full of vectors: x, L, S and so on Classically, a (three-diensional) vector is defined by its properties under
More informationSpine Fin Efficiency A Three Sided Pyramidal Fin of Equilateral Triangular Cross-Sectional Area
Proceedings of the 006 WSEAS/IASME International Conference on Heat and Mass Transfer, Miai, Florida, USA, January 18-0, 006 (pp13-18) Spine Fin Efficiency A Three Sided Pyraidal Fin of Equilateral Triangular
More informationOBJECTIVES INTRODUCTION
M7 Chapter 3 Section 1 OBJECTIVES Suarize data using easures of central tendency, such as the ean, edian, ode, and idrange. Describe data using the easures of variation, such as the range, variance, and
More informationPhysics 215 Winter The Density Matrix
Physics 215 Winter 2018 The Density Matrix The quantu space of states is a Hilbert space H. Any state vector ψ H is a pure state. Since any linear cobination of eleents of H are also an eleent of H, it
More informationQualitative Modelling of Time Series Using Self-Organizing Maps: Application to Animal Science
Proceedings of the 6th WSEAS International Conference on Applied Coputer Science, Tenerife, Canary Islands, Spain, Deceber 16-18, 2006 183 Qualitative Modelling of Tie Series Using Self-Organizing Maps:
More informationRESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS
BIT Nuerical Matheatics 43: 459 466, 2003. 2003 Kluwer Acadeic Publishers. Printed in The Netherlands 459 RESTARTED FULL ORTHOGONALIZATION METHOD FOR SHIFTED LINEAR SYSTEMS V. SIMONCINI Dipartiento di
More informationOptimal Resource Allocation in Multicast Device-to-Device Communications Underlaying LTE Networks
1 Optial Resource Allocation in Multicast Device-to-Device Counications Underlaying LTE Networks Hadi Meshgi 1, Dongei Zhao 1 and Rong Zheng 2 1 Departent of Electrical and Coputer Engineering, McMaster
More informationGenetic Algorithm Search for Stent Design Improvements
Genetic Algorith Search for Stent Design Iproveents K. Tesch, M.A. Atherton & M.W. Collins, South Bank University, London, UK Abstract This paper presents an optiisation process for finding iproved stent
More informationA note on the realignment criterion
A note on the realignent criterion Chi-Kwong Li 1, Yiu-Tung Poon and Nung-Sing Sze 3 1 Departent of Matheatics, College of Willia & Mary, Williasburg, VA 3185, USA Departent of Matheatics, Iowa State University,
More informationConsistent Multiclass Algorithms for Complex Performance Measures. Supplementary Material
Consistent Multiclass Algoriths for Coplex Perforance Measures Suppleentary Material Notations. Let λ be the base easure over n given by the unifor rando variable (say U over n. Hence, for all easurable
More informationA Simplified Analytical Approach for Efficiency Evaluation of the Weaving Machines with Automatic Filling Repair
Proceedings of the 6th SEAS International Conference on Siulation, Modelling and Optiization, Lisbon, Portugal, Septeber -4, 006 0 A Siplified Analytical Approach for Efficiency Evaluation of the eaving
More informationGEE ESTIMATORS IN MIXTURE MODEL WITH VARYING CONCENTRATIONS
ACTA UIVERSITATIS LODZIESIS FOLIA OECOOMICA 3(3142015 http://dx.doi.org/10.18778/0208-6018.314.03 Olesii Doronin *, Rostislav Maiboroda ** GEE ESTIMATORS I MIXTURE MODEL WITH VARYIG COCETRATIOS Abstract.
More informationCourse Notes for EE227C (Spring 2018): Convex Optimization and Approximation
Course Notes for EE7C (Spring 018: Convex Optiization and Approxiation Instructor: Moritz Hardt Eail: hardt+ee7c@berkeley.edu Graduate Instructor: Max Sichowitz Eail: sichow+ee7c@berkeley.edu October 15,
More informationFast Montgomery-like Square Root Computation over GF(2 m ) for All Trinomials
Fast Montgoery-like Square Root Coputation over GF( ) for All Trinoials Yin Li a, Yu Zhang a, a Departent of Coputer Science and Technology, Xinyang Noral University, Henan, P.R.China Abstract This letter
More informationStochastic Subgradient Methods
Stochastic Subgradient Methods Lingjie Weng Yutian Chen Bren School of Inforation and Coputer Science University of California, Irvine {wengl, yutianc}@ics.uci.edu Abstract Stochastic subgradient ethods
More informationOptimal Jamming Over Additive Noise: Vector Source-Channel Case
Fifty-first Annual Allerton Conference Allerton House, UIUC, Illinois, USA October 2-3, 2013 Optial Jaing Over Additive Noise: Vector Source-Channel Case Erah Akyol and Kenneth Rose Abstract This paper
More information3.3 Variational Characterization of Singular Values
3.3. Variational Characterization of Singular Values 61 3.3 Variational Characterization of Singular Values Since the singular values are square roots of the eigenvalues of the Heritian atrices A A and
More informationDetection and Estimation Theory
ESE 54 Detection and Estiation Theory Joseph A. O Sullivan Sauel C. Sachs Professor Electronic Systes and Signals Research Laboratory Electrical and Systes Engineering Washington University 11 Urbauer
More informationRecursive Algebraic Frisch Scheme: a Particle-Based Approach
Recursive Algebraic Frisch Schee: a Particle-Based Approach Stefano Massaroli Renato Myagusuku Federico Califano Claudio Melchiorri Atsushi Yaashita Hajie Asaa Departent of Precision Engineering, The University
More information