如何使用..命令

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")