小编sbm*_*sbm的帖子

我有一个名为data1的数据数组,其中包含从0到超过一千的值.我只想要一个直方图和一个从0到10的值的KDE.因此我写道:

sns.distplot(data1, kde=True, hist=True, hist_kws={"range": [0,10]})
plt.show()

然而,我得到的是所有值的直方图(直到2000年代).

问题 我想通过我的主管从未发表的论文中实现算法,作为其中的一部分,我需要使用本文中给出的一些规则来构建协方差矩阵C. 我来自Matlab并希望借此机会最终学习Python,因此我的问题是:我如何以最有效(快速)的方式在Python(包括numpy,scipy)中做到这一点？

子问题1:

• 选项1:我使用2 for循环,循环遍历所有行和所有列.我认为这是最糟糕的事情.
• 选项2:使用列表推导,我构建了一个欧几里德对列表,然后迭代该列表.这就是我现在正在做的事情.

有没有更好的方法？

子问题2

• 选项1:我迭代矩阵中的所有元素.
• 选项2:我只在下三角部分(无对角线)上迭代,然后添加转置(因为协方差矩阵是对称的),然后添加对角线.

我相信子问题1是不费脑筋但我不知道子问题2.我可能也应该说我正在处理的矩阵可能是2*10 ^ 4 x 2*10 ^ 4.

谢谢!

编辑 我不想给出实际的协方差矩阵,但由于人们想要一个例子,假设我们想要构造一个称为"布朗桥"的随机过程的协方差矩阵.它的结构由下式给出:

cov(Xs,Xt)= min {s,t} - st

因为我们说s,t∈{1,...,100}.你会如何建造它？

设定：

• Python 3.5.3 | Continuum Analytics，Inc. | （默认值，2017年3月6日，12:15:08）
• Mac OSX 10.13.1

问题：

我已经下载了以下R脚本https://github.com/daleroberts/heston/blob/master/heston.r ，我使用RPy2软件包通过Python反复调用了其中的一个函数。现在，对于我输入到R函数中的某些输入，R返回以下错误：

rpy2.rinterface.RRuntimeError：积分错误（PIntegrand，Lower = 0，upper = Inf，lambda，vbar，eta，：检测到舍入错误

如何在Python中捕获此RuntimeError？