https://github.com/tensorly/tensorly
Tip revision: 23c6d135ae81f12205d4d57c22652628ecc2adf1 authored by Jean Kossaifi on 06 January 2017, 13:34:08 UTC
Better github ribbon (especially for mobiles)
Better github ribbon (especially for mobiles)
Tip revision: 23c6d13
_higher_order_moment.py
from ._khatri_rao import khatri_rao
# Author: Jean Kossaifi
# License: BSD 3 clause
def higher_order_moment(matrix, order=3):
"""Higher order moment of a matrix of observations
Computes the higher moment of order `order` of `matrix`
Each row of `matrix` represents a samples
(i.e. an observation)
Parameters
----------
matrix : 2D-array
array of shape (n_samples, n_features)
i.e. each row is a sample
order : int, optional
order of the moment to compute
"""
matrix = matrix - matrix.mean(axis=0)
n_features = matrix.shape[-1]
t = khatri_rao([matrix.T] * order).mean(axis=1)
return t.reshape([n_features] * order)