譯殊襯 6 語
襯 蟲.
> #讌螻磯逢覯
> k = (factorial(6)/(factorial(1)*factorial(2)*factorial(3)))
> k * (1/6)^1 * (2/6)^2 * (3/6)^3
[1] 0.1388889
> #dmultinom() 伎
> dmultinom(c(1,2,3), prob=c(1/6, 2/6, 3/6))
[1] 0.1388889
>
企 螻殊殊 , 譴, , 覿蟆 4譬襯襦 覿襯.
觜企 , 襦 20螳襯 豬 蟆曙 , 譴, , 覿蟆 螳螳 5,8,4,3 螳 襯 蟲.
> #讌螻磯逢覯
> n <- factorial(20)
> x1 <- factorial(5)
> x2 <- factorial(8)
> x3 <- factorial(4)
> x4 <- factorial(3)
> p1 <- 0.25^5
> p2 <- 0.40^8
> p3 <- 0.25^4
> p4 <- 0.10^3
> (n / (x1*x2*x3*x4)) * (p1*p2*p3*p4) #蟯碁ゼ 讌 ｌ伎 .
[1] 0.008729721
> #dmultinom() 伎
> dmultinom(c(5,8,4,3), prob=c(0.25, 0.40, 0.25, 0.10))
[1] 0.008729721
>
4 螳 譯殊襯 襴曙朱 讌 , 1 2螳, 2 1螳 襯 (狩覿)
> p <- factorial(4) / (factorial(2)*factorial(1)*factorial(1))* (1/6)^2 * (1/6)^1 * (1/3)^1
> p * 2
[1] 0.03703704