new RowMatrix(rows, nRowsopt, nColsopt)
Represents a row-oriented distributed Matrix with no meaningful row indices.
Parameters:
Name | Type | Attributes | Description |
---|---|---|---|
rows |
module:eclairjs.RDD | stored as an RDD[Vector] | |
nRows |
number |
<optional> |
number of rows. A non-positive value means unknown, and then the number of rows will be determined by the number of records in the RDD `rows`. |
nCols |
number |
<optional> |
number of columns. A non-positive value means unknown, and then the number of columns will be determined by the size of the first row. |
Example
var RowMatrix = require('eclairjs/mllib/linalg/distributed/RowMatrix');
var Vectors = require("'eclairjs/mllib/linalg/Vectors");;
var rowsList = [Vectors.dense([1.12, 2.05, 3.12]), Vectors.dense([5.56, 6.28, 8.94]), Vectors.dense([10.2, 8.0, 20.5])];
var rows = sc.parallelize(rowsList);
var mat = new RowMatrix(rows);
Extends
Methods
columnSimilarities(thresholdopt) → {CoordinateMatrix}
Compute similarities between columns of this matrix using a sampling approach.
The threshold parameter is a trade-off knob between estimate quality and computational cost.
Setting a threshold of 0 guarantees deterministic correct results, but comes at exactly
the same cost as the brute-force approach. Setting the threshold to positive values
incurs strictly less computational cost than the brute-force approach, however the
similarities computed will be estimates.
The sampling guarantees relative-error correctness for those pairs of columns that have
similarity greater than the given similarity threshold.
To describe the guarantee, we set some notation:
Let A be the smallest in magnitude non-zero element of this matrix.
Let B be the largest in magnitude non-zero element of this matrix.
Let L be the maximum number of non-zeros per row.
For example, for {0,1} matrices: A=B=1.
Another example, for the Netflix matrix: A=1, B=5
For those column pairs that are above the threshold,
the computed similarity is correct to within 20% relative error with probability
at least 1 - (0.981)^10/B^
The shuffle size is bounded by the *smaller* of the following two expressions:
O(n log(n) L / (threshold * A))
O(m L^2^)
The latter is the cost of the brute-force approach, so for non-zero thresholds,
the cost is always cheaper than the brute-force approach.
Parameters:
Name | Type | Attributes | Description |
---|---|---|---|
threshold |
number |
<optional> |
Set to 0 for deterministic guaranteed correctness. Similarities above this threshold are estimated with the cost vs estimate quality trade-off described above. between columns of this matrix. |
Returns:
An n x n sparse upper-triangular matrix of cosine similarities
- Type
- CoordinateMatrix
computeCovariance() → {module:eclairjs/mllib/linalg.Matrix}
Computes the covariance matrix, treating each row as an observation. Note that this cannot
be computed on matrices with more than 65535 columns.
Returns:
a local dense matrix of size n x n
computeGramianMatrix() → {module:eclairjs/mllib/linalg.Matrix}
Computes the Gramian matrix `A^T A`. Note that this cannot be computed on matrices with
more than 65535 columns.
Returns:
computePrincipalComponents(k) → {module:eclairjs/mllib/linalg.Matrix}
Computes the top k principal components.
Rows correspond to observations and columns correspond to variables.
The principal components are stored a local matrix of size n-by-k.
Each column corresponds for one principal component,
and the columns are in descending order of component variance.
The row data do not need to be "centered" first; it is not necessary for
the mean of each column to be 0.
Note that this cannot be computed on matrices with more than 65535 columns.
Parameters:
Name | Type | Description |
---|---|---|
k |
integer | number of top principal components. |
Returns:
a matrix of size n-by-k, whose columns are principal components
computeSVD(k, computeU, rCond) → {module:eclairjs/mllib/linalg.SingularValueDecomposition}
Computes singular value decomposition of this matrix. Denote this matrix by A (m x n). This
will compute matrices U, S, V such that A ~= U * S * V', where S contains the leading k
singular values, U and V contain the corresponding singular vectors.
At most k largest non-zero singular values and associated vectors are returned. If there are k
such values, then the dimensions of the return will be:
- U is a RowMatrix of size m x k that satisfies U' * U = eye(k),
- s is a Vector of size k, holding the singular values in descending order,
- V is a Matrix of size n x k that satisfies V' * V = eye(k).
We assume n is smaller than m, though this is not strictly required.
The singular values and the right singular vectors are derived
from the eigenvalues and the eigenvectors of the Gramian matrix A' * A. U, the matrix
storing the right singular vectors, is computed via matrix multiplication as
U = A * (V * S^-1^), if requested by user. The actual method to use is determined
automatically based on the cost:
- If n is small (n < 100) or k is large compared with n (k > n / 2), we compute
the Gramian matrix first and then compute its top eigenvalues and eigenvectors locally
on the driver. This requires a single pass with O(n^2^) storage on each executor and
on the driver, and O(n^2^ k) time on the driver.
- Otherwise, we compute (A' * A) * v in a distributive way and send it to ARPACK's DSAUPD to
compute (A' * A)'s top eigenvalues and eigenvectors on the driver node. This requires O(k)
passes, O(n) storage on each executor, and O(n k) storage on the driver.
Several internal parameters are set to default values. The reciprocal condition number rCond
is set to 1e-9. All singular values smaller than rCond * sigma(0) are treated as zeros, where
sigma(0) is the largest singular value. The maximum number of Arnoldi update iterations for
ARPACK is set to 300 or k * 3, whichever is larger. The numerical tolerance for ARPACK's
eigen-decomposition is set to 1e-10.
Parameters:
Name | Type | Description |
---|---|---|
k |
integer | number of leading singular values to keep (0 < k <= n). It might return less than k if there are numerically zero singular values or there are not enough Ritz values converged before the maximum number of Arnoldi update iterations is reached (in case that matrix A is ill-conditioned). |
computeU |
boolean | whether to compute U |
rCond |
float | the reciprocal condition number. All singular values smaller than rCond * sigma(0) are treated as zero, where sigma(0) is the largest singular value. |
Returns:
SingularValueDecomposition(U, s, V). U = null if computeU = false.
multiply(B) → {module:eclairjs/mllib/linalg/distributed.RowMatrix}
Multiply this matrix by a local matrix on the right.
Parameters:
Name | Type | Description |
---|---|---|
B |
module:eclairjs/mllib/linalg.Matrix | a local matrix whose number of rows must match the number of columns of this matrix which preserves partitioning |
Returns:
a [[org.apache.spark.mllib.linalg.distributed.RowMatrix]] representing the product,
numCols() → {integer}
Gets or computes the number of rows.
- Inherited From:
- Source:
Returns:
- Type
- integer
numRows() → {integer}
Gets or computes the number of rows.
- Inherited From:
- Source:
Returns:
- Type
- integer
rows() → {module:eclairjs.RDD}
Gets RDD of Vectors
Returns:
RDD of Vectors
- Type
- module:eclairjs.RDD
tallSkinnyQR(computeQ) → {module:eclairjs/mllib/linalg.QRDecomposotion}
Compute QR decomposition for RowMatrix. The implementation is designed to optimize the QR
decomposition (factorization) for the RowMatrix of a tall and skinny shape.
Reference:
Paul G. Constantine, David F. Gleich. "Tall and skinny QR factorizations in MapReduce
architectures" ([[http://dx.doi.org/10.1145/1996092.1996103]])
Parameters:
Name | Type | Description |
---|---|---|
computeQ |
boolean | whether to computeQ |
Returns:
QRDecomposition(Q, R), Q = null if computeQ = false.
- Type
- module:eclairjs/mllib/linalg.QRDecomposotion