matlab - 如何返回 Akaike Criteria 选择的 VAR(2) 模拟中选择的最佳滞后的百分比频率

Question

如果我在模拟 3 维 VAR(2) 模型之前放置一个计数器循环，并且在计算之后假设 AIC 和滞后 3 被选为最佳拟合滞后，那么我很难计算最佳滞后的百分比频率（或其他单词在 100 次模拟中选择了多少次这个最佳滞后）？

Answer 1

在您的循环中，将当前模拟的滞后保存在向量中。例如：

sims_count = 100;
lags = NaN(sims_count,1);

for i = 1:sims_count
    % your simulation here...
    lags(i) = simulation_lag;
end

然后，一旦迭代结束并且您的lags向量已被填充，为了计算与最佳滞后匹配的滞后频率（等于3），只需使用以下代码：

best_lag_freq = sum(lags == 3) / sims_count;

一个实用的独立示例：

lags = randi(3,100,1);
best_lag_freq = sum(lags == 3) / 100;