Comment puis-je obtenir plus de chiffres pour une valeur de p?

Quelqu'un connais un moyen facile d'obtenir Stata pour afficher plus de trois chiffres pour la p-valeur lors de l'exécution d'une régression Tobit?

Normalement Stata rapports que la valeur de p est .001 ou .065, mais j'aimerais voir plus de chiffres, par exemple, .0011123 ou .065320.

Pour être clair, je ne veux pas (nécessairement) de modifier la façon dont les données sont produites dans la régression de la table.

Je veux seulement être en mesure d'obtenir Stata pour afficher plus de chiffres pour ceux qui les valeurs de p qui m'intéressent.

  • Avez-vous vraiment besoin de plus de 3 chiffres? Je suppose que la plupart des gens seraient rapport p-valeurs à 2 chiffres pour lesquels une sortie de 3 chiffres est plus que suffisant. Pourriez-vous préciser pourquoi vous avez besoin de plus de chiffres?
  • Corrections de Bonferroni avec des modèles de grande taille (plusieurs centaines de variables) est l'une des raisons.
  • C'est un très bon point.
  • Je le soupçonne de Gunter observation a été faite dans un contexte différent, probablement celui où quelqu'un essayait de overinterpret quelque chose comme une p-valeur de 10^-100 ou à peu près. Le fait est que l'interprétation et l'utilisation relativement faibles valeurs de p (10^-4 à 10^-6) est un enjeu important: il a été l'objet de réglementations nationales et fédérales litige aux états-unis, par exemple, il n'est donc pas quelque chose à être rejeté avec une désinvolte commentaire. Je vous remercie pour la prise clair que vous offre le devis, au moins partiellement, dans la plaisanterie!
  • Mes excuses les plus sincères, le commentaire est totalement à la blague. Je l'ai enlevé car il n'appartient pas à une discussion sérieuse de ce genre. En fait je ne savais pas que la p-valeur faible est un gros problème en NOUS. Merci de me laisser savoir. Je me demande ce que gagnent-ils gain de confiance p-valeur faible en dehors de l'acquisition de ces connaissances "quelque chose ne va pas".
  • Il apparaît dans les règlements régissant la surveillance de l'environnement. Les programmes de surveillance peut entraîner dans la conduite des dizaines de milliers de tests simultanément, la création d'une nécessité d'éviter des résultats faussement positifs.
  • Merci beaucoup ! .. J'ai vraiment la valeur de la connaissance que j'ai de gain à partir de ce forum!
  • Merci à tous pour les commentaires!

InformationsquelleAutor | 2010-10-27
  1. 8

    Suivi de l'tobie commande avec

    est tab, p(%12.10g)

    (par exemple). Cela devrait fonctionner, même dans de très anciennes versions du logiciel Stata. Un peu moins facile est d'écrire votre propre sortie de la procédure.

    • Bonne idée. J'ajouterais que -est - onglet est une abréviation (compris par Stata) pour la commande intégrée -estimations de la table. Qui a été introduite avec la version de logiciel Stata 8.0 en 2003.
    • OP ici. Merci pour votre réponse, fonctionne parfaitement! 🙂
    InformationsquelleAutor whuber
  2. 6

    Stata 11.1 introduit un ensemble pformat commande qui spécifie le format de sortie de p-valeurs dans les tableaux de coefficients. (Je ne sais pas à propos de STATA je crains que je pense qui a été abandonné peu de temps dans les années 1980).

    Par le chemin, vous seriez probablement mieux de poser de telles complètement Stata questions spécifiques sur Statalist plutôt qu'ici.

    • Merci pour votre réponse! 🙂
    • Pour éviter la nécessité pour la recherche de formats - set pformat %5.4f, perm va augmenter la valeur affichée dans le tableau de sortie de 3 à 4 décimales
    InformationsquelleAutor onestop
  3. 2

    Un grand nombre de fois, vous pouvez obtenir la plus grande précision si vous connaissez votre valeur de p par son nom interne. J'ai l'habitude de type return list ou ereturn list après presque chaque commande que je vais sérieusement à utiliser, puis saisir les résultats qui peuvent ressembler à des e(p) ou r(p) ou e(p_chi2) ou quel que soit le scalaire qui contient la valeur de p peut être.

    • Chaque fois que possible, je le fais aussi. Cependant, les valeurs de p sont pas toujours des "laissés derrière" après l'estimation de commandes (ereturn list). Une belle manière de faire les p-valeurs de ces commandes est décrite ici (stata-journal.com/sjpdf.html?articlenum=st0137) et utilise une base mata fonctions. Cependant, la set pformat et set tab solutions sont beaucoup plus pratiques.
    InformationsquelleAutor StasK
  4. 1

    Après une régression tobit, vous pouvez utiliser le test commande pour obtenir la valeur de p à partir de l'hypothèse nulle x1=0: 

    sysuse auto
tobit weight trunk length headroom, ll(1500)
test trunk

    Le résultat retourné dans r(p)

    return list
    InformationsquelleAutor Keith
  5. 1

    En utilisant le premier exemple de la tobit commande fichier d'aide:

    . sysuse auto, clear
. generate wgt = weight /1000

. tobit mpg wgt, ll(17)

Tobit regression                                Number of obs     =         74
                                                LR chi2(1)        =      72.85
                                                Prob > chi2       =     0.0000
Log likelihood = -164.25438                     Pseudo R2         =     0.1815

------------------------------------------------------------------------------
         mpg |      Coef.   Std. Err.      t    P>|t|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
         wgt |   -6.87305   .7002559    -9.82   0.000    -8.268658   -5.477442
       _cons |   41.49856    2.05838    20.16   0.000     37.39621     45.6009
-------------+----------------------------------------------------------------
      /sigma |   3.845701   .3663309                      3.115605    4.575797
------------------------------------------------------------------------------
            18  left-censored observations at mpg <= 17
            56     uncensored observations
             0 right-censored observations

    Vous pouvez facilement obtenir toutes les p-valeur dans les résultats retournés dans r():

    . matrix list r(table)

r(table)[9,3]
             model:      model:      sigma:
               wgt       _cons       _cons
     b  -6.8730504   41.498557   3.8457011
    se   .70025591   2.0583803   .36633085
     t  -9.8150552   20.160782          .b
pvalue   5.610e-15   1.471e-31          .b
    ll  -8.2686584   37.396211   3.1156048
    ul  -5.4774424   45.600903   4.5757975
    df          73          73          73
  crit   1.9929971   1.9929971   1.9929971
 eform           0           0           0

    Et ensuite le format en conséquence:

    . matrix results = r(table)

. display %18.17f results[4,1]
0.00000000000000561

    Type help format de Stata l'invite de commande pour plus d'informations.

    InformationsquelleAutor Pearly Spencer