r - 在 R 中写一个双和

Question

我是 R 的新手，正在尝试编写一个带有两个输入 p 和 q 的双和函数。例如，要计算 WolframAlpha 上的总和，我使用

Sum[p^25 * (1-q)^j * Sum[(25 choose r) * ((j-1) choose (r-1)) * ((q * (1-p)) / (p * (1 - q)))^r,{r,1,25}],{j,0,23}]

Here's an example为 p = .6和 q = .5
这应该很容易用 R 编写，但我不知道如何这样做。 (我意识到选择函数变成 choose(n,k) ，但事实证明总和是难以处理的。)

谢谢你的帮助。

Answer 1

这是一种方法:

p <- 0.6
q <- 0.5
qp <- q*(1-p)/(p*(1-q))
f <- function(r,j) choose(25,r)*choose(j-1,r-1)*qp^r
g <- function(j) p^25*(1-q)^j*sum(sapply(1:25,f,j))
sum(sapply(1:23,g))
# [1] 0.721113

您可以使用匿名函数代替 f和 g并将其全部放在一行上，但 IMO 这一点不那么令人困惑。