[Numpy-discussion] Producing a Histogram When Bins Are Known
josef.pktd at gmail.com
josef.pktd at gmail.com
Fri Nov 27 22:01:22 EST 2009
On Fri, Nov 27, 2009 at 9:44 PM, Wayne Watson
<sierra_mtnview at sbcglobal.net> wrote:
> 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
(values is the data array, ... no idea about scientificstat.standardDeviation)
Sebastian's example is only part of a larger script that defines many
of the variables and functions that are used.
If you are not yet familiar with these examples, maybe you look at the
self contained examples in the matplotlib docs. At least that's what I
do when I only have a rough idea about what graph I want to do but
don't know how to do it with matplotlib. I usually just copy a likely
looking candidate and change it until it (almost) produces what I
want.
For example look at histogram examples in
http://matplotlib.sourceforge.net/examples/index.html
Josef
> 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
>>
>>
>>>> ax = fig.add_subplot(1,1,1)
>>>> 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
>>>>>
>>>>> ax = fig.add_subplot(1,1,1)
>>>>> 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
>>>>>
>>>>>
>>>> --
>>>> Wayne Watson (Watson Adventures, Prop., Nevada City, CA)
>>>>
>>>> (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
>>
>>
>
> --
> Wayne Watson (Watson Adventures, Prop., Nevada City, CA)
>
> (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
>
More information about the NumPy-Discussion
mailing list