如何在Python中计算派生量中的错误

Question

我有两个量x和y以及它们的协方差矩阵cov(x,y)，我想计算导出量z=1/(x-y)中的误差。是否有一个包来计算z和sigma(xz)和sigma(zy)的平均值

事先非常感谢

Answer 1

假定y和x服从正态分布，导出的量一般不会服从正态分布。您可以用数字来估计误差分布：

import scipy.stats
import numpy as np
import matplotlib.pyplot as plt

# sample from multivariate normal
x = scipy.stats.multivariate_normal(mean=[0,1], cov=[[1,0.5],[0.5,2]]).rvs(10000)
z = 1 / (x[0] - x[1])
print(z.mean())
print(z.var())
# Create kde of the distribution
kde = scipy.stats.gaussian_kde(z)
grid = np.linspace(z.min(), z.max(), 1000)
plt.plot(grid, kde.evaluate(grid))
plt.show()