Struct petal_decomposition::Pca

pub struct Pca<A> where
    A: Scalar,  { /* private fields */ }

Principal component analysis.

This reduces the dimensionality of the input data using Singular Value Decomposition (SVD). The data is centered for each feature before applying SVD.

Examples

use petal_decomposition::PcaBuilder;

let x = ndarray::arr2(&[[0_f64, 0_f64], [1_f64, 1_f64], [2_f64, 2_f64]]);
let y = PcaBuilder::new(1).build().fit_transform(&x).unwrap();  // [-2_f64.sqrt(), 0_f64, 2_f64.sqrt()]
assert!((y[(0, 0)].abs() - 2_f64.sqrt()).abs() < 1e-8);
assert!(y[(1, 0)].abs() < 1e-8);
assert!((y[(2, 0)].abs() - 2_f64.sqrt()).abs() < 1e-8);

Implementations

impl<A> Pca<A> where
    A: Scalar, 

pub fn new(n_components: usize) -> Self

Creates a PCA model with the given number of components.

impl<A> Pca<A> where
    A: FromPrimitive + Lapack,
    A::Real: ScalarOperand, 

pub fn components(&self) -> &Array2<A>

Returns the principal axes in feature space.

pub fn mean(&self) -> &Array1<A>

Returns the per-feature empirical mean.

pub fn n_components(&self) -> usize

Returns the number of components.

pub fn singular_values(&self) -> &Array1<A::Real>

Returns sigular values.

pub fn explained_variance_ratio(&self) -> Array1<A::Real>

Returns the ratio of explained variance for each component.

pub fn fit<S>(
    &mut self,
    input: &ArrayBase<S, Ix2>
) -> Result<(), DecompositionError> where
    S: Data<Elem = A>, 

Fits the model with input.

Errors

• DecompositionError::InvalidInput if any of the dimensions of input is less than the number of components, or the layout of input is incompatible with LAPACK.
• DecompositionError::LinalgError if the underlying Singular Vector Decomposition routine fails.

pub fn transform<S>(
    &self,
    input: &ArrayBase<S, Ix2>
) -> Result<Array2<A>, DecompositionError> where
    S: Data<Elem = A>, 

Applies dimensionality reduction to input.

Errors

• DecompositionError::InvalidInput if the number of features in input does not match that of the training data.

pub fn fit_transform<S>(
    &mut self,
    input: &ArrayBase<S, Ix2>
) -> Result<Array2<A>, DecompositionError> where
    S: Data<Elem = A>, 

Fits the model with input and apply the dimensionality reduction on input.

This is equivalent to calling both fit and transform for the same input, but more efficient.

Errors

• DecompositionError::InvalidInput if any of the dimensions of input is less than the number of components, or the layout of input is incompatible with LAPACK.
• DecompositionError::LinalgError if the underlying Singular Vector Decomposition routine fails.

pub fn inverse_transform<S>(
    &self,
    input: &ArrayBase<S, Ix2>
) -> Result<Array2<A>, DecompositionError> where
    S: Data<Elem = A>, 

Transforms data back to its original space.

Errors

Returns DecompositionError::InvalidInput if the number of rows of input is different from that of the training data, or the number of columns of input is different from the number of components.