1 / 28

Design of Experiments ------Morris Method

Design of Experiments ------Morris Method. Yaomin Jin 01-03-2002. Outline of the presentation. Introduction of screening technique Morris method Examples Conclusions. Model. Input factors x 1 x 2 x 3 … x k. Output y=f(x 1 ,x 2 ,…,x k ). Which factor is important?.

alanreed
Download Presentation

Design of Experiments ------Morris Method

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Design of Experiments------Morris Method Yaomin Jin 01-03-2002

  2. Outline of the presentation • Introduction of screening technique • Morris method • Examples • Conclusions

  3. Model Input factors x1 x2 x3 … xk Output y=f(x1,x2,…,xk) Which factor is important? Screening technique • large-scale models • requirement of considerable computer time for each run • depend on a large number of input variables

  4. Morris method(1991) • OAT(one factor at a time) • the baseline changes at each step wanders in the input factors space • Estimate the main effect of a factor by computing r number of local measures at x1,x2,…xr in the input space then take average.

  5. Elementary effects • Reduce the dependence on the specific point that a local experiment has. • Determine which factor have: • negligible effects • linear and additive effects • non-linear and interaction effects

  6. Elementary effects • k dimensional factors vector x for the simulation model has components xi that have p-values in the set {0, 1/(p-1),…,1} • The region of the experiment is a k dimensional p level grid. In practical applications, the values sampled in are subsequently rescaled to generate the actual values of the simulation factors. • Δ=1/(p-1).

  7. Elementary effects of i-th factor at given point x where x is any value in selected such that the perturb point is still in . A finite distribution Fi of elementary effects for the i-th input factor is obtained by sampling x from . The number of elements of each Fi is

  8. Economy of Morris method • In the simplest form, the total computational effect required for a random sample of r values from each distribution Fi is n=2rk runs. Each elementary effect requires the evaluation of y twice. • The simplest form of Morris design has an economy rk/2rk=1/2.

  9. Economical design Based on the construction of a matrix B* with rows that represent input vectors x, for which the correspondingexperiment provides k elementary effects (one for each input factor) from k+1 runs. Economy of the design is k/(k+1). assume that p is even and , each of the elementary effects for the i-th input factor has equal probability of being selected. The key idea is: Base value x* is randomly chosen from the vector x, each component xi being sampled from the set One or more of the k components of x* are increased by such that vector x(1) still in

  10. The estimated elementary effect of the i-th component of x(1) (if the i-th component of x(1) has been changed by )if x(1) increased by Δ;if x(1) decreased by Δ. Let x(2) be the new vector , select a third vector x(3) such that differs from x(2) only one component j: Economical design (continue)

  11. Economical design (continue) else Repeat the upper step get the k+1 input vectors x(1),x(2),…,x(k+1) , any component i of x* is selected at least once to be increased by . To estimate one elementary effect for each factor.

  12. Economical design (continue) The rows of orientation matrix B* are the vectors describe above. This provides a single elementary effect per factor. To build B*, Restrict attention: a. p is even; b. Firstly, selection of sampling matrixB with elements that are 0 or 1, such that every column there are two rows of that differ in only one element.

  13. Economical design (continue) In particular, B may be chosen to be a strictly lower triangular matrix of 1, consider B’ given by

  14. Economical design (continue) Jk+1,1 is a matrix of 1. as a design matrix(i.e. each row a value for x) x* is randomly chosen base value of x. B’ could be used as a design matrix. Each element is randomly assigned a value from with equal probability. Since it would provide k elementary effects; one effect each input factor, with a computation cost of k+1 runs. However the problem is that the k elementary effects B’ produces would not be randomly selected.

  15. Economical design (continue) A randomised version of the design matrix is given by where • D* is diagonal matrix in which each diagonal element is either +1 or -1 • P* is random permutation matrix, in which each column contains one element equal to 1 and all the others equal to 0, and no two columns have 1’s in the same position • B* provides one elementary effect per factor that is randomly selected.

  16. Suppose that p=4, k=4 and , that is, four factors that may have values in the set {0,1/3,2/3, 1}. Then B5*4 is given by Example

  17. Example (continue) and the randomly generated x*, D* and P* happen to be

  18. Example (continue)

  19. Example (continue)

  20. Elementary effects To estimate the mean and variance of the distribution Fi(i=1,…,k), take a random sample of r elements; that is sample r mutually independent orientation matrices. Since each orientation matrix provides one elementary effect for every factor, the r matrices together provide r×k dimensional samples, one for each Fi(i=1,…,k). We use the classic estimate for every factor’s mean and standard deviation.

  21. Elementary effects The characterization of the distribution Fi through its mean and standard deviation gives useful information about the influence on the output; • a high mean indicates a factor with an important overall influence on the output, • a high standard deviation indicates either a factor interacting with other factors or a factor whose effect is non-linear.

  22. The lines constituting a wedge, are described by ; where is the standard deviation of the mean elementary effect. If the parameter has coordinates below the wedge, i.e. , this is a strong indication that the mean elementary effect of the parameter is non-zero. A location of the parameter coordinates above the wedge indicates that interaction effects with other parameters or non-linear effects are dominant. Standard of importance

  23. Example 1

  24. Result of example 2 (p=4, Δ=2/3 and r=4)

  25. Example 2 where wi=2(xi-0.5) except for i=3, wi=2(1.1xi/(xi+0.1)-0.5), otherwise. Coefficients of relatively large value were assigned as

  26. Result of example 2 (from graph) • Input 1-5 are clearly separated from the cluster of remaining outputs, which have means and standard deviations close to 0. • In particular, inputs 4,5 have mean elementary effects that are substantially different from 0 while having small deviations. • Consider both means and standard deviations together, we conclude that the first 5 inputs are important, and that of these the first three appear to have effects that involve either curvature or interactions. • This is coincide with the model.

  27. Conclusions • Economical for models with a large numbers of parameters • Does not dependent on any assumptions about the relationship between parameters and outputs • Results are easily explained in a graph • Drawback is not consider the dependencies between parameters(such as interactions)

  28. Question?

More Related