如何查找函数内部正在发生的事情ARMAtoMA？我得到了以下结果。
ARMAtoMA
#function (ar = numeric(), ma = numeric(), lag.max)
#.Call(C_ARMAtoMA, as.double(ar), as.double(ma), as.integer(lag.max))
#<bytecode: 0x000000001a5f5700>
#<environment: namespace:stats>
您可以在R的源文件中找到它：
R-<version>/src/library/stats/src/pacf.c
R-<version>/src/library/stats/是因为它来自stats包装。R-<version>/src/library/stats/src因为它是C代码。.c文件R-<version>/src/library/stats/src。如果在这些文件中搜索“ ARMAtoMA”，则可以找到pacf.c。我在这里为您制作一份副本：
SEXP ARMAtoMA(SEXP ar, SEXP ma, SEXP lag_max)
{
int i, j, p = LENGTH(ar), q = LENGTH(ma), m = asInteger(lag_max);
double *phi = REAL(ar), *theta = REAL(ma), *psi, tmp;
SEXP res;
if(m <= 0 || m == NA_INTEGER)
error(_("invalid value of lag.max"));
PROTECT(res = allocVector(REALSXP, m));
psi = REAL(res);
for(i = 0; i < m; i++) {
tmp = (i < q) ? theta[i] : 0.0;
for(j = 0; j < min(i+1, p); j++)
tmp += phi[j] * ((i-j-1 >= 0) ? psi[i-j-1] : 1.0);
psi[i] = tmp;
}
UNPROTECT(1);
return res;
}
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句