|
Документ взят из кэша поисковой машины. Адрес
оригинального документа
: http://num-anal.srcc.msu.ru/lib_na/cat/sf_htm_p/sf56r_p.htm
Дата изменения: Fri Nov 6 13:43:42 2015 Дата индексирования: Sun Apr 10 03:30:38 2016 Кодировка: Windows-1251 |
|
Текст подпрограммы и версий sf56r_p.zip |
Тексты тестовых примеров tsf56r_p.zip |
Вычисление обратной функции t - распределения вероятностей Стьюдента.
По заданному значению вероятности P SF56R вычисляет такое значение X, что абсолютное значение случайной величины, подчиненной закону t - распределения вероятностей Стьюдента со степенью свободы DF (DF ≥ 1), больше X.
G.W.Hill, Algorithm 396, Student's t - quantilies, CACM, 13 (10), 1970.
function SF56R(P :Real; DF :Real; var IERR :Integer): Real;
Параметры
| P - | заданное значение вероятности (тип: вещественный); |
| DF - | заданное значение степени свободы t - распределения вероятностей Стьюдента (тип: вещественный); |
| IERR - | целая переменная, служащая для сообщения об ошибках, обнаруженных в ходе работы функции; при этом: |
| IERR=65 - | когда заданная степень свободы t - распределения Стьюдента меньше 1; значение функции полагается равным 3.4E38; |
| IERR=66 - | когда заданное значение вероятности не принадлежит интервалу (0, 1); значение функции полагается равным 3.4E38. |
Версии: нет
Вызываемые подпрограммы
| SF51R - | вычисление обратной функции нормального pаспределения вероятностей. |
| UTSF14 - | подпрограмма выдачи диагностических сообщений в ходе работы функции SF56R. |
Замечания по использованию: нет
Unit tsf56r_p;
interface
uses
SysUtils, Math, { Delphi }
Lstruct, Lfunc, UtRes_p, SF56R_p;
function tsf56r: String;
implementation
function tsf56r: String;
var
IERR :Integer;
DF,P,X :Real;
begin
Result := ''; { результат функции }
DF := 1.0;
P := 0.7999537599;
X := SF56R(P,DF,IERR);
Result := Result + Format('%s',[' P=']);
Result := Result + Format('%20.16f ',[P]);
Result := Result + Format('%s',[' DF']);
Result := Result + Format('%20.16f ',[DF]);
Result := Result + Format('%s',[' X=']);
Result := Result + Format('%20.16f ',[X]) + #$0D#$0A;
Result := Result + Format('%s',[' IERR=']);
Result := Result + Format('%5d ',[IERR]) + #$0D#$0A;
UtRes('tsf56r',Result); { вывод результатов в файл tsf56r.res }
exit;
end;
end.
Результаты: X = 0.325000000004 , IERR = 0