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 <[EMAIL PROTECTED]> > To: [EMAIL PROTECTED] > > 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
gauss.ps
Description: PostScript document
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
