首页 文章
  • 0 votes
     answers
     views

    从scipy CSR矩阵索引到numpy数组的最有效方法?

    我有一个形状为 (4000, 3) 的numpy ndarray X ,其中 X 中的每个样本都是一个3D坐标(x,y,z) . 我有一个形状 (4000, 4000) 的scipy csr矩阵 nn_rad_csr ,这是从 sklearn.neighbors.radius_neighbors_graph(X, 0.01, include_self=True) 生成的最近邻居图 . nn_ra...
  • 1 votes
     answers
     views

    限制scikit-learn或statsmodels中线性回归的调整参数的界限

    是否可以限制scikit-learn或statsmodels中线性回归的调整参数的界限,例如:在statsmodels.regression.linear_model.OLS或sklearn.linear_model.LinearRegression? http://statsmodels.sourceforge.net/devel/generated/statsmodels.regression...
  • 1 votes
     answers
     views

    使用sci-kit中的训练/测试数据学习曲线而不是交叉验证

    我将我的训练和测试数据分开(从不同的CSV加载到不同的pandas数据帧中)我想用这个训练和测试数据绘制学习曲线,而不是使用交叉验证训练和测试训练集自身生成的测试数据(似乎成为learn_curve工作的常用方式) . 似乎scikit希望您的测试和训练数据存在于同一个Dataframe中,但这样分类器也会学习测试数据,这不是我想要的 . 我怎样才能解决这个问题?我是sci-kit的新手 .
  • 6 votes
     answers
     views

    python - 如何在数据矩阵中计算nans的相关矩阵

    当数据中存在NaN时,我找不到计算包含多于两个变量的观测值的数组的相关系数矩阵的函数 . 有一些函数可以为变量对执行此操作(或者只使用~is.nan()来掩盖数组) . 但是通过循环遍历大量变量来使用这些函数,计算每对变量的相关性可能非常耗时 . 所以我自己尝试并很快意识到这样做的复杂性是协方差正确归一化的问题 . 我对你如何做的意见非常感兴趣 . 这是代码: def nancorr(X,nanf...
  • 2 votes
     answers
     views

    加速Python中的一对多关联计算

    我想计算一个向量和Python中每行数组之间的Pearson相关系数(假设为numpy和/或scipy) . 由于实际数据阵列的大小和存储器约束,将不可能使用标准相关矩阵计算功能 . 这是我天真的实现: import numpy as np import scipy.stats as sps np.random.seed(0) def correlateOneWithMany(one, man...
  • 9 votes
     answers
     views

    在Python中,如何计算两个数据阵列之间的相关性和统计显着性?

    我有两个同样长的数据数组的数据集,或者我可以创建一个两项条目数组,我想计算数据所代表的相关性和统计显着性(可能紧密相关,或者可能有无统计学意义的相关性) . 我用Python编程并安装了scipy和numpy . 我查看并找到Calculating Pearson correlation and significance in Python,但这似乎想要操纵数据,因此它落入指定范围 . 我假设,有...
  • 15 votes
     answers
     views

    Pearson相关系数双尾p值意义[关闭]

    从我使用的sciPy库:scipy.stats.stats导入pearsonr来计算两个数组的相关系数,我得到一个值:(0.80751532276005755,0.19248467723994242) . 我以为我会在-1到1的范围内得到一个值,所以我不确定如何解释这两个结果 . 这是我的两个数组: x = [50,500,1500,2500] y = [17, 6, 6, 194] 我做了: ...
  • 0 votes
     answers
     views

    皮尔逊相关零值

    我正在计算Pearson Correlation . 最后我得到了如下结果(correlation1) . 我想知道为什么我的所有第二个系数都有0.0作为相关的结果1 . 有谁能解释一下?而且,我的相关代码运行缓慢 . 我怎么能快点? 结果(样本):(0.52543523179249552, 0.0), (0.52543905756911169, 0.0), (0.5254419657220660...
  • 1 votes
     answers
     views

    Python:分类数据的排序顺序相关性

    我对编程和统计有点新,所以如果它正式不正确,请帮助我改进这个问题 . 我有很多参数和一些我在MonteCarlo模拟中生成的结果向量 . 现在我想测试每个参数对结果的影响 . 我已经有了一个与Kendall的Tau合作的剧本 . 现在我想与Spearman和Pearson rho进行比较 . 一个例子: from scipy.stats import spearmanr, kendalltau, ...
  • 89 votes
     answers
     views

    Numpy:用向量元素划分每一行

    假设我有一个numpy数组: data = np.array([[1,1,1],[2,2,2],[3,3,3]]) 我有一个相应的“矢量:” vector = np.array([1,2,3]) 如何对每行 data 进行减法或除法,结果如下: sub_result = [[0,0,0], [0,0,0], [0,0,0]] div_result = [[1,1,1], [1,1,1], [1...
  • 1 votes
     answers
     views

    具有快速访问的稀疏矩阵

    在使用大型SciPy CSR稀疏矩阵时,我注意到切割矩阵从矩阵中获取单行非常慢,因为它似乎可以复制 . 有没有办法制作一个稀疏矩阵,它接受现有行的引用而不是复制它,也许有一个比CSR矩阵更合适的实现? 我的实现需要的是快速查找元素和行以及快速查找向量的所有非零索引 . 我永远不需要以任何方式更改矩阵或对矩阵执行其他操作 .
  • 2 votes
     answers
     views

    返回scipy稀疏矩阵的对角元素

    给定 scipy.sparse.csr.csr_matrix ,是否有快速返回对角线元素的方法? 我想这样做的原因是计算 inv(D) A ,其中 D 是对角线矩阵,其对角线条目与 A 一致( A 是我的稀疏矩阵,保证在对角线上具有非零值) .
  • 3 votes
     answers
     views

    找到具有numpy / sympy的线性系统的整数解(diophantine)

    我需要用numpy或sympy来解决一个线性不定方程组的系统 . 有没有办法约束numpy的linalg.solve / linalg.lstsq方法只返回整数解? (可能不是,但我想应该问) 我查看了Sympy的丢番图解算器,它似乎不适用于解决整个系统 我正在努力解决的问题是: P1(X) + P2(Y) = TargetPro F1(X) + F2(Y) = TargetFat C1(X) +...
  • 5 votes
     answers
     views

    将满足特定条件的scipy.sparse矩阵行设置为零

    我想知道用稀疏矩阵替换不满足某个条件的行的最佳方法是什么 . 例如(我使用普通数组进行说明): 我想用一行零替换总和大于10的每一行 a = np.array([[0,0,0,1,1], [1,2,0,0,0], [6,7,4,1,0], # sum > 10 [0,1,1,0,1], ...
  • 0 votes
     answers
     views

    在python中对角元素非零的稀疏对称矩阵生成

    我必须生成给定维数n * n的稀疏对称矩阵,其中所有对角元素都是非零的 . 并且非对角线元素可以期望可能的k个非零值 . 如果k = 3,并且n = 4,则我们应该具有大小为4 * 4的对称稀疏矩阵 . 并且所有对角元素都是非零的 . 对于第一行,可以有3个非零值,所有其他值都为零 . 我怎样才能实现它?
  • 5 votes
     answers
     views

    稀疏矩阵的元素划分,忽略0/0

    我有两个稀疏矩阵E和D,它们在同一个地方有非零项 . 现在我希望将 E/D 作为稀疏矩阵,仅在D非零时定义 . 例如,请使用以下代码: import numpy as np import scipy E_full = np.matrix([[1.4536000e-02, 0.0000000e+00, 0.0000000e+00, 1.7914321e+00, 2.6854320e-01, 4.1...
  • 12 votes
     answers
     views

    将sympy表达式转换为numpy数组的函数

    我有一个同意写的ODE系统: from sympy.parsing.sympy_parser import parse_expr xs = symbols('x1 x2') ks = symbols('k1 k2') strs = ['-k1 * x1**2 + k2 * x2', 'k1 * x1**2 - k2 * x2'] syms = [parse_expr(item) for item...
  • 2 votes
     answers
     views

    在numpy或pytorch中自动获得对角矩阵条纹

    我需要获得矩阵的对角条纹(这里不确定术语,对角矩阵条纹似乎最好描述它) . 比方说,我有一个大小为KxN的矩阵,其中K和N是任意大小,K> N.说,我有一个矩阵: [[ 0 1 2] [ 3 4 5] [ 6 7 8] [ 9 10 11]] 从中我需要提取对角条纹,在这种情况下,是通过截断原始条带创建的矩阵MxV大小: [[ 0 x x] [ 3 4 x] ...
  • 0 votes
     answers
     views

    Scipy的优化与Multiprocessing不兼容?

    在尝试使用Scipy的优化算法来最小化在子过程中计算其值的函数时,我发现基于梯度的算法(到目前为止流域购物和L-BFGS-B)在optimize.py的第562行遇到以下错误: grad [k] =(f(*((xk d,)args)) - f0)/ d [k]TypeError:不支持的操作数类型 - :'NoneType'和'NoneType' 以下是生成此错误的代码的简单示例: impo...
  • 0 votes
     answers
     views

    Numpy - AttributeError:'Zero' object没有属性'exp'

    我在解决运行时断开的东西之间的差异时遇到了麻烦,但在python控制台中使用完全相同的数据和操作,让它工作正常 . # f_err - currently has value 1.11819388872025 # l_scales - currently a numpy array [1.17840183376334 1.13456764589809] sq_euc_dists = self.se...
  • 2 votes
     answers
     views

    scipy.integrate.solve_ivp不清楚如何解决0 = F(t,y(t),y'(t))形式的隐式ODE

    目前,我确实使用assimulos solver suite来求解形式为0 = F(t,y(t),y'(t))的隐式微分方程 . 我想使用一个带有原生scipy安装的解算器,并查看scipy.integrate.solve_ivp的文档 . 文档说解算器解决了以下形式的系统: 解决ODE系统的初始值问题 . 给定初始值dy / dt = f(t,y)y(t0)= y0时,该函数对常微分方程组进行...
  • 2 votes
     answers
     views

    Scipy最小化,fmin,minimalsq类型问题(设置数组元素与序列),不合适

    我在使用scipy.optimize.fmin和scipy.optimize.minimize函数时遇到了麻烦 . 我检查并确认传递给函数的所有参数都是numpy.array类型,以及错误函数的返回值 . 此外,carreau函数返回标量值 . 一些额外参数(例如大小)的原因是:我需要使用给定模型(Carreau)来拟合数据 . 数据是在不同的温度下获得的,这些温度用移位因子(也由模型拟合)校正,...
  • 3 votes
     answers
     views

    BFGS无法融合

    我正在研究的模型是一个多项logit选择模型 . 这是一个非常具体的数据集,因此其他现有的MNLogit库不适合我的数据 . 所以基本上,它是一个非常复杂的函数,它接受11个参数并返回对数似然值 . 然后我需要找到可以使用scipy.optimize.minimize最小化对数似然的最佳参数值 . 以下是我遇到的不同方法的问题: 'Nelder-Mead':效果很好,总能给我正确答案 . 但是,它...
  • 0 votes
     answers
     views

    使用args传递单个参数时,Scipy.optimize.minimize给出错误

    给出以下值 a = 100.0 b = 50.0 c = 200.0 以下代码有效(它没有做任何有趣的事) def excessRevenue1(taxH,blah1,blah2): return taxH**2 + blah1 + blah2 print optimize.minimize( excessRevenue1,c,args=(a,b) ).x 但是下面的代码产生错...
  • 26 votes
     answers
     views

    sigmoidal回归与scipy,numpy,python等

    我有两个变量(x和y)彼此之间有一些S形关系,我需要找到某种预测方程,这将使我能够在给定任何x值的情况下预测y的值 . 我的预测方程需要显示两个变量之间的某种S形关系 . 因此,我不能满足于产生线的线性回归方程 . 我需要看到两个变量图的右侧和左侧出现的斜率的逐渐曲线变化 . 我在googling曲线回归和python之后开始使用numpy.polyfit,但是如果你运行下面的代码,这给了我可怕的...
  • 4 votes
     answers
     views

    最小的置信区间适合scipy python

    如何在python中计算最小二乘拟合(scipy.optimize.leastsq)的置信区间?
  • 1 votes
     answers
     views

    Scipy LeastSq错误栏

    我使用来自SciPy的LeastSq将实验光谱拟合到理论预期 . 当然有与实验值相关的错误 . 如何将这些提供给LeastSq或者我需要一个不同的例程?我在文档中找不到任何内容 .
  • 1 votes
     answers
     views

    python scipy leastsq适合复杂的数字

    我有一个复数的数据集,我希望能够找到最适合数据的参数 . 您可以使用pyp中的scipy实现的leastsq来复杂数字的数据吗? 例如,我的代码是这样的: import cmath from scipy.optimize import leastsq def residuals(p,y,x): L,Rs,R1,C=p denominator=1+(x**2)*(C**2...
  • 2 votes
     answers
     views

    Pyinstaller与scipy.signal ImportError:DLL加载失败

    使用PyInstaller 3.3,Python 3.4我编译了一个需要使用scipy.signal.periodogram()函数的程序,它导致.exe在它启动的瞬间崩溃 . 我一直在搜索堆栈溢出和其他网站,试图包含正确的DLL,但无济于事 . 请帮我弄清楚我是否需要在我的.spec文件中包含一些特定的DLL或包含一些其他的隐藏导入! 这是我的PyInstaller命令: pyinstaller...
  • 2 votes
     answers
     views

    Python Scipy,插值数组到图像

    我的基本脚本有问题,加载图像,然后将其转换为数组,通过插值缩放,并应将缩放图像保存到文件,但最后一行不起作用 . 我不知道如何解决它,我知道它与缩放数组有关 . 我会很高兴得到任何帮助 . 脚本代码: import numpy as np import scipy.ndimage as ndimage from scipy import misc import matplotlib.pyplot ...

热门问题