New directories and files introduced.

load ../../psData_Merged_n184 # load data file

ana_expert      # psuade> analysis expert mode on

printlevel 3    # psuade> printlevel set from 0 to 3

rscheck         # Check for best data-fitted surrogate model

10               # Select Gaussian Process Model (GPM) Tong implementation as response surface method

1               # Quantity of interest no

y               # Perform cross-cross_validation

120             # Enter the number of groups to validate : (2 - 120) 120 fold cross-validation

n               # Random selection of leave-out groups ? (y or n)

sys mv RSFA_CV_err.m R2_GPM10_184F_RSFA_CV_err.m # Rename the CV error file before next one

sys mv RSFA_training_err.m R2_GPM10_RSFA_training_err.m # Rename the training error file before next one

rs2

256             # Grid resolution ? (32 - 256)

10               # Available response surface tools:  select Gaussian Process Model

6               # Enter the input for x axis (1 - 7) :

5               # Enter the input for y axis (1 - 7), not 1 : 2

y               # Set other nominal values automatically ? (y or n)

1               # Enter output number (1 - 1) :

n               #  Ymin and Ymax found = 2.037723e-01 3.421112e-01. Set lower threshold ? (y or n) :

n               # Set upper threshold ? (y or n) : n

sys mv matlabrs2.m R2_GPM10_rs2_x6x5.m    # Rename the response surface constructed

rs2

256             # Grid resolution ? (32 - 256)

10               # Available response surface tools:  select Gaussian Process Model

6               # Enter the input for x axis (1 - 7) :

3               # Enter the input for y axis (1 - 7), not 1 : 2

y               # Set other nominal values automatically ? (y or n)

1               # Enter output number (1 - 3) :

n               #  Ymin and Ymax found = 2.037723e-01 3.421112e-01. Set lower threshold ? (y or n) :

n               # Set upper threshold ? (y or n) : n

sys mv matlabrs2.m R2_GPM10_rs2_x6x3.m    # Rename the response surface constructed

sys ls -alrt R2_GPM*  # Check most one created after mv

q

% This file contains total order Sobol' indices

% with error bars coming from bootstrapping.

% to select the most important ones to display,

% set sortFlag = 1 and set nn to be the number

% of inputs to display.

sortFlag = 0;

nn = 6;

Means = [

  1.5794846448076320e-02

  6.4653198143395704e-03

  2.6574227265898410e-02

  1.8567679906782156e-02

  4.8209734080710644e-02

  9.3791789323028563e-01

];

Stds = [

  7.1169553996481014e-03

  5.5558263793413446e-03

  8.9121938307062231e-03

  9.4909574446616109e-03

  1.1372238912996293e-02

  1.3435554885050963e-01

];

  Str = {'t1:P0','t2:beta','t3:StatWeight','t4:ep_g','t5:VelfacCoeff','x1:InitConc'};

hold off

if (sortFlag == 1)

  [Means, I2] = sort(Means,'descend');

  Stds = Stds(I2);

  I2 = I2(1:nn);

  Means = Means(1:nn);

  Stds = Stds(1:nn);

  Str  = Str(I2);

end

ymin = min(Means-Stds);

if ymin < 0 

    ymin = 0;

end;

ymax = max(Means+Stds);

h2 = 0.05 * (ymax - ymin);

bar(Means,0.8);

for ii = 1:nn

   if (ii == 1)

hold on

   end;

   XX = [ii ii];

   YY = [Means(ii)-Stds(ii) Means(ii)+Stds(ii)];

   if YY(1) < 0 

      YY(1) = 0;

   end;

   plot(XX,YY,'-ko','LineWidth',3.0,'MarkerEdgeColor','k','MarkerFaceColor','g','MarkerSize',12)

end;

set(gca,'linewidth',2)

set(gca,'fontweight','bold')

set(gca,'fontsize',12)

grid on

box on

axis([0  nn+1 ymin ymax])

set(gca,'XTickLabel',[]);

th=text(1:nn, repmat(ymin-0.05*(ymax-ymin),nn,1),Str,'HorizontalAlignment','left','rotation',90);

set(th, 'fontsize', 12)

set(th, 'fontweight', 'bold')

title('Total Order Sobol Indices (with bootstrap)','FontWeight','bold','FontSize',12)

ylabel('Total Order Sobol Index (Normalized)','FontWeight','bold','FontSize',12)

hold off