Fwd: Re: <OT> unix math function: norm

[Joel Hammer]

I ran the data both on Excel and with my own bash script. They get close
results.
I have attached my data.txt and a better ps file, with the labels better
spaced out.

Joel

from the On Tue, May 28, 2002 at 12:07:53AM +0200, Klaus-Peter Schrage wrote:
> linux-sxs.org seems to be unreachable, so I'll send this one directly to you.
> Klaus
> 
> ----------  Weitergeleitete Nachricht  ----------
> 
> Subject: Re: <OT> unix math function: norm
> Date: Mon, 27 May 2002 19:06:26 +0200
> From: Klaus-Peter Schrage <kpschrage at gmx.de>
> To: linux-users at linux-sxs.org
> 
> Am Montag, 27. Mai 2002 00:07 schrieb Joel Hammer:
> 
> Ok, I can see now what you want - thought about replying off list, but I have
> seen more devious discussions on this list -;)
> 
> Although I don't have the data underlying your bar graph, one can see by mere
> visual inspection of your plot (gauss.ps) that the variance (.475 ...) resp.
> standard deviation (.689..) are not properly calculated, which yields the
> peak in your normal curve: the smaller the variance or deviation, the higher
> the peak. I guess the variance to be something around 5 or 6.
> Can you send me your raw data? I'll try to check the calculations.
> Klaus
> 
> BTW: I have checked now
> http://csep1.phy.ornl.gov/mc/node19.html
> You are definitely right, their formula is definitely wrong.
> 
> > What I would like to do is make a Gaussian normal curve that will
> > superimpose itself over bar graphs showing a population distribution.
> > The idea is to give an immediate visual impression of how far from
> > "normal" the population data is given the population mean and std dev.
> > I haven't had success with this. I can't seen to get it right.  What I
> > see is a much higher peak of my normal curve than what I see in my
> > data. I have attached a plot in fact. Here is the plot file for this.
> > (I use a big bash script to see this stuff up for gnuplot.)
> >
> > set key left Left
> > u=13.35000000000000000000
> > var=.47548245614035087719
> > display_v=.475
> > display_u=13.350
> > set label 1 "mean = 13.350 " at 15.2,.13205 right
> > set label 2 "std dev = .689 " at 15.2,.13205*.95 right
> > set label 3 "std error mean = .052 " at 15.2 ,.13205*.90 right
> > set label 4 "count = 172 " at 15.2,.13205*.85 right
> > std=.68955235924500387223
> > count=172
> > stderrormean=.0027
> > set ylabel "Result Result Frequency"
> > set xlabel " "
> > f(x)=exp(-((x-u)**2/(2*var)))/(sqrt(2*pi*var))
> > plot "/tmp/plot_data_bar" using 2:1 notitle with boxes , f(x)
> 
> -------------------------------------------------------
> 
> -- 
> Klaus-Peter Schrage
> Fridtjof-Nansen-Str. 21
> D-38108 Braunschweig
-------------- next part --------------
A non-text attachment was scrubbed...
Name: gauss.ps
Type: application/postscript
Size: 14190 bytes
Desc: not available
Url : http://mail.linux-sxs.org/pipermail/linux-users/attachments/20020527/92144af1/gauss.ps
-------------- next part --------------
11.2
11.3
11.4
11.7
12.0
12.2
12.2
12.3
12.3
12.3
12.4
12.4
12.4
12.4
12.4
12.4
12.5
12.5
12.5
12.5
12.5
12.5
12.5
12.5
12.6
12.6
12.6
12.6
12.6
12.7
12.7
12.7
12.7
12.7
12.7
12.7
12.7
12.8
12.8
12.9
12.9
12.9
12.9
12.9
12.9
12.9
12.9
12.9
12.9
12.9
13.0
13.0
13.0
13.0
13.0
13.0
13.1
13.1
13.1
13.1
13.1
13.1
13.2
13.2
13.2
13.2
13.2
13.2
13.2
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.3
13.4
13.4
13.4
13.4
13.4
13.4
13.4
13.4
13.4
13.5
13.5
13.5
13.5
13.5
13.5
13.5
13.5
13.5
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.6
13.7
13.7
13.7
13.7
13.7
13.7
13.7
13.7
13.8
13.8
13.8
13.8
13.8
13.8
13.8
13.8
13.8
13.8
13.8
13.8
13.8
13.9
13.9
13.9
13.9
13.9
13.9
14.0
14.0
14.1
14.1
14.1
14.1
14.2
14.2
14.2
14.2
14.3
14.3
14.4
14.4
14.4
14.5
14.5
14.5
14.5
14.5
14.6
14.6
14.6
14.6
14.7
14.7
14.9
15.3