Effacer les filtres
Effacer les filtres

Looping outputs URGENT please help

1 vue (au cours des 30 derniers jours)
Robert
Robert le 26 Fév 2015
Commenté : dpb le 26 Fév 2015
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 commentaires
Robert
Robert le 26 Fév 2015
Thanks for your help! I didn't delete the question, not sure what you mean...
dpb
dpb le 26 Fév 2015
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.

Connectez-vous pour commenter.

Réponse acceptée

dpb
dpb le 26 Fév 2015
...
[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.

Plus de réponses (0)

Catégories

En savoir plus sur Matrices and Arrays dans Help Center et File Exchange

Community Treasure Hunt

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

Start Hunting!

Translated by