Как рассчитать значение p, если у меня есть t-статистика и d.f. (в Perl)? - PullRequest
4 голосов
/ 05 мая 2009

Я написал Perl-скрипт, который выполняет множество t-тестов с одним образцом. Я получаю тысячи t-статистики с их степенями свободы (df). Мне нужно обновить скрипт, чтобы он также возвращал их p-значения (их слишком много, чтобы искать их вручную в таблице). Есть ли какая-то формула, которую я могу использовать для этого с t-статистикой и d.f в качестве входных данных? Я надеюсь, что кто-то может помочь мне с этим, большое спасибо заранее!

* 1003 A.A. *

Ответы [ 3 ]

7 голосов
/ 05 мая 2009

Использование Статистика :: Распределения кажется довольно простым:

print Statistics::Distributions::tprob($dof,$tstat);
6 голосов
/ 05 мая 2009

Поиск в MetaCPAN обнаруживает следующее:

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

3 голосов
/ 05 мая 2009

Если вы проводите двухсторонний тест, тогда ваше значение p = 2 * P (T> t), где t - это ваша расчетная статистика теста. По сути, вам нужен способ смоделировать T-dist для интеграции (T-dist от t до INFINITY). Вот демо: http://www.stat.tamu.edu/~west/applets/tdemo.html

Я не знаком с Perl и его библиотеками, но, надеюсь, это поможет вам начать. Вы можете написать элементарный интегратор и проверить некоторые значения, чтобы убедиться, что он достаточно точен.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...