R语言实现ARFIMA，源码和dll文件

R语言实现ARFIMA，源码和dll文件

/**************************************************************************************
durlevsim: adapted （with permission） from DLSim of the ltsa package.
input parameters:
	error - vector of length n innovation sequence
	nR - length of time series
	c - vector of length n containing the autocovariances
	EPSL - machine epsilon in double precision
	
	output parameters:
	z - vector of length n simulated time series
	ifault -  0 OK ;  1 c is not p.d.    
***************************************************************************************/

#include 

void durlevsim（double *z double *error int *nR double *c double *EPSL int *fault）  
{
	double sum;
	int ijkn;
	double EPS;

	
	*fault = 0;
	n = *nR;
	EPS = *EPSL;
	if （n < 1） *fault = 1;
	double v[n]phi[n]phiki[n]phikj[n]; 
	
	v[0] = c[0];
	z[0] = error[0]*sqrt（v[0]）;
	
	if （c[0] <= EPS） *fault = 1;
	phi[0] = c[1] / c[0];
	phiki[0] = phi[0];
	v[1] = v[0] * （1.0 - phi[0]*phi[0]）;
	z[1] = sqrt（v[1]）*error[1] + phi[0] * z[0];
	
	if （v[1] <= EPS） *fault = 1;
	
	for （k = 2; k < n; k++）
	{
		sum = 0.0;
		for （i = 1; i < k; i++）
			sum += phiki[i - 1] * c[k - i];	         
		phi[k - 1] = （c[k] - sum） / v[k - 1]; 
		for （j = 1; j< k; j++）
		{
			phikj[j - 1] = phiki[j - 1] - phi[k - 1] * phiki[k - j - 1];
		}
		phikj[k - 1] = phi[k - 1];
		sum = 0.0;
		for （j = 1; j<=k; j++）
		{
			sum += phikj[j - 1] * z[k - j];
			phiki[j - 1] = phikj[j - 1];
		}
		v[k] = v[k - 1] * （1.0 - phi[k - 1]*phi[k - 1]）;
		if （v[k] <= EPS） *fault = 1;
		error[k] *= sqrt（v[k]）; 
		z[k] = error[k] + sum;
		}
	
	return;    
}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2015-11-17 20:53  arfima\
     文件         437  2015-08-23 04:31  arfima\CITATION
     文件         673  2015-08-23 04:31  arfima\DEscriptION
     文件        2889  2015-08-23 04:31  arfima\INDEX
     文件        1299  2015-08-23 04:31  arfima\MD5
     目录           0  2015-11-17 20:45  arfima\meta\
     文件        1570  2015-08-23 04:31  arfima\meta\Rd.rds
     文件         117  2015-08-23 04:31  arfima\meta\data.rds
     文件        1426  2015-08-23 04:31  arfima\meta\hsearch.rds
     文件         554  2015-08-23 04:31  arfima\meta\links.rds
     文件         539  2015-08-23 04:31  arfima\meta\nsInfo.rds
     文件         824  2015-08-23 04:31  arfima\meta\package.rds
     文件         714  2015-08-23 04:31  arfima\NAMESPACE
     目录           0  2015-11-17 20:45  arfima\R\
     文件        1056  2015-08-23 04:31  arfima\R\arfima
     文件       68735  2015-08-23 04:31  arfima\R\arfima.rdb
     文件        1309  2015-08-23 04:31  arfima\R\arfima.rdx
     文件      207018  2015-11-17 20:53  arfima\arfima.pdf
     目录           0  2015-11-17 20:45  arfima\data\
     文件        1958  2015-08-23 04:31  arfima\data\SeriesJ.RData
     文件         872  2015-08-23 04:31  arfima\data\tmpyr.RData
     目录           0  2015-11-17 20:45  arfima\help\
     文件         881  2015-08-23 04:31  arfima\help\AnIndex
     文件         344  2015-08-23 04:31  arfima\help\aliases.rds
     文件       95862  2015-08-23 04:31  arfima\help\arfima.rdb
     文件         879  2015-08-23 04:31  arfima\help\arfima.rdx
     文件         417  2015-08-23 04:31  arfima\help\paths.rds
     目录           0  2015-11-17 20:45  arfima\html\
     文件        6187  2015-08-23 04:31  arfima\html\00Index.html
     文件        1319  2015-08-23 04:31  arfima\html\R.css
     目录           0  2015-11-17 20:47  arfima\inst\
............此处省略14个文件信息