Como faço para calcular um valor p se eu tiver a estatística t e d.f. (Em Perl)?

Question

Eu escrevi um script Perl que executa muitos testes t de uma amostragem. Eu recebo milhares de estatísticas t com seus graus de liberdade (df). Eu preciso atualizar o script para também retornar seus valores p (há muitos para procurá-los manualmente em uma tabela). Existe algum tipo de fórmula que posso usar para isso com a estatística t e d.f como entrada? Espero que alguém possa me ajudar com isso, muito obrigado antecipadamente!

A.A.

Answer 1

Estatísticas :: Distribuições parece bastante simples:

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

Answer 2

A busca de MetaCPAN revela o seguinte:

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

Answer 3

Se você está fazendo um teste bicaudal, então o seu valor-p = 2 * P (T> t), onde t é o seu teste estatístico calculado. Então, basicamente, você precisa encontrar uma maneira de modelar o T-dist para integrar (T-dist a partir de t para o infinito). Aqui está uma demonstração: http://www.stat.tamu.edu/ ~ oeste / applets / tdemo.html

Eu não estou familiarizado com Perl e suas bibliotecas, mas espero que isto é uma introdução. Você pode escrever um integrador rudimentar e verificar alguns valores para se certificar de que é suficientemente preciso.