function out_star=boot_reg(rs)
global c cx q0 A n

E=A*rs;
SSE=sum(E.^2);
s2hat=SSE/n;
delta_star=cx'*rs;
T_star0=sqrt(n)*delta_star/sqrt(s2hat*q0);
out_star=[T_star0];

return