Statalist The Stata Listserver


[Date Prev][Date Next][Thread Prev][Thread Next][Date index][Thread index]

st: InverseGamma with beta<0


From   "Carlo Lazzaro" <[email protected]>
To   <[email protected]>
Subject   st: InverseGamma with beta<0
Date   Fri, 15 Jun 2007 20:20:13 +0200



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/



© Copyright 1996–2024 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   What's new   |   Site index