Re: Gaussian PDF Overlay

• To: mathgroup at smc.vnet.net
• Subject: [mg19283] Re: [mg19277] Gaussian PDF Overlay
• From: BobHanlon at aol.com
• Date: Sat, 14 Aug 1999 01:45:11 -0400
• Sender: owner-wri-mathgroup at wolfram.com

```Lawrence,

Use a GeneralizedBarChart:

Needs["Statistics`ContinuousDistributions`"]
Needs["Statistics`DataManipulation`"];
Needs["Graphics`Graphics`"];

dist = NormalDistribution[mu = 4., sigma = 2.];

xmin = mu - 3sigma; xmax = mu + 3sigma;

plt1 = Plot[PDF[dist, x], {x, xmin, xmax}, PlotStyle -> AbsoluteThickness[2],
DisplayFunction -> Identity];

nbrSamples = 1000;
data = RandomArray[dist, nbrSamples];
nbrBins = 10;
step = (xmax - xmin)/nbrBins;
freq = BinCounts[data, {xmin, xmax, step}]/nbrSamples;
midPoints = Table[xmin + step(k - 1/2), {k, nbrBins}];

plt2 = GeneralizedBarChart[Transpose[{midPoints, freq, Table[1, {nbrBins}]}],
DisplayFunction -> Identity];

Show[{plt2, plt1}, DisplayFunction -> \$DisplayFunction];

Bob Hanlon

In a message dated 8/12/99 5:35:54 AM, lwalker701 at earthlink.net writes:

>I've tried to overlay a continuous gaussian probability
>density curve onto a bar chart of similar distribution.  The
>values for the barchart were generated by bin counting the
>data using the BinCounts function.  The data is a list
>generated by Random[NormalDistribution[]].
>
>At first I simply tried to overlay the plot and barchart
>using Show.  The x-axis scale of the curve was different
>from the
>barchart's scale.  So I had to manually change the position
>and rescale the curve until it overlayed the barchart using
>the program below.  I didn't like the fact that the x-scale
>shows 0 to 40 instead of -4 to 4.
>
>Is there a better way of doing this?
>

```

• Prev by Date: Re: Subscripts, Doh!!!