擅长:python、mysql、java
<p>您可以使用<a href="http://docs.scipy.org/doc/scipy/reference/generated/scipy.stats.binned_statistic.html" rel="nofollow">^{<cd1>}</a>轻松完成。在</p>
<pre><code>import scipy.stats
import numpy
import matplotlib.pyplot as plt
x = numpy.random.rand(10000)
y = x + scipy.stats.norm(0, 0.2).rvs(10000)
means_result = scipy.stats.binned_statistic(x, [y, y**2], bins=50, range=(0,1), statistic='mean')
means, means2 = means_result.statistic
standard_deviations = numpy.sqrt(means2 - means**2)
bin_edges = means_result.bin_edges
bin_centers = (bin_edges[:-1] + bin_edges[1:])/2.
plt.errorbar(x=bin_centers, y=means, yerr=standard_deviations, linestyle='none', marker='.')
</code></pre>