¿Cómo calculo un valor p si tengo el estadístico t y d.f. (en Perl)?

Question

He escrito un script en Perl que realiza muchas pruebas t de una muestra. Obtengo miles de estadísticas con sus grados de libertad (df). Necesito actualizar el script para devolver también sus valores p (hay demasiados para buscarlos manualmente en una tabla). ¿Hay algún tipo de fórmula que pueda usar para esto con el estadístico t y d.f como entrada? Espero que alguien pueda ayudarme con esto, muchas gracias de antemano!

A.A.

Answer 1

Usar Statistics :: Distributions parece bastante sencillo:

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

Answer 2

Una búsqueda de MetaCPAN revela lo siguiente:

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

Answer 3

Si está realizando una prueba de dos colas, entonces su valor p = 2 * P (T > t), donde t es su estadística de prueba calculada. Básicamente, necesita una forma de modelar el T-dist para poder integrarlo (T-dist desde t hasta INFINITY). Aquí hay una demostración: http://www.stat.tamu.edu/ ~ west / applets / tdemo.html

No estoy familiarizado con Perl y sus bibliotecas, pero espero que esto te ayude a comenzar. Puede escribir un integrador rudimentario y verificar algunos valores para asegurarse de que sea lo suficientemente preciso.