clean up get_ar_ec_matrices
parent
397a850ce5
commit
1929a07d71
|
@ -138,25 +138,15 @@ for i = 1:length(rhsvars)
|
||||||
ndiffs = ndiffs - 1;
|
ndiffs = ndiffs - 1;
|
||||||
end
|
end
|
||||||
for k = 0:ndiffs
|
for k = 0:ndiffs
|
||||||
if mod(k, 2) == 0
|
|
||||||
oo_.var.(var_model_name).ar(i, rhsvars{i}.arRhsIdxs(j), lag + k) = ...
|
oo_.var.(var_model_name).ar(i, rhsvars{i}.arRhsIdxs(j), lag + k) = ...
|
||||||
oo_.var.(var_model_name).ar(i, rhsvars{i}.arRhsIdxs(j), lag + k) + nchoosek(ndiffs,k) * g1(i, g1col);
|
oo_.var.(var_model_name).ar(i, rhsvars{i}.arRhsIdxs(j), lag + k) + (-1)^k * nchoosek(ndiffs,k) * g1(i, g1col);
|
||||||
else
|
|
||||||
oo_.var.(var_model_name).ar(i, rhsvars{i}.arRhsIdxs(j), lag + k) = ...
|
|
||||||
oo_.var.(var_model_name).ar(i, rhsvars{i}.arRhsIdxs(j), lag + k) - nchoosek(ndiffs,k) * g1(i, g1col);
|
|
||||||
end
|
|
||||||
end
|
end
|
||||||
elseif rhsvars{i}.ecRhsIdxs(j) > 0
|
elseif rhsvars{i}.ecRhsIdxs(j) > 0
|
||||||
% Fill EC
|
% Fill EC
|
||||||
[lag, ndiffs] = findLagForVar(var, -rhsvars{i}.lags(j), 0, ecRhsVars);
|
[lag, ndiffs] = findLagForVar(var, -rhsvars{i}.lags(j), 0, ecRhsVars);
|
||||||
for k = 0:ndiffs
|
for k = 0:ndiffs
|
||||||
if mod(k, 2) == 0
|
|
||||||
oo_.var.(var_model_name).ec(i, rhsvars{i}.ecRhsIdxs(j), lag + k) = ...
|
oo_.var.(var_model_name).ec(i, rhsvars{i}.ecRhsIdxs(j), lag + k) = ...
|
||||||
oo_.var.(var_model_name).ec(i, rhsvars{i}.ecRhsIdxs(j), lag + k) + nchoosek(ndiffs,k) * g1(i, g1col);
|
oo_.var.(var_model_name).ec(i, rhsvars{i}.ecRhsIdxs(j), lag + k) + (-1)^k * nchoosek(ndiffs,k) * g1(i, g1col);
|
||||||
else
|
|
||||||
oo_.var.(var_model_name).ec(i, rhsvars{i}.ecRhsIdxs(j), lag + k) = ...
|
|
||||||
oo_.var.(var_model_name).ec(i, rhsvars{i}.ecRhsIdxs(j), lag + k) - nchoosek(ndiffs,k) * g1(i, g1col);
|
|
||||||
end
|
|
||||||
end
|
end
|
||||||
else
|
else
|
||||||
error('Shouldn''t arrive here');
|
error('Shouldn''t arrive here');
|
||||||
|
|
Loading…
Reference in New Issue