📄 qstatlb.m
字号:
function [qstat, msig, autorho] = qstatlb(zz, p,k);
% Leung-Box Q Statistic for Regression Residual
% Input: zz data, p: order of lags, k: number of regressors
% Output: qstat, msig, autocorrelations
zzlag = mylag(zz, p);
zzz = [zz(p+1:end,:) zzlag];
junk = cov(zzz);
zzz = zz(p+1,end,:);
zzsq = zzz .^ 2;
sumzzsq = sum(zzsq);
autorho = junk(1,2:p+1) ./ junk(1,1);
rhosq = autorho .^2;
[tt,cc] = size(zz);
for j = 1:p, ttjj(j) = tt - j; end;
rhosqj = rhosq ./ ttjj;
qstat = tt * (tt+2) * sum(rhosqj);
msig = 1 - chi2cdf(qstat, p - k);
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -