【发布时间】：2013-08-22 00:06:54
【问题描述】：
假设我有两个不同的参数向量用作函数的输入。此外，假设两个参数向量的长度不同。有没有办法输出该函数的所有可能值？我知道如果我使用两个不同长度的参数向量，那么较短的参数向量只会重复，这样就行不通了。我可以“手动”解决这个问题，如下所示。但是，我想找到一种更有效的方式来计算所有可能的组合。一个例子如下：
假设我正在使用 dbinom() 函数，该函数的输入为 x（样本中“成功”的次数）、n（样本中的观察次数）和 p（每个“成功”的概率X）。 n 保持在 20 不变；但是，我的 x 从 7,8,9,...,20（下面的“x7”）或 0,1（下面的“x1”）变化。另外，我想在 p 的不同值下评估 dbinom()，特别是从 0 到 1 以 0.1 为增量（下面的“p”）。可以看到三个参数向量x7、x1、p的长度都不同，分别是14、2、11。
> p<-seq(from=0,to=1,by=.1)
> x7<-seq.int(7,20)
> x1<-c(0,1)
> n<-20
我可以通过使用 dbinom() 中的一个向量（x7/x2 或 p）来评估每个组合，然后为剩余的参数选择一个值。正如您在下面看到的，我使用了向量 x7 或 x2，然后“手动”将 p 更改为等于 0,.1,.2,...,1。
> sum(dbinom(x7,n,.1))
[1] 0.002386089
> sum(dbinom(x7,n,.1))+sum(dbinom(x1,n,.1))
[1] 0.3941331
> sum(dbinom(x7,n,.2))+sum(dbinom(x1,n,.2))
[1] 0.1558678
> sum(dbinom(x7,n,.3))+sum(dbinom(x1,n,.3))
[1] 0.3996274
> sum(dbinom(x7,n,.4))+sum(dbinom(x1,n,.4))
[1] 0.7505134
> sum(dbinom(x7,n,.5))+sum(dbinom(x1,n,.5))
[1] 0.9423609
> sum(dbinom(x7,n,.6))+sum(dbinom(x1,n,.6))
[1] 0.9935345
> sum(dbinom(x7,n,.7))+sum(dbinom(x1,n,.7))
[1] 0.999739
> sum(dbinom(x7,n,.8))+sum(dbinom(x1,n,.8))
[1] 0.9999982
> sum(dbinom(x7,n,.9))+sum(dbinom(x1,n,.9))
[1] 1
> sum(dbinom(x7,n,1))+sum(dbinom(x1,n,1))
[1] 1
基本上，我想知道是否有办法让 R 使用单行输入或其他更有效的方法来改变参数 p 从 0.3941331,0.1558678,...,1 打印上述所有总和无需简单地复制和更改每一行的 p。
*我是 Stackoverflow 的新手，所以如果我没有按照惯例提出我的问题，我提前道歉。
【问题讨论】：
-
你想用
sum(dbinom(....))实现什么你肯定想要pbinom？
标签： r vector parameters