<P>For a real data column X value&nbsp;ranged between (56.4521,32317.9) with missing values, I need to overlay 2 plots: histogram &amp; fitted gamma density.</P>
<P>I use following to generate histogram.</P>
<P>xbk_seq(50,33000,by=100)</P>
<P>hist(x,breaks=xbk)</P>
<P>But I don't know how to get "fitted gamma density"?</P>
<P>In SAS proc capability,&nbsp;I got Shape=2.59, Scale=3481).</P>
<P>But when I do</P>
<P>plot(dgamma(x, shape=2.59, rate=1,scale=3481)), it gives me a flat densed line close to 0, which does not look correct.</P>
<P>Would anybody help me on this?</P>
<P>Thanks!</P>
<P>Shelton</P><p><br><hr size=1>Do you Yahoo!?<br>
<a href="http://faith.yahoo.com">Faith Hill</a> - Exclusive Performances, Videos, & more<br>
<a href="http://faith.yahoo.com">faith.yahoo.com</a>