Não conheço nenhuma distribuição multimodal.
Por que todas as distribuições conhecidas são unimodais? Existe alguma distribuição "famosa" que tenha mais de um modo?
Obviamente, as misturas de distribuições geralmente são multimodais, mas eu gostaria de saber se existem distribuições "sem mistura" que tenham mais de um modo.
distributions
mode
Miroslav Sabo
fonte
fonte
Respostas:
A primeira parte da pergunta é respondida em comentários à pergunta: muita "de marca" distribuições são multimodal, tal como qualquer Beta de distribuição com um < 1 e b < 1(a,b) a<1 b<1 . Vamos passar para a segunda parte da pergunta.
Todas as distribuições discretas são claramente misturas (de átomos, que são unimodais).
Mostrarei que a maioria das distribuições contínuas também são misturas de distribuições unimodais. A intuição por trás disso é simples: podemos "lixar" os solavancos de um gráfico irregular de um PDF, um por um, até que o gráfico esteja na horizontal. Os solavancos se tornam os componentes da mistura, cada um dos quais é obviamente unimodal.
Consequentemente, exceto talvez por algumas distribuições incomuns cujos PDFs são altamente descontínuos, a resposta para a pergunta é "nenhuma": todas as distribuições multimodais absolutamente contínuas, discretas ou uma combinação dessas duas são misturas de distribuições unimodais.
Considere distribuições contínuas cujos PDFs fF f sejam contínuos (essas são as distribuições "absolutamente contínuas"). (A continuidade não é muito uma limitação; pode ser relaxada ainda mais com uma análise mais cuidadosa, assumindo apenas que os pontos de descontinuidade são discretos.)
Para lidar com "platôs" de valores constantes que podem ocorrer, defina um "modo" como um intervalo (que pode ser um ponto único em que x l = x u ), de modo quem=[xl,xu] xl=xu
tem um valor constante em m , digamos y .f m, y
não é constante em nenhum intervalo que contenha estritamente m .f m
Existe um número positivo tal que o valor máximo de f atingido em [ x l - ϵ , x u + ϵ ] é igual a y .ϵ f [xl−ϵ,xu+ϵ] y
Seja qualquer modo de f . Como f é contínuo, existem intervalos [ x ′ l , x ′ u ] contendo m para os quais f não diminui em [ x ′ l , x l ] (que é um intervalo adequado, não apenas um ponto) e não aumenta em [ x u , x ′ u ]m=[xl,xu] f f [x′l,x′u] m f [x′l,xl] [xu,x′u] (que também é um intervalo adequado). Vamos ser o infinimum de todos esses valores e x ' u o supremo de todos esses valores.x′l x′u
Essa construção definiu uma "corcunda" no gráfico de estende de x ′ l a x ′ u . Vamos y ser o maior de f ( x ' G ) e f ( x ' u ) . Por construção, o conjunto de pontos x em [ x ′ l , x ′ u ] para o qual f ( x ) ≥ y é um intervalo adequado m ′f x′l x′u y f(x′l) f(x′u) x [x′l,x′u] f(x)≥y m′ contendo estritamente (porque contém todo o conjunto de [ x ′ l , x l ] ou [ x u , x ′ u ] ).m [x′l,xl] [xu,x′u]
Nesta ilustração de um PDF multimodal, um modo é identificado por um ponto vermelho no eixo horizontal. A extensão horizontal da parte vermelha do preenchimento é o intervalo m ′ : é a base da corcunda determinada pelo modo m . A base dessa corcunda está na altura y = 0,16 . O PDF original é a soma do preenchimento vermelho e do preenchimento azul. Observe que o preenchimento azul possui apenas um modo próximo a 2 ; o modo original em [ 0 , 0 ] foi removido.m=[0,0] m′ m y≈0.16 2 [0,0]
escrevendo para o comprimento de m ′ , defina|m′| m′
e
quando e f m ( x ) = 0, caso contrário. (Isso torna f m uma função contínua, aliás.) O numerador é a quantidade pela qual f se eleva acima de ye o denominador p m é a área entre o gráfico de f e y . Assim, f m não é negativo e tem a área total 1 : é o PDF de uma distribuição de probabilidade. Por construção, tem um modo único m .x∈m′ fm(x)=0 fm f y pm f y fm 1 m
Também por construção, a função
é um PDF fornecido . (Obviamente, se p m = 1 , não resta nada de f , que deve ter sido unimodal para começar.) Além disso, ele não possui modos no intervalo m ′ (onde é constante, e é por isso que a cuidadosa definição anterior de necessário como intervalo). Além disso,pm<1 pm=1 f, m′
is a mixture of the unimodal PDFfm and the PDF f′m .
of unimodal distributions, QED.
fonte
By unimodal, I think the OP plainly means that there is just one interior mode (i.e. excluding corner solutions). The question is thus really asking ...
i.e. why do most brand name distributions look something like this:
... plus or minus some skewness or some discontinuities? When the question is posed thus, the Beta distribution would not be a valid counter example.
It appears the OP's conjecture has some validity: most common brand name distributions do not allow for more than one interior mode. There may be theoretical reasons for this. For example, any distribution that is a member of the Pearson family (which includes the Beta) will necessarily be (interior) unimodal, as a consequence of the parent differential eqn that defines the entire family. And the Pearson family nests most of the best-known brand names.
Nevertheless, here are some brand name counter examples ...
Counter example
One brand-name counter-example is theSinc2 distribution with pdf:
defined on the real line. Here is a plot of theSinc2 pdf:
We could also perhaps add the family of cardiod and distributions related to this class ... with pdf plots such as:
The family of reflected brand name distributions would also perhaps be possible brand name contenders (though, these might be thought of as a 'cheat solution' ... but they are still brand names) such as the Reflected Weibull shown here:
fonte
That you mightn't think of any doesn't mean there aren't any.
I can name "known" distributions that aren't unimodal.
For example, a Beta distribution withα and β both <1 .
http://en.wikipedia.org/wiki/Beta_distribution
also see
http://en.wikipedia.org/wiki/U-quadratic_distribution
(This isn't a special case of the beta distribution, in spite of the comment that says it is. The two families have some overlap, however.)
Mixture distributions are certainly known, and many of those are multimodal.
fonte
Alpha-skew-normal distribution (Elal-Olivero 2010) has a PDF:
whereφ is the PDF of a standard Gaussian.
For|α|>1.34 the distribution is bimodal. Examplary plot for μ=1,σ=0.5,a=2 :
fonte