hotelling t2 matlab,Hotelling的T ^ 2在python中得分

62 阅读 0 评论 41 点赞

我是靠谱客的博主繁荣皮卡丘，最近开发中收集的这篇文章主要介绍hotelling t2 matlab,Hotelling的T ^ 2在python中得分，觉得挺不错的，现在分享给大家，希望可以做个参考。

matplotlib的

PCA类不包括Hotelling T2计算,但只需几行代码即可完成.以下代码包括计算每个点的T2值的函数. __main__脚本将PCA应用于

Matlab’s pca documentation中使用的相同示例,因此您可以验证函数是否生成与Matlab相同的值.

from __future__ import print_function, division

import numpy as np

from matplotlib.mlab import PCA

def hotelling_tsquared(pc):

"""`pc` should be the object returned by matplotlib.mlab.PCA()."""

x = pc.a.T

cov = pc.Wt.T.dot(np.diag(pc.s)).dot(pc.Wt) / (x.shape[1] - 1)

w = np.linalg.solve(cov, x)

t2 = (x * w).sum(axis=0)

return t2

if __name__ == "__main__":

hald_text = """Y X1 X2 X3 X4

78.5 7 26 6 60

74.3 1 29 15 52

104.3 11 56 8 20

87.6 11 31 8 47

95.9 7 52 6 33

109.2 11 55 9 22

102.7 3 71 17 6

72.5 1 31 22 44

93.1 2 54 18 2

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。

本图文内容来源于网友提供，作为学习参考使用，或来自网络收集整理，版权属于原作者所有。