Comment calculer une p-valeur si j'ai les statistiques t et d.f. (en Perl)?

Question

J'ai écrit un script Perl qui effectue de nombreux tests t à un échantillon. Je reçois des milliers de statistiques t avec leurs degrés de liberté (df). J'ai besoin de mettre à jour le script pour qu'il retourne aussi ses p-valeurs (il y en a trop pour les chercher manuellement dans un tableau). Y at-il une sorte de formule que je peux utiliser pour cela avec la statistique t et d.f en entrée? J'espère que quelqu'un pourra m'aider pour cela, merci d'avance!

A.A.

Answer 1

L'utilisation de Statistics :: Distributions semble assez simple:

print Statistics::Distributions::tprob($dof,$tstat);

Answer 2

Une recherche sur MetaCPAN révèle les éléments suivants:

https://metacpan.org/pod/Statistics::Distributions

Answer 3

Si vous effectuez un test bilatéral, votre valeur p = 2 * P (T > t), où t correspond à votre statistique de test calculée. Donc, essentiellement, vous avez besoin d’un moyen de modéliser le T-dist afin de l’intégrer (T-dist de t à INFINITY). Voici une démo: http://www.stat.tamu.edu/ ~ west / applets / tdemo.html

Je ne connais pas bien Perl ni ses bibliothèques, mais j'espère que cela vous aidera à démarrer. Vous pouvez écrire un intégrateur rudimentaire et vérifier certaines valeurs pour vous assurer qu’elles sont suffisamment précises.