在MATLAB中生成异方差级数？

Question

我正在尝试在MATLAB中生成异方差错误项。

我找到了一个很好的链接来帮助here，我遇到的问题是在MATLAB中复制它。

这是我的尝试：

n = [(1:100) (1:100)]';
sigma2 = n.^(1.3);
t = size(n,1);
for i = 1:200
    eps(i) = normrnd(0, sqrt(sigma2(1)));
end
eps = eps'
h = archtest(eps)

但是，针对异质性的测试表明我仍然没有异方差数据，任何人都可以看到我要去哪里。

Answer 1

好像您将sigma2固定在循环内的第一个值中。

替换

eps(i) = normrnd(0, sqrt(sigma2(1)));

与

eps(i) = normrnd(0, sqrt(sigma2(i)));