Quel est le minimum de sur toutes les distributions unimodales continues sur un intervalle borné ?

8

Toutes les distributions sur un intervalle borné satisfont:[0,1]

σ2μ(1μ)

où est la moyenne et la variance.μσ2

Supposons maintenant que la distribution soit unimodale, en ce sens qu'elle a au plus un maximum local. Quelle est la valeur minimale que le rapport suivant peut avoir:

μ(1μ)σ2?
Becko
la source
... votre première équation implique que le rapport ne peut pas être inférieur à 1. Demandez-vous quelle distribution la rend égale à 1?
user603
Jetez un oeil à un Bernoulli avec . Il est assez courant que les solutions à ces types de problèmes extrêmes soient discrètes et ne portent que sur quelques points. Vous semblez avoir fait plusieurs publications de type "bookwork". Est-ce que ce travail est pour un sujet? (p)μ=p
Glen_b -Reinstate Monica
@Glen_b La question demande cependant une distribution unimodale, ce qui n'est pas le cas d'une version tachée d'un Bernoulli.
Dougal
La distribution uniforme sur donne une valeur de 3. Les distributions bêta donnent et ne sont unimodales que si , , donc aussi 3 (quand il est également uniforme). J'ai essayé plusieurs autres familles de distribution nommées ( d'ici ) et n'ai jamais obtenu une valeur meilleure que 3. J'ai également commencé à l'écrire comme un problème d'optimisation en faisant une interpolation linéaire entre les points, mais cela ressemblait à un problème d'optimisation difficile, et je me suis arrêté avant en fait le codage et l'essayer. [0,1]α+β+1α>1β>1
Dougal
1
Interrogé simultanément sur math.SE où il a déjà reçu deux réponses (dont une a été supprimée par l'auteur de la réponse à cause de l'impolitesse perçue de l'OP).
Dilip Sarwate

Réponses:

6

Un minimum n'existe pas. Cependant, un infimum fait. Il résulte du fait que

Le supremum de la variance des distributions unimodales définies sur [0,1] avoir méchant μ est μ(23μ)/3 (0μ1/2) ou (1μ)(3μ1)/3 (1/2μ1).

Le supremum est en fait atteint par une distribution qui - bien qu'elle n'ait pas de fonction de densité - peut encore (dans un sens généralisé) être considérée comme "unimodale"; il aura un atome0 (quand μ<1/2) ou un atome 1 (quand μ>1/2) mais sinon être uniforme.


Je vais esquisser l'argument. La question nous demande d'optimiser une fonction linéaire

Lx2:D[0,1]R

soumis à diverses contraintes d’égalité et d’inégalité, où D[0,1] est l'ensemble des mesures (signées) sur l'intervalle [0,1]. Pour différenciableF:[0,1]R et g:[0,1]R toute fonction continue, définir

Lg[F]=01g(x)dF(x),

et étendre L à tous D[0,1] par continuité.

Les contraintes d'égalité sont

L1[F]=1

et

Lx[F]=μ.

Les contraintes d'inégalité sont que

f(x)0

et il existe λ[0,1] (un "mode") tel que pour tous 0xyλ et tout λyx1,

f(x)f(y).

Ces contraintes déterminent un domaine convexeXD[0,1] par-dessus lequel Lx2 doit être optimisé.

Comme pour tout programme linéaire dans un espace de dimension finie, les extrêmes de Lg sera atteint aux sommets de X. Ce sont évidemment les mesures, absolument continues par rapport à la mesure de Lebesgue, qui sont constantes par morceaux , parce que les sommets sont là où presque toutes les inégalités deviennent des égalités: et la plupart de ces inégalités sont associées à l'unimodalité deF (comportement de queue non croissant).

Afin de satisfaire les deux contraintes d'égalité, nous devons faire une seule coupure dans le graphique de f, disons à un certain nombre 0<λ<1. Laisser la valeur constante sur l'intervalle[0,λ) être a et la valeur constante sur (λ,1] être b, un calcul simple basé sur les contraintes d'égalité donne

a=1+λ2μλ, b=2μλ1λ.

Figure 1: Tracé d'un $ f _ {(\ lambda, \ mu)} $ typique.

Cette figure dit tout: elle représente graphiquement la fonction de distribution localement constante de la moyenneμ avec au plus une seule pause à λ. (L'intrigue def(λ,μ) pour μ>1/2 ressemble à l'inversion de celui-ci.)

La valeur de Lx2 à de telles mesures (que je désignerai f(λ,μ), la densité d'une distribution F(λ,μ)) est tout aussi facilement calculé

Lx2[f(λ,μ)]=13(2μ+(2μ1)λ).

Cette expression est linéaire dansλ, ce qui implique qu'il est maximisé à 0 (quand μ<1/2), 1 (quand μ>1/2), ou à n'importe quelle valeur (lorsque μ=1/2). Cependant, sauf lorsqueμ=1/2, les valeurs limites des mesures f(λ,μ) ne sont plus continus: la distribution correspondante F=limλ0F(λ,μ) ou F=limλ1F(λ,μ) a une discontinuité de saut à 0 ou 1 (mais pas les deux).

Figure 2: Graphique de $ F $ optimal pour $ \ mu = 2/5 $.

Cette figure représente la valeur optimale F pour une moyenne de μ2/5.

Quoi qu'il en soit, la valeur optimale est

σμ2=supλLx2[f(λ,μ)]=13μ(23μ).

Par conséquent, l'infimum de μ(1μ)/σ2 pour 0μ<1/2 est

μ(1μ)/σμ2=33μ23μ,

avec une expression comparable lorsque 1/2<μ1 (obtenu en remplaçant μ par 1μ).

Figure 3: Graphique de l'infimum en fonction de $ \ mu $.

Ce chiffre trace le supremum μ(1μ)/σμ2 contre μ.

whuber
la source
1
Je pense que c'est une belle réponse. Est-ce basé sur un manuel ou sur du papier? Existe-t-il une référence avec plus de résultats comme celui-ci?
Becko
1
@becko Merci. J'aimerais pouvoir aider, mais c'est une solution originale. Je ne sais pas par où commencer pour chercher d'autres résultats, car je ne suis pas spécialiste des inégalités de répartition.
whuber