Files
playbook/docs/tsl/codegen/dotnet/runtime_support/pending_migrated.md
T

55 KiB
Raw Blame History

运行时支撑函数 - 未细分函数

函数

clusterKmeans(samples, class_count, options)

执行 K-means 聚类。

参数 类型 说明
samples array 样本矩阵
class_count any 聚类数量
options array 参数表;常见键包括 maxitertoldistanceout

返回:array

clusterWard(samples)

执行 Ward 聚类。

参数 类型 说明
samples array 样本矩阵

返回:array

eigenJacobi(matrix, eigen_values, eigen_vectors)

用 Jacobi 方法计算矩阵特征值和特征向量。

参数 类型 说明
matrix array 待分解矩阵
eigen_values any 接收特征值数组
eigen_vectors any 接收特征向量矩阵

返回:void

ozProg(f, a, b)

求解 0/1 规划问题。

参数 类型 说明
f array 目标函数线性系数
a array 线性不等式约束左侧矩阵
b array 线性不等式约束右侧向量

返回:array

示例

f:=array(-3,2,-5);
A:=array((1,2,-1),(1,4,1),(1,1,0),(4,0,1));
B:=array(2,4,3,6);
OZProg(f,A,B);
//结果:x1=1,x2=0,x3=1,f(x)=-8

regressRidge(y_values, x_values[, k_values][, trace_mode])

执行岭回归,或按 k_values 输出岭迹。

参数 类型 说明
y_values array 因变量样本
x_values array 自变量样本矩阵
k_values array 岭参数序列;省略时自动选择
trace_mode bool 是否按岭迹模式输出

返回:array

regressWls(y_values, x_values, weights, residuals)

执行加权最小二乘回归并输出残差。

参数 类型 说明
y_values array 因变量样本
x_values array 自变量样本
weights array 权重序列
residuals any 接收残差序列

返回:array

invokeinarray(obj, member_name, member_value, params)

调用类的成员方法或对类的成员变量进行赋值操作。支持对重载方法的调用。基本功能与invoke一样,区别在于设置成员变量的值与调用成员方法的参数值通过数组进行传递。另外,支持命名参数的调用。

参数 类型 说明
obj object 实例对象
member_name string 字符串(小写),成员(方法或属性)的名称
member_value bool object
params array 数组,若是成员属性则数组传递的是设置的成员属性值,若是成员方法则表示成员方法的参数,成员方法没有参数时须初始化为空数组。命名参数调用时,下标字符串作为参数名。

返回:array

frameconcat(arrdic, type)

数据表-拼接

参数 类型 说明
arrdic array Array of number,数组的数组,需要拼接的数组集合
type any Usrdefine,数组拼接方式,默认为0行并=0;列并=1;非完全矩阵列并=2

返回:array

示例

a := ones(3,3);
b := ones(2,3);
frameConCat(array(a,b));

randgam(a, b, row, col)

产生伽马分布的随机数

参数 类型 说明
a float 形状参数,大于0的实数
b float 尺度参数,大于0的实数
row array 行数或者行名字符串数组
col array 列数或者列名字符串数组

返回:any

示例

randgam(2,2,100,array("a","b"));

sfBetapdf(x, a, b, v)

贝塔分布的密度函数,输入随机变量,返回随机变量对应的密度函数值

参数 类型 说明
x array :随机变量,可以使用实数,也可以是一个实数数组
a float beta分布参数1,实数,大于0
b float beta分布参数2,实数,大于0
v array 变参返回,该函数也支持直接返回,类型和x一样

返回:any

hypothSnuTest(x, u, d, tail, alpha)

单个正态总体的均值检验,分为标准差未知和已知两个情况,返回结果是一个数组,第一个数为统计值,二为p值,三为原假设的值(0表示拒绝,1表示接受) 模型表达: 原假设:tail取值both:;left:;right: 检验统计量:分为两种情况 1)标准差已知时,即d有取值,为U检验: 2)标准差未知时,即d取nil,为T检验: 其中:x为样本均值,为检验均值,为正态总体标准差,n为样本容量,s为样本标准差

参数 类型 说明
x array 样本序列,为一维数组类型
u float 假设检验均值,为实数类型
d float 正态总体的标准差,为实数类型,若未知,则输入nil
tail string 假设检验的三种类型,字符串,取值为both、left和right
alpha float 显著性水平,实数或整数,在01之间

返回:array

hypothTnuTest(x, uu, d, tail, alpha, isequal)

两个正态总体均值差的检验,分为标准差未知和已知两个情况,返回的结果是一个数组,第一个数为统计值,二为p值,三为是否接受原假设 模型表达: 原假设:tail取值both:;left:;right: 检验统计量: 标准差已知时,即d数组有取值,为u检验: 标准差未知时,即d取(nil,nil),为t检验:, 其中:为样本x的均值,为样本y的均值,,为正态总体x,y的方差,分别为样本x和y的容量,为联合样本标准差

参数 类型 说明
x array 样本序列矩阵,为二维数组类型,第一行为样本1,第二行为样本2
uu float 假设检验均值差,为实数类型
d float 正态总体的标准差数组,为一维数组类型,若未知,则输入(nil,nil)
tail string 假设检验的三种类型,字符串,取值为both、left和right
alpha float 显著性水平,实数或整数,在01之间
isequal bool object

返回:array

hypothSndTest(x, d2, tail, alpha)

单个正态总体方差的卡方检验,返回的结果是一个数组,第一个数为统计值,二为p值,三为原假设的值(0表示拒绝,1表示接受) 模型表达: 原假设:tail取值both:;left:;right: 检验统计量: 其中:为假设检验的方差,n为样本容量,s为样本标准差

参数 类型 说明
x array 样本序列,为一维数组类型
d2 float 假设检验的方差,为实数类型
tail string 假设检验类型
alpha float 显著性水平,实数或整数,在01之间

返回:array

hypothTndTest(x, tail, alpha)

两个正态总体方差比的F检验,返回的结果是一个数组,第一个数为统计值,二为p值,三为原假设的值(0表示拒绝,1表示接受) 模型表达: 原假设:tail取值both:;left:;right: 检验统计量: 其中:Sx为x的样本标准差,Sy为y的样本标准差

参数 类型 说明
x array 样本序列矩阵,第一行为样本1,第二行为样本2,为二维数组类型
tail string 假设检验类型
alpha float 显著性水平,实数或整数,在01之间

返回:array

hypothLsuTest(mx, s, n, u, tail, alpha)

单个或两个总体大样本的均值检验,即大样本U检验,返回的结果是一个数组,第一个数为统计值,二为p值,三为原假设的值(0表示拒绝,1表示接受) 模型表达 单个总体大样本的均值原假设:tail取值both:;left:;right 两个总体大样本的均值假设:tail取值both:;left;right: 检验统计量: 或 其中:x为样本x的均值,y为样本y的均值,为正态总体x方差,为正态总本y方差,m,n分别为样本x和y的容量s

参数 类型 说明
mx array 样本均值数组,第一行为样本1,第二行为样本2,为一维数组类型
s any 样本方差数组,第一行为样本1,第二行为样本2,为一维数组类型
n array 样本容量,第一行为样本1,第二行为样本2,为一维数组类型
u float 单个总体时,为假设检验均值,两个总体时,为均值差
tail string 假设检验类型
alpha float 显著性水平,实数或整数,在01之间

返回:array

anovaSingle(y, alpha)

在1-alpha置信水平下,在假设服从正态分布,且方差齐性,样本之间独立的基础上,判断各水平的均值是否显著相等 模型表达: , 其中:效应平方和:,误差平方和: 总离差

参数 类型 说明
y array 数据矩阵,为二维数组类型,每一行表示一个水平
alpha float 显著性水平,实数或整数,在01之间

返回:array

anovaTwice(y, alpha)

在1-alpha的置信水平下,进行双因子方差分析,假设每一种情况都服从一个正态分布,且方差相等,数据独立,判断各水平的均值是否显著相等 模型表达: 有A,B两个因素,因素A有r个水平,B有s个水平,对于A,B每一个组合(Ai,Bj)都进行t次试验,得到结果: 因素B 因素A B1 B2 … Bs A1 X111, X112,…,X11t X121, X122,…,X12t … X1s1, X1s2,…,X1st A2 X211, X212,…, X21t X221, X222,…, X22t … X2s1, X2s2,…, X2st … … … … Ar Xr11, Xr12,…, Xr1t Xr21, Xr22,…, Xr2t … Xrs1, Xrs2,…, Xrst 设:且相互独立 ,其中,是水平Ai的效应,是水平Bj的效应,称为水平Ai和Bj的交互效应。, 原假设:,, 现构造了3个F统计量来分别检验假设,原理略

参数 类型 说明
y array 数据矩阵,为二维数组类型,每rep行,每一列表示一个交叉项
alpha float 显著性水平,实数或整数,在01之间;Rep:每一个交叉项,做rep次试验,rep=1时,没有交叉效应为0的假设

返回:array

anovaBartlett(y, alpha)

在1-alpha的置信水平下,Bartlett方法检验方差齐性;B-Stat为统计值,P-Value为p值,df为自由度,Hypothesis为是否接受原假设 模型表达: Bartlett检验思想来源于几何平均数小于或等于算术平均数 原假设:是水平个数 ,其中 修正的Bartlett检验:

参数 类型 说明
y array 数据矩阵,为二维数组类型,每一行表示一个水平
alpha float 显著性水平,实数或整数,在01之间

返回:array

anovaLevene(y, alpha)

在1-Alpha的置信水平下,Levene方法检验方差齐性,和Bartlett检验不同的是,Levene检验不要求总体具有正态性。返回一个数组,F-Stat为统计值,P-Value为p值,df1df2为自由度,Hypothesis为是否接受原假设 模型表达: 原假设:是水平个数 ,其中是观察数目,是水平数目

参数 类型 说明
y array 数据矩阵,为二维数组类型,每一行表示一个水平
alpha float 显著性水平,实数或整数,在01之间

返回:array

hypothSignTest(x, pi, q0, tail, alpha)

非参数检验-符号检验。对单样本进行符号检验 可以按照实际需求选择双侧或单侧检验 其中,检验变量服从二项分布

参数 类型 说明
x array 样本序列,为一维数字数组类型
pi float 分位点,为实数类型
q0 float 检验值,即检验分位点的值,为实数类型
tail string 假设检验的三种类型,both(双侧检验),left(左侧检验),right(右侧检验),为字符串类型
alpha float 显著性水平,在0-1之间,为实数类型

返回:array

hypothPSignTest(x, pi, q0, tail, alpha)

非参数检验-配对样本符号检验,配对样本的符号检验 要求样本矩阵中两个样本是配对出现的 可以按照实际需求选择双侧或单侧检验

参数 类型 说明
x array 一维数字数组,样本矩阵
pi float 实数,分位点
q0 float 实数,检验值
tail string 字符串,双侧检验 or 左侧检验 or 右侧检验
alpha float 实数,显著性水平

返回:array

correlationMatrix(a1, t)

相关系数矩阵求解

参数 类型 说明
a1 array
t float

返回:array

示例

a := rand(4,3);
CorrelationMatrix(a,1);

covariance(a)

计算协方差矩阵

参数 类型 说明
a array

返回:any

示例

a := rand(4,3);
Covariance(A);

corrPartial(x, y, z, type)

偏相关系数,返回x和y在控制变量z下的偏相关系数,行为样本,列为指标

参数 类型 说明
x array 变量一,二维实数数组
y array 变量二,二维实数数组
z array 控制变量,二维实数数组
type any 1:pearson偏相关;2:Spearman偏相关;tail:假设检验的三种类型,字符串,取值为both、left和right,缺省是只输出相关系数

返回:any

factor(x, var_value)

对样本数据进行因子分析,分别采用主成分分析法估计因子载荷矩阵(因子载荷矩阵不唯一),方差最大法做载荷矩阵的正交旋转矩阵,回归方法计算因子得分

参数 类型 说明
x array 二维维实数数组,表示样本矩阵;
var_value float 实数,表示方差贡献率,控制因子个数;;返回结果:ssRet["Correlation Matrix",0]:样本数据的相关系数矩阵;Ret["kmo",0]:KMO检验值;0.9<KMO: 非常适合;0.8<KMO<0.9:适合;0.7<KMO<0.8:一般; 0.6<KMO<0.7:适合度较低;KMO<0.6:适合度很低;ret["Total Variance Explained",0]:特征值、方差贡献率、累计方差贡献率等信息;ret["Component Matrix",0]:因子载荷矩阵;ret["Communalities",0]:共同变量度,即载荷矩阵行的平方和;ret["thogonal matrix",0]:正交变换后的载荷矩阵;ret["Rotation matrix",0]:旋转矩阵;ret["Factor score",0]:因子得分

返回:array

示例

stockn:=getbk("上证A股");
vo:=BackUpSystemParameters();
SetSysParam(pn_date(),today());
sto:=array();
for nI:=0 to length(stockn)-1 do
begin
SetSysParam(pn_stock(),stockn[nI]);
if IsStockGoMarket(inttodate(20011231)) then
sto[length(sto)]:=stockn[nI];
end
st:=array2str(sto,";");
A1:= Query("",st,True,"","代码",DefaultStockID(),
"主营业务利润",report(46008,20101231),
"利润总额",report(46024,20101231),
"净利润",report(46033,20101231),
"加权每股收益",report(42003,20101231),
"每股净资产",report(42006,20101231),
"加权净资产收益率",report(42013,20101231),
"资产总计",report(44059,20101231),
"股本",report(44099,20101231));
RestoreSystemParameters(vo);
a:=zeros(mrows(a1),mcols(a1));
a[:,:]:=a1[:,:];
stock:=a[:,0];
data:=a[:,1:8];
//取得纯数据矩阵
br:=factor(data,0.85); //
br;

histc(a, b)

端点定位频数,判断序列a 中各个元素在b中各个段位中的频数

参数 类型 说明
a array 一维数组
b array 一维数组

返回:array

示例

a:=rand(1000)*5;
b:=array(0,1,2,3,4,5);
HISTC(a,b);

regressDominFormula(y, x)

相对重要性分析-R2公式分解法 功能描述: 1、相对重要性分析 2、算法:根据含常数项的多元回归的r2的分解公式,来计算 不同解释变量的相对贡献度 3、公式法: 未考虑到不同解释变量之间的相互关系, 贡献度可能为负数

参数 类型 说明
y array Array of number,因变量
x array Array of number,自变量

返回:array

regressMls(y, x, u)

线性回归方程的最小二乘法参数估计(包含常数项),返回回归方程的系数,常数项排在第一项

参数 类型 说明
y array 被解释变量序列,为一维数组类型
x array 解释变量矩阵,为二维数组类型,每一列为一个解释变量
u array 残差序列,为一维数组类型,变参

返回:array

regressCmls(y, x, u, constant)

线性回归方程的最小二乘法参数估计(可选择是否包含常数项),返回回归方程的系数,如果有常数项则排在第一项

参数 类型 说明
y array 被解释变量序列,为一维数组类型
x array 解释变量矩阵,为二维数组类型,每一列为一个解释变量
u array 残差序列,为一维数组类型,变参
constant bool 是否包含常数项,为布尔类型,缺省为true

返回:array

regressDwTest(u)

德宾-沃森检验:检验随机误差项是否存在一阶自相关。结果返回DW统计量 DW统计量:(一般认为值在2附近,则随机误差项不存在自相关)

参数 类型 说明
u array 方程残差,为一维数组类型;

返回:any

regressTTestRwls(x, a, u, k, alpha, constant, weight, restrictiveness)

T检验,在alpha的置信水平下,n元线性回归变量系数的显著性检验 原假设:回归方程系数不显著。带约束的加权最小二乘 RWLS 模型,对回归系数进行显著性检验, 适用于以下

参数 类型 说明
x array 一维数字数组,自变量矩阵,注意:含常数项的回归不需要给全1列
a array 一维数字数组,方程系数;有常数项时,a 为 常数项系数 union 自变量系数;无常数项时,a 为 自变量系数
u array 一维数字数组,残差序列
k any 整型,自变量的个数;k = length(a)-constant;有常数项回归时,k := mcols(x);;无常数项回归,k := mcols(x)-1;
alpha float 实数,显著性水平
constant bool Bool,是否为含常数项回归
weight array 数组,加权回归的权重,默认等权重
restrictiveness array 数组,线性约束,只需要给对应于回归系数个数长度的一维数组,默认无约束;例如:有5个系数,分别为x0,x1,x2,x3,x4;约束为: x1+x2+3x4=3;则所给的参数为:array(0,1,1,0,3)

返回:array

regressionOlsNw(y, x, constant, alpha, al)

NW调整的OLS回归,此调整用于当残差或者因变量(两者等价)存在异方差性并且自相关性时,利用OLS回归后对回归系数的协方差矩阵进行NW调整计算,得到每个系数White调整后的标准误,再进行T检验。T检验部分相关公式详见Regress_TTest_NW,回归系数的OLS解如下β*=(X'X)-1X'y 其余部分定义与普通最小二乘OLS一致

参数 类型 说明
y array 一维数字数组, 因变量
x array 二维数字数组, 自变量,不需要包含全1列
constant bool object
alpha float 实数,显著性水平,应为0到1之间的值,默认0.05
al any 整数,自相关最大滞后阶数,可缺省,默认用Newey and West (1994)?提出的方法计算,至少为1

返回:array

leastSquaresCurveFitting(f, y)

最小二乘法进行回归

参数 类型 说明
f array 二维数字数组类型,系数矩阵
y array 一维数字数组类型,观察数据值

返回:array

multipleRegression(x, y, call_back)

多元回归

参数 类型 说明
x array 二维数字数组类型,观察数据点
y array 一维数字数组类型,观察数据值
call_back string 字符串类型,计算系数矩阵的回调

返回:array

示例

x:=array((0.00,0.00),(0.00,1.00),(1.00,0.00),(1.00,1.00));
y:=array(1.00,2.00,2.00,3.00);
CallBack:='MultipleRegression_CallBack';
MultipleRegression(x,y,CallBack);
//结果:

singleUnLineralRegressionCallBack(x)

单元非线性回归的回调函数

参数 类型 说明
x array 一维数字数组类型

返回:array

eastrq(b, c, q)

Eastrq 相关函数。

参数 类型 说明
b array 一维数字数组类型
c array 一维数字数组类型
q array 二维数字数组类型

返回:any

rotatedComponentMatrix(cm, eigen_value, eigen_vector)

该函数已废弃

参数 类型 说明
cm array 二维数字数组,相关系数矩阵
eigen_value array 一维数字数组,特征根
eigen_vector array 二维数字数组,特征向量

返回:any

示例

//设置样本矩阵
zt:=array((0.7592,1.7593,21.9279,48.5626,55.1679),
(0.0490,0.4457,2.2255,82.2658,192.4348),
(-2.2313,-1.3096,-59.8883,3.9922,-80.2419),
(0.1364,0.539,4.3910,55.1438,46.0612),
(0.7233,1.9467,15.6127,111.3018,58.3969));
zT_s:=Standarize_II(zT);  //标准化
mt_Transposition(zT_s,Z);  //转置
R:=CorrelationMatrix(Z);  //相关系数矩阵
Eigen_Jacobi(R,EigenValue,EigenVector); //特征值&特征向量
Cm:=ComponentMatrix(EigenValue,EigenVector,MaxEigenValue); //载荷矩阵
RotatedComponentMatrix(Cm);
//结果:

getEigenValue(a, eigen_value)

通过赋值得到EigenValue

参数 类型 说明
a array 二维数字数组类型
eigen_value array 一维数字数组类型

返回:any

timeDiff(y, n)

序列差分

参数 类型 说明
y array 需要做差分的数组数据
n any 差分的阶数,0阶不做处理

返回:any

timeAutoCov(y, l)

自协方差,求样本序列跟滞后L阶序列的协方差,即自协方差 自协方差计算公式如下: 其中:yt是样本序列,yL是滞后L阶的序列。

参数 类型 说明
y array 样本序列,为一维数组类型;
l any 滞后阶数,为整数类型;

返回:any

timeAcf(y)

自相关系数,计算随机序列与所有滞后阶数的相关系数。

参数 类型 说明
y array 数组,样本序列;;p:整数,返回样本自相关系数序列的长度,默认值为3,允许缺省;

返回:array

timeDanielTest(y, alpha)

在1-alpha置信水平下,进行随机时间序列的Daniel平稳性检验。返回结果为一个数组,第一个为q值(q>0:序列有上升趋势,q<0:序列有下降趋势),第二个为T统计量,第三个为p值,第四个为是否接受原假设。原假设为x是平稳序列,所以H为1表示接受原假设,0表示拒绝原假设 其中:qs是Spearman相关系数,Rt为秩统计量,n为样本容量

参数 类型 说明
y array 随机时间序列,为一维数组类型
alpha float 显著性水平,实数类型,取值在0-1之间

返回:array

arma(series, p, q)

自回归移动平均模型,ARMA(p,q)模型,采用逆函数去估计自回归系数和移动平均系数,返回结果有自回归系数、移动平均系数、残差方差等信息 ARMA模型: 其中:为原序列零均值化的新序列,是自回归方程系数,为移动平均系数,是零均值、方差为的平稳白噪声

参数 类型 说明
series array 样本序列,为一维数组类型,注意要零均值化
p any 自回归阶数,整数类型
q any 移动平均阶数,整数类型;返回结果:Ret["AR"]:参数估计结果,,第一项为1Ret["MA"]:参数估计结果,,第一项为1Ret["T-test"]:AR模型参数的T检验;Ret["LossFunction"]Ret["FPE"]:预测最大误差;Ret["AIC"]:模型的AIC信息准则;Ret["SC"]:模型的SC信息准则;Ret["e"]:残差信息;Ret["Inverted AR Roots"]:AR特征方程的特征根;Ret["Inverted MA Roots"]:MA特征方程的特征根;Ret["R_squared"]:AR模型的可决系数;Ret["Adjusted_R_squared "]:AR模型的修正可决系数;Ret["Sum_squared_resid"]:AR模型的离差平方和;Ret["Durbin Watson stat"]:DW统计量;Ret["Log Likelihood"]:对数极大似然值;Ret["S.E. of regression"]:回归标准差

返回:array

示例

elps:=randnorm(0,1,200);
y:=array();
y[0]:=0;
for i:=1 to 199 do
y[i]:=0.8*y[i-1]+elps[i]-0.4*elps[i-1];
ARMA(y,1,1);

timeArch(r_series, q, l)

自回归条件异方差模型,本模型的ARCH模型如下: sss 可把异方差方程进行整合,得到 用最小二乘法去估计方程的系数了,得到系数的相合估计。 或者用最大似然估计去估计参数,得到精度估计,最大似然估计的目标函数值为:

参数 类型 说明
r_series array 资产收益率,为一维数组类型;
q any 自回归阶数,为整数类型;
l any 预测步长,为整数类型

返回:array

timeGarch(rseries, model)

广义自回归条件异方差模型,包含四种模型:GARCH(1,1)、GJR(1,1)、IGARCH(1,1)和EGARCH(1,1),本模型采用了最大似然估计法去估计模型的系数,可做一步预测。 GARCH(1,1) GJR(1,1) IGARCH(1,1) EGARCH(1,1) 最大似然估计的目标函数值为:

参数 类型 说明
rseries array 资产收益率,为一维数组类型;
model string 模型类别,为字符串类型;

返回:any

timesCointergrationTest(y, x, constant, alpha_type, style, criterion)

E-G两步法协整检验变量的协整关系,E-G法基本步骤如下: 1)在因变量与自变量存在同阶单整关系的条件,进行最小二乘法回归; 2)把回归后的残差进行平稳性检验,如果残差序列是平稳的,则两者存在协整关系,否则不存在。

参数 类型 说明
y array 被解释变量,为一维数组类型
x array 解释变量,为一维数组类型
constant bool 回归方程是否含有常数项,为布尔类型
alpha_type float 显著性水平,为实数类型,一般在0-1之间
style any ADF法检验单整阶数时是否存在常数或趋势项,整数类型,无常数无趋势使用0、有常数无趋势1、有常数有趋势2
criterion string 滞后阶数判断信息准则,字符串类型,"SBIC""AIC"

返回:array

timesEcm(y, x, resid, alpha, q, p)

两变量间的误差修正模型,准确的刻画了两序列间的短期波动特征 误差修正模型是时间序列间协整关系的主要表现形式,准确的刻画了序列间的长期均衡关系和短期波动特征。设一阶自回归分布滞后模型为: 其中,是均值为零,方差为的白噪声,进行适当的整理有: 其中若记 称为误差修正项。

参数 类型 说明
y array 被解释变量,为一维数组类型
x array 解释变量,为一维数组类型
resid array 协整检验时的残差,没有输入,使用最小二乘得到
alpha float 显著性水平
q any 整数类型,被解释变量的滞后项
p any 整数类型,解释变量的滞后项

返回:any

timesGranger(y, x, q, alpha)

Granger因果检验,原假设:X不是Y的Granger原因,Hypothesis为1表示接受原假设,0表示拒绝原假设 原假设: 分别对上述两个回归方程做回归,提取残差平方和XESS跟YESS,则F统计量如下计算:

参数 类型 说明
y array 被解释变量,为一维数组类型
x array 解释变量,为一维数组类型
q any 滞后阶数,为整数类型
alpha float 显著性水平,为实数类型,一般在0-1之间

返回:array

timeARp(y, p, l)

建议使用AR模型替换,自回归模型,AR(p)模型,用Yule-Walker方程去估计p阶自回归模型的系数,返回结果有自回归系数、残差方差、预测值及预测波动。估计过程如下 AR(p)模型: 如有常数项,则 ,为自协方差 其中:为原序列零均值化的新序列,为原序列的均值,是自相关系数,是自回归方程系数,是零均值,方差为的平稳白噪声

参数 类型 说明
y array 样本序列,为一维数组类型;
p any 自回归阶数;
l any 预测步长;

返回:array

timeMAq(y, q, l)

建议使用ARMA模型替换,移动平均模型,MA(q)模型,用条件最小二乘法去估计q阶移动平均模型的系数,返回结果有移动平均系数、残差方差、预测值及预测波动。 估计过程如下: MA(q)模型: 具有逆转形式: 用算子形式比较,可得: s 条件最小二乘法使得下列残差平方和最小: 可将上述转化为矩阵形式,先求得的取值,在利用上上式计算移动平均系数 如有常数,则常数项为原序列的均值。 ,为自协方差 其中:yt为原序列零均值化的新序列,μ为原序列的均值,θq为移动平均系数,εt是零均值,方差为σ2的平稳白噪声

参数 类型 说明
y array 样本序列,为一维数组类型;
q any 移动平均阶数;
l any 预测步长;

返回:array

timeArma(y, p, q, l)

建议使用ARMA模型替换,自回归移动平均模型,ARMA(p,q)模型,采用逆函数去估计自回归系数和移动平均系数,返回结果有自回归系数、移动平均系数、残差方差、预测值及预测波动。 估计过程如下: ,为自协方差 其中:yt为原序列零均值化的新序列,μ为原序列的均值,ρp是自相关系数,φq是自回归方程系数,θq为移动平均系数,εt是零均值、方差为σ2的平稳白噪声

参数 类型 说明
y array 样本序列,为一维数组类型;
p any 自回归阶数;
q any 移动平均阶数
l any 预测步长,整数

返回:array

timeInvFun(ap, mq)

中间函数,模型逆函数,ARMA模型的逆函数,Green函数的递推公式如下: 式中: ,,是自回归方程系数,为移动平均系数,

参数 类型 说明
ap array 自相关系数,为一维数组类型;
mq array 移动平均系数,为一维数组类型;

返回:array

timeAnalyse(y, model, p, q, l)

ARMA模型分析,时间序列的ARMA模型分析,可以自动的选择q与p

参数 类型 说明
y array 样本序列,为一维数组类型;
model string 模型类别为AR、MA、ARMA,为字符串类型;
p any 自回归阶数,取0则自动选择阶数,为整数类型;
q any 移动平均阶数,取0则自动选择阶数,为整数类型;
l any 预测步长,为整型类型;

返回:any

timeChoseP(y, model, f)

判断AR模型最佳的阶数,采用信息准则函数判断模型的最佳阶数 其中:T为样本容量,K为待估计参数

参数 类型 说明
y array 样本序列,为一维数组类型;
model string 模型类别为AR、MA、ARMA,为字符串类型;
f string 判断准则为SBIC和AIC,为字符串类型

返回:any

timesCointergration(y, x, alpha)

建议使用Times_Cointergration_test函数替换本函数,协整检验—E-G两步法检验变量的协整关系,E-G法如下: 1)在因变量与自变量存在同阶单整关系的条件,进行最小二乘法回归; 2)把回归后的残差进行平稳性检验,如果残差序列是平稳的,则两者存在协整关系,否则不存在。

参数 类型 说明
y array 因变量,为一维数组类型;
x array 自变量矩阵,为二维数组类型;
alpha float 1-alpha为置信水平,一般为5%,也有1%和10%,数值越小,判断越严格;

返回:string

timeArchTest(u, q, alpha)

时间序列异方差的ARCH检验,原假设都为没有异方差,即不存在ARCH效应,Hypothesis为1表示接受原假设,0表示拒绝原假设 检验统计量为: 对方程进行估计,得到F统计量和T×R2统计量并进行检验;

参数 类型 说明
u array 残差序列,为一维数组类型
q any 滞后阶数,为整数类型
alpha float 1-alpha为置信水平,一般为5%,也有1%和10%,数值越小,判断越严格

返回:array

timesJohansen(x, p, k)

进行Johansen协整检验,第一列的迹、最大特征值统计量的原假设是:没有协整关系;第二列的迹、最大特征值统计量的原假设是最多只有一个协整关系;第三列的迹、最大特征值统计量的原假设是最多只有两个协整关系……

参数 类型 说明
x array 输入多元时间序列,每一列为一指标,每行为一时间 (nobs*m)
p any 是否含时间趋势及常数项p = -1,没有时间趋势以及常数项p =0,含常数项p =1,含常数项以及一次时间项p >1,含常数项以及高次时间项
k any 滞后阶数

返回:array

progQ(h, f, a, b, aeq, beq, lb, ub, x0, options)

二次规划问题求解(采用作用集方法) Prog_Q(H,f,A,b) 求解不等式约束 Prog_Q(H,f,A,b,Aeq,beq) 求解不等式约束、等式约束 Prog_Q(H,f,A,b,Aeq,beq,LB,UB) 求解不等式约束、等式约束、上下界限约束 Prog_Q(H,f,A,b,Aeq,beq,LB,UB,x0,options) 模型表达: min= `x* H* x /2 + sum(f*x) st. A:*x<=b Aeq :* x = Beq; lb<x<ub;

参数 类型 说明
h array 二次型,正定二次型二维数字数组,不可缺省
f array 线性目标,一维数字数组,不可缺省
a array 线性约束,二维数字数组,不可缺省
b array 线性约束,一维数字数组,不可缺省
aeq array 线性约束,二维数字数组,缺省为空
beq array 线性约束,一维数字数组,缺省为空
lb array 下界约束,一维数字数组,缺省为空
ub array 上界约束,一维数字数组,缺省为空
x0 array 初始值,一维数组
options array 选项设置,数据表类型,缺省为array('maxiter': 400 , //最大迭代次数'tolcon': 1.0e-6, //约束容忍误差'tolfun' : 1.0e-6 , //目标最小误差'tolx' : 1.0e-6 //变量最小变化)

返回:array

nonLpFminbnd(fun, a, b, methods)

一维最小值搜索,求解模型如下: min fun(x)

参数 类型 说明
fun string 目标函数(函数名字符串)
a float 下界,实数
b float 上界,实数
methods any 方法选择0 :黄金分割法+抛物线法(默认)1 : 盲人收索法2 :黄金分割法3 :斐波那契法

返回:array

nonLpRosenbrock(fun, x0)

(无约束最小值求解建议统一使用NonLP_Fminsearch函数)Rosenbrock求无约束多维极值问题,返回最优点和最优值。

参数 类型 说明
fun string 目标函数表达式,为字符串类型;
x0 array 初值向量,为一维数组类型

返回:array

ozProg(f, a, b)

(线性规划,包括整数规划,01规划建议统一使用prog_L函数),0-1规划,解决如下0-1整数规划模型问题 返回该0-1规划模型的最小值及其对应变量的取值,其中可以转化为和;

参数 类型 说明
f array 一维实数数组,表示0-1规划模型的目标函数系数;
a array 二维实数数组,表示0-1规划模型的带“”约束方程系数矩阵,如果是带“”,则将该方程两边同乘以“-1”;
b array 一维实数数组,表示0-1规划模型带“”约束方程的约束值;

返回:array

示例

f:=array(-3,2,-5);
A:=array((1,2,-1),(1,4,1),(1,1,0),(4,0,1));
B:=array(2,4,3,6);
OZProg(f,A,B);
//结果:x1=1,x2=0,x3=1,f(x)=-8

hJMethodPgm(call_back)

H_J算法测试

参数 类型 说明
call_back string 字符串类型,计算目标值的回调

返回:array

outPutGrowValue2(v1, v2)

增长

参数 类型 说明
v1 float 实型,期末值
v2 float 实型,期初值

返回:float

示例

//返回13与14差值
OutPutGrowValue2(13,14);
//结果:-1

outPutGrowRatio(v1, v2)

增长率(%)

参数 类型 说明
v1 float 实型,期末值
v2 float 实型,期初值

返回:float

示例

//返回13到14的增长率(%)
OutPutGrowRatio(14,13);
//结果:7.69

outPutGrowRatio2(v1, v2)

增长率(%)

参数 类型 说明
v1 float 实型,期末值
v2 float 实型,期初值

返回:float

示例

//返回13到14的增长率(%)
OutPutGrowRatio2(14,13);
//结果:7.69

logined()

判定是否已经登录成功。只读属性。

返回:any

fUrl()

字符串,TS-OPS服务注册地址 类型:读写

返回:string

fHeader()

字符串,TS-OPS服务注册时需要指定的头部信息,比如TS-OPS-CALL等。 类型:读写

返回:string

fOnMessage()

函数指针,消息处理过程,指向一个函数。可通过Findfunction("函数名")产生 类型:读写

返回:any

create(url, header, onmessage, u, p, s, reconn)

构造方法

参数 类型 说明
url string TS-OPS服务注册地址
header string TS-OPS服务注册时需要指定的头部信息,比如TS-OPS-CALL等
onmessage any 消息处理过程,指向一个函数。可通过Findfunction("函数名")产生
u string 用户名
p string 用户密码
s string 租户Session Key;使用特殊租户(调用功能列表为*)替代用户时指定
reconn bool object

返回:any

timeSeriesSplit()

时间序列数据处理;将时间序列数据:时间、数据、(股票)处理成为(时间股票数据)的格式。注意,该函数主要用于客户端结果处理,不提供给用户程序编辑。

返回:array

reportSeriesSplit(table_name)

报告期数据处理。注意,该函数主要用于客户端结果处理,不提供给用户程序编辑。

参数 类型 说明
table_name string 字符串,用户数据名称

返回:any