Looping outputs URGENT please help

1 回表示 (過去 30 日間)
Robert
Robert 2015 年 2 月 26 日
コメント済み: dpb 2015 年 2 月 26 日
I have the following code. As you can see it sorts some data, then performs a stepwise regression with the data. I need the output for each loop printed, preferably all on the same sheet (B is the only output I want). You can see where I have made an attempt under %Get Output. I am a total matlab newbie and horrible at it so as much detail as possible would be appreciated. Thank you in advance for any help!
%Sort the Data [rows, cols]=size(A03); for i = 1:rows
for j = 1:cols
y=zeros(8,1)
x=zeros(8,3)
for m=1:8
y(m)=NDVI((m-1)*rows+i,j)
x(m,1)=PET((m-1)*rows+i,j)
x(m,2)=P((m-1)*rows+i,j)
x(m,3)=Amp((m-1)*rows+i,j)
end
%Run the stepwise Regression
if sum(sum(isnan(x)),2)==0 && sum(isnan(y))==0
xx=x(:,1:3); yy=y(:,1);
[B,SE,PVAL,INMODEL,STATS,NEXTSTEP,HISTORY]= ...
stepwisefit(xx,yy,'penter',.05);
%Get Outputs
B = zeros(size(A03));
[B(1:1, 1:3, :)] = stepwisefit(xx,yy,'penter',.05);
end
end
end
  3 件のコメント
Robert
Robert 2015 年 2 月 26 日
Thanks for your help! I didn't delete the question, not sure what you mean...
dpb
dpb 2015 年 2 月 26 日
Hmmm....somebody did; the same question w/o the "URGENT" label disappeared between the time I first responded but was gone by time tried to Submit...if it wasn't you who did it, sorry.

サインインしてコメントする。

採用された回答

dpb
dpb 2015 年 2 月 26 日
...
[B,SE,PVAL,INMODEL,STATS,NEXTSTEP,HISTORY]= ...
stepwisefit(xx,yy,'penter',.05);
%Get Outputs
B = zeros(size(A03));
[B(1:1, 1:3, :)] = stepwisefit(xx,yy,'penter',.05);
...
Above you computed the coefficients and then immediately wiped them out and then do the work all over again...what's the point in that??? :)
If all you're wanting is to output the coefficient array, then
[B,SE,PVAL,INMODEL,STATS,NEXTSTEP,HISTORY]= ...
stepwisefit(xx,yy,'penter',.05);
dlmwrite('file.txt',B,'delimiter',' ','precision','%14.8f','-append')
Adjust format to suit and obviously use a better file name and it would be good to also add some sort of identifier as to case number or something to identify who's who in the zoo of models.

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeDates and Time についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by