Package lintools. July 30, 2018
|
|
- Avis Gordon
- 5 years ago
- Views:
Transcription
1 Package lintools July 30, 2018 Maintainer Mark van der Loo License GPL-3 Title Manipulation of Linear Systems of (in)equalities LazyData no Type Package LazyLoad yes Variale elimination (Gaussian elimination, Fourier-Motzkin elimination), Moore-Penrose pseudoinverse, reduction to reduced row echelon form, value sustitution, projecting a vector on the convex polytope descried y a system of (in)equations, simplify systems y removing spurious columns and rows and collapse implied equalities, test if a matrix is totally unimodular, compute variale ranges implied y linear (in)equalities. Version URL BugReports Imports utils Suggests testthat, knitr VignetteBuilder knitr RoxygenNote NeedsCompilation yes uthor Mark van der Loo [aut, cre], Edwin de Jonge [aut] Repository CRN Date/Pulication :10:05 UTC R topics documented: lock_index compact
2 2 lock_index echelon eliminate is_feasile is_totally_unimodular lintools normalize pinv project ranges sparse_constraints sparse_project sust_value Index 18 lock_index Find independent locks of equations. Find independent locks of equations. lock_index(, eps = 1e-08) eps [numeric] Matrix [numeric] Coefficients with asolute value < eps are treated as zero. list containing numeric vectors, each vector indexing an independent lock of rows in the system x <=. <- matrix(c( 1,0,2,0,0, 3,0,4,0,0, 0,5,0,6,7, 0,8,0,0,9 ),yrow=true,nrow=4) <- rep(0,4) i <- lock_index() lapply(i,function(ii) compact([ii,,drop=flse],=[ii])$)
3 compact 3 compact Remove spurious variales and restrictions system of linear (in)equations can e compactified y removing zero-rows and zero-columns (=variales). Such rows and columns may arise after sustitution (see sust_value) or eliminaton of a variale (see eliminate). compact(,, x = NULL, neq = nrow(), nleq = 0, eps = 1e-08, remove_columns = TRUE, remove_rows = TRUE, deduplicate = TRUE, implied_equations = TRUE) x neq nleq eps [numeric] matrix [numeric] vector [numeric] vector [numeric] The first neq rows in and are treated as linear equalities. [numeric] The nleq rows after neq are treated as inequations of the form a.x<=. ll remaining rows are treated as strict inequations of the form a.x<. [numeric] nything with asolute value < eps is considered zero. remove_columns [logical] Toggle remove spurious columns from and variales from x remove_rows [logical] Toggle remove spurious rows deduplicate [logical] Toggle remove duplicate rows implied_equations [logical] replace cases of a.x<= and a.x>= with a.x==. list with the following elements. : The compactified version of input : The compactified version of input x: The compactified version of input x neq: numer of equations in new system nleq: numer of inequations of the form a.x<= in the new system cols_removed: [logical] indicates what elements of x (columns of ) have een removed It is assumend that the system of equations is in normalized form (see link{normalize}).
4 4 echelon echelon Reduced row echelon form Transform the equalities in a system of linear (in)equations or Reduced Row Echelon form (RRE) echelon(,, neq = nrow(), nleq = 0, eps = 1e-08) neq nleq eps [numeric] matrix [numeric] vector [numeric] The first neq rows of, are treated as equations. [numeric] The nleq rows after neq are treated as inequations of the form a.x<=. ll remaining rows are treated as strict inequations of the form a.x<. [numeric] s of magnitude less than eps are considered zero (for the purpose of handling machine rounding errors). list with the following components: : the matrix with equalities transformed to RRE form. : the constant vector corresponding to neq: the numer of equalities in the resulting system. nleq: the numer of inequalities of the form a.x <=. This will only e passed to the output. The parameters, and neq descrie a system of the form x<=, where the first neq rows are equalities. The equalities are transformed to RRE form. system of equations is in reduced row echelon form when ll zero rows are elow the nonzero rows For every row, the leading coefficient (first nonzero from the left) is always right of the leading coefficient of the row aove it. The leading coefficient equals 1, and is the only nonzero coefficient in its column.
5 eliminate 5 echelon( = matrix(c( 1,3,1, 2,7,3, 1,5,3, 1,2,0), yrow=true, nrow=4), = c(4,-9,1,8), neq=4 ) eliminate Eliminate a variale from a set of edit rules Eliminating a variale amounts to deriving all (non-redundant) linear (in)equations not containing that variale. Geometrically, it can e interpreted as a projection of the solution space (vectors satisfying all equations) along the eliminated variale s axis. eliminate(,, neq = nrow(), nleq = 0, variale, H = NULL, h = 0, eps = 1e-08) neq nleq variale H h eps [numeric] Matrix [numeric] vector [numeric] The first neq rows in and are treated as linear equalities. [numeric] The nleq rows after neq are treated as inequations of the form a.x<=. ll remaining rows are treated as strict inequations of the form a.x<. [numeric logical character] Index in columns of, representing the variale to eliminate. [numeric] (optional) Matrix indicating how linear inequalities have een derived. [numeric] (optional) numer indicating how many variales have een eliminated from the original system using Fourier-Motzkin elimination. [numeric] Coefficients with asolute value <= eps are treated as zero.
6 6 is_feasile list with the folowing components : the corresponding to the system with variales eliminated. : the constant vector corresponding to the resulting system neq: the numer of equations H: The memory matrix storing how each row was derived h: The numer of variales eliminated from the original system. For equalities Gaussian elimination is applied. If inequalities are involved, Fourier-Motzkin elimination is used. In principle, FM-elimination can generate a large numer of redundant inequations, especially when applied recursively. Redundancies can e recognized y recording how new inequations have een derived from the original set. This is stored in the H matrix when multiple variales are to e eliminated (Kohler, 1967). References D.. Kohler (1967) Projections of convex polyhedral sets, Operational Research Center Report, ORC 67-29, University of California, Berkely. H.P. Williams (1986) Fourier s method of linear programming and its dual. merican Mathematical Monthly 93, pp # Example from Williams (1986) <- matrix(c( 4, -5, -3, 1, -1, 1, -1, 0, 1, 1, 2, 0, -1, 0, 0, 0, 0, -1, 0, 0, 0, 0, -1, 0),yrow=TRUE,nrow=6) <- c(0,2,3,0,0,0) L <- eliminate(=, =, neq=0, nleq=6, variale=1) is_feasile Check feasiility of a system of linear (in)equations Check feasiility of a system of linear (in)equations
7 is_totally_unimodular 7 is_feasile(,, neq = nrow(), nleq = 0, eps = 1e-08, method = "elimination") neq nleq eps method [numeric] matrix [numeric] vector [numeric] The first neq rows in and are treated as linear equalities. [numeric] The nleq rows after neq are treated as inequations of the form a.x<=. ll remaining rows are treated as strict inequations of the form a.x<. [numeric] solute values < eps are treated as zero. [character] t the moment, only the elimination method is implemented. # n infeasile system: # x + y == 0 # x > 0 # y > 0 <- matrix(c(1,1,1,0,0,1),yrow=true,nrow=3) <- rep(0,3) is_feasile(=,=,neq=1,nleq=0) # feasile system: # x + y == 0 # x >= 0 # y >= 0 <- matrix(c(1,1,1,0,0,1),yrow=true,nrow=3) <- rep(0,3) is_feasile(=,=,neq=1,nleq=2) is_totally_unimodular Test for total unimodularity of a matrix. Check wether a matrix is totally unimodular. is_totally_unimodular() n oject of class matrix.
8 8 is_totally_unimodular matrix for which the determinant of every square sumatrix equals 1, 0 or 1 is called totally unimodular. This function tests wether a matrix with coefficients in { 1, 0, 1} is totally unimodular. It tries to reduce the matrix using the reduction method descried in Scholtus (2008). Next, a test ased on Heller and Tompkins (1956) or Raghavachari is performed. logical References Heller I and Tompkins CB (1956). n extension of a theorem of Danttzig s In kuhn HW and Tucker W (eds.), pp Princeton University Press. Raghavachari M (1976). constructive method to recognize the total unimodularity of a matrix. _Zeitschrift fur operations research_, *20*, pp Scholtus S (2008). lgorithms for correcting some ovious inconsistencies and rounding errors in usiness survey data. Technical Report 08015, Netherlands. # Totally unimodular matrix, reduces to nothing <- matrix(c( 1,1,0,0,0, -1,0,0,1,0, 0,0,01,1,0, 0,0,0,-1,1),nrow=5) is_totally_unimodular() # Totally unimodular matrix, y Heller-Tompson criterium <- matrix(c( 1,1,0,0, 0,0,1,1, 1,0,1,0, 0,1,0,1),nrow=4) is_totally_unimodular() # Totally unimodular matrix, y Raghavachani recursive criterium <- matrix(c( 1,1,1, 1,1,0, 1,0,1)) is_totally_unimodular()
9 lintools 9 lintools Tools for manipulating linear systems of (in)equations Tools for manipulating linear systems of (in)equations This package offers a asic and consistent interface to a numer of operations on linear systems of (in)equations not availale in ase R. Except for the projection on the convex polytope, operations are currently supported for dense matrices only. The following operations are implemented. Split matrices in independent locks Remove spurious rows and columns from a system of (in)equations Rewrite equalities in reduced row echelon form Eliminate variales through Gaussian or Fourier-Motzkin elimination Determine the feasiility of a system of linear (in)equations Compute Moore-Penrose Pseudoinverse Project a vector onto the convec polytope descried y a set of linear (in)equations Simplify a system y sustituting values Most functions assume a system of (in)equations to e stored in a standard form. The normalize function can ring any system of equations to this form. normalize Bring a system of (in)equalities in a standard form Bring a system of (in)equalities in a standard form normalize(,, operators, unit = 0) [numeric] Matrix [numeric] vector operators [character] operators in {<,<=,==,>=,>}. unit [numeric] (nonnegative) Your unit of measurement. This is used to replace strict inequations of the form a.x < with a.x <= -unit. Typically, unit is related to the units in which your data is measured. If unit is 0, inequations are not replaced.
10 10 pinv list with the folowing components : the corresponding to the normalized sytem. : the constant vector corresponding to the normalized system neq: the numer of equations nleq: the numer of non-strict inequations (<=) order: the index vector used to permute the original rows of. For this package, a set of equations is in normal form when The first neq rows represent linear equalities. The next nleq rows represent inequalities of the form a.x <= ll other rows are strict inequalities of the form a.x < If unit>0, the strict inequalities a.x < are replaced with inequations of the form a.x <= -unit, where unit represents the precision of measurement. <- matrix(1:12,nrow=4) <- 1:4 ops <- c("<=","==","==","<") normalize(,,ops) normalize(,,ops,unit=0.1) pinv Moore-Penrose pseudoinverse Compute the pseudoinverse of a matrix using the SVD-construction pinv(, eps = 1e-08) eps [numeric] matrix [numeric] tolerance for determining zero singular values
11 project 11 The Moore-Penrose pseudoinverse (sometimes called the generalized inverse) + of a matrix has the property that + + =. It can e constructed as follows. Compute the singular value decomposition = UDV T Replace diagonal elements in D of which the asolute values are larger than some limit eps with their reciprocal values Compute + = UDV T References S Lipshutz and M Lipson (2009) Linear lgera. In: Schuam s outlines. McGraw-Hill <- matrix(c( 1, 1, -1, 2, 2, 2, -1, 3, -1, -1, 2, -3 ),yrow=true,nrow=3) # multiply y 55 to retrieve whole numers pinv() * 55 project Project a vector on the order of the region defined y a set of linear (in)equality restrictions. Compute a vector, closest to x in the Euclidean sense, satisfying a set of linear (in)equality restrictions. project(x,,, neq = length(), w = rep(1, length(x)), eps = 0.01, maxiter = 1000L) x neq w eps maxiter [numeric] Vector that needs to satisfy the linear restrictions. [matrix] Coefficient matrix for linear restrictions. [numeric] Right hand side of linear restrictions. [numeric] The first neq rows in and are treated as linear equalities. The others as Linear inequalities of the form x <=. [numeric] Optional weight vector of the same length as x. Must e positive. The maximum allowed deviation from the constraints (see details). maximum numer of iterations
12 12 project list with the following entries: x: the adjusted vector status: Exit status: 0: success 1: could not allocate enough memory (space for approximately 2(m + n) doules is necessary). 2: divergence detected (set of restrictions may e contradictory) 3: maximum numer of iterations reached eps: The tolerance achieved after optimizing (see ). iterations: The numer of iterations performed. duration: the time it took to compute the adjusted vector ojective: The (weighted) Euclidean distance etween the initial and the adjusted vector The tolerance eps is defined as the maximum asolute value of the difference vector x for equalities. For inequalities, the difference vector is set to zero when it s value is lesser than zero (i.e. when the restriction is satisfied). The algorithm iterates until either the tolerance is met, the numer of allowed iterations is exceeded or divergence is detected. See lso sparse_project # the system # x + y = 10 # -x <= 0 # ==> x > 0 # -y <= 0 # ==> y > 0 # <- matrix(c( 1,1, -1,0, 0,-1), yrow=true, nrow=3 ) <- c(10,0,0) # x and y will e adjusted y the same amount project(x=c(4,5), =, =, neq=1) # One of the inequalies violated project(x=c(-1,5), =, =, neq=1) # Weighted distances: 'heavy' variales change less
13 ranges 13 project(x=c(4,5), =, =, neq=1, w=c(100,1)) # if w=1/x0, the ratio etween coefficients of x0 stay the same (to first order) x0 <- c(x=4,y=5) x1 <- project(x=x0,=,=,neq=1,w=1/x0) x0[1]/x0[2] x1$x[1] / x1$x[2] ranges Derive variale ranges from linear restrictions Gaussian and/or Fourier-Motzkin elimination is used to derive upper and lower limits implied y a system of (in)equations. ranges(,, neq = nrow(), nleq = 0, eps = 1e-08) neq nleq eps [numeric] Matrix [numeric] vector [numeric] The first neq rows in and are treated as linear equalities. [numeric] The nleq rows after neq are treated as inequations of the form a.x<=. ll remaining rows are treated as strict inequations of the form a.x<. [numeric] Coefficients with asolute value <= eps are treated as zero. using Fourier-Motzkin elimination. sparse_constraints Generate sparse set of constraints. Generate a constraint set to e used y sparse_project Read sparse constraints from a data.frame Print sparse_constraints oject
14 14 sparse_constraints sparse_constraints(oject,...) sparseconstraints(oject,...) ## S3 method for class 'data.frame' sparse_constraints(oject,, neq = length(), ase = 1L, sorted = FLSE,...) ## S3 method for class 'sparse_constraints' print(x, range = 1L:10L,...) oject R oject to e translated to sparse_constraints format.... options to e passed to other methods Constant vector neq The first new equations are interpreted as equality constraints, the rest as <= ase are the indices in oject[,1:2] ase 0 or ase 1? sorted x Note range is oject sorted y the first column? an oject of class sparse_constraints integer vector stating which constraints to print Oject of class sparse_constraints (see details). s of version , sparseconstraints is deprecated. Use sparse_constraints instead. The sparse_constraints ojects holds coefficients of and of the system x in sparse format, outside of R s memory. It can e reused to find solutions for vectors to adjust. In R, it is a reference oject. In particular, it is meaningless to Copy the oject. You only will only generate a pointer to physically the same oject. Save the oject. The physical oject is destroyed when R closes, or when R s garage collector cleans up a removed sparse_constraints oject. The $project method Once a sparse_constraints oject sc is created, you can reuse it to optimize several vectors y calling sc$project() with the following parameters: x: [numeric] the vector to e optimized
15 sparse_project 15 w: [numeric] the weight vector (of length(x)). By default all weights equal 1. eps: [numeric] desired tolerance. By default 10 2 maxiter: [integer] maximum numer of iterations. By default The return value of $spa is the same as that of sparse_project. See lso sparse_project, project # The following system of constraints, stored in # row-column-coefficient format # # x1 + x8 == 950, # x3 + x4 == 950, # x6 + x7 == x8, # x4 > 0 # <- data.frame( row = c( 1, 1, 2, 2, 3, 3, 3, 4), col = c( 1, 2, 3, 4, 2, 5, 6, 4), coef = c(-1,-1,-1,-1, 1,-1,-1,-1) ) <- c(-950, -950, 0,0) sc <- sparse_constraints(,, neq=3) # djust the 0-vector minimally so all constraints are met: sc$project(x=rep(0,8)) # Use the same oject to adjust the 100*1-vector sc$project(x=rep(100,8)) # use the same oject to adjust the 0-vector, ut with different weights sc$project(x=rep(0,8),w=1:8) sparse_project Successive projections with sparsely defined restrictions Compute a vector, closest to x satisfying a set of linear (in)equality restrictions.
16 16 sparse_project sparse_project(x,,, neq = length(), w = rep(1, length(x)), eps = 0.01, maxiter = 1000L,...) x [numeric] Vector to optimize, starting point. [data.frame] Coeffiencient matrix in [row,column,coefficient] format. [numeric] Constant vector of the system x neq [integer] Numer of equalities w [numeric] weight vector of same length of x eps maximally allowed tolerance maxiter maximally allowed numer of iterations.... extra parameters passed to sparse_constraints list with the following entries: x: the adjusted vector status: Exit status: 0: success 1: could not allocate enough memory (space for approximately 2(m + n) doules is necessary). 2: divergence detected (set of restrictions may e contradictory) 3: maximum numer of iterations reached eps: The tolerance achieved after optimizing (see ). iterations: The numer of iterations performed. duration: the time it took to compute the adjusted vector ojective: The (weighted) Euclidean distance etween the initial and the adjusted vector The tolerance eps is defined as the maximum asolute value of the difference vector x for equalities. For inequalities, the difference vector is set to zero when it s value is lesser than zero (i.e. when the restriction is satisfied). The algorithm iterates until either the tolerance is met, the numer of allowed iterations is exceeded or divergence is detected. See lso project, sparse_constraints
17 sust_value 17 # the system # x + y = 10 # -x <= 0 # ==> x > 0 # -y <= 0 # ==> y > 0 # Defined in the row-column-coefficient form: <- data.frame( row = c(1,1,2,3), col = c(1,2,1,2), coef= c(1,1,-1,-1) ) <- c(10,0,0) sparse_project(x=c(4,5),=,=) sust_value Sustitute a value in a system of linear (in)equations Sustitute a value in a system of linear (in)equations sust_value(,, variales, values, remove_columns = FLSE) [numeric] matrix [numeric] vector variales [numeric logical character] vector of column indices in values [numeric] vecor of values to sustitute. remove_columns [logical] Remove spurious columns when sustituting? list with the following components: : the corresponding to the simplified sytem. : the constant vector corresponding to the new system system of the form x <= can e simplified if one or more of the x[i] values is fixed.
18 Index lock_index, 2 compact, 3 echelon, 4 eliminate, 3, 5 is_feasile, 6 is_totally_unimodular, 7 lintools, 9 lintools-package (lintools), 9 matrix, 7 normalize, 9, 9 pinv, 10 print.sparse_constraints (sparse_constraints), 13 project, 11, 15, 16 ranges, 13 sparse_constraints, 13, 16 sparse_project, 12, 13, 15, 15 sparseconstraints (sparse_constraints), 13 sust_value, 3, 17 18
Package deductive. June 2, 2017
Package deductive June 2, 2017 Maintainer Mark van der Loo License GPL-3 Title Data Correction and Imputation Using Deductive Methods LazyData no Type Package LazyLoad yes Attempt
More informationPackage ECctmc. May 1, 2018
Type Package Package ECctmc May 1, 2018 Title Simulation from Endpoint-Conditioned Continuous Time Markov Chains Version 0.2.5 Date 2018-04-30 URL https://github.com/fintzij/ecctmc BugReports https://github.com/fintzij/ecctmc/issues
More informationPackage lumberjack. R topics documented: July 20, 2018
Package lumberjack July 20, 2018 Maintainer Mark van der Loo License GPL-3 Title Track Changes in Data LazyData no Type Package LazyLoad yes A function composition ('pipe') operator
More informationLegal and impossible dependences
Transformations and Dependences 1 operations, column Fourier-Motzkin elimination us use these tools to determine (i) legality of permutation and Let generation of transformed code. (ii) Recall: Polyhedral
More informationPackage widyr. August 14, 2017
Type Package Title Widen, Process, then Re-Tidy Data Version 0.1.0 Package widyr August 14, 2017 Encapsulates the pattern of untidying data into a wide matrix, performing some processing, then turning
More informationPackage ECOSolveR. February 18, 2018
Package ECOSolveR Type Package Title Embedded Conic Solver in R Version 0.4 Date 2018-02-16 VignetteBuilder knitr SystemRequirements GNU make February 18, 2018 URL https://github.com/bnaras/ecosolver BugReports
More informationDM545 Linear and Integer Programming. Lecture 2. The Simplex Method. Marco Chiarandini
DM545 Linear and Integer Programming Lecture 2 The Marco Chiarandini Department of Mathematics & Computer Science University of Southern Denmark Outline 1. 2. 3. 4. Standard Form Basic Feasible Solutions
More informationx = 12 x = 12 1x = 16
2.2 - The Inverse of a Matrix We've seen how to add matrices, multiply them by scalars, subtract them, and multiply one matrix by another. The question naturally arises: Can we divide one matrix by another?
More informationAutomatic editing of numerical data with R
Automatic editing of numerical data with R Mark van der Loo, Edwin De Jonge, and Sander Scholtus Statistics Netherlands, dpt. of Methodology and Quality, PO box 24500, 2490 HA The Hague, The Netherlands
More informationRecognition, SVD, and PCA
Recognition, SVD, and PCA Recognition Suppose you want to find a face in an image One possibility: look for something that looks sort of like a face (oval, dark band near top, dark band near bottom) Another
More informationPackage mixsqp. November 14, 2018
Encoding UTF-8 Type Package Version 0.1-79 Date 2018-11-05 Package mixsqp November 14, 2018 Title Sequential Quadratic Programming for Fast Maximum-Likelihood Estimation of Mixture Proportions URL https://github.com/stephenslab/mixsqp
More informationTheory of Integers. CS389L: Automated Logical Reasoning. Lecture 13: The Omega Test. Overview of Techniques. Geometric Description
Theory of ntegers This lecture: Decision procedure for qff theory of integers CS389L: Automated Logical Reasoning Lecture 13: The Omega Test şıl Dillig As in previous two lectures, we ll consider T Z formulas
More informationMATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix.
MATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix. Row echelon form A matrix is said to be in the row echelon form if the leading entries shift to the
More informationPackage fastdummies. January 8, 2018
Type Package Package fastdummies January 8, 2018 Title Fast Creation of Dummy (Binary) Columns and Rows from Categorical Variables Version 1.0.0 Description Creates dummy columns from columns that have
More informationEXTENSION. a 1 b 1 c 1 d 1. Rows l a 2 b 2 c 2 d 2. a 3 x b 3 y c 3 z d 3. This system can be written in an abbreviated form as
EXTENSION Using Matrix Row Operations to Solve Systems The elimination method used to solve systems introduced in the previous section can be streamlined into a systematic method by using matrices (singular:
More informationIntroduction to Machine Learning Spring 2018 Note Sparsity and LASSO. 1.1 Sparsity for SVMs
CS 189 Introduction to Machine Learning Spring 2018 Note 21 1 Sparsity and LASSO 1.1 Sparsity for SVMs Recall the oective function of the soft-margin SVM prolem: w,ξ 1 2 w 2 + C Note that if a point x
More informationArray Dependence Analysis as Integer Constraints. Array Dependence Analysis Example. Array Dependence Analysis as Integer Constraints, cont
Theory of Integers CS389L: Automated Logical Reasoning Omega Test Işıl Dillig Earlier, we talked aout the theory of integers T Z Signature of T Z : Σ Z : {..., 2, 1, 0, 1, 2,..., 3, 2, 2, 3,..., +,, =,
More informationSolving Systems Using Row Operations 1 Name
The three usual methods of solving a system of equations are graphing, elimination, and substitution. While these methods are excellent, they can be difficult to use when dealing with three or more variables.
More informationNumerical Linear Algebra
Numerical Linear Algebra Probably the simplest kind of problem. Occurs in many contexts, often as part of larger problem. Symbolic manipulation packages can do linear algebra "analytically" (e.g. Mathematica,
More informationPackage svd. R topics documented: September 26, 2017
Package svd September 26, 2017 Type Package Imports methods Suggests testthat (>= 0.7) Title Interfaces to Various State-of-Art SVD and Eigensolvers Version 0.4.1 Author Anton Korobeynikov [aut, cre],
More informationPackage projector. February 27, 2018
Package projector February 27, 2018 Title Project Dense Vectors Representation of Texts on a 2D Plan Version 0.0.2 Date 2018-02-27 Maintainer Michaël Benesty Display dense vector representation
More informationPackage quadprogxt. February 4, 2018
Package quadprogxt February 4, 2018 Title Quadratic Programming with Absolute Value Constraints Version 0.0.2 Description Extends the quadprog package to solve quadratic programs with absolute value constraints
More informationPackage filematrix. R topics documented: February 27, Type Package
Type Package Package filematrix February 27, 2018 Title File-Backed Matrix Class with Convenient Read and Write Access Version 1.3 Date 2018-02-26 Description Interface for working with large matrices
More informationTHEORY OF LINEAR AND INTEGER PROGRAMMING
THEORY OF LINEAR AND INTEGER PROGRAMMING ALEXANDER SCHRIJVER Centrum voor Wiskunde en Informatica, Amsterdam A Wiley-Inter science Publication JOHN WILEY & SONS^ Chichester New York Weinheim Brisbane Singapore
More informationFor example, the system. 22 may be represented by the augmented matrix
Matrix Solutions to Linear Systems A matrix is a rectangular array of elements. o An array is a systematic arrangement of numbers or symbols in rows and columns. Matrices (the plural of matrix) may be
More informationCS6015 / LARP ACK : Linear Algebra and Its Applications - Gilbert Strang
Solving and CS6015 / LARP 2018 ACK : Linear Algebra and Its Applications - Gilbert Strang Introduction Chapter 1 concentrated on square invertible matrices. There was one solution to Ax = b and it was
More informationPackage PUlasso. April 7, 2018
Type Package Package PUlasso April 7, 2018 Title High-Dimensional Variable Selection with Presence-Only Data Version 3.1 Date 2018-4-4 Efficient algorithm for solving PU (Positive and Unlabelled) problem
More informationLARP / 2018 ACK : 1. Linear Algebra and Its Applications - Gilbert Strang 2. Autar Kaw, Transforming Numerical Methods Education for STEM Graduates
Triangular Factors and Row Exchanges LARP / 28 ACK :. Linear Algebra and Its Applications - Gilbert Strang 2. Autar Kaw, Transforming Numerical Methods Education for STEM Graduates Then there were three
More informationPackage nmslibr. April 14, 2018
Type Package Title Non Metric Space (Approximate) Library Version 1.0.1 Date 2018-04-14 Package nmslibr April 14, 2018 Maintainer Lampros Mouselimis BugReports https://github.com/mlampros/nmslibr/issues
More informationCurriculum Map: Mathematics
Curriculum Map: Mathematics Course: Honors Advanced Precalculus and Trigonometry Grade(s): 11-12 Unit 1: Functions and Their Graphs This chapter will develop a more complete, thorough understanding of
More informationPackage kdtools. April 26, 2018
Type Package Package kdtools April 26, 2018 Title Tools for Working with Multidimensional Data Version 0.3.1 Provides various tools for working with multidimensional data in R and C++, including etremely
More informationPackage geojsonsf. R topics documented: January 11, Type Package Title GeoJSON to Simple Feature Converter Version 1.3.
Type Package Title GeoJSON to Simple Feature Converter Version 1.3.0 Date 2019-01-11 Package geojsonsf January 11, 2019 Converts Between GeoJSON and simple feature objects. License GPL-3 Encoding UTF-8
More informationPackage jdx. R topics documented: January 9, Type Package Title 'Java' Data Exchange for 'R' and 'rjava'
Type Package Title 'Java' Data Exchange for 'R' and 'rjava' Package jdx January 9, 2018 Description Simplifies and extends data exchange between 'R' and 'Java'. Version 0.1.0 License GPL (>= 2 BSD_3_clause
More informationPackage canvasxpress
Version 1.18.2 Package canvasxpress Title Visualization Package for CanvasXpress in R January 19, 2018 Enables creation of visualizations using the CanvasXpress framework in R. CanvasXpress is a standalone
More information10/26/ Solving Systems of Linear Equations Using Matrices. Objectives. Matrices
6.1 Solving Systems of Linear Equations Using Matrices Objectives Write the augmented matrix for a linear system. Perform matrix row operations. Use matrices and Gaussian elimination to solve systems.
More informationFinite Math Linear Programming 1 May / 7
Linear Programming Finite Math 1 May 2017 Finite Math Linear Programming 1 May 2017 1 / 7 General Description of Linear Programming Finite Math Linear Programming 1 May 2017 2 / 7 General Description of
More informationPackage BigVAR. R topics documented: April 3, Type Package
Type Package Package BigVAR April 3, 2017 Title Dimension Reduction Methods for Multivariate Time Series Version 1.0.2 Date 2017-03-24 Estimates VAR and VARX models with structured Lasso Penalties. Depends
More informationPackage geneslope. October 26, 2016
Type Package Package geneslope October 26, 2016 Title Genome-Wide Association Study with SLOPE Version 0.37.0 Date 2016-10-26 Genome-wide association study (GWAS) performed with SLOPE, short for Sorted
More informationPackage elmnnrcpp. July 21, 2018
Type Package Title The Extreme Learning Machine Algorithm Version 1.0.1 Date 2018-07-21 Package elmnnrcpp July 21, 2018 BugReports https://github.com/mlampros/elmnnrcpp/issues URL https://github.com/mlampros/elmnnrcpp
More informationPackage RSpectra. May 23, 2018
Type Package Package RSpectra May 23, 2018 Title Solvers for Large-Scale Eigenvalue and SVD Problems Version 0.13-1 Date 2018-05-22 Description R interface to the 'Spectra' library
More informationPackage cattonum. R topics documented: May 2, Type Package Version Title Encode Categorical Features
Type Package Version 0.0.2 Title Encode Categorical Features Package cattonum May 2, 2018 Functions for dummy encoding, frequency encoding, label encoding, leave-one-out encoding, mean encoding, median
More informationhp calculators hp 39g+ & hp 39g/40g Using Matrices How are matrices stored? How do I solve a system of equations? Quick and easy roots of a polynomial
hp calculators hp 39g+ Using Matrices Using Matrices The purpose of this section of the tutorial is to cover the essentials of matrix manipulation, particularly in solving simultaneous equations. How are
More informationPackage raker. October 10, 2017
Title Easy Spatial Microsimulation (Raking) in R Version 0.2.1 Date 2017-10-10 Package raker October 10, 2017 Functions for performing spatial microsimulation ('raking') in R. Depends R (>= 3.4.0) License
More informationFinite Math - J-term Homework. Section Inverse of a Square Matrix
Section.5-77, 78, 79, 80 Finite Math - J-term 017 Lecture Notes - 1/19/017 Homework Section.6-9, 1, 1, 15, 17, 18, 1, 6, 9, 3, 37, 39, 1,, 5, 6, 55 Section 5.1-9, 11, 1, 13, 1, 17, 9, 30 Section.5 - Inverse
More informationPackage NNLM. June 18, 2018
Type Package Package NNLM June 18, 2018 Title Fast and Versatile Non-Negative Matrix Factorization Description This is a package for Non-Negative Linear Models (NNLM). It implements fast sequential coordinate
More informationExercise Set Decide whether each matrix below is an elementary matrix. (a) (b) (c) (d) Answer:
Understand the relationships between statements that are equivalent to the invertibility of a square matrix (Theorem 1.5.3). Use the inversion algorithm to find the inverse of an invertible matrix. Express
More informationPackage frequencyconnectedness
Type Package Package frequencyconnectedness September 24, 2018 Title Spectral Decomposition of Connectedness Measures Version 0.2.1 Date 2018-09-24 Accompanies a paper (Barunik, Krehlik (2018) )
More informationLinear and Integer Programming :Algorithms in the Real World. Related Optimization Problems. How important is optimization?
Linear and Integer Programming 15-853:Algorithms in the Real World Linear and Integer Programming I Introduction Geometric Interpretation Simplex Method Linear or Integer programming maximize z = c T x
More informationPackage errorlocate. R topics documented: March 30, Type Package Title Locate Errors with Validation Rules Version 0.1.3
Type Package Title Locate Errors with Validation Rules Version 0.1.3 Package errorlocate March 30, 2018 Errors in data can be located and removed using validation rules from package 'validate'. License
More informationPackage ramcmc. May 29, 2018
Title Robust Adaptive Metropolis Algorithm Version 0.1.0-1 Date 2018-05-29 Author Jouni Helske Package ramcmc May 29, 2018 Maintainer Jouni Helske Function for adapting the shape
More informationPackage Cubist. December 2, 2017
Type Package Package Cubist December 2, 2017 Title Rule- And Instance-Based Regression Modeling Version 0.2.1 Date 2017-12-01 Maintainer Max Kuhn Description Regression modeling using
More informationPackage nngeo. September 29, 2018
Type Package Title k-nearest Neighbor Join for Spatial Data Version 0.2.2 Package nngeo September 29, 2018 K-nearest neighbor search for projected and non-projected 'sf' spatial layers. Nearest neighbor
More informationChap5 The Theory of the Simplex Method
College of Management, NCTU Operation Research I Fall, Chap The Theory of the Simplex Method Terminology Constraint oundary equation For any constraint (functional and nonnegativity), replace its,, sign
More informationPackage filling. December 11, 2017
Type Package Package filling December 11, 2017 Title Matrix Completion, Imputation, and Inpainting Methods Version 0.1.0 Filling in the missing entries of a partially observed data is one of fundamental
More informationPackage svalues. July 15, 2018
Type Package Package svalues July 15, 2018 Title Measures of the Sturdiness of Regression Coefficients Version 0.1.6 Author Carlos Cinelli Maintainer Carlos Cinelli Implements
More informationSome Advanced Topics in Linear Programming
Some Advanced Topics in Linear Programming Matthew J. Saltzman July 2, 995 Connections with Algebra and Geometry In this section, we will explore how some of the ideas in linear programming, duality theory,
More informationPackage ompr. November 18, 2017
Type Package Package ompr November 18, 2017 Title Model and Solve Mixed Integer Linear Programs Version 0.7.0 Model mixed integer linear programs in an algebraic way directly in R. The is solver-independent
More informationPackage tropicalsparse
Type Package Title Sparse Tropical Algebra Version 0.1.0 Date 2018-06-25 Package tropicalsparse July 1, 2018 Some of the basic tropical algebra functionality is provided for sparse matrices by applying
More informationPackage Rtsne. April 14, 2017
Type Package Package Rtsne April 14, 2017 Title T-Distributed Stochastic Neighbor Embedding using a Barnes-Hut Implementation Version 0.13 Description An R wrapper around the fast T-distributed Stochastic
More informationInteger Programming Theory
Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x
More informationPackage mgc. April 13, 2018
Type Package Title Multiscale Graph Correlation Version 1.0.1 Date 2018-04-12 Package mgc April 13, 2018 Maintainer Multiscale Graph Correlation (MGC) is a framework developed by Shen
More informationPackage MTLR. March 9, 2019
Type Package Package MTLR March 9, 2019 Title Survival Prediction with Multi-Task Logistic Regression Version 0.2.0 Author Humza Haider Maintainer Humza Haider URL https://github.com/haiderstats/mtlr
More informationMultiple View Geometry in Computer Vision
Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 Structure Computation Lecture 18 March 22, 2005 2 3D Reconstruction The goal of 3D reconstruction
More informationPackage svd. R topics documented: March 24, Type Package. Title Interfaces to various state-of-art SVD and eigensolvers. Version 0.3.
Package svd March 24, 2013 Type Package Title Interfaces to various state-of-art SVD and eigensolvers Version 0.3.1-1 Date 2013-03-24 Author Anton Korobeynikov Maintainer Anton Korobeynikov
More informationInterlude: Solving systems of Equations
Interlude: Solving systems of Equations Solving Ax = b What happens to x under Ax? The singular value decomposition Rotation matrices Singular matrices Condition number Null space Solving Ax = 0 under
More informationPackage spark. July 21, 2017
Title 'Sparklines' in the 'R' Terminal Version 2.0.0 Author Gábor Csárdi Package spark July 21, 2017 Maintainer Gábor Csárdi A 'sparkline' is a line chart, without axes and labels.
More informationEE613 Machine Learning for Engineers LINEAR REGRESSION. Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov.
EE613 Machine Learning for Engineers LINEAR REGRESSION Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov. 4, 2015 1 Outline Multivariate ordinary least squares Singular value
More informationPackage validara. October 19, 2017
Type Package Title Validate Brazilian Administrative Registers Version 0.1.1 Package validara October 19, 2017 Maintainer Gustavo Coelho Contains functions to validate administrative
More informationPackage modmarg. R topics documented:
Package modmarg February 1, 2018 Title Calculating Marginal Effects and Levels with Errors Version 0.9.2 Calculate predicted levels and marginal effects, using the delta method to calculate standard errors.
More informationPackage EvolutionaryGames
Type Package Package EvolutionaryGames November 27, 2017 Title Important Concepts of Evolutionary Game Theory Version 0.1.0 Maintainer Jochen Staudacher A comprehensive
More informationPackage sparsehessianfd
Type Package Package sparsehessianfd Title Numerical Estimation of Sparse Hessians Version 0.3.3.3 Date 2018-03-26 Maintainer Michael Braun March 27, 2018 URL http://www.smu.edu/cox/departments/facultydirectory/braunmichael
More informationPackage semisup. March 10, Version Title Semi-Supervised Mixture Model
Version 1.7.1 Title Semi-Supervised Mixture Model Package semisup March 10, 2019 Description Useful for detecting SNPs with interactive effects on a quantitative trait. This R packages moves away from
More informationPackage eply. April 6, 2018
Type Package Title Apply a Function Over Expressions Version 0.1.2 Package eply April 6, 2018 Evaluate a function over a data frame of expressions. License GPL-3 Depends R (>= 3.0.0) Imports magrittr,
More informationOptimization of Design. Lecturer:Dung-An Wang Lecture 8
Optimization of Design Lecturer:Dung-An Wang Lecture 8 Lecture outline Reading: Ch8 of text Today s lecture 2 8.1 LINEAR FUNCTIONS Cost Function Constraints 3 8.2 The standard LP problem Only equality
More informationPackage bisect. April 16, 2018
Package bisect April 16, 2018 Title Estimating Cell Type Composition from Methylation Sequencing Data Version 0.9.0 Maintainer Eyal Fisher Author Eyal Fisher [aut, cre] An implementation
More informationEE613 Machine Learning for Engineers LINEAR REGRESSION. Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov.
EE613 Machine Learning for Engineers LINEAR REGRESSION Sylvain Calinon Robot Learning & Interaction Group Idiap Research Institute Nov. 9, 2017 1 Outline Multivariate ordinary least squares Matlab code:
More informationPackage nlgeocoder. October 8, 2018
Title Geocoding for the Netherlands Version 0.1.3 Package nlgeocoder October 8, 2018 R interface to the open location server API of 'Publieke Diensten Op de Kaart' (). It offers geocoding,
More informationPackage balance. October 12, 2018
Title Visualize Balances of Compositional Data Version 0.1.6 URL http://github.com/tpq/balance Package balance October 12, 2018 BugReports http://github.com/tpq/balance/issues Balances have become a cornerstone
More informationPackage RODBCext. July 31, 2017
Version 0.3.1 Package RODBCext July 31, 2017 Title Parameterized Queries Extension for RODBC An extension for RODBC package adding support for parameterized queries. URL https://github.com/zozlak/rodbcext
More informationPackage pwrab. R topics documented: June 6, Type Package Title Power Analysis for AB Testing Version 0.1.0
Type Package Title Power Analysis for AB Testing Version 0.1.0 Package pwrab June 6, 2017 Maintainer William Cha Power analysis for AB testing. The calculations are based
More informationPackage PCADSC. April 19, 2017
Type Package Package PCADSC April 19, 2017 Title Tools for Principal Component Analysis-Based Data Structure Comparisons Version 0.8.0 A suite of non-parametric, visual tools for assessing differences
More informationSolving Systems of Equations Using Matrices With the TI-83 or TI-84
Solving Systems of Equations Using Matrices With the TI-83 or TI-84 Dimensions of a matrix: The dimensions of a matrix are the number of rows by the number of columns in the matrix. rows x columns *rows
More informationPackage rnn. R topics documented: June 21, Title Recurrent Neural Network Version 0.8.1
Title Recurrent Neural Network Version 0.8.1 Package rnn June 21, 2018 Implementation of a Recurrent Neural Network in R. Depends R (>= 3.2.2) License GPL-3 LazyData true RoxygenNote 6.0.1 URL http://qua.st/rnn,
More informationPackage assertr. R topics documented: February 23, Type Package
Type Package Package assertr February 23, 2018 Title Assertive Programming for R Analysis Pipelines Version 2.5 Provides functionality to assert conditions that have to be met so that errors in data used
More informationPackage bayesdp. July 10, 2018
Type Package Package bayesdp July 10, 2018 Title Tools for the Bayesian Discount Prior Function Version 1.3.2 Date 2018-07-10 Depends R (>= 3.2.3), ggplot2, survival, methods Functions for data augmentation
More informationPackage pdfsearch. July 10, 2018
Type Package Version 0.2.3 License MIT + file LICENSE Title Search Tools for PDF Files Package pdfsearch July 10, 2018 Includes functions for keyword search of pdf files. There is also a wrapper that includes
More informationPackage ADMM. May 29, 2018
Type Package Package ADMM May 29, 2018 Title Algorithms using Alternating Direction Method of Multipliers Version 0.3.0 Provides algorithms to solve popular optimization problems in statistics such as
More informationOn Unbounded Tolerable Solution Sets
Reliable Computing (2005) 11: 425 432 DOI: 10.1007/s11155-005-0049-9 c Springer 2005 On Unbounded Tolerable Solution Sets IRENE A. SHARAYA Institute of Computational Technologies, 6, Acad. Lavrentiev av.,
More informationPackage datasets.load
Title Interface for Loading Datasets Version 0.1.0 Package datasets.load December 14, 2016 Visual interface for loading datasets in RStudio from insted (unloaded) s. Depends R (>= 3.0.0) Imports shiny,
More informationPackage glmnetutils. August 1, 2017
Type Package Version 1.1 Title Utilities for 'Glmnet' Package glmnetutils August 1, 2017 Description Provides a formula interface for the 'glmnet' package for elasticnet regression, a method for cross-validating
More informationLARGE SCALE LINEAR AND INTEGER OPTIMIZATION: A UNIFIED APPROACH
LARGE SCALE LINEAR AND INTEGER OPTIMIZATION: A UNIFIED APPROACH Richard Kipp Martin Graduate School of Business University of Chicago % Kluwer Academic Publishers Boston/Dordrecht/London CONTENTS Preface
More informationPackage splithalf. March 17, 2018
Type Package Package splithalf March 17, 2018 Title Calculate Task Split Half Reliability Estimates Version 0.3.1 Maintainer Sam Parsons A series of functions to calculate the
More informationMEMORY EFFICIENT WDR (WAVELET DIFFERENCE REDUCTION) using INVERSE OF ECHELON FORM by EQUATION SOLVING
Available Online at www.ijcsmc.com International Journal of Computer Science and Mobile Computing A Monthly Journal of Computer Science and Information Technology IJCSMC Vol. 3 Issue. 7 July 2014 pg.512
More informationPackage simr. April 30, 2018
Type Package Package simr April 30, 2018 Title Power Analysis for Generalised Linear Mixed Models by Simulation Calculate power for generalised linear mixed models, using simulation. Designed to work with
More informationCS354 Computer Graphics Rotations and Quaternions
Slide Credit: Don Fussell CS354 Computer Graphics Rotations and Quaternions Qixing Huang April 4th 2018 Orientation Position and Orientation The position of an object can be represented as a translation
More informationPackage colf. October 9, 2017
Type Package Package colf October 9, 2017 Title Constrained Optimization on Linear Function Version 0.1.3 URL https://github.com/lyzander/colf BugReports https://github.com/lyzander/colf/issues Depends
More informationObjectives and Homework List
MAC 1140 Objectives and Homework List Each objective covered in MAC1140 is listed below. Along with each objective is the homework list used with MyMathLab (MML) and a list to use with the text (if you
More informationArtificial Intelligence
Artificial Intelligence Combinatorial Optimization G. Guérard Department of Nouvelles Energies Ecole Supérieur d Ingénieurs Léonard de Vinci Lecture 1 GG A.I. 1/34 Outline 1 Motivation 2 Geometric resolution
More informationVector Algebra Transformations. Lecture 4
Vector Algebra Transformations Lecture 4 Cornell CS4620 Fall 2008 Lecture 4 2008 Steve Marschner 1 Geometry A part of mathematics concerned with questions of size, shape, and relative positions of figures
More informationPackage lhs. R topics documented: January 4, 2018
Package lhs January 4, 2018 Version 0.16 Date 2017-12-23 Title Latin Hypercube Samples Author [aut, cre] Maintainer Depends R (>= 3.3.0) Suggests RUnit Provides a number of methods
More information