Matrix Functions for Teaching and Learning Linear Algebra and Multivariate Statistics, http://friendly.github.io/matlib/
Version 0.9.7
These functions are mainly for tutorial purposes in teaching & learning matrix algebra ideas using R. In some cases, functions are provided for concepts or computations available elsewhere in R, but where the name is not obvious, e.g., R()
for the rank of a matrix, or tr()
for matrix trace.
In other cases, we provide cover functions to show or demonstrate an algorithm in more detail, sometimes providing a verbose =
argument to print the details of computations, e.g., Det()
for a matrix determinant, Inverse()
for a matrix inverse, using gaussianElimination()
to show the steps.
In addition, a collection of functions are provided for drawing vector diagrams in 2D and 3D, illustrating various concepts of linear algebra more concretely than has been available before. For example,
showEqn(A, b)
shows the matrix equations Ax = b in text or LaTeX form, while plotEqn(A, b)
and plotEqn3d(A, b)
plots those equations in 2D or 3D space.
vectors()
, vectors3d()
plot geometric vector diagrams in 2D or 3D, with other functions to draw angles and arcs.
regvec3d()
calculates and plot vectors representing a bivariate regression model, lm(y ~ x1 + x2)
Get the released version from CRAN:
The development version can be installed to your R library directly from this repo via:
if (!require(remotes)) install.packages("remotes")
remotes::install_github("friendly/matlib", build_vignettes = TRUE)
This installs the package from the source and creates the package vignettes, so you will need to have R Tools installed on your system. R Tools for Windows takes you to the download page for Windows. R Tools for Mac OS X has the required programs for Mac OS X.
The functions that draw 3D graphs use the rgl package. On macOS, rgl requires that XQuartz be installed. After installing XQuartz, it’s necessary either to log out of and back into your macOS account or to reboot your Mac.
The functions in this package are grouped under the following topics
tr()
- trace of a matrixR()
- rank of a matrixJ()
- constant vector, matrix or arraylen()
- Euclidean length of a vector or columns of a matrixvec()
- vectorize a matrixProj(y, X)
- projection of vector y on columns of matrix X
mpower(A, p)
- matrix powers for a square symmetric matrixxprod(...)
- vector cross-productminor()
- Minor of A[i,j]cofactor()
- Cofactor of A[i,j]rowMinors()
- Row minors of A[i,]rowCofactors()
- Row cofactors of A[i,]Det()
- Determinants by elimination or eigenvaluesrowadd()
- Add multiples of rows to other rowsrowmult()
- Multiply rows by constantsrowswap()
- Interchange two rows of a matrixshowEqn(A, b)
- show matrices (A, b) as linear equationsplotEqn(A, b)
, plotEqn3d(A, b)
- plot matrices (A, b) as linear equationsverbose=TRUE
argument to show the intermediate steps and a fractions=TRUE
argument to show results using MASS::fractions()
.gaussianElimination(A, B)
- reduces (A,B) to (I,A−1B)
Inverse(X)
, inv()
- uses gaussianElimination
to find the inverse of X, X−1
echelon(X)
- uses gaussianElimination
to find the reduced echelon form of XGinv(X)
- uses gaussianElimination
to find the generalized inverse of XLU(X)
- LU decomposition of a matrix Xcholesky(X)
- calculates a Cholesky square root of a matrixswp()
- matrix sweep operatorEigen()
- eigenvalues and eigenvectorsSVD()
- singular value decomposition, $mathbf{X = U D V}$powerMethod()
- find the dominant eigenvector using the power methodshowEig()
- draw eigenvectors on a 2D scatterplot with a dataEllipseMoorePenrose()
- illustrates how the Moore-Penrose inverse can be calculated using SVD()
arrows3d()
- draw nice 3D arrowscorner()
, arc()
- draw a corner or arc showing the angle between two vectors in 2D/3DpointOnLine()
- position of a point along a linevectors()
, vectors3d()
- plot geometric vector diagrams in 2D/3Dregvec3d()
- calculate and plot vectors representing a bivariate regression model, lm(y ~ x1 + x2)
in mean-deviation form.A small collection of vignettes is now available. Use browseVignettes("matlib")
to explore them.
Vignette | Title |
---|---|
det-ex1 | Properties of determinants |
det-ex2 | Evaluation of determinants |
inv-ex1 | Inverse of a matrix |
inv-ex2 | Matrix inversion by elementary row operations |
ginv | Generalized inverse |
eigen-ex1 | Eigenvalues and Eigenvectors: Properties |
eigen-ex2 | Eigenvalues: Spectral Decomposition |
linear-equations | Solving Linear Equations |
gramreg | Gram-Schmidt Orthogonalization and Regression |
data-beta | Vector Spaces of Least Squares and Linear Equations |
See also:
Fox & Friendly, Visualizing Simultaneous Linear Equations, Geometric Vectors, and Least-Squares Regression with the matlib Package for R, June 2016, useR! Conference, Stanford.
Ivan Savov, Linear algebra explained in four pages