Shortcuts

torch.cov

torch.cov(input, *, correction=1, fweights=None, aweights=None) 张量

估计由input矩阵给出的变量的协方差矩阵,其中行是变量,列是观测值。

协方差矩阵是一个方阵,给出每对变量之间的协方差。对角线包含每个变量的方差(变量与其自身的协方差)。根据定义,如果 input 表示单个变量(标量或一维),则返回其方差。

变量 xxyy 的样本协方差由以下公式给出:

cov(x,y)=i=1N(xixˉ)(yiyˉ)max(0, N  δN)\text{cov}(x,y) = \frac{\sum^{N}_{i = 1}(x_{i} - \bar{x})(y_{i} - \bar{y})}{\max(0,~N~-~\delta N)}

其中 xˉ\bar{x}yˉ\bar{y} 分别是 xxyy 的简单平均值,并且 δN\delta Ncorrection

如果提供了fweights和/或aweights,则计算加权协方差,其公式为:

covw(x,y)=i=1Nwi(xiμx)(yiμy)max(0, i=1Nwi  i=1Nwiaii=1Nwi δN)\text{cov}_w(x,y) = \frac{\sum^{N}_{i = 1}w_i(x_{i} - \mu_x^*)(y_{i} - \mu_y^*)} {\max(0,~\sum^{N}_{i = 1}w_i~-~\frac{\sum^{N}_{i = 1}w_ia_i}{\sum^{N}_{i = 1}w_i}~\delta N)}

其中 ww 表示 fweightsaweights(为简洁起见,fa)基于提供的哪一个,或者如果两者都提供,则为 w=f×aw = f \times a,并且 μx=i=1Nwixii=1Nwi\mu_x^* = \frac{\sum^{N}_{i = 1}w_ix_{i} }{\sum^{N}_{i = 1}w_i} 是变量的加权平均值。如果没有提供,f 和/或 a 可以被视为适当大小的 1\mathbb{1} 向量。

Parameters

输入 (张量) – 一个包含多个变量和观测值的二维矩阵,或表示单个变量的标量或一维向量。

Keyword Arguments
  • 校正 (int, 可选) – 样本大小与样本自由度之间的差异。 默认为贝塞尔校正,correction = 1,即使指定了fweightsaweights,也会返回无偏估计。correction = 0 将返回简单平均值。默认为1

  • fweights (tensor, 可选) – 一个标量或1D张量,表示观测向量的频率,表示每个观测值应重复的次数。其numel必须等于input的列数。必须具有整数数据类型。如果为None,则忽略。默认为None

  • aweights (tensor, 可选) – 一个标量或一维观测向量权重数组。 这些相对权重通常对于被认为是“重要”的观测值较大,而对于被认为是较不“重要”的观测值较小。其数量必须等于input的列数。 必须具有浮点型数据类型。如果为None,则忽略。默认为None

Returns

(Tensor) 变量的协方差矩阵。

另请参阅

torch.corrcoef() 标准化协方差矩阵。

Example::
>>> x = torch.tensor([[0, 2], [1, 1], [2, 0]]).T
>>> x
tensor([[0, 1, 2],
        [2, 1, 0]])
>>> torch.cov(x)
tensor([[ 1., -1.],
        [-1.,  1.]])
>>> torch.cov(x, correction=0)
tensor([[ 0.6667, -0.6667],
        [-0.6667,  0.6667]])
>>> fw = torch.randint(1, 10, (3,))
>>> fw
tensor([1, 6, 9])
>>> aw = torch.rand(3)
>>> aw
tensor([0.4282, 0.0255, 0.4144])
>>> torch.cov(x, fweights=fw, aweights=aw)
tensor([[ 0.4169, -0.4169],
        [-0.4169,  0.4169]])