R: каков вектор квантилей в функции плотности dvmnorm

library(mvtnorm)
dmvnorm(x, mean = rep(0, p), sigma = diag(p), log = FALSE)

dmvnorm обеспечивает функцию плотности для многомерного нормального распределения. Что именно представляет первый параметр x? В документации написано «вектор или матрица квантилей. Если x - матрица, каждая строка считается квантилем».

> dmvnorm(x=c(0,0), mean=c(1,1))
[1] 0.0585

Вот пример кода на странице справки. В этом случае вы генерируете вероятность наличия квантиля 0 при нормальном распределении со средним значением 1 и sd 1 (при условии, что это значение по умолчанию). Поскольку это многомерная функция нормальной плотности и был передан вектор квантилей (0, 0), почему выходные данные не являются вектором вероятностей?


r statistics probability-density
person Adrian    schedule 05.06.2017    source источник

Ответы (1)


arrow_upward
1
arrow_downward

Просто возьмем в качестве примера двумерный нормальный (X1, X2), передав x = (0, 0), вы получите P(X1 = 0, X2 = 0), которое является единственным значением. Почему вы рассчитываете получить вектор?

Если вам нужен вектор, вам нужно передать матрицу. Например, x = cbind(c(0,1), c(0,1)) дает

P(X1 = 0, X2 = 0)
P(X1 = 1, X2 = 1)

В этой ситуации каждая строка матрицы обрабатывается параллельно.

person Zheyuan Li    schedule 05.06.2017