从文档中,rbinom
有三个参数,n
,size
和prob
。 this question的答案指出,如果prob
是一个向量,它将被回收直到达到n
。我想找到那个代码!
R源立即调用C函数:
> rbinom
function (n, size, prob)
.Call(C_rbinom, n, size, prob)
<bytecode: 0x113526220>
<environment: namespace:stats>
在Github镜像中搜索R gets me the this rbinom
C function,但它只需要两个参数。 n
在哪里? prob
回收到n
的长度在哪里?我在这个函数定义中找不到它。
double rbinom(double nin, double pp)
{
...
}
搜索all C files containing rbinom
,我找到is this one的唯一其他相关的看似C文件,但我不明白。这是回收发生的地方,还是我缺少的另一个功能定义?
在堆栈溢出的有用指导下,我设法满意地回答了这个问题。我不确定细节,但是这里有:
> rbinom
叫这行
rbinom <- function(n, size, prob) .Call(C_rbinom, n, size, prob)
位于/src/library/stats/R/distn.R
的149号线
这反过来打电话
static R_INLINE SEXP random2(SEXP sn, SEXP sa, SEXP sb, ran2 fn, SEXPTYPE type)
与sn=n
,sa=size
,sb=prob
和fn=rbinom
。
在该函数中,您想要的代码在第185行:qazxsw poi,它显示rbinom被称为回收size和prob参数。
random2和rbinom之间的连接是使用与函数random2相同的文件中的预处理器建立的。
rx = fn(ra[i % na], rb[i % nb]);
DEFRAND2_INT(rbinom)