Dear Statalisters,
I beg Your pardon for posting a possible basic thread.
I have fitted a Gamma distribution with Stata 9/SE from a dataset containing
savings (ie: cost difference all < 0) accrued to patients who underwent
healthcare Programme A (50,000 patients) instead of healthcare Programme B
(50,000 patients)
Parameter alfa was>0 but beta was less <0.
Then I have drawn 50,000 random samples via InvGamma:
gen InvGamma=beta*invgammap(alfa, uniform())
InvGamma results confirmed that healthcare Programme A was always less
costly than Healthcare programme B.
However, being Gamma constrained on 0 to 1 interval, I really do not know
whether the results of InvGamma are or not reliable.
Any hints would be highly appreciated.
Thanks a lot,
Carlo
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
*
* For searches and help try:
* http://www.stata.com/support/faqs/res/findit.html
* http://www.stata.com/support/statalist/faq
* http://www.ats.ucla.edu/stat/stata/
