Ordonner les statistiques (par exemple, au minimum) d'une collection infinie de variables khi-deux?

C'est ma première fois ici, alors faites-moi savoir si je peux clarifier ma question de quelque manière que ce soit (y compris le formatage, les balises, etc.). (Et j'espère pouvoir éditer plus tard!) J'ai essayé de trouver des références, et j'ai essayé de me résoudre en utilisant l'induction, mais j'ai échoué dans les deux cas.

J'essaie de simplifier une distribution qui semble se réduire à une statistique d'ordre d'un ensemble infiniment dénombrable de variables aléatoires indépendantes avec différents degrés de liberté; en particulier, quelle est la distribution de la e plus petite valeur parmi les ? $\chi^2$ $m$ $\chi^2_2,\chi^2_4,\chi^2_6,\chi^2_8,\ldots$

Je serais intéressé par le cas particulier $m=1$ : quelle est la distribution du minimum de (indépendant) $\chi^2_2,\chi^2_4,\chi^2_6,\ldots$ ?

Pour le cas du minimum, j'ai pu écrire la fonction de distribution cumulative (CDF) comme un produit infini, mais je ne peux pas le simplifier davantage. J'ai utilisé le fait que le CDF de $\chi^2_{2m}$ est

F_{2 m} (x) = γ (m, x / 2) / Γ (m) = γ (m, x / 2) / (m - 1)! = 1 - e^{- x / 2} \sum_{k = 0}^{m - 1} x^{k} / (2^{k} k!) .

$F_{2m}(x)=\gamma(m,x/2)/\Gamma(m)=\gamma(m,x/2)/(m-1)!=1-e^{-x/2}\sum_{k=0}^{m-1}x^k/(2^k k!).$ (Avec

m = 1

$m=1$ , cela confirme le deuxième commentaire ci-dessous sur l'équivalence avec une distribution exponentielle avec l'espérance 2.) Le CDF du minimum peut alors être écrit comme

F_{m i n} (x) = 1 - (1 - F_{2} (x)) (1 - F_{4} (x)) \dots = 1 - \prod_{m = 1}^{\infty} (1 - F_{2 m} (x))

$F_{min}(x) = 1-(1-F_2(x))(1-F_4(x))\ldots = 1-\prod_{m=1}^\infty (1-F_{2m}(x))$

= 1 - \prod_{m = 1}^{\infty} (e^{- x / 2} \sum_{k = 0}^{m - 1} \frac{x^{k}}{2^{k} k!}) .

$= 1- \prod_{m=1}^\infty \left(e^{-x/2}\sum_{k=0}^{m-1}\frac{x^k}{2^k k!}\right).$ Le premier terme du produit est simplement

e^{- x / 2}

$e^{-x/2}$ , et le "dernier" terme est

e^{- x / 2} \sum_{k = 0}^{\infty} x^{k} / (2^{k} k!) = 1

$e^{-x/2}\sum_{k=0}^\infty x^k/(2^k k!)=1$ . Mais je ne sais pas (si possible?) Le simplifier à partir de là. Ou peut-être qu'une approche totalement différente est meilleure.

Un autre rappel potentiellement utile: $\chi^2_2$ est identique à une distribution exponentielle avec l'espérance 2, et $\chi^2_4$ est la somme de deux de ces exponentielles, etc.

Si quelqu'un est curieux, j'essaie de simplifier le théorème 1 dans cet article pour le cas de la régression sur une constante ( $x_i=1$ pour tout $i$ ). (J'ai $\chi^2$ au lieu des distributions $\Gamma$ car j'ai multiplié par $2\kappa$ .)

distributions chi-squared exponential order-statistics minimum David M Kaplan
la source

Est- ce que cela répond à votre question?

mpiktas

@mpiktas: merci pour la suggestion. C'est similaire, sauf qu'au lieu d'exponentielles avec des paramètres de taux différents, j'ai des chi-carrés avec différents degrés de liberté (et un nombre infini d'entre eux, pas fini). Et tandis que est une exponentielle, ne le sont pas; ce sont des sommes d'exponentielles, mais les sommes d'exponentielles ne sont pas elles-mêmes exponentielles. (Et idéalement, j'espère une statistique d'ordre général, bien que le min soit un bon début.)

χ_{2}^{2}

$\chi^2_2$

χ_{4}^{2}, χ_{6}^{2}, \dots

$\chi^2_4,\chi^2_6,\ldots$

David M Kaplan

Je doute qu'il y ait un formulaire fermé pour cela. Il a cependant une curieuse caractérisation: lorsque sont iid Poisson ( ) varie, , alors est la chance que tous les .

X_{k}

$X_k$

λ / 2

$\lambda/2$

k = 1, 2, \dots

$k=1,2,\ldots$

1 - F_{m i n} (λ)

$1-F_{min}(\lambda)$

X_{k} \leq k

$X_k \le k$

whuber

@whuber: Ce n'est peut-être pas tout à fait aussi curieux en termes de processus de Poisson, qui est la formulation avec laquelle je jouais. Soit des variables aléatoires iid , avec le processus de Poisson correspondant de taux . Soit , , , etc. Ensuite, les sont indépendants et par la propriété stationnaire des incréments indépendants d'un processus de Poisson, avoir .

T_{1}, T_{2}, \dots

$T_1, T_2, \ldots$

E x p (1 / 2)

$\mathrm{Exp}(1/2)$

N (t) := sup {n : \sum_{i = 1}^{n} T_{i} \leq t}

$N(t) := \sup\{n: \sum_{i=1}^n T_i \leq t\}$

1 / 2

$1/2$

U_{1} = T_{1}

$U_1 = T_1$

U_{2} = T_{2} + T_{3}

$U_2 = T_2 + T_3$

U_{3} = T_{4} + T_{5} + T_{6}

$U_3 = T_4 + T_5 + T_6$

U_{i} \sim χ_{2 i}^{2}

$U_i\sim\chi_{2i}^2$

P (U_{i} \geq t) = P (N (t) \leq i)

$\mathbb{P}(U_i \geq t) = \mathbb{P}( N(t) \leq i)$

Cardinal

@Cardinal Bien sûr: c'est une bonne façon de le voir. La curiosité n'est pas dans la relation entre Poissons et Gammas; cela réside dans la description de l'événement lui-même!

whuber

Réponses:

Les zéros du produit infini seront l'union des zéros des termes. Le calcul jusqu'au 20e trimestre montre la tendance générale:

tracé de zéros complexes

Ce tracé des zéros dans le plan complexe distingue les contributions des termes individuels dans le produit au moyen de symboles différents: à chaque étape, les courbes apparentes sont étendues davantage et une nouvelle courbe est lancée encore plus à gauche.

La complexité de cette image démontre qu'il n'existe pas de solution de forme fermée en termes de fonctions bien connues de l'analyse supérieure (telles que les gammas, les thétas, les fonctions hypergéométriques, etc., ainsi que les fonctions élémentaires, comme enquêté dans un texte classique comme Whittaker & Watson ).

Ainsi, le problème pourrait être posé de manière plus fructueuse un peu différemment : que devez-vous savoir sur les distributions des statistiques de commande? Des estimations de leurs fonctions caractéristiques? Des moments de faible ordre? Approximations aux quantiles? Autre chose?

whuber
la source

Pourquoi les zéros du produit sont-ils importants? Je sens que je manque quelque chose de trivial.

mpiktas

@mp Les zéros et les pôles montrent quelque chose sur la complexité de la fonction. Les fonctions rationnelles en ont un nombre fini. Les fonctions élémentaires ont généralement une ligne de zéros, comme à , intégrale, pour ; les fonctions "transcendantales" typiques ont des motifs de zéros légèrement plus complexes, comme par exemple tous les entiers non positifs (réciproques de la fonction Gamma) ou sur un réseau de points (fonctions thêta et fonctions elliptiques). Le schéma compliqué présenté ici suggère qu'il sera difficile, voire impossible, d'exprimer le CDF en termes de ces fonctions familières.

2 i π n

$2i\pi n$

n

$n$

\exp ()

$\exp()$

whuber

@whuber (1/2), merci! Je ne connaissais pas les différentes classes de fonctions ayant ces différents motifs de zéros dans le plan complexe; cela semble très utile, et votre graphique semble répondre à ma question (telle que posée).

David M Kaplan

@whuber (2/2), il s'agissait de vérifier un cas particulier de distribution (compliquée) d'un estimateur donné dans un autre article. Ils ont utilisé l'existence de la distribution pour justifier l'utilisation du bootstrap; mon conseiller m'a suggéré d'essayer d'approcher la distribution. Il semble que leur distribution pourrait être désactivée pour ce cas spécial (où je sais ce que cela devrait être), donc je vérifierai avec mon conseiller après sa date limite de subvention; mais potentiellement, je serais en train de prendre une expansion d'ordre supérieur de la ième ordre stat (divisé par ) que , dans un cadre plus compliqué. Publiera à nouveau si oui; Merci encore!

m

$m$

m

$m$

m \to \infty

$m\to\infty$

David M Kaplan

quelle est la distribution du minimum de (indépendant) ? $\chi^2_2,\chi^2_4,\chi^2_6,\ldots$

Toutes mes excuses pour votre arrivée avec 6 ans de retard. Même si le PO est probablement passé à d'autres problèmes, la question reste fraîche et j'ai pensé que je pourrais suggérer une approche différente.

On nous donne où où avec le pdf : $(X_1, X_2, X_3, \dots)$ $X_i \sim \text{Chisquared}(v_i)$ $v_i= 2i$ $f_i(x_i)$

Voici un tracé du du pdf correspondant , à mesure que la taille de l'échantillon augmente, pour : $f_i(x_i)$ $i = 1 \text{ to } 8$

Nous nous intéressons à la distribution de . $\text{min}(X_1, X_2, X_3, \dots)$

Chaque fois que nous ajoutons un terme supplémentaire, le pdf du dernier terme marginal ajouté se déplace de plus en plus vers la droite, de sorte que l'effet de l'ajout de termes de plus en plus devient non seulement de moins en moins pertinent, mais après quelques termes seulement , devient presque négligeable - sur l'échantillon minimum. Cela signifie, en effet, que seul un très petit nombre de termes est susceptible d'avoir une importance réelle ... et l'ajout de termes supplémentaires (ou la présence d'un nombre infini de termes) est largement hors de propos pour le problème minimum de l'échantillon.

Tester

Pour tester cela, j'ai calculé le pdf de à 1 terme, 2 termes, 3 termes, 4 termes, 5 termes, 6 termes, 7 termes, 8 termes, à 9 termes et à 10 termes. Pour ce faire, j'ai utilisé la fonction de mathStatica , l'instruisant ici pour calculer le pdf de l'échantillon minimum (la statistique d'ordre ) dans un échantillon de taille , et où le paramètre (à la place d'être fixé) est : $\text{min}(X_1, X_2, X_3, \dots)$ OrderStatNonIdentical $1^{\text{st}}$ $j$ $i$ $v_i$

Cela devient un peu compliqué à mesure que le nombre de termes augmente ... mais j'ai montré la sortie pour 1 terme (1ère ligne), 2 termes (deuxième ligne), 3 termes (3ème ligne) et 4 termes ci-dessus.

Le diagramme suivant compare le pdf de l'échantillon minimum à 1 terme (bleu), 2 termes (orange), 3 termes et 10 termes (rouge). Notez la similitude des résultats avec seulement 3 termes contre 10 termes:

Le diagramme suivant compare 5 termes (bleu) et 10 termes (orange) - les graphiques sont tellement similaires qu'ils s'oblitèrent et on ne peut même pas voir la différence:

En d'autres termes, l'augmentation du nombre de termes de 5 à 10 n'a pratiquement aucun impact visuel perceptible sur la distribution de l'échantillon minimum.

Approximation semi-logistique

Enfin, une excellente approximation simple du pdf de l'échantillon min est la distribution semi-logistique avec pdf:

g (x) = \frac{2 e^{- x}}{{(e^{- x} + 1)}^{2}} for x > 0

$g(x) = \frac{2 e^{-x}}{\left(e^{-x}+1\right)^2} \quad \text{ for } x>0$

Le diagramme suivant compare la solution exacte à 10 termes (qui ne se distingue pas de 5 ou 20 termes) et l'approximation semi-logistique (en pointillés):

L'augmentation à 20 termes ne fait aucune différence perceptible.

Wolfies
la source