uqlab

% LR
InputOpts.Marginals(1).Name = 'LR';
InputOpts.Marginals(1).Type = 'Uniform';
InputOpts.Marginals(1).Moments  = [3 20];

% E2
InputOpts.Marginals(2).Name = 'RT';
InputOpts.Marginals(2).Type = 'Uniform';
InputOpts.Marginals(2).Moments  = [20.1 30];

% E3
InputOpts.Marginals(3).Name = 'A1';
InputOpts.Marginals(3).Type = 'Uniform';
InputOpts.Marginals(3).Moments  = [0 90];

InputOpts.Marginals(4).Name = 'A2';
InputOpts.Marginals(4).Type = 'Uniform';
InputOpts.Marginals(4).Moments  = [0 90];

InputOpts.Marginals(5).Name = 'E1';
InputOpts.Marginals(5).Type = 'Gaussian';
InputOpts.Marginals(5).Moments = [102000, 10200];
% InputOpts.Marginals(5).Cov = 0.1;
myInput = uq_createInput(InputOpts);

MetaOpts.Type = 'Metamodel';
MetaOpts.MetaType = 'PCE';

MetaOpts.ExpDesign.X = X;
MetaOpts.ExpDesign.Y = Y;

MetaOpts.Degree = 2:10;

MetaOpts.ValidationSet.X = Xval;
MetaOpts.ValidationSet.Y = Yval;

myPCE = uq_createModel(MetaOpts);

uq_print(myPCE)

YPCE = uq_evalModel(myPCE,Xval);
