全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
22426 8
2010-12-07
悬赏 140 个论坛币 未解决
R语言实现LM检验
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2010-12-8 00:03:15
arch {vars}R Documentation
ARCH-LM testDescriptionThis function computes univariate and multivariate ARCH-LM tests for aVAR(p).
Usagearch(x, lags.single = 16, lags.multi = 5, multivariate.only = TRUE)
Arguments
xObject of class ‘varest’; generated byVAR(), or an object of class ‘vec2var’; generated byvec2var().
lags.singleAn integer specifying the lags to be used for theunivariate ARCH statistics.
lags.multiAn integer specifying the lags to be used for themultivariate ARCH statistic.
multivariate.onlyIf TRUE (the default), onlythe multivariate test statistic is computed.
DetailsThe multivariate ARCH-LM test is based on the following regression(the univariate test can be considered as special case of theexhibtion below and is skipped):

vech(hat{u}_t hat{u}_t') = β_0 + B_1vech(hat{u}_{t-1} hat{u}_{t-1}') + ... + B_qvech(hat{u}_{t-q} hat{u}_{t-q}' + v_t)

whereby v_t assigns a spherical error process andvech is the column-stacking operator for symmetric matricesthat stacks the columns from the main diagonal on downwards. Thedimension of β_0 is frac{1}{2}K(K +1) and forthe coefficient matrices B_i with i=1, ..., q,frac{1}{2}K(K +1) times frac{1}{2}K(K +1). The nullhypothesis is: H_0 := B_1 = B_2 = ... = B_q = 0 and thealternative is: H_1: B_1 neq 0 or B_2 neq 0 or ... B_q  neq0.The test statistic is:

VARCH_{LM}(q) = frac{1}{2}T K (K + 1)R_m^2 quad ,

with

R_m^2 = 1 - frac{2}{K(K+1)}tr(hat{Omega} hat{Omega}_0^{-1})quad ,

and hat{Omega} assigns the covariance matrix of the abovedefined regression model. This test statistic is distributed aschi^2(qK^2(K+1)^2/4).
ValueA list with class attribute ‘varcheck’ holding thefollowing elements:

residA matrix with the residuals of the VAR.
arch.uniA list with objects of class ‘htest’containing the univariate ARCH-LM tests per equation. This elementis only returned if multivariate.only = FALSE is set.
arch.mulAn object with class attribute ‘htest’containing the multivariate ARCH-LM statistic.
Author(s)Bernhard Pfaff
ReferencesDoornik, J. A. and D. F. Hendry (1997), Modelling DynamicSystems Using PcFiml 9.0 for Windows, International ThomsonBusiness Press, London.
Engle, R. F. (1982), Autoregressive conditional heteroscedasticitywith estimates of the variance of United Kingdom inflation,Econometrica, 50: 987-1007.
Hamilton, J. (1994), Time Series Analysis, PrincetonUniversity Press, Princeton.
L
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-12-8 09:24:51
这个。。。楼上的基本就是抄help吧,呵呵
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2010-12-8 21:28:09
是啊,不过不是基本抄,简直就是完全的复制粘贴.希望能有帮助,钱是不要的.
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2015-9-15 22:10:33
用专门做金融时间序列的FinTS包,其中ArchTest()函数就是用来做Arch LM检验的
例子:
dat <- ArchTest(x, lags=12, demean = FALSE)
参数解释:
x                numeric vector  #数值型向量
lags                positive integer number of lags #滞后阶数
demean        logical: If TRUE, remove the mean before computing the test statistic. #是否做中心化处理
可以参考蔡瑞胸的金融时间序列分析一书哦~
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-3-20 08:52:13
Sunshine飞鸟 发表于 2015-9-15 22:10
用专门做金融时间序列的FinTS包,其中ArchTest()函数就是用来做Arch LM检验的
例子:
dat
那你说的,对面板数据怎么检验啊?解释变量有四五个呢
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群