Talk:Generalized normal distribution

From Wikipedia, the free encyclopedia
Jump to: navigation, search
WikiProject Statistics (Rated Start-class, Low-importance)
WikiProject icon

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-Class article 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.
WikiProject Mathematics (Rated Start-class, Low-importance)
WikiProject Mathematics
This article is within the scope of WikiProject Mathematics, a collaborative effort to improve the coverage of Mathematics on Wikipedia. If you would like to participate, please visit the project page, where you can join the discussion and see a list of open tasks.
Mathematics rating:
Start Class
Low Importance
 Field:  Probability and statistics

Fat Error[edit]

The PDF formula presented for the generalized normal distribution version 1, aka exponential power distribution aka generalized error distribution, is NOT a generalized NORMAL distribution, but most probably a generalized Laplace distribution. The PDF formula presented cannot reproduce the unmodified normal distribution! alpha=1 leads to sd= 0.707 (Normal distribution: sd=1) and alpha= sqrt(2) leads to density(0)= 0.2829 (Normal Distribution: density(0)= 0.3989). In case of doubt, I can prove the R code to reproduce my calculations. — Preceding unsigned comment added by Consuli74 (talkcontribs) 09:15, 29 June 2016 (UTC)

Merger proposal[edit]

This is the same distribution as Generalized normal distribution. While it shouldn't be too hard to merge the two, the main question is which name to use? Any comments? -3mta3 (talk) 12:46, 2 March 2009 (UTC)

One thing that needs thinking about is how to handle other distributions which are also called "generalised normal". For example, documentaion for R includes a generalised normal distribution which is not the one referred to here (see eg. [1] ). Nether of these distributions is what is called either the skew-generalised normal or skew-normal distribution. Also I note that the ISI glossary has "Kapteyn's univariate distribution" as an alternative name for a generalised normal (not clear which) [2], but I guess this name should be avoided. Melcombe (talk) 14:31, 2 March 2009 (UTC)

I agree. I added this page recently because I didn't find this distribution (under either name) on the List of probability distributions page. Since "Gaussian distribution" redirects to "Normal distribution", I propose that we merge these two pages under "generalized normal distribution" with a redirect from "generalized Gaussian distribution." Then we can add a comment in the text about the "Kapteyn" name. I don't know what to do about the other generalized normal distribution. Maybe we could have them on the same page with two copies of the probability distribution template. Skbkekas (talk) 16:39, 2 March 2009 (UTC)

On further investigation, the generalized normal distribution referred to in the R documentation cited above ([3]) does not appear to include the normal distribution as a special case (also, the literature reference in the R code is to a 1990 paper of Hosking's, but this paper does not discuss anything like the distribution in the R code, and I didn't find any use of the term "generalized normal" in Hosking's other papers on JSTOR). In that sense, it is "generalized" in the same way that the lognormal, inverse normal, and half-normal distributions are (i.e. derived from the normal via a transformation). I think we can clarify that "generalized here" means a parameteric family that includes the normal distribution as a special case. This includes skew-normal, which already has a page that we can link to. Skbkekas (talk) 19:22, 2 March 2009 (UTC)

The above was moved from Talk:Generalized Gaussian distribution. Melcombe (talk) 10:34, 3 March 2009 (UTC)

Actually the generalized normal of Hosking does include the normal distribution as a special case, as well as the three-parameter log-normal of both positive and negative skewness. In one sense it might really be only be a form of reparameterisation of the three-parameter lognormal distribution, but the point here is that it appears in the literature under the name "generalized normal". Melcombe (talk) 10:46, 3 March 2009 (UTC)
Given that Hosking's version is in most cases just a log-normal, it might be better overall to have a separate article for that parameterisation but calling it something like "three-parameter lognormal" or "shifted lognormal", or even "generalised log-normal". I think putting several diffeent families of distributions in a single article is too confusing. The diagrams are an excellent contribution though. Melcombe (talk) 09:45, 5 March 2009 (UTC)

Melcombe is correct that these are both generalizations of the normal distribution in the same sense. I have included a discussion of both of them on this page, calling them "version 1" and "version 2" for lack of better terms. I don't have any references for "version 2" except the R documentation, so if someone could add the appropriate reference to Hosking's book or paper that would be great. Skbkekas (talk) 15:37, 5 March 2009 (UTC)

I have added a reference for this. However, it seems that by the time of that book they had decided to changec the name to just "lognormal distribution" and it appears only under that name in the book. However, it is certainly the same distribution as Hosking originally called the "generalized normal" . The earliest ref is "Hosking, J. R. M.: 1986, The theory of probability weighted moments. IBM Research Report, RC12210" but this is not readiliy accessible. The "generalized normal" terminology has been used by others based on this earlier report: for example, from 1998, , and, from 2007, . Melcombe (talk) 12:34, 6 March 2009 (UTC)

Multivariate Version[edit]

We need a multivariate generalization for this function. Currently, only the univariate version is given (where x and alpha are scalars). I was almost able to figure out the multivariate version (where x is a vector and alpha is a matrix), but I couldn't figure out the scale factor (to ensure unit variance).

Almon.David.Ing (talk) 17:22, 17 June 2009 (UTC)

Version 1 Questions[edit]

Does the claim of continuous derivatives under "Parameter Estimation" refer to only to the \beta? The LaPlace distribution (\beta=1) has no derivative in \mu at zero, which contradicts the claim of floor(\beta)=1 continuous derivatives in the text. Similarly, I think that the loglikelihood is infinitely smooth is \alpha. Some clarification would be nice. Also, I don't understand the CDF plot; since all of the exponential power densities as illustrated have \mu=0, the CDF of the \beta=0.5 case should be 0.5 when x=0 (like all of the others)...? Actually, the CDF of the \beta=1 case is wrong too; there must be a bug in the integrating function used for the plots. I looked at the python code and couldn't find it, unless scipy's gammainc function is buggy which would be odd. (talk) 04:36, 19 March 2012 (UTC)

The CDF of version 1 seems to have a value other than zero at x = -\inf which is odd. It works out correct if the \Gamma(1/\beta) is not present in the denominator of the second term. — Preceding unsigned comment added by (talk) 06:13, 28 March 2012 (UTC)

ah, okay. python must use a normalized incomplete gamma. i would redo the graphs but i have a mac and right now i don't feel like jumping through the hoops to install numpy and scipy.

The CDF plot is indeed wrong (as pointed out, since the PDFs are even CDFs should all equal 0.5 at x = 0); it seems like it just shows the integral of the PDFs as shown in the figure above, i.e. on the interval [-3, 3]. Would anyone who's good with Python be able to fix this? — Preceding unsigned comment added by (talk) 08:19, 25 July 2013 (UTC)

Fixed the CDFs. The plot is visually accurate when simply re-running the original code in a modern python setup. jugander (t) 00:44, 28 April 2015 (UTC)

Assessment comment[edit]

The comment(s) below were originally left at Talk:Generalized normal distribution/Comments, and are posted here for posterity. Following several discussions in past years, these subpages are now deprecated. The comments may be irrelevant or outdated; if so, please feel free to remove this section.

This article addresses the univariate function (where x and alpha are univariate). It could be expanded to address the multivariate function (where x is a vector and alpha is expressed as a matrix).

I was almost able to derive this, but I am not sure about the scale factor (to ensure unit area).

Almon.David.Ing (talk) 17:12, 17 June 2009 (UTC)

Last edited at 17:12, 17 June 2009 (UTC). Substituted at 03:11, 3 May 2016 (UTC)

External links modified[edit]

Hello fellow Wikipedians,

I have just modified one external link on Generalized normal distribution. Please take a moment to review my edit. If you have any questions, or need the bot to ignore the links, or the page altogether, please visit this simple FaQ for additional information. I made the following changes:

When you have finished reviewing my changes, you may follow the instructions on the template below to fix any issues with the URLs.

You may set the |checked=, on this template, to true or failed to let other editors know you reviewed the change. If you find any errors, please use the tools below to fix them or call an editor by setting |needhelp= to your help request.

  • If you have discovered URLs which were erroneously considered dead by the bot, you can report them with this tool.
  • If you found an error with any archives or the URLs themselves, you can fix them with this tool.

If you are unable to use these tools, you may set |needhelp=<your help request> on this template to request help from an experienced user. Please include details about your problem, to help other editors.

Cheers.—InternetArchiveBot (Report bug) 04:27, 9 January 2017 (UTC)

Version 2 issues[edit]

I've been trying to use the PDF that's listed for version 2, but it doesn't seem to work. In order to better understand the PDF I tried to follow the links (11 in particular since it's online, but it's broken). I was able to Google the documentation in question, but it does not give details on what the actual function in R does (or what formula it's following). The details listed for 'version' 2 are pretty sparse and I've spent well over an hour at this point trying to make the function work as listed (and to track down more details on it, short of trying to get the book listed through an interlibrary loan).

The PDF listed suggests that it is proportional to Φ(y)/x, where Φ(y)= 1/sqrt(2) exp(-y^2). However, it then suggests that y (when κ≠0) is -1/κ...where the negative sign would be counteracted by the square...?

The log function (I assume base 10 log?) leads to imaginary values in the exponential, which leads to large oscillations/asymptotes in the final plot.

Anyway, I'm not sure what the right answer is, but Version 2 section could use a bit of clarification. This seemed like the best place to bring the subject up. — Preceding unsigned comment added by Plasma geek (talkcontribs) 16:36, 29 June 2017 (UTC)