Comment obtenir plus de chiffres décimaux de la valeur p du khi-deux dans la table bidirectionnelle

Question

J'utilise la commande cchi2 pour afficher la contribution de chaque cellule au khi-deux de Pearson dans une table à deux entrées dans Stata. La valeur p de chaque cellule est affichée sous la forme d'une décimale, par exemple, la valeur p est 0,0, mais j'aimerais voir plus de chiffres, par ex. 0,052 ou 0,050.

Est-il possible de définir les chiffres de la valeur décimale p?

Answer 1

Ce manque un bon exemple reproductible avec un appel de données, nous pouvons comprendre et donne en fait pas de code que ce soit. Il semble également confus en termes de statistiques et de statistiques.

Il existe une option cchi2 pour la commande tabulate lorsqu'elle est utilisée avec deux variables. cchi2 est pas une commande séparée; il donne la contribution au chi-carré et a plus de sens lorsqu'il est combiné avec l'option chi2, par ex.

. sysuse auto, clear 
(1978 Automobile Data) 

. tab foreign rep78, chi2 cchi2 

+-------------------+ 
| Key    | 
|-------------------| 
|  frequency  | 
| chi2 contribution | 
+-------------------+ 

      |     Repair Record 1978 
    Car type |   1   2   3   4   5 |  Total 
-----------+-------------------------------------------------------+---------- 
    Domestic |   2   8   27   9   2 |  48 
      |  0.3  1.1  1.8  1.0  4.2 |  8.3 
-----------+-------------------------------------------------------+---------- 
    Foreign |   0   0   3   9   9 |  21 
      |  0.6  2.4  4.1  2.3  9.5 |  19.0 
-----------+-------------------------------------------------------+---------- 
    Total |   2   8   30   18   11 |  69 
      |  0.9  3.5  5.9  3.3  13.7 |  27.3 

      Pearson chi2(4) = 27.2640 Pr = 0.000 

Ici, nous pouvons voir la contribution au chi-carré; dans ce cas, 9.5 de la statistique 27.264 chi-carré provient de la cellule en bas à droite. Ce ne sont pas des valeurs P individuelles; il n'y a qu'une seule valeur P, pour toute la table, donnée ici entre 0,000 et 3 dp.

Une façon d'obtenir plus de détails est avec la commande tabchi téléchargeable avec ssc install tab_chi. Ici l'option pearson donne les résidus de Pearson les plus utiles, (observés - attendus)/racine carrée de expected, qui sont les racines carrées signées de la contribution au khi-carré.

. tabchi foreign rep78, pearson 

      observed frequency 
      expected frequency 
      Pearson residual 

-------------------------------------------------- 
      |   Repair Record 1978   
Car type |  1  2  3  4  5 
----------+--------------------------------------- 
Domestic |  2  8  27  9  2 
      | 1.391 5.565 20.870 12.522 7.652 
      | 0.516 1.032 1.342 -0.995 -2.043 
      | 
    Foreign |  0  0  3  9  9 
      | 0.609 2.435 9.130 5.478 3.348 
      | -0.780 -1.560 -2.029 1.505 3.089 
-------------------------------------------------- 

4 cells with expected frequency < 5 
1 cell with expected frequency < 1 

      Pearson chi2(4) = 27.2640 Pr = 0.000 
likelihood-ratio chi2(4) = 29.9121 Pr = 0.000 

Taper return list après cette commande donne plus décimales pour la valeur P:

. ret li 

scalars: 
        r(N) = 69 
        r(r) = 2 
        r(c) = 5 
       r(chi2) = 27.26396103896104 
        r(p) = .0000175796084266 

Dans votre cas, nous pouvons utiliser vos fréquences affichées pour obtenir les résidus Pearson de la commande tabchii dans le même paquet.

. tabchii 3 10 2 \ 11 54 10, pearson 

      observed frequency 
      expected frequency 
      Pearson residual 

---------------------------------- 
      |   col   
     row |  1  2  3 
----------+----------------------- 
     1 |  3  10  2 
      | 2.333 10.667 2.000 
      | 0.436 -0.204 0.000 
      | 
     2 |  11  54  10 
      | 11.667 53.333 10.000 
      | -0.195 0.091 0.000 
---------------------------------- 

2 cells with expected frequency < 5 

     Pearson chi2(2) = 0.2786 Pr = 0.870 
likelihood-ratio chi2(2) = 0.2643 Pr = 0.876 

Dans votre cas, comme la statistique totale du chi carré se trouve être inférieur à 1, puis toutes les contributions, les places de (observée - attendue)/prévu, sont également tous à moins de 1, mais (un) qui n'est pas vrai en général (b) ils ne sont pas des valeurs P (il n'y a pas de test séparé effectué dans chaque cellule).

Les deux tabchi et tabchii n'ont également cont options équivalentes à l'option cchi2 de tabulate. Vous pouvez également régler le nombre de décimales affichées à l'aide des options tabdisp, comme indiqué.

Answer 2

Il existe un moyen facile de le faire sans utiliser d'autre logiciel que votre navigateur. Cette calculatrice statistique en ligne peut fournir une valeur de p pour un maximum de 12 chiffres après la virgule: https://www.icalcu.com/stat/chisqtest.html

Pour la première question, il suffit de coller les chiffres ci-dessous:

3 10 2 
11 54 10 

et vous obtiendrez une valeur p de 0.869979427395

Pour le deuxième jeu de données, il suffit de coller les chiffres ci-dessous:

2  8  27  9  2 
0  0  3  9  9 

et vous obtiendrez une valeur p de ,000017579608.