如何在测量误差的数据样条上使用scipy quad时传播错误？

Question

我有一个带有N个点的数据集，我使用scipy.integrate.quad来拟合样条并进行整合。我想使用N个相关的测量误差来对最终的积分值进行误差估计。

我最初尝试使用不确定性包但x +/- stddev对象不适用于scipy。

List<String> mutable = new ArrayList<>(Arrays.asList("a", "b", "c"));
ImmutableList<String> immutable = Lists.immutable.withAll(mutable);
List<String> mutableAgain = immutable.toList();

有任何想法吗？

Answer 1

合成数据集。您的数据点有错误。现在生成1000个数据集，每个点从正常分布中绘制，该分布以测量点为中心，此时由errror给出标准偏差。适合每个数据集。整合。重复。现在你有1000个积分值。计算这些值的均值和标准值。