R Multinom 多項式分布

說明

生成多項分布的隨機數向量並計算多項概率。

用法

rmultinom(n, size, prob)
dmultinom(x, size = NULL, prob, log = FALSE)

參數

細節

如果 x 是 K 分量向量，則 dmultinom(x, prob) 是概率

其中C是“多項式係數”C = N! / (x_1! \cdots x_K!)和N = \sum_{j=1}^K x_j.
根據定義，每個組件X_j二項分布為Bin(size, prob[j])為了j = 1, \ldots, K.

rmultinom() 算法從 Bin(n_j,P_j) 中依次提取二項式 X_j，其中 n_1 = N (N := size )、P_1 = \pi_1(\pi 是 prob 縮放至和 1)，對於 j \ge 2 ，遞歸地，n_j = N - \sum_{k=1}^{j-1} X_k 和 P_j = \pi_j / (1 - \sum_{k=1}^{j-1} \pi_k) 。

值

對於 rmultinom() ，整數 K \times n 矩陣，其中每列都是根據所需多項式定律生成的隨機向量，因此求和為 size 。雖然轉置結果一開始看起來更自然，但由於按列存儲，返回的矩陣更有效。

注意

dmultinom 目前根本沒有矢量化，也沒有 C 接口 (API)；將來可能會對此進行修改。

例子

rmultinom(10, size = 12, prob = c(0.1,0.2,0.8))

pr <- c(1,3,6,10) # normalization not necessary for generation
rmultinom(10, 20, prob = pr)

## all possible outcomes of Multinom(N = 3, K = 3)
X <- t(as.matrix(expand.grid(0:3, 0:3))); X <- X[, colSums(X) <= 3]
X <- rbind(X, 3:3 - colSums(X)); dimnames(X) <- list(letters[1:3], NULL)
X
round(apply(X, 2, function(x) dmultinom(x, prob = c(1,2,5))), 3)

也可以看看

Distributions 用於標準發行版，包括 dbinom，這是概念上的特殊情況。