[Numpy-discussion] Producing a Histogram When Bins Are Known

Wayne Watson sierra_mtnview at sbcglobal.net
Fri Nov 27 21:44:44 EST 2009

Joseph,
That got it by the fig problem but there is yet another one. value is
not defined on the very long line:
range = ...
Wayne

josef.pktd at gmail.com wrote:
> On Fri, Nov 27, 2009 at 9:05 PM, Sebastian <sebas0 at gmail.com> wrote:
>
> ...
> you need to create a figure, before you can use it
>
> fig = pylab.figure()
>
> Josef
>
>
>>> pylab.title(r'\Large  BCG NO radio distribution $\rm{TITLE}$')
>>> n, bins, patches = pylab.hist(values, bins=math.sqrt(len(values)),
>>>
>>> range=(numpy.mean(values)-3*scientificstat.standardDeviation(values),numpy.mean(values)+3*scientificstat.standardDeviation(values)),
>>> normed=1, facecolor='y', alpha=0.5)
>>> ax.set_xlabel(r'\Large$\rm{values}$')
>>> ax.set_ylabel(r'\Large Delatavalue/Value')
>>>
>>>
>>> gausx=numpy.arange(numpy.mean(Value)-3*scientificstat.standardDeviation(Value),numpy.mean(Value)+3*scientificstat.standardDeviation(bpty_plt),0.1)
>>>
>>> gaus=normpdf(gausx,numpy.mean(Value),scientificstat.standardDeviation(Value))
>>> pylab.plot(gausx,gaus, color='red', lw=2)
>>> ax.set_xlim(-1.5, 1.5)
>>> ax.grid(True)
>>>
>>> Sebastian wrote:
>>>
>>>> Did you try using the parameter range?
>>>> I do something like this.
>>>> regards
>>>>
>>>>     pylab.title(r'\Large  BCG NO radio distribution $\rm{TITLE}$')
>>>>     n, bins, patches = pylab.hist(values, bins=math.sqrt(len(values)),
>>>>
>>>> range=(numpy.mean(values)-3*scientificstat.standardDeviation(values),numpy.mean(values)+3*scientificstat.standardDeviation(values)),
>>>>     normed=1, facecolor='y', alpha=0.5)
>>>>     ax.set_xlabel(r'\Large$\rm{values}$')
>>>>     ax.set_ylabel(r'\Large Delatavalue/Value')
>>>>
>>>>
>>>> gausx=numpy.arange(numpy.mean(Value)-3*scientificstat.standardDeviation(Value),numpy.mean(Value)+3*scientificstat.standardDeviation(bpty_plt),0.1)
>>>>
>>>> gaus=normpdf(gausx,numpy.mean(Value),scientificstat.standardDeviation(Value))
>>>>     pylab.plot(gausx,gaus, color='red', lw=2)
>>>>     ax.set_xlim(-1.5, 1.5)
>>>>     ax.grid(True)
>>>>
>>>>
>>>> On Fri, Nov 27, 2009 at 4:38 PM, Christopher Barker
>>>> <Chris.Barker at noaa.gov <mailto:Chris.Barker at noaa.gov>> wrote:
>>>>
>>>>     josef.pktd at gmail.com <mailto:josef.pktd at gmail.com> wrote:
>>>>     > On Fri, Nov 27, 2009 at 12:57 PM, Skipper Seabold
>>>>     <jsseabold at gmail.com <mailto:jsseabold at gmail.com>> wrote:
>>>>
>>>>     >>  This kind of info might be useful to other newcomers
>>>>     >> somewhere...  <http://www.scipy.org/History_of_SciPy>?  Thoughts
>>>> on
>>>>     >> posting this on the wiki here?
>>>>     >
>>>>     > I also agree. It will improve with the newly redesigned website
>>>>     for scipy.org <http://scipy.org>
>>>>     > However, I cannot find the link right now for the development
>>>>     version of
>>>>     > the new website.
>>>>
>>>>     Feel free to crib whatever you want from my post for that -- or
>>>>     suggest
>>>>     a place for me to put it, and I'll do it. I'm just not sure where it
>>>>     should go at this point.
>>>>
>>>>     -Chris
>>>>
>>>>
>>>>     --
>>>>     Christopher Barker, Ph.D.
>>>>     Oceanographer
>>>>
>>>>     Emergency Response Division
>>>>     NOAA/NOS/OR&R            (206) 526-6959   voice
>>>>     7600 Sand Point Way NE   (206) 526-6329   fax
>>>>     Seattle, WA  98115       (206) 526-6317   main reception
>>>>
>>>>     Chris.Barker at noaa.gov <mailto:Chris.Barker at noaa.gov>
>>>>     _______________________________________________
>>>>     NumPy-Discussion mailing list
>>>>     NumPy-Discussion at scipy.org <mailto:NumPy-Discussion at scipy.org>
>>>>     http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>>>
>>>>
>>>> ------------------------------------------------------------------------
>>>>
>>>> _______________________________________________
>>>> NumPy-Discussion mailing list
>>>> NumPy-Discussion at scipy.org
>>>> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>>>
>>>>
>>> --
>>>
>>>             (121.015 Deg. W, 39.262 Deg. N) GMT-8 hr std. time)
>>>              Obz Site:  39° 15' 7" N, 121° 2' 32" W, 2700 feet
>>>
>>>                   350 350 350 350 350 350 350 350 350 350
>>>                     Make the number famous. See 350.org
>>>            The major event has passed, but keep the number alive.
>>>
>>>                    Web Page: <www.speckledwithstars.net/>
>>>
>>> _______________________________________________
>>> NumPy-Discussion mailing list
>>> NumPy-Discussion at scipy.org
>>> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>>
>> _______________________________________________
>> NumPy-Discussion mailing list
>> NumPy-Discussion at scipy.org
>> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>>
>>
>>
> _______________________________________________
> NumPy-Discussion mailing list
> NumPy-Discussion at scipy.org
> http://mail.scipy.org/mailman/listinfo/numpy-discussion
>
>

--

(121.015 Deg. W, 39.262 Deg. N) GMT-8 hr std. time)
Obz Site:  39° 15' 7" N, 121° 2' 32" W, 2700 feet

350 350 350 350 350 350 350 350 350 350
Make the number famous. See 350.org
The major event has passed, but keep the number alive.

Web Page: <www.speckledwithstars.net/>