mvtcdf

Многомерная кумулятивная функция распределения t

Синтаксис

y = mvtcdf(X,C,DF)
y = mvtcdf(xl,xu,C,DF)
[y,err] = mvtcdf(...)
[...] = mvntdf(...,options)

Описание

y = mvtcdf(X,C,DF) возвращает интегральную вероятность многомерного распределения t параметрами корреляции C и степени свободы DF, оцененный в каждой строке X. Строки n-by-d матричный X соответствуйте наблюдениям или точкам, и столбцы соответствуют переменным или координатам. y n- 1 вектор.

C симметричное, определенное положительное, d-by-d матрица, обычно корреляционная матрица. Если его диагональные элементы не 1, mvtcdf шкалы C к форме корреляции. mvtcdf не повторно масштабирует X. DF скаляр или вектор с элементами n.

Многомерная интегральная вероятность t в X задан как вероятность что случайный векторный T, распределенный как многомерный t, будет находиться в пределах полубесконечного прямоугольника с верхними пределами, заданными X, т.е. Pr{T(1)X(1),T(2)X(2),...T(d)X(d)}.

y = mvtcdf(xl,xu,C,DF) возвращает многомерную интегральную вероятность t, оцененную по прямоугольнику с нижними и верхними пределами, заданными xl и xu, соответственно.

[y,err] = mvtcdf(...) возвращает оценку ошибки в y. Для двумерных и trivariate распределений, mvtcdf использует адаптивную квадратуру на преобразовании плотности t, на основе методов, разработанных Genz, как описано в ссылках. Допуском абсолютной погрешности по умолчанию к этим случаям является 1e-8. Для четырех или больше размерностей, mvtcdf использует алгоритм интегрирования квази-Монте-Карло на основе методов, разработанных Genz и Bretz, как описано в ссылках. Допуском абсолютной погрешности по умолчанию к этим случаям является 1e-4.

[...] = mvntdf(...,options) указывает, что параметры управления для численного интегрирования использовались для расчета y. Этот аргумент может быть создан вызовом statset. Выбор statset параметры:

  • 'TolFun' — Максимальный допуск абсолютной погрешности. Значением по умолчанию является 1e-8 когда d <4, или 1e-4 когда d ≥ 4.

  • 'MaxFunEvals' — Максимальное количество оценок подынтегрального выражения позволило когда d ≥ 4. Значением по умолчанию является 1e7. 'MaxFunEvals' проигнорирован когда d <4.

  • 'Display' — Level of display выводится. Выбором является 'off' (значение по умолчанию), 'iter', и 'final'Отображение проигнорирован когда d <4.

Примеры

свернуть все

Вычислите cdf многомерного t распределения параметрами корреляции C = [1 .4; .4 1] и 2 степени свободы.

C = [1 .4; .4 1];
df = 2;
[X1,X2] = meshgrid(linspace(-2,2,25)',linspace(-2,2,25)');
X = [X1(:) X2(:)];
p = mvtcdf(X,C,df);

Постройте cdf.

figure;
surf(X1,X2,reshape(p,25,25));

Ссылки

[1] Genz, A. “Численный расчет Прямоугольных Двумерных и Нормальных и t Вероятностей Trivariate”. Статистика и Вычисление. Издание 14, № 3, 2004, стр 251–260.

[2] Genz, A. и Ф. Брец. “Численный расчет Многомерных t Вероятностей с Приложением к Расчету мощности Нескольких Контрастов”. Журнал Статистического Расчета и Симуляции. Издание 63, 1999, стр 361–378.

[3] Genz, A. и Ф. Брец. “Сравнение Методов для Расчета Многомерных t Вероятностей”. Журнал Вычислительной и Графической Статистики. Издание 11, № 4, 2002, стр 950–971.

Введен в R2006a

Для просмотра документации необходимо авторизоваться на сайте