Quel est le CDF à deux échantillons de

9

J'essaie de comprendre comment obtenir des valeurs de pour le test unilatéral de Kolmogorov-Smirnov , et j'ai du mal à trouver des CDF pour et dans le cas à deux échantillons. Ce qui suit est cité à quelques endroits comme le CDF pour dans un cas à un échantillon:pDn1,n2+Dn1,n2Dn+

pn+(x)=P(Dn+x|H0)=xj=0n(1x)(nj)(jn+x)j1(1xjn)nj

De plus, whuber sez il y a une formulation légèrement différente de ce CDF à un échantillon (je remplace x pour t dans sa citation pour la cohérence avec ma notation ici):

En utilisant la transformée intégrale de probabilité, Donald Knuth dérive leur distribution (commune) sur p. 57 et exercice 17 du TAoCP Volume 2. Je cite:

(Dn+xn)=xnnckx(nk)(kx)k(x+nk)nk1

Cela s'appliquerait aux hypothèses unilatérales dans le cas d'un échantillon, telles que: H 0F(x)F00 , où F(x) est le CDF empirique de x , et F0 est un CDF.

Je pense que le x dans ce cas est la valeur de Dn+ dans son échantillon, et que n(1x) est le plus grand entier de nnx . (Est-ce correct?)

Mais quel est le CDF pour (ou ) quand on a deux échantillons? Par exemple, lorsque H pour les CDF empiriques de et ? Comment obtenir ?Dn1,n2+Dn1,n20FA(x)FB(x)0ABpn1,n2+

Alexis
la source
1
Tout comme un pointeur pour quiconque cherche à répondre à cette question - ma réponse à la question précédente d'Alexis (qui est liée à la question ci-dessus) a des liens vers plusieurs références avec une discussion de l'histoire, chacune avec un certain nombre de références pertinentes. Vous pouvez vérifier ces documents et leur liste de références.
Glen_b -Reinstate Monica
@Glen_b Merci! J'apprécie vraiment votre excellente réponse à mon autre question, et j'ai suivi les ressources citées, mais je n'ai eu aucune traction sur le CDF pour là-bas, et plutôt que d'enliser les commentaires, je pensais que j'ouvrirais simplement une nouvelle requête . Des références supplémentaires sont les bienvenues, si vous en connaissez qui fonctionneront pour cela. D+
Alexis
Alexis: aucune critique n'a été voulue par mon commentaire; votre choix d'ouvrir une nouvelle question était tout à fait juste (à mon avis). Je voulais juste épargner aux gens un peu de travail pour retrouver certaines des références pertinentes - je me suis dit que tout le monde ne pourrait pas nécessairement suivre votre lien vers l'autre question, et cela pourrait ne pas arriver aux personnes qui ont fait ces liens dans mon la réponse contenait des références dont ils voudraient peut-être avoir connaissance.
Glen_b -Reinstate Monica

Réponses:

6

Ok, je vais essayer ça. Les informations critiques sont les bienvenues.

À la page 192, Gibbons et Chakraborti (1992), citant Hodges, 1958, commencent par un CDF à petit échantillon (exact?) Pour le test bilatéral (j'échange leur notation et pour et , respectivement):m,ndn1,n2x

P(Dn1,n2x)=1P(Dn1,n2x)=1A(n1,n2)(n1+n2n1)

Où est produit par une énumération de chemins (augmentant de façon monotone dans et ) de l'origine au point travers un graphique avec — en remplaçant par —les valeurs de l' axe x et de l' axe y sont et . Les chemins doivent en outre obéir à la contrainte de rester à l'intérieur des limites (où est la valeur de la statistique de test de Kolmogorov-Smirnov):A(n1,n2)n1n2(n1,n2)Sm(x)Fn1(x)n1F1(x)n2F2(x)x

n2n1±(n1+n2)x(n1+n2n1)

Ci-dessous, leur image Figure 3.2 fournit un exemple pour , avec 12 de ces chemins:A(3,4)

Figure 3.2 de la page 193 de Gibbons et Chakraborti (1992) Inférence statistique non paramétrique.

Gibbons et Chakaborti continuent en disant que la valeur unilatérale est obtenue en utilisant cette même méthode graphique, mais avec seulement la borne inférieure pour , et seulement la partie supérieure pour .pDn1,n2+Dn1,n2

Ces approches à petit échantillon impliquent des algorithmes d'énumération de chemin et / ou des relations de récurrence, ce qui rend sans aucun doute souhaitable des calculs asymptotiques. Gibbons et Chakraborti notent également les CDF limites lorsque et approchent de l'infini, de :n1n2Dn1,n2

limn1,n2P(n1n2n1+n2Dn1,n2x)=12i=1(1)i1e2i2x2

Et ils donnent le CDF limite de (ou ) comme:Dn1,n2+Dn1,n2

limn1,n2P(n1n2n1+n2Dn1,n2+x)=1e2x2

Parce que et sont strictement non négatifs, le CDF ne peut prendre que des valeurs non nulles sur :D+D[0,)

CDF de $ D ^ {+} $ (ou $ D ^ {-} $)


Références
Gibbons, JD et Chakraborti, S. (1992). Inférence statistique non paramétrique . Marcel Decker, Inc., 3e édition, édition révisée et augmentée.

Hodges, JL (1958). La probabilité de signification du test à deux échantillons de Smirnov. Arkiv för matematik . 3 (5): 469--486.

Alexis
la source
1
Le cdf réel existe partout, mais pour le cdf sera nul; la forme fonctionnelle que vous avez donnée ne s'applique qu'à (cela se prête à un raisonnement simple; qu'est-ce que ?(,0)x0P(D+<0)
Glen_b -Reinstate Monica