Re: mathematica fit data to inverse gamma distribution
- To: mathgroup at smc.vnet.net
- Subject: [mg117404] Re: [mg117354] mathematica fit data to inverse gamma distribution
- From: Alexei Boulbitch <alexei.boulbitch at iee.lu>
- Date: Thu, 17 Mar 2011 06:33:26 -0500 (EST)
Hi, Alex, put "ProbabilityDensity" instead of "PDF". In other words, try this: data1 = RandomVariate[InverseGammaDistribution[0.5, 3], 10^4]; Show[Histogram[data1, {0, 20, 1}, "ProbabilityDensity"], Plot[PDF[TruncatedDistribution[{0, 20}, InverseGammaDistribution[0.5, 3]], x], {x, 0, 20}, PlotStyle -> Thick]] Have fun, Alexei Hi, I am trying to obtain same figure that is provided in Mathematica's website related to fitting an inverse gamma to a dataset. Here is the code I enter: data1 = RandomVariate[InverseGammaDistribution[0.5, 3], 10^4]; Show[ Histogram[data1, {0, 20, 1}, "PDF"], Plot[PDF[ TruncatedDistribution[{0, 20}, InverseGammaDistribution[0.5, 3]], x], {x, 0, 20}, PlotStyle -> Thick]] and here is the error I receive: Histogram::hspec: The height specification PDF is not one of the named height functions or a function that can be used to compute the heights for each bin.>> Show::"gcomb" : StyleBox[Could not combine the graphics objects in ..... I would appreciate anyone's help. Regards Alex -- Alexei Boulbitch, Dr. habil. Senior Scientist Material Development IEE S.A. ZAE Weiergewan 11, rue Edmond Reuter L-5326 CONTERN Luxembourg Tel: +352 2454 2566 Fax: +352 2454 3566 Mobile: +49 (0) 151 52 40 66 44 e-mail: alexei.boulbitch at iee.lu www.iee.lu -- This e-mail may contain trade secrets or privileged, undisclosed or otherwise confidential information. If you are not the intended recipient and have received this e-mail in error, you are hereby notified that any review, copying or distribution of it is strictly prohibited. Please inform us immediately and destroy the original transmittal from your system. Thank you for your co-operation.