help for qinvgamma 

Quantile-quantile plot for data versus fitted inverse gamma distribution

qinvgamma varname [if exp] [in range] [weight] [, grid param(# #) show(condition) generate(newvar) graph_options ]

where graph_options are

rlopts(cline_options) plot(plot) (Stata 8 only) addplot(plot) (Stata 9 only) scatter_options twoway_options

fweights and aweights are allowed; see help weights.


qinvgamma plots the quantiles of varname against the quantiles of a two-parameter inverse gamma distribution with probability density function for variable x > 0, shape parameter a > 0 and scale parameter b > 0 of (b^a / Gamma(a)) x^(-a - 1) exp(-b / x). By default, maximum likelihood estimation is carried out, using invgammafit, which should be installed separately.


In the majority of cases, qinvgamma will be used to fit an inverse gamma distribution on the fly and to assess that fit. In some cases, it may be of interest to compare data with an inverse gamma distribution with known or hypothesised parameters, which may be specified using the param() option. In all cases it is important not only to specify any if or in restrictions, but also to specify relevant weights, which will (unless constant) affect the configuration of the plot.

Stata programmers will find appended to the code various Mata functions for the inverse gamma that may be useful in other circumstances.


grid adds grid lines at the .05, .10, .25, .50, .75, .90, and .95 quantiles.

param() may be used to supply parameter values (namely, shape and scale parameters) directly for use in estimation of fitted quantiles. a and b should be provided as separate values in precisely that order.

show() may be used to specify that you wish to restrict the graph according to some condition, say looking at one tail of the distribution only. Note that if and in should not be used for this purpose.

generate() specifies the name of a new variable to hold values of fitted inverse gamma quantiles.

graph_options are rlopts(cline_options), addplot(plot) (Stata 9 only), plot(plot) (Stata 8 only), scatter_options, and twoway_options.

rlopts(cline_options) affect the rendition of the reference line; see help cline_options.

addplot(addplot) provides a way to add other plots to the generated graph; see help addplot_option.

plot(plot) provides a way to add other plots to the generated graph; see help plot_option.

scatter_options affect the rendition of the plotted points; see help scatter.

twoway_options are any of the options documented in help twoway_options excluding by(). These include options for titling the graph (see help title_options) and options for saving the graph to disk (see help saving_option).


. qinvgamma spacing

. qinvgamma spacing, param(20 3)

. qinvgamma spacing, gen(spacing_g)


Nicholas J. Cox, Durham University, U.K. n.j.cox@durham.ac.uk


Evans, M., Hastings, N. and Peacock, B. 2000. Statistical distributions. New York: John Wiley.

Kleiber, C. and Kotz, S. 2003. Statistical size distributions in economics and actuarial sciences. Hoboken, NJ: John Wiley.

Also see

Online: help for invgammafit, pinvgamma (if installed), graph, diagplots