How can I repeat the following text 5000 times ?

조회 수: 1 (최근 30일)
Fox
Fox 2015년 12월 22일
편집: goerk 2015년 12월 23일
Hello, I make some wild bootstrap here and I know want to repeat the following text(formulas)as a whole 5000 times:
% Estimate the p-values by bootstrapping the residuals.
e=[e1,e2];
% resample the residuals with replacement in pairs
res=datasample(e,743,'Replace',true);
% create a Random Variable which follows the Rademacher distribution
rd = ((rand(1,743)<.5)*2 - 1)'; %
% multiply the res with a random variable of v Rademacher distribution
resrd1=res(:,1).*rd;
resrd2=res(:,2).*rd;
% Time Series of Predictors x(t)and r(t)
% computing unconditional mean a2/(1-b), a2=b(1), b=b(2)and build
% time series x(t)* = x2
b(1)/(1-b(2))
x2=b(1)+ b(2)*(b(1)/(1-b(2)))+resrd2;
% build time series for r(t)*=rt
rt=a+resrd1;
% make regressions with the before constructed x2 and rt
% getting with OLS ß*=beta
x3=[ones(size(rt)),x2]; % to have an intercept in the regression
betas=regress(rt,x3);
How can I do this in Matlab ?
  댓글 수: 1
Fox
Fox 2015년 12월 22일
How can I save the 5000 different betas ?

댓글을 달려면 로그인하십시오.

채택된 답변

goerk
goerk 2015년 12월 22일
Put the code in a scriptfile (e.g. doThisCode.m) and call it from anotherone with.
for i=1:5000;
doThisCode
end
This will lead to the same result as if you press F5 5000 times.
  댓글 수: 2
Fox
Fox 2015년 12월 22일
Thanks. Do you know how I can save my 5000 solutions for the betas in a matrix ?
goerk
goerk 2015년 12월 23일
편집: goerk 2015년 12월 23일
for a scalar value of beta without preallocation:
beta_vec = [];
for i = 1:5000
doThisCode
beta_vec = [beta_vec beta];
end

댓글을 달려면 로그인하십시오.

추가 답변 (0개)

카테고리

Help CenterFile Exchange에서 Resampling Techniques에 대해 자세히 알아보기

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by