scipy相关内容

Hotelling 在 Python 中的 T^2 分数

我在 python 中使用 matplotlib 对数据集应用了 pca.但是，matplotlib 不像 Matlab 那样提供 t 平方分数.有没有办法像 Matlab 一样计算 Hotelling 的 T^2 分数? 谢谢. 解决方案 matplotlib 的 PCA 类不包括 Hotelling T2 计算，但只需几行代码即可完成.以下代码包含一个函数，用于计算每个点的 T2 ..
发布时间：2022-01-07 23:36:15 Python

用于查找对的压缩矩阵函数

对于一组观察: [a1,a2,a3,a4,a5] 它们的成对距离 d=[[0,a12,a13,a14,a15][a21,0,a23,a24,a25][a31,a32,0,a34,a35][a41,a42,a43,0,a45][a51,a52,a53,a54,0]] 以压缩矩阵形式给出(上面的上三角，从 scipy.spatial.distance.pdist 计算): c=[a12,a1 ..
发布时间：2022-01-07 23:32:03 Python

如何在python中实现概率分布的Conflation?

我在网上查找了将多个连续概率分布组合成一个连续概率分布的方法.这种方法叫做Conflation，该方法可以在下面的文章中找到:An Optimal合并不同实验数据的方法.在这篇文章中，我发现最好执行 Conflation 而不是平均来组合分布. 从我从文章中了解到的是，方程通过将来自多个概率分布的每个概率密度值相乘除以来自连续分布的多个概率分布的每个概率密度值的乘积的积分来执行，而对于离散分 ..
发布时间：2022-01-07 23:27:11 Python

离散数据拟合:负二项式、泊松、几何分布

在 scipy 中，不支持使用数据拟合离散分布.我知道有很多关于这个的主题. 例如，如果我有一个如下所示的数组: x = [2,3,4,5,6,7,0,1,1,0,1,8,10,9,1,1,1,0,0] 我无法申请此数组: from scipy.stats import nbinom参数 = nbinom.fit(x) 但是我想问你最新的，有没有什么办法适合这三个离散分布，然 ..
发布时间：2022-01-07 23:24:43 Python

如何在 ipython notebook 中显示图形

尝试在 SymPy 中进行绘图 - 根据这个视频，我写了: from sympy.plotting 导入图，plot_parametrice = sin(2*sin(x**3))plot(e, (x, 0, 5)); 但是在评估那个单元格之后我没有得到任何输出?没有错误或任何东西，它只是不显示任何东西. 另一个测试: from sympy import *从 sympy.plot ..
发布时间：2021-12-31 23:44:15 Python

在 DNN 训练结束时返回逆 Hessian 矩阵，并在输入中使用偏导数

使用 Keras 和 Tensorflow 作为后端，我构建了一个 DNN，它以恒星光谱作为输入(7213 个数据点)并输出三个恒星参数(温度、重力和金属丰度).网络在我的测试集上训练得很好，预测也很好，但为了让结果在科学上有用，我需要能够估计我的错误.这样做的第一步是获得逆 Hessian 矩阵，这似乎无法仅使用 Keras.因此，我试图用 scipy 创建一个解决方法，使用 scipy.opt ..
发布时间：2021-12-31 17:06:16 Python

交换 numpy 数组中的两个值.

是否有比以下代码更有效的方法来交换 numpy 一维数组的两个值? input_seq = arange(64)ix1 = randint(len(input_seq))ixs2 = randint(len(input_seq))temp = input_seq[ix2]输入序列[ix2] = 输入序列[ix1]input_seq[ix1] = temp 解决方案 您可以使用元组解包.元组解 ..
发布时间：2021-12-31 12:30:40 Python

在 Python 中读取 scipy/numpy 中的 csv 文件

我在 python 中读取由制表符分隔的 csv 文件时遇到问题.我使用以下功能: def csv2array(filename, skiprows=0, delimiter='\t', raw_header=False, missing=None, with_header=True):"""将文件名解析为数组.返回数组和其他标题行.默认情况下，将标题行解析为字典，假设参数是数字，使用'pars ..
发布时间：2021-12-31 12:30:30 Python

不均匀间隔点之间的直观插值

我有以下图表，我想使用 Python 和 Matplotlib 将其数字化为高质量的出版物级图表: 我使用数字化仪程序从 3 个数据集之一中抓取了一些样本: x_data = np.array([1、1.2371，1.6809,2.89151，5.13304，9.23238，])y_data = np.array([0.0688824,0.0490012,0.0332843，0.02358 ..
发布时间：2021-12-31 12:30:24 Python

带有屏蔽数据的 Scipy 插值?

我正在尝试插入一个包含屏蔽数据的二维数组.我使用了一些 SciPy 模块的可用方法，包括 interp2d、bisplrep/bisplev 以及 RectBivariateSpline.作为附加信息，我的数据是一个规则数组，这意味着网格具有相同的维度(在本例中为 1ºX1º). 话虽如此，有没有什么办法可以用Python在数组中插入避免屏蔽数据?我还是新手，使用 Python 和 NumP ..
发布时间：2021-12-31 12:30:16 Python

有效地从矩阵中减去向量 (Scipy)

我有一个存储为 scipy.sparse.csc_matrix 的大矩阵，我想从大矩阵的每一列中减去一个列向量.当您执行标准化/标准化等操作时，这是一项非常常见的任务，但我似乎无法找到有效执行此操作的正确方法. 这是一个演示示例: # mat 是一个 3x3 矩阵mat = scipy.sparse.csc_matrix([[1, 2, 3],[2, 3, 4],[3, 4, 5]])#v ..
发布时间：2021-12-31 12:30:08 Python

在 Python 中集成返回数组的函数

我有很多数据需要整合，我想找到一种只用矩阵来完成所有操作的方法，并且愿意在准确性上妥协以提高性能.我的想法是这样的: 导入numpy进口scipya = np.array([1,2,3])定义函数(x):返回 x**2 + xdef func2(x):全球 a返回 a*x定义被积函数(x):返回 func(x)*func2(x)积分 = 四边形(被积函数，0, 1) 所以我试图整合来自 int ..
发布时间：2021-12-31 12:29:58 Python

Python - 在 2D numpy 数组中查找特定值的最大区域的有效方法

我有一个 2D numpy 数组，其中有些值为零，有些值不是.我试图找到一种有效的方法来找到数组中最大的零块(通过返回零的数量，以及对中心位置的粗略了解) 例如在这个数组中，我想找到以 (3,4) 为中心的 9 块: [[ 1, 1, 1, 0, 0 ],[ 1, 0, 1, 1, 0 ],[ 1, 1, 1, 1, 1 ],[ 1, 1, 0, 0, 0 ],[ 1, 1, 0, 0, ..
发布时间：2021-12-31 12:29:49 Python

分配问题，一个 NumPy 函数?

由于 赋值问题 可以以单个矩阵的形式提出，我是想知道 NumPy 是否有一个函数来解决这样的矩阵.到目前为止，我没有找到.也许你们中的一个人知道 NumPy/SciPy 是否有分配问题解决功能? 编辑:与此同时，我在 http://software.clapper.org/munkres/.我仍然认为 NumPy/SciPy 实现会更快，对吗? 解决方案 不，NumPy 不包含这样的 ..
发布时间：2021-12-31 12:29:34 Python

替代 scipy.misc.imresize()

我想使用仍然使用 scipy.misc.imresize() 的旧脚本，该脚本不仅已弃用，而且已完全从 scipy 中删除.相反，开发人员建议使用 numpy.array(Image.fromarray(arr).resize()) 或 skimage.transform.resize(). 不再工作的确切代码行是这样的: new_image = scipy.misc.imresize(o ..
发布时间：2021-12-31 12:29:24 Python