问一下:
如果我用bootstrap方法做ols回归?模拟一下回归结果怎么搞?
自己搞了一个,不对。
各位大侠看看我哪里出错了?
具体步骤如下:
close all;
clear all;
%% 拟建立一个线性方程
% y=0.8x
%% 生成标准正态分布随机数
rn=randn(100,1);
x=0.1:0.1:10;x=x';
y=0.8*x+rn;
%% ols
beta0=inv(x'*x)*x'*y;
%% bootstrap: bootstrap error term
beta=inv(x'*x)*x'*y;
e=y-beta*x;
for i=1:20000
eb=bootrsp(e);
yb=beta*x+eb;
beta1(i)=inv(x'*x)*x'*yb;
beta=beta1(i);
end
%%求其均值作为拟合的结果
beta1_hat=mean(beta1);
是方法错了?还是程序错了?