eic-opticks/ana/normal.py

0001 #!/usr/bin/env python
0002 #
0003 # Copyright (c) 2019 Opticks Team. All Rights Reserved.
0004 #
0005 # This file is part of Opticks
0006 # (see https://bitbucket.org/simoncblyth/opticks).
0007 #
0008 # Licensed under the Apache License, Version 2.0 (the "License");
0009 # you may not use this file except in compliance with the License.
0010 # You may obtain a copy of the License at
0011 #
0012 #   http://www.apache.org/licenses/LICENSE-2.0
0013 #
0014 # Unless required by applicable law or agreed to in writing, software
0015 # distributed under the License is distributed on an "AS IS" BASIS,
0016 # WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
0017 # See the License for the specific language governing permissions and
0018 # limitations under the License.
0019 #
0020
0021
0022
0023 import numpy as np
0024
0025 # Sample from a normal distribution using numpy's random number generator
0026 samples = np.random.normal(size=10000)
0027
0028 # Compute a histogram of the sample
0029 bins = np.linspace(-5, 5, 30)
0030 histogram, bins = np.histogram(samples, bins=bins, normed=True)
0031
0032 bin_centers = 0.5*(bins[1:] + bins[:-1])
0033
0034 # Compute the PDF on the bin centers from scipy distribution object
0035 #from scipy import stats
0036 #pdf = stats.norm.pdf(bin_centers)
0037
0038 from matplotlib import pyplot as plt
0039 plt.figure(figsize=(6, 4))
0040 plt.plot(bin_centers, histogram, label="Histogram of samples")
0041 #plt.plot(bin_centers, pdf, label="PDF")
0042 plt.legend()
0043 plt.show()
0044
0045