Unit: Multicriteria decision analysis Learning goals I. General - - PowerPoint PPT Presentation

unit multicriteria decision analysis learning goals
SMART_READER_LITE
LIVE PREVIEW

Unit: Multicriteria decision analysis Learning goals I. General - - PowerPoint PPT Presentation

Unit: Multicriteria decision analysis Learning goals I. General classification of strategies in Multi-objective Optimization based on the time of decision maker interaction. II. Ways to structure the decision making process III. What are


slide-1
SLIDE 1

Unit: Multicriteria decision analysis

slide-2
SLIDE 2

Learning goals

I. General classification of strategies in Multi-objective Optimization based on the time of decision maker interaction. II. Ways to structure the decision making process III. What are utility functions? How can we construct rankings?

  • IV. Pareto optima and related definitions

V. Interpretation and Visualization of Pareto fronts

slide-3
SLIDE 3

Multicriteria Decision Analysis

Kaiza Miettinen, Finnish Mathematician

Miettinen, Kaisa. Nonlinear multiobjective optimization. Springer, 1999.

slide-4
SLIDE 4

A priori multicriteria decision making

slide-5
SLIDE 5

Questions in MCDA

slide-6
SLIDE 6

Utility Function, Indifference Curves

better worse indifferent U (f1,f2)  2 U(f1,f2) 3 U(f1,f2) 4 f1 f2 f1 f2

slide-7
SLIDE 7

Multi-Attribute Utility Theory (MAUT)

Keeney, R. L. and Raiffa, H. (1976). Decisions with Multiple Objectives: Preferences and Value Tradeoffs. Wiley,New York. Reprinted, Cambridge Univ. Press, New York (1993). Keeney, R. L. (2009). Value-focused thinking: A path to creative decision making. Harvard Univ. Press.

slide-8
SLIDE 8

Two ways to define parameters of utility function

slide-9
SLIDE 9

Ordinal regression (example)

slide-10
SLIDE 10

Ordinal regression (solution)

slide-11
SLIDE 11

Ordinal regression

slide-12
SLIDE 12

Further Reading

  • Greco, Salvatore, Roman Słowiński, José Rui Figueira, and Vincent
  • Mousseau. "Robust ordinal regression." In Trends in multiple criteria

decision analysis, pp. 241-283. Springer US, 2010.

  • Greco, Salvatore, Vincent Mousseau, and Roman Słowiński.

"Ordinal regression revisited: multiple criteria ranking using a set of additive value functions." European Journal of Operational Research 191, no. 2 (2008): 416-436.

slide-13
SLIDE 13

Multi-Attribute Utility Theory (MAUT)

Keeney, R.L. (1992). Value-focused thinking—A Path to Creative Decision making. Harvard University Press.

slide-14
SLIDE 14

Client-Theory by Kahneman and Tversky

FKahneman: ‘Thinking fast and slow’ Penguin press. 2014.

slide-15
SLIDE 15

Desirability functions (DFs)

1.0

100 km/h 200 km/h 0 km/h

not acceptable fully satisfied

1000 Euro 2000 Euro 0 Euro

not acceptable totally satisfied

10 l/100km 15 l/100km 5 l/100km

not acceptable totally satisfied

1.0 1.0

v.d. Kuijl, Emmerich, Li: A robust multi-objective resource allocation scheme, Concurrency & Computation 22 (3), (2009)

slide-16
SLIDE 16

Desirability functions (DFs)

  • Harrington, E. C. "The desirability function." Industrial quality

control 21, no. 10 (1965): 494-498. (DFs type Harrington)

  • Derringer, G., & Suich, R. (1980). Simultaneous optimization of

several response variables. Journal of quality technology, 12(4), 214-219. (DFs Type Derringer-Suich)

  • Van der Kuijl, Emmerich, M., & Li, H. (2010). A robust multi‐objective

resource allocation scheme incorporating uncertainty and service

  • differentiation. Concurrency and Computation: Practice and

Experience, 22(3), 314-328. (DFs Type v.d.Kuijl et al.)

  • Trautmann, Heike, and Claus Weihs. Pareto-Optimality and

Desirability Indices. Universität Dortmund, 2004.

slide-17
SLIDE 17

Harrington vs. Derringer Suich type of DFs

n1 n3 n2 n1  n2  n3

slide-18
SLIDE 18

Harrington vs. Derringer Suich type of DFs (2)

n1 n3 n2 n1  n2  n3

slide-19
SLIDE 19

A posteriori multicriteria decision making

slide-20
SLIDE 20

Multicriteria Problems: Car Example

Speed Cost Your objectives: Cost  min, Speed  max Length Add constraint: Only red cars !

slide-21
SLIDE 21

Decision vs. Objective Space, pre-images

Speed(x)max Cost(x)min

.

X={Beetle, CV2, Ferrari}

Decision space Objective space

  • Def.:The space of candidate solutions is

called the decision space, for instance 𝑌 ={𝐶𝑓𝑓𝑢𝑚𝑓, 𝐷𝑊2, 𝐺𝑓𝑠𝑠𝑏𝑠𝑗} in the figure.

  • Def.: The space of objective function

values (vectors) is called the objective space, for instance ℝ2 in the figure.

  • Def.: For a point in the objective space

the corresponding point(s) in the decision space are called their pre- image(s). In the figure, Beetle and CV2 are pre-images of the green point.

  • Remark: Two different points in the

decision space (e.g. Beetle and CV2) can map to the same point in the

  • bjective space, but two points in the
  • bjective space have never the same

preimage in the decision space.

slide-22
SLIDE 22

Incomparable and indifferent solutions

Speed(x)max Cost(x)min

.

  • Def.: Given two solutions and some

criterion functions, the solutions are said to be incomparable, if and only if 1. the first solution is better than the second solution in one or more criterion function value and 2. the second solution is better than the first alternative in one or more

  • ther criterion function values.
  • Def.: Given two alternatives and some

criterion functions, the alternatives are said to be indifferent with respect to each other, if and only if they share exactly the same criterion function values. Remark: In both cases additional preference information is required to decide which solution is best.

incomparable indifferent

X={Beetle, CV2, Ferrari}

Decision space Objective space

slide-23
SLIDE 23

Pareto dominance

Speed(x)max Cost(x)min

.

  • Def.: Given two alternative solutions and

some objective functions, the first solution ist said to Pareto dominate the second solution, if and only if 1. the first solution is better or equal in all objective function values, and 2. the first solution is better in at least

  • ne objective function value.

Francis Y. Edgeworth Irish Economist 1845-1926 Vilfredo Pareto Italian Economist 1848-1923

slide-24
SLIDE 24

Pareto Dominance, Pareto front

Given: A decision space 𝑌 comprising all (feasible) decision alternatives, a number of criterion functions 𝑔𝑗 :𝑌 →ℝ , 𝑗 = 1, …, 𝑛 Def.: A decision alternative x in X dominates a solution x’ in X, iff it is not worse in each objective function value, and better in at least one objective function value. Def.: If a solution 𝑦∈𝑌 is not dominated by any other solution in X , then it is called Edgeworth-Pareto optimal (or Pareto optimal) (in X). Def.: The set of all Pareto optimal solutions in X is called efficient set. Def.: The set of all Pareto optimal function vectors of solutions in the efficient set is called the Pareto front.

Speed(x)max Cost(x)min X={Beetle, Limosine, BMW, Ferrari}

Efficient set

XE={Beetle, BMW, Ferrari}

Pareto front

= {(Speed(Beetle), Cost(Beetle))T, (Speed(BMW), Cost(BMW))T, (Speed(Ferrari), Cost(Ferrari))T}

slide-25
SLIDE 25

Dominating Subspace

Fundamental Concepts: Dominance diagram

f1 (min) f2(min)

Dominated Subspace Space of incomparable solutions Space of incomparable solutions Reference Solution y

What about the points on the dark red lines? How would this diagram look for maximization of f1 and minimization of f2? extended to infinity extended to - infinity

slide-26
SLIDE 26

Construction of Pareto front in 2-D

Geometrical construction for separating dominated and non dominated points in 2-D: 1. Indicate for each point the dominated subspace by shading 2. The covered subspace consists of dominated points within the set, that is points that are dominated by at least

  • ne other point

3. The outer corner points on the lower left boundary form the Pareto front of the point set.

f1 min f2min

Non-dominated solutions Pareto front Pareto dominated solutions

slide-27
SLIDE 27

Construction for 2-D continuous functions

x1 x2 x3 Search space (decision space) f2 Objective space (solution space) Image set f(X) Pareto front

f

f2 Efficient set

  • For f=(f1, ...,fm) the image set

f(X) is defined as the set of all (f1(x), ..., fm(x)) for xX.

  • The non-dominated solutions

in f(X) are located at the lower left boundary and form the Pareto front.

  • Note, that for unbounded or

non-closed sets f(X) the Pareto front does not always exist.

  • If it exists, the Pareto front of a

m-objective problem has at most m-1 dimensions.

  • The set of all preimages of

points in the Pareto front is the efficient set.

slide-28
SLIDE 28

Interpreting and visualizing Pareto fronts

slide-29
SLIDE 29

Thinking about trade-off, knee points

f2min

Trade-off curve (Pareto front)

f1min

e.g. Emissions e.g. Cost

f2 Image under f Pareto front

Moving from x1 to x2 is a balanced trade-off. x1 x2 x3 x4 x1 and x4 are ideal solutions Moving from x2 to x1 is a unbalanced tradeoff Region of good compromise Solutions (knee point region)

f1

slide-30
SLIDE 30

De Kruijf, Niek, et al. "Topological design of structures and composite materials with multiobjectives" Int. Journal of Solids and Structures 44.22 (2007) Deb et al.:. "Innovization: Innovating design principles through optimization" GECCO. ACM, 2006

Innovization – Design principles from multiobjective optimization

  • By looking at changing

designs across the Pareto front, designers can study design principle

  • How does the design change

when moving across the Pareto front?

  • Innovization: Finding design

principles by multicriteria

  • ptimization
slide-31
SLIDE 31

Construction of Pareto front for 3-D point set

  • A single in the objective space

point dominates a 3-D cuboid

  • the point is the upper corner
  • the lower corner is

(- ,-,- )T

  • Dominated cuboids can be

drawn in a perspectivic plot.

  • The boundary between

dominated and non-dominated space is called attainment surface

  • The points that belong to the

Pareto front are located at its

  • uter corners.
slide-32
SLIDE 32

Figure: Approximation to 3-D Pareto Front Here, minimization is the goal.

slide-33
SLIDE 33

Example Pareto front: embedded systems design

  • cost
  • latency

CPU utilization Configuration

  • f circuit is the same

but parameters change Li, R., Etemaadi, R., Emmerich, M. T., & Chaudron, M. R. (2011, June). An evolutionary multiobjective optimization approach to component-based software architecture design. In Evolutionary Computation (CEC), 2011 IEEE Congress on (pp. 432-439). IEEE.

slide-34
SLIDE 34

Figure: 4-D Pareto Front visualization. All functions to be maximized. VisuWei tool. http://natcomp.liacs.nl f4=0.00

slide-35
SLIDE 35

Parallel Coordinates Diagrams (PCDs) are a common way to visualize solutions with many (typically > 3) attributes (multivariate data) (1) Each criterion is represented by a vertical axis. (2) A solution is represented by means

  • f a polyline connecting points (criteria

values) of that solutions. Normalize the range of criteria, that is minimum, maximum of criterion in data set determine scale of its axis. Choose sign of the criteria can be chosen such that all criteria are to be minimized.

N-Dimensions: Parallel Coordinates Diagram

  • 2
  • 3
  • 3.8
  • 4

700 800 900 1000 2.4 2.2 2.5 2.6 3.0 3.2 low medium high Notebook3 Notebook2

Noise min Weightmin Pricemin

  • GHz min

Notebook1 Example: Notebook comparison with three notebooks

slide-36
SLIDE 36
  • 2
  • 3
  • 3.8
  • 4

700 800 900 1000 2.4 2.2 2.5 2.6 3.0 3.2 low medium high dill ifrog lightspeed powernote shark

Noise min Weightmin Pricemin

  • GHz min

N-Dimensions: Parallel Coordinates Diagram

Find pairs where A dominates B? How to determine indifference and incomparability? Use definition of Pareto dominance!

Normalize: All to be minimized Example: Comparing notebooks Ordinal scale

slide-37
SLIDE 37

Brushing is an interactive tool that can be used to explore sensitivity to constraints (1) Indicate ranges by gray area using, for instance, a computer mouse. (2) Highlighted solutions (polylines) are solutions that fall into the specified ranges. In the freeware tool XMDV the Parallel Coordinates diagram allows brushing, see example (gray area indicates selected ranges). Example (left): Filtering cars from Europe and in certain ranges of the continuous output variables.

Interactive Parallel Coordinates - Brushing tool

Brushing tool of XMDV Parallel Coordinates Example: Car data.

Europe USA Japan

slide-38
SLIDE 38

Star plots are a similar idea to visualize multivariate data. (1) Axis originate from the same point and are spread in an equi-angular way. (2) The scaling of the axis is from minimum to maximum. (3) Solutions are represented by closed polygons. Remarks:

  • Shapes are easier

perceived and remembered (so called Star-glyphs).

  • Brushing, selection and

grouping tools more difficult

Star Plots or Radar Plots

slide-39
SLIDE 39

Scatterplot Matrix - Example

  • The scatter plot matrix consists
  • f scatter plots for all values of

fi,fj with i < j and i  {1, …, m}, j={1, …,m}

  • Only the upper diagonal matrix

needs to be depicted

  • Points on the diagonal are not

interesting, but often the diagonal is used to plot total frequency of values

  • Scatterplots serves to

investigate the correlation between two objectives functions

X f1 f2 f3 x1 2 3 2 x2 1 1 2 x3 2 1 3

f1 f2 f2 f3 f1 f3

slide-40
SLIDE 40

Multivariate Data Visualization (MDV)

  • MDV tools
  • XMDV
  • Ggobi
  • Spreadsheet tools
  • e.g.MS Excel, Open Office
  • Scientific visualization

and programming tools

  • MATLAB/Octave
  • SciLab
  • Python/Matplotlib
  • R: Parallel Coordinates

https://www.safaribooksonline.com/blog/2014/03 /31/mastering-parallel-coordinate-charts-r/

slide-41
SLIDE 41

Scatter plot matrix

  • Criteria scatter plot matrix shows all pairwise scatter plots.
  • (Anti)correlation between criteria is revealed in this plot:
  • Positive slope  Criteria are positively correlated plot
  • Negative slope  Criteria are anticorrelated (conflicting)
slide-42
SLIDE 42

*Progressive and Interactive MCDA Methods

Interactive decision making tools work with questions to the decision maker, typically s/he is asked to do pairwise comparisons“ Preference elicitation: Deducting a utility function from comparisons and questionnaire data

  • Prometheus: Construction of a utility function from questionnaires
  • ELECTRE: Partial orders; resolve inconsistencies in ranking
  • AHP: Pairwise comparison, hierarchy of criteria => ranking
  • NIMBUS: Interactive industrial optimization and decision making

tool developed by the finnish group of Miettinen

  • Robust Ordinal Regression: Constructs possibility space of utility

function that are consistent with pairwise comparisons

  • NEMO: Combines robust ordinal regression with heuristic
  • ptimization
slide-43
SLIDE 43

Summary: Take home messages

1. In multioobjective optimization a priori, a posteriori, and progressive methods are distinguished, depending on when the DM interacts. 2. Multicriteria decision analysis structures decision process 3. Utilitiy functions capture user preferences: Linear weighting, Keeney Raiffa, Derringer-Suich type and Harrington Desirability functions. 4. In multicriteria optimization and decision analysis, two solutions can be incomparable, indifferent, or one solution dominates the other. 5. Pareto dominance and (Edgeworth-)Pareto optimality characterizes Pareto optima in multiobjective optimization. 6. Pareto fronts can be used to reason about trade-offs and also to find new design principles. They reveal the nature of the conflict(s). 7. Pareto fronts can be interpreted as trade-off curves (2-D) or as a trade-off surfaces in 3-D. For >= 4-D: Parallel coordinates, Radar plots, Scatter plots.