Skip to content

Commit

Permalink
More eficient plotting (max number of data)
Browse files Browse the repository at this point in the history
  • Loading branch information
cultpenguin committed May 11, 2009
1 parent bea2487 commit 5f11d6c
Showing 1 changed file with 7 additions and 7 deletions.
14 changes: 7 additions & 7 deletions krig_optim_mcmc.m
Original file line number Diff line number Diff line change
Expand Up @@ -219,10 +219,10 @@
t_old_plot=now;
subplot(2,1,1)
plotyy(1:nacc,L_acc,1:nacc,-be_acc);

nn=size(par2,1);nmax=400;ndd=ceil(nn/nmax);ii=[ndd:ndd:nn];
if size(par2,2)==1
% ONLY PLOT NMAX DATA
nn=size(par2,1);nmax=400;ndd=ceil(nn/nmax);ii=[ndd:ndd:nn];


subplot(2,3,4)
%plot(par2(:,1),L_acc,'k.')
Expand Down Expand Up @@ -256,15 +256,15 @@
drawnow;
elseif size(par2,2)==2
subplot(2,3,4)
scatter(par2(:,1),par2(:,2),22,L_acc,'filled')
scatter(par2(ii,1),par2(ii,2),22,exp(L_acc(ii)),'filled')
xlabel('Range 1');ylabel('Range 2');title('Likelihood')
%colorbar
subplot(2,3,5)
scatter(par2(:,1),par2(:,2),22,-be_acc,'filled')
xlabel('Range 1');ylabel('Range 2');title('-be')
%%subplot(2,3,5)
%%scatter(par2(:,1),par2(:,2),22,-be_acc,'filled')
%%xlabel('Range 1');ylabel('Range 2');title('-be')
%colorbar
subplot(2,3,6)
scatter3(par2(:,1),par2(:,2),nugfrac_acc,20,L_acc,'filled');
scatter3(par2(ii,1),par2(ii,2),nugfrac_acc(ii),20,exp(L_acc(ii)),'filled');
xlabel('Range 1');ylabel('Range 2');zlabel('Nugget Fraction');title('Likelihood')
drawnow;
elseif size(par2,2)>2
Expand Down

0 comments on commit 5f11d6c

Please sign in to comment.