如何使用..命令

Question

這是我所做的。我需要5個隨機數，所以第一個命令是好的。之後，我想了解概率公式，但是我無法正確寫入。 我需要找到五個隨機數中每一個的概率。之後，我需要找到1中的所有N個數的概率89，做這樣一個在這裏的圖形..

Answer 1

你在你的函數有一個額外的括號，應

b <- function(n){ 
    x<- prod(1:365) 
    P <- ((x)/(365^n)*(prod(1:365-n))) 
    Pn <- 1-P 
} 

Answer 2

prod(1:365)（例如，365！）評估爲Inf，所以這不起作用。相反，使用此：

b <- function(n) 1 - prod(365:(365-(n-1))/365) 

這需要兩個事實的優點：

1. 365 /（365-N）！ = prod（365,364，...，365-（n-1）），即只有n個項的乘積。
2. 其中每個術語除以365（通過分母中的365^n）。

所以結果與prod（1,364/365,363/365，...，（365-（n-1））/ 365））相同。這個產品是可計算的。

但還有一個額外的問題：在plot(...)中使用此函數時，它必須是「向量化的」，也就是說它必須將一個向量作爲參數並返回一個長度與結果相同的向量。要做到這一點，使用Vectorize(...)功能R.

b <- Vectorize(function(n) 1 - prod(365:(365-(n-1))/365)) 

現在你可以使用plot(...)

plot(1:89,b(1:89), type="l")