# Talk:Inverse Gaussian distribution

WikiProject Statistics (Rated Start-class, Low-importance)

This article is within the scope of the WikiProject Statistics, a collaborative effort to improve the coverage of statistics on Wikipedia. If you would like to participate, please visit the project page or join the discussion.

Start  This article has been rated as Start-Class on the quality scale.
Low  This article has been rated as Low-importance on the importance scale.

## graphics

it would be illustrative to have a graphical represaentation of the distribution

The last change (November 4) to the Brownian motion section was made by me... reference is "The Inverse Gaussian Distribution: Theory, methodology and applications" p29.

Deavik 23:39, 4 November 2007 (UTC)

## Parameterizations

There are several parameterizations of the inverse Gaussian distribution, including one which makes the relationship between it and a Brownian motion with drift more explicit (in terms of the drift parameter v and the variance parameter ${\displaystyle \sigma ^{2}}$). These other parameterizations should be at least mentioned. I am not aware of any good reason to present the given one as "canonical"; though if such a reason exists, it, too, should be presented. Cheers, Eliezg 05:04, 6 November 2007 (UTC)

## Error in likelihood function.

I think the likelihood function should be changed to:

${\displaystyle L(\mu ,\lambda )=\left({\frac {\lambda }{2\pi }}\right)^{\frac {n}{2}}\left(\prod _{i=1}^{n}{\frac {w_{i}}{X_{i}^{3}}}\right)^{\frac {1}{2}}\exp \left({\frac {\lambda }{\mu }}-{\frac {\lambda }{2\mu ^{2}}}\sum _{i=1}^{n}w_{i}X_{i}-{\frac {\lambda }{2}}\sum _{i=1}^{n}w_{i}{\frac {1}{X_{i}}}\right).}$

Without this additional term, you can not solve the first order condition for the MLE's, as they are given in this article.

Ryantg (talk) 23:14, 15 March 2009 (UTC)

I agree - fixed Batman50 (talk) 15:17, 4 May 2010 (UTC)

the article mentioned the name is misleading, and the inverse gaussian distribution is not the distribution of ${\displaystyle Y={\frac {1}{X}}}$ (X is normal), then what is the name of the distribution of Y???? It will be great if we say something about that. Jackzhp (talk) 14:11, 28 July 2009 (UTC)

If X is normal, the distribution of 1/X doesn't seem to have a special name, but its pdf is straightforward to work out from probability density function#Dependent variables and change of variables. I can't find any references but from a quick bit of experimental mathematics i think it doesn't have any moments if the mean of X is zero and none apart from the mean otherwise, and it's oddly-behaved around zero in general, although if the coefficient of variation of X is large you may not notice in practice. See also Box-Cox distribution for the case when X is truncated at zero. As i don't have any references, it's rather messy and doesn't seem to have much theoretical or practical importance, i'm not going to add anything to the article myself. Qwfp (talk) 18:01, 25 May 2010 (UTC)

## Error in the labels of the figure

Apart from formatting differences, a figure identical to the one shown can be generated in R using

plot(x, dinvgauss(x, 1, 1), type="l", xaxs="i", yaxs="i", xlab="", ylab="", col=1); lines(x,dinvgauss(x,1,0.2),col=2); lines(x,dinvgauss(x,1,3),col=3); lines(x,dinvgauss(x,3,1),col=4); lines(x,dinvgauss(x,3,0.2),col=5)

However, according to the R manual (e.g., [1]), the second parameter of the function dinvgauss corresponds to 1/lambda, instead of lambda. Hence, the figure's labels should be replaced. Alternatively, use the following code to generate an interesting sequence of densities:

x<-seq(0,3,0.01); plot(x, dinvgauss(x, 1, 1/4), type="l", xaxs="i", yaxs="i", xlab="", ylab="", col=1,ylim=c(0,2.5)); lines(x,dinvgauss(x, 1, 1/2), type="l", col=2); lines(x,dinvgauss(x, 1, 1), type="l", col=3); lines(x,dinvgauss(x, 1, 2), type="l", col=4); lines(x,dinvgauss(x, 1, 4), type="l", col=5); lines(x,dinvgauss(x, 1, 8), type="l", col=6); lines(x,dinvgauss(x, 1, 16), type="l", col=7); lines(x,dinvgauss(x, 1, 32), type="l", col=8)

Btw, it would be useful to add code snippets that generate the graphs. —Preceding unsigned comment added by Szepi (talkcontribs) 21:16, 6 February 2010 (UTC)

# Clarifying the initial description

I think one of the introductory paragraphs is likely to be misconstrued, specifically the clause, "while the Gaussian describes the distribution of distance at fixed time in Brownian motion . . . ." It is so commonplace to generalize Brownian motion to higher dimensions that I think it is confusing to claim the process's distance is Gaussian (since in higher dimensions, it isn't). I think this article should not assume that readers be familiar with Brownian motion. If we ditched the word "distance" and replaced it with "level," I think we would not only clear up the ambiguity, but the wording would nicely complement the following sentence: Gaussian describes level at a certain time; inverse Gaussian describes time to hit a certain level. Disadvantage: "level" isn't very precise, although my gut tells me it is fairly intuitive. LandruBek (talk) 22:07, 6 April 2010 (UTC)

# If there is no drift...

The article should, I think, also mention the special case of this distribution as μ tends to infinity. This is how I was first introduced to the IG distribution -- it describes the hitting time of a drift-less Wiener process at level λ². The derivation is an elementary one from the CDF of a Gaussian using a simple change of variables, and thus I think readers are likely to expect to see this case mentioned explicitly. LandruBek (talk) 22:07, 6 April 2010 (UTC)

## Weights?

What are the ${\displaystyle w_{i}}$ terms occurring in some of the equations? Weights?

Any information on the quantiles? --178.2.135.21 (talk) 15:39, 17 August 2013 (UTC)