主要内容

lbqtest

Ljung-Box Q-test残差自相关

描述

例子

h= lbqtest (res)返回拒绝的决定h进行一个Ljung-Box Q-test自相关的残留系列res

例子

(h,pValue,统计,cValue)= lbqtest (res)还返回p价值pValue、测试数据统计,临界值cValue的测试。

例子

StatTbl= lbqtest (资源描述)返回的表StatTbl包含变量的测试结果、统计数据和设置进行Ljung-Box Q-test残差自相关的输入表的最后一个变量或时间表资源描述。选择一个不同的变量资源描述测试中,使用DataVariable名称-值参数。

例子

(___)= lbqtest (___,名称=值)指定选项使用一个或多个名称参数除了任何输入参数组合在以前的语法。lbqtest返回输出参数组合对应的输入参数。

一些选项测试的数量进行控制。当下列条件适用lbqtest进行多个测试:

  • lbqtest将每个测试视为独立于所有其他测试。

  • 如果您指定res,所有的输出都是向量。

  • 如果您指定资源描述的每一行StatTbl包含相应的测试的结果。

例如,lbqtest(资源描述,DataVariable = " ResidualGDP ",α= 0.025,滞后= (1 - 4))进行两次测试,在0.025的显著性水平,存在残差自相关的变量ResidualGDP表的资源描述。第一个测试包括1滞后的检验统计量,第二个测试包括4滞后。

例子

全部折叠

测试剩余的时间序列自相关的使用默认选项lbqtest。输入时间序列数据作为数字向量。

加载马克/英镑汇率数据集。

负载Data_MarkPound

数据是一个时间序列向量的日常马克/英镑双边现货汇率。

情节。

情节(数据)标题(“\男朋友马克/英镑双边即期汇率”)ylabel (“即期汇率”)包含(“自1984年1月2日,工作日”)

图包含一个坐标轴对象。坐标轴对象与标题空白D e u t s c h m r k / B r我t s h双边即期汇率英镑,包含工作日自1月2日,1984年,ylabel即期汇率包含一个类型的对象。

该系列出现不稳定。

稳定系列,现货汇率转换为回报。

回报= price2ret(数据);情节(回报)标题(“\男朋友马克/英镑双边即期汇率”)ylabel (“回归”)包含(“自1984年1月3日,工作日”)

图包含一个坐标轴对象。坐标轴对象与标题空白D e u t s c h m r k / B r我t s h双边即期汇率英镑,包含工作日自1月3日,1984年,ylabel返回包含一个类型的对象。

返回系列的计算偏差的意思。

残差=回报-意味着(回报);

在0.05水平的意义,测试剩余系列自相关使用默认选项Ljung-Box Q-test。

h = lbqtest(残差)
h =逻辑0

结果h= 0表示拒绝零假设的存在证据不足没有通过残差自相关20滞后。

加载马克/英镑汇率数据集。

负载Data_MarkPound

预处理的数据遵循这个过程:

  1. 稳定的系列计算每日的回报。

  2. 计算偏离均值回归。

回报= price2ret(数据);残差=回报-意味着(回报);

测试剩余系列从1到20滞后显著相关。回归测试的决定, p 值,测试数据,和重要价值。

[h, pValue,统计,cValue] = lbqtest(残差)
h =逻辑0
pValue = 0.1131
统计= 27.8445
cValue = 31.4104

测试一个时间序列,这是表中的一个变量,对残差自相关使用的默认选项lbqtest

加载股指数据集Data_EquityIdx。预处理每日纳斯达克收盘价格通过执行以下操作:

  1. 转换价格系列通过使用一系列回报率price2ret

  2. 代表了残差波动在一个恒定水平的系列定心的回报系列。

表中存储剩余系列与其它数据。因为price-to-return转换减少了样本大小的头系列,替换丢失的残余

负载Data_EquityIdxret = 100 * price2ret (DataTable.NASDAQ);res = ret -意味着(ret);数据表。Residuals_NASDAQ =[南;res];DataTable.Properties.VariableNames{结束}
ans = ' Residuals_NASDAQ '

剩余系列表中的最后一个变量。

开展Ljung-Box Q-test剩余序列在5%的显著性水平上通过提供整个数据集lbqtest

StatTbl = lbqtest(数据表)
StatTbl =表1×7h pValue stat cValue滞后α景深_____ __________交___ _____ ___测试1真正的2.8182 e-11 92.395 - 31.41 20 0.05 20

lbqtest返回表中的测试结果和设置StatTbl,变量对应于测试结果(h,pValue,统计,cValue)和设置(滞后,α,景深),行对应于个人测试(在这种情况下,lbqtest进行一个测试)。

h = 1pValue = 2.82 e-11拒绝零假设的证据,表明至少有一个显著的自相关滞后1到20日在纳斯达克返回剩余系列是强大的。

默认情况下,lbqtest测试表中最后一个变量。选择一个变量从一个输入表测试,设置DataVariable选择。

加载马克/英镑汇率数据集。

负载Data_MarkPound

转换价格的回报。

回报= price2ret(数据);

计算返回系列的偏差。

res =回报-意味着(回报);

测试的假设残余系列不是autocorrelated,使用默认的滞后。

h1 = lbqtest (res)
h1 =逻辑0

h1 = 0表明没有足够证据拒绝零假设的剩余工资回报不是autocorrelated。

测试的假设有显著的拱效应,使用默认的滞后[3]。

h2 = lbqtest (> ^ 2)
h2 =逻辑1

h2 = 1表明有明显的拱效应的残差的回报。

测试剩余的异方差性archtest。指定一个替代拱( l )模型, l 20. 的一致性,h2

h3 = archtest (res,滞后= 20)
h3 =逻辑1

h3 = 1表明没有残余异方差性应该拒绝零假设的支持一个拱( l )模型, l 20. 。这个结果是一致的h2

为自相关进行多个Ljung-Box Q-tests指定几个滞后测试统计。数据集是一个时间序列的连续57天overshorts从一个地下油箱在科罗拉多州[2]。也就是说,当前overshort ( y t )代表测量的燃料量的准确性:

  • 坦克在最后的一天 t

  • 坦克在最后的一天 t - - - - - - 1

  • 交付给坦克天 t

  • 天卖 t

加载数据集。

负载Data_OvershortT =身高(数据表);图绘制(DataTable.OSHORT)标题(“每日汽油Overshorts”)

图包含一个坐标轴对象。坐标轴对象与标题每日汽油Overshorts包含一个类型的对象。

lbqtest适合一系列以一个恒定的意思。因为这个系列似乎在不断波动的意思是,你不需要稳定。

计算偏差的意思。

数据表。Residuals_OSHORT =数据表。OSHORT——意味着(DataTable.OSHORT);

评估是否autocorrelated残差。包括5、10和15落后检验统计量,并调整每个测试的显著性水平0.05 / 3

StatTbl = lbqtest (DataTable, DataVariable =“Residuals_OSHORT”,滞后=(5 10 15),α= 0.05 / 3)
StatTbl =表3×7h pValue stat cValue滞后α景深__ ________ ___ _____ _____交测试1真的0.016667 0.0016465 19.36 13.839 5 5测试2真的0.016667 0.00068328 30.599 21.707 10 10测试3真的0.016667 0.001281 36.964 28.88 15 15

StatTbl包含单独测试的结果为每个指定的延迟进行。每个测试拒绝零假设在0.0167水平的意义。

推断出残差估计ARIMA模型,评估是否残差自相关使用lbqtest

加载澳大利亚消费者价格指数(CPI)数据集。时间序列(消费者价格指数)是日志季度消费者价格指数从1972年到1991年。删除趋势系列通过第一个区别。

负载Data_JAustraliancpi = DataTable.PAU;T =长度(cpi);3 = diff (cpi);dt = datetime(日期、ConvertFrom =“datenum”);图绘制(dt (2: T), 3)标题(“澳大利亚CPI差”)包含(“年”)ylabel (“消费者价格指数增长率”)轴

图包含一个坐标轴对象。坐标轴对象与标题差澳大利亚CPI包含,ylabel CPI增长包含一个类型的对象。

差系列似乎静止。

适合的AR(1)模型系列,然后推断出残差估计模型。

Mdl = arima (1,0,0);EstMdl =估计(Mdl, 3);
ARIMA(1,0,0)模型(高斯分布):价值StandardError TStatistic PValue _____ _________________ __________ __________常数0.015564 0.0028766 5.4106 6.2808 e-08 AR{1}方差0.29646 0.11048 2.6834 0.0072876 0.0001038 - 1.1932 8.6994 - 3.3362 e-05 e-18
res =推断(EstMdl, 3);stdRes = res /√(EstMdl.Variance);%标准化残差

评估是否残差由进行Ljung-Box autocorrelated Q-test。标准化残差源自估计模型(EstMdl包含参数)。当使用这样的残差,执行以下操作:

  • 调整自由度(景深检验统计量的分布占估计参数。

  • 滞后的数量设置为包含在测试统计。

  • 当你计算估计参数,跳过常数和方差参数。

滞后= 10;景深= - 1落后;%一个自回归参数[h, pValue] = lbqtest (stdRes滞后=滞后,景深=景深)
h =逻辑1
pValue = 0.0119

pValue = 0.0119表明残差有明显的自相关的至少一个滞后落后1到5,在5%的水平。

输入参数

全部折叠

剩余系列,指定为一个数值向量。的每个元素res对应于一个观察。

通常情况下,res包含(标准化)的残差模型适合观测时间序列。

数据类型:

时间序列数据,指定为一个表或时间表。每一行的资源描述是一个观察。

指定单个剩余系列(变量)来测试使用DataVariable论点。所选变量必须是一个数字。

请注意

指定失踪观察使用。的lbqtest处理缺失值的函数失踪的完全随机

名称-值参数

指定可选的双参数作为Name1 = Value1,…,以=家,在那里的名字参数名称和吗价值相应的价值。名称-值参数必须出现在其他参数,但对的顺序无关紧要。

R2021a之前,用逗号来分隔每一个名称和值,并附上的名字在报价。

例子:lbqtest(资源描述,DataVariable = " ResidualGDP ",α= 0.025,滞后= (1 - 4))进行两次测试,在0.025的显著性水平,存在残差自相关的变量ResidualGDP表的资源描述。第一个测试包括1滞后的检验统计量,第二个测试包括4滞后。

数量的滞后l包括检验统计量,小于指定为一个正整数T或一个向量的正整数,T是有效的样本大小(常数的值输入系列)。

lbqtest为每个元素进行单独测试滞后

例子:滞后= [1 - 4]进行两个测试。第一个测试只包括第一滞后的AR模型残差平方,和第二个测试包括第一到第四滞后。

数据类型:

假设检验的显著性水平,指定为一个数字标量的间隔(0,1)或数字向量的值。

lbqtest为每个值进行单独测试α

例子:α= (0.01 - 0.05)使用一个水平的意义0.01在第一个测试中,然后使用一个水平的意义0.05第二个测试。

数据类型:

自由度的检验统计量的渐近卡方分布在零假设下,指定为一个正整数或向量的正整数。

lbqtest为每个值进行单独测试景深

如果景深是一个整数,那么它必须小于或等于滞后。否则,的每个元素景深必须小于或等于相应的元素滞后

例子:景深= 15指定15自由度检验统计量的分布。

数据类型:

变量资源描述测试,指定字符串标量或字符向量包含一个变量名Tbl.Properties.VariableNames,或者一个整数或逻辑向量代表指数的一个名字。

例子:DataVariable = " ResidualGDP "

例子:DataVariable =(假真假假)DataVariable = 2测试第二个表变量。

数据类型:|逻辑|字符|字符串

输出参数

全部折叠

测试拒绝的决定,作为一个逻辑返回标量或矢量长度相等数量的测试。lbqtest返回h当你提供输入res

  • 的值1表明拒绝零假设没有残差自相关的选择。

  • 的值0显示失败拒绝零假设无残差自相关。

检验统计量p值,作为数字返回标量或矢量长度等于测试的数量。lbqtest返回pValue当你提供输入res

测试统计数据,作为数字返回标量或矢量长度相等数量的测试。lbqtest返回统计当你提供输入res

测试关键值,由α,作为一个数字返回标量或矢量长度相等数量的测试。lbqtest返回cValue当你提供输入res

测试总结,作为一个表返回的变量输出h,pValue,统计,cValue为每个测试,行。lbqtest返回StatTbl当你提供输入资源描述

StatTbl包含指定的测试设置的变量滞后,α,景深

更多关于

全部折叠

Ljung-Box Q-Test

的Ljung-Box Q-test是一个“混合”测试,评估一系列的零假设残差的展品没有固定数量的滞后的自相关l(见滞后),对一些自相关系数的选择ρ(k),k= 1,…,l是零。

测试数据

= T ( T + 2 ) k = 1 l ( ρ ( k ) 2 ( T k ) ) ,

在哪里T是样品的尺寸,l是自相关的数量落后,ρ(k)在滞后,样本自相关吗k。零假设下的渐近分布是卡方l的自由度。

失踪的完全随机

观察一个随机变量失踪的完全随机如果一个观察失踪的趋势是独立随机变量和所有其他观测失踪的趋势。

提示

如果你获得输入剩余系列数据拟合模型,减少了自由度景深估计系数的数量,不包括常数。例如,如果您获得输入残差拟合ARMA (p,)模型,设置景深=lp,在那里l的值是滞后

算法

  • 的值滞后论点l影响的力量测试。

    • 如果l太小了,测试没有检测到高阶的自我。

    • 如果l测试失去权力太大,当一个显著相关滞后被无关紧要的相关性在其他落后淘汰。

    • 盒子,詹金斯,Reinsel建议默认设置滞后=分钟(20 t - 1)[1]

    • -蔡引用模拟证据显示更好的动力性能进行测试l大约是日志(T)[5]

  • lbqtest不直接测试序列自相关以外的依赖关系。然而,你可以用它来识别条件异方差(ARCH效应)通过测试方残差[4]

    恩格尔的测试评估的意义直接拱效应。有关详细信息,请参见archtest

引用

[1]盒子,乔治·e·P。,Gwilym M. Jenkins, and Gregory C. Reinsel.时间序列分析:预测与控制。第三。恩格尔伍德悬崖,新泽西:普伦蒂斯霍尔,1994年。

[2]Brockwell, p . j . r·a·戴维斯。介绍了时间序列和预测。第二版。纽约,纽约州:施普林格,2002年。

[3]Gourieroux C。拱模型和金融应用程序。纽约:斯普林格出版社,1997年版。

[4]麦克劳德,人工智能和w·k·李。“诊断检查ARMA时间序列模型使用Squared-Residual自我。”Journal of Time Series Analysis. Vol. 4, 1983, pp. 269–273.

[5]-蔡,r S。金融时间序列的分析。第二版,霍博肯,台北:约翰·威利& Sons Inc ., 2005年。

版本历史

之前介绍过的R2006a