Published On August 5, 2024
Journal Issue LJER Volume 24 Issue 5

Decisions with Uncertainty using Deterministic Analysis

Dr. Acosta-Flores Jos Jes s
Dr. Acosta-Flores Jos Jes s
Article Fingerprint
Research ID 404V6

IntelliPaper

Abstract

The purpose of this article is to present a simple algorithm for the analysis of decisions with multiple objectives whose measures of effectiveness are random variables. This paper discusses the possibility of using the probability functions: uniform, normal, exponential, Cauchy, Chi-square, Erlang, Gamma, and Laplace. The algorithm, based on the concepts and methodology of Decision Theory, guides the analyst so that he can interact with the decision- maker. First, the analyst asks the decision-maker to define his objectives in the problem he is going to analyze, as well as the measures of effectiveness to evaluate its achievement. Then, he asks questions to the decision-maker to determine his type of behavior: aversion, proneness or neutrality to risk, for each measure of effectiveness. Next, for each of them, calculate its utility function. At that point, he asks the decision-maker to specify the alternatives to be analyzed with their probability function, their range, their mean, and standard deviation.

Explore Digital Article Text

I. INTRODUCCIÓN

En muchos proyectos se tienen incertidumbre y varios objetivos, por lo que se requiere tener varias medidas de efectividad, al menos una para cada objetivo, y una función utilidad conjunta que evalúe las consecuencias de las decisiones posibles mediante sus medidas de efectividad. Con el fin de tener esa función utilidad hay necesidad de verificar diversos supuestos como la independencia preferencial, la independencia en utilidad y la independencia aditiva, (Keeney and Raiffa, 1976). Dependiendo del cumplimiento de esos supuestos se tienen funciones utilidad de diferentes tipos, que requieren respuestas a preguntas que pueden ser difíciles de contestar para un decisor.

Acosta (2019) desarrolla un algoritmo donde transforma el problema aleatorio en uno determinista y luego aplica el método ProAct para determinar la mejor decisión.

Para transformar el problema aleatorio en uno determinista, Acosta calcula para cada alternativa de decisión i y cada medida de efectividad j, el equivalente bajo certeza, , considerando donde es la función densidad de probabilidad de la medida de efectividad j para la alternativa i; y es la utilidad del decisor sobre la medida de efectividad j.

En el primer artículo (Acosta, 2019), se consideran la función densidad de probabilidad uniforme y la función utilidad exponencial.

En un segundo artículo (Acosta, 2021) utiliza además de la uniforme, las funciones normal y triangular.

Cuando se ha transformado el problema en determinista, se emplea el método ProAct (Hammond et al, 2000) que alterna los análisis de dominancia con los de permutas compensatorias.

En el análisis de dominancia se eliminan todas las opciones dominadas. Se considera que una opción A domina a B cuando todas las medidas de efectividad en A son mejores o iguales que las de B y al menos una de ellas es estrictamente mejor.

Para las permutas compensatorias se le pide al decisor que considere las variaciones en las medidas de efectividad que se deben hacer para que todas las alternativas tengan igual una de esas medidas. En cuanto se logra ello se elimina esa medida de efectividad que es igual para todas las alternativas.

Así, alternando dominancia (para eliminar alternativas) y permutas compensatorias (para eliminar medidas de efectividad) se llega, después de varias iteraciones, a alternativas evaluadas con una sola medida de efectividad. La que tenga el mejor valor en esa medida de efectividad será la mejor opción.

Se presenta a continuación la transformación del problema aleatorio en uno determinista considerando además de la uniforme y la normal, las funciones densidad de probabilidad Cauchy, Chi-cuadrada, Erlang, Gamma y Laplace.

II. TRANSFORMACIÓN DEL PROBLEMA ALEATORIO EN UNO DETERMINISTA

Para esta transformación se calcula, como se mencionó anteriormente, para cada alternativa de decisión i y cada medida de efectividad j, el equivalente bajo certeza, , considerando donde es la función densidad de probabilidad de la medida de efectividad j cuando se toma la alternativa i; es la utilidad del decisor sobre la medida de efectividad j.

Usando la integral anterior, se obtienen las fórmulas que permiten el cálculo de EBC dependiendo de la función densidad de probabilidad, el tipo de comportamiento y si la función utilidad es monotónica creciente o decreciente.

Keeney y Raiffa (1976) establecen que cuando las preferencias son crecientes y se tiene aversión constante al riesgo y para propensión constante al riesgo ; cuando las preferencias son decrecientes y se tiene aversión constante al riesgo y para propensión constante al riesgo

Por ejemplo, para determinar la fórmula para calcular EBC si las preferencias son crecientes, con propensión constante al riesgo y una función densidad de probabilidad de Laplace,

\[u(EBC) = \int_{x_0=-\infty}^{\infty} f_0(x_0)u(x_0)dx_0 = \int_{x_0=-\infty}^{\infty} \frac{a}{2} e^{-a|x_0-b|} e^{c x_0} dx_0\]

Como se tiene el valor absoluto \left|x{0}-b\right|, entonces \left|x{0}-b\right|=\left{x{0}-b\quad si x{0}\geq b\quad b-x{0}\quad si x{0}<b\right}

\[Luego u(EB C) = \int_{x_0=-\infty}^b e^{-a(b-x_0)}e^{(a+c)x_0}dx_0 + \int_{x_0=b}^\infty \frac{a}{2}e^{-a(x_0-b)}e^{cx_0}dx_0 = \\= \frac{ae^{-ab}}{2}\int_{x_0=-\infty}^b e^{(a+c)x_0}dx_0 + \frac{ae^{ab}}{2}\int_{x_0=b}^\infty e^{(-a+c)x_0}dx_0 = \\= \frac{ae^{-ab}}{2(a+c)}e^{(a+c)b} + \frac{ae^{ab}}{2(a-c)}e^{(c-a)b} = \frac{a^2e^{bc}}{a^2-c^2}\]

Como entonces empleando el logaritmo natural en ambos lados de la ecuación y despejando EBC, queda

La fórmula anterior se colocó en su lugar en la Tabla 2.

De la misma manera se desarrollaron las demás fórmulas que se muestran en las Tablas 1 y 2, donde el de la función uniforme está en el artículo de Acosta (2019) y el de la función normal en el de 2021.

Tabla 1: EBC u(x) creciente y aversión constante al riesgo o u(x) decreciente y propensión

Función de probabilidadparámetrosrangoEBC
Uniforme, $\frac{1}{b-a}$ a, b, media = $\frac{a+b}{2}$ $a \leq x_0 \leq b$ $\frac{-1}{c}Ln\left|\frac{e^{-cb}-e^{-ca}}{c(a-b)}\right|$
Normal $\frac{1}{\sqrt{2\pi}\sigma}e^{\frac{-(x_0-m)^2}{2\sigma^2}}$ Media = m,varianza = $\sigma^2$ $-\infty < x_0 < \infty$ $m - \frac{c\sigma^2}{2}$
Exponencial, $\lambda e^{-\lambda x_0}$ Media = $\frac{1}{\lambda}$ $0 \leq x_0$ $\frac{-1}{c}Ln\frac{\lambda}{c+\lambda}$
Cauchy, $\frac{1}{\pi} \frac{a}{a^{2} + (x_{0} - b)^{2}}$ Media = b, $a > 0, -\infty < b <$ $-\infty < x_{0} < \infty$ b + a
Chi-cuadrada, $\frac{x_{0}^{\left(\frac{n}{2}-1\right)} e^{-\frac{x_{0}}{2}}}{2^{\frac{n}{2}\left(\frac{n}{2}-1\right)!}}$ Media = n $x_{0} > 0$ $\frac{-n}{2c}Ln\left(\frac{1}{1+2c}\right)$
Erlang, $\frac{a^{n}x_{0}^{n-1}e^{-ax_{0}}}{(n-1)!}$ Media = $\frac{n}{a}, a > 0$ $x_{0} > 0$ $\frac{-n}{c}Ln\left(\frac{a}{c+a}\right)$
Gamma, $\frac{x_{0}^{a}e^{\frac{-x_{0}}{b}}}{a!b^{a+1}}$ Media = $(a + 1)b, a > -1, b > 0$ $x_{0} > 0$ $\frac{-1}{c}Ln\left(\frac{-1}{(1+bc)^{a+1}}\right)$
Laplace, $\frac{a}{2}e^{-a|x_{0} - b|}$ Media = b, a > 0, $-\infty < b < \infty$ $-\infty < x_{0} < \infty$ $-\frac{1}{c}Ln\left(\frac{a^{2}e^{-bc}}{a^{2}-c^{2}}\right)$

Tabla 2: EBC u(x) creciente y propensión constante al riesgo o u(x) decreciente y aversión

Función de probabilidadparámetrosrangoEBC
Uniforme, $\frac{1}{b-a}$ a, b $a \leq x_0 \leq b$ $\frac{1}{c}Ln\left|\frac{e^{cb}-e^{ca}}{c(b-a)}\right|$
Normal $\frac{1}{\sqrt{2\pi}\sigma}e^{\frac{-(x_0-m)^2}{2\sigma^2}}$ Media = m, varianza = $\sigma^2$ $-\infty < x_0 < \infty$ $m + \frac{c\sigma^2}{2}$
Exponencial, $\lambda e^{-\lambda x_0}$ Media = $\frac{1}{\lambda}$ $0 \leq x_0$ $\frac{1}{c}Ln\frac{\lambda}{|c-\lambda|}$
Cauchy, $\frac{1}{\pi}\frac{a}{a^2+(x_0-b)^2}$ Media = b, $a > 0, -\infty < b < \infty$ $-\infty < x_0 < \infty$ b-a
Chi-cuadrada, $\frac{x_0^{\left(\frac{n}{2}-1\right)}e^{-\frac{x_0}{2}}}{2^{\frac{n}{2}}\left(\frac{n}{2}-1\right)!}$ Media = n $x_0 > 0$ $\frac{n}{2c}Ln\left(\frac{1}{2c-1}\right)$
Erlang, $\frac{a^n x_0^{n-1}e^{-ax_0}}{(n-1)!}$ Media = $\frac{n}{a}, a > 0$ $x_0 > 0$ $\frac{n}{c}Ln\left(\frac{a}{c-a}\right)$
Gamma, $\frac{x_0^a e^{-\frac{x_0}{b}}}{a!b^{a+1}}$ Media = (a + 1)b, $a > -1, b > 0$ $x_0 > 0$ $\frac{1}{c}Ln\left(\frac{1}{(bc-1)^{a+1}}\right)$
Laplace, $\frac{a}{2}e^{-a|x_0-b|}$ Media = b, a > 0, $-\infty < b < \infty$ $-\infty < x_0 < \infty$ $\frac{1}{c}Ln\left(\frac{a^2e^{bc}}{a^2-c^2}\right)$

III. MÉTODO PROACT MODIFICADO

El algoritmo propuesto en este artículo consiste en transformar el problema aleatorio en uno determinista y después aplicar el método ProAct modificándolo, mediante la utilización de una función valor de dos atributos en las permutas compensatorias en vez de preguntar directamente al decisor, como sugieren los autores del método. Si la función valor de dos atributos cumple con la condición de intercambios correspondientes es de tipo aditivo, es decir, .

La condición de intercambios correspondientes (Keeney y Raiffa, 1976) consiste en lo siguiente:

Considere cuatro puntos

\[A: (x_{1}, y_{1}), B: (x_{1}, y_{2}), C: (x_{2}, y_{1}) \text{y} D: (x_{2}, y_{2})\]

Suponga que se cumple lo siguiente:

  1. En un aumento b en Y compensa una disminución a en X;

  2. En un aumento c en Y compensa una disminución a en X;

  3. En un aumento b en Y compensa una disminución d en X;

Si en un aumento c en Y compensa una disminución d en X y esto se cumple independientemente de las cantidades , se dice que se satisface la condición de intercambios correspondiente.

Keeney y Raiffa (1976) mencionan incluso que la función valor puede utilizarse, por facilidad de análisis, aún cuando sus resultados sean aproximaciones a los reales

También, como los tres axiomas de funciones valor son los tres primeros de los siete de las funciones utilidad (De Neufville, 1990) es posible concluir que todas las funciones utilidad son funciones valor y que no todas las funciones valor son funciones utilidad.

\[\text{Luego}, v\big(x_{0},y_{0}\big) = k_{1}u_{x}\big(x_{0}\big) + k_{2}u_{y}\big(y_{0}\big)\]

A continuación, se presenta el algoritmo donde interactúan el decisor y el analista.

Algoritmo

Parte uno: Definición de objetivos, medidas de efectividad y función utilidad

Paso 1-1 El analista le pide al decisor que defina el horizonte de planeación, los objetivos, las medidas de efectividad y sus rangos, estableciendo si las medidas son monotónicas crecientes o decrecientes.

Paso 1-2 El analista para cada medida de efectividad forma la lotería L cuyas consecuencias son los extremos del rango proporcionado en el paso anterior, le asigna probabilidades de 0.5 a cada una y le pregunta al decisor el equivalente bajo certeza (EBC) de cada lotería.

Paso 1-3 El analista calcula el valor esperado (VE) de cada lotería y determina la función utilidad, u(x) de cada medida de efectividad. Como establecen Keeney y Raiffa (1976) si las preferencias son crecientes (por ejemplo, ganancias) entonces, si EBC = VE, neutralidad al riesgo, u(x) = x; si EBC < VE, aversión al riesgo, ; si EBC > VE, propensión al riesgo, ; si las preferencias son decrecientes (por ejemplo, costos) entonces, si EBC = VE, neutralidad al riesgo, u(x) = -x; si EBC > VE, aversión al riesgo ; si EBC < VE, propensión al riesgo,

Paso 1-4 El analista, para cada medida de efectividad donde el decisor tiene aversión o propensión al riesgo, determina el valor de c > 0, mediante la ecuación u(EBC) = u(L). Para ello puede utilizar el algoritmo Newton-Raphson, ; donde es la primera derivada de f con respecto a c

Parte dos: Estipulación de alternativas

Paso 2-2 El analista calcula el EBC de cada alternativa y medida de efectividad, utilizando las tablas 1 y 2 cuando el decisor tiene aversión o propensión. En el caso de neutralidad EBC = media de la función densidad de probabilidad. El analista presenta los EBC calculados en una Tabla cuyos renglones son las alternativas y las columnas las medidas de efectividad.

Parte tres: Selección de la mejor alternativa.

Paso 3-2 El analista ve las alternativas que no han sido eliminadas. Si sólo queda una, esa es la mejor opción y termina el algoritmo. Si no es así, habrá que continuar con el Paso 3-3.

Paso 3-3 El analista le hará una pregunta al decisor para determinar la función valor para dos medidas de efectividad. Sean X, Y dos medidas de efectividad, donde son los peores valores y son los mejores. El analista le pregunta la cantidad de x donde . Como hay indiferencia

Si se cumple la condición de intercambios correspondientes la función valor es de tipo aditivo, o sea, . Como todas las funciones utilidad son funciones valor, se emplean las funciones determinadas en el Paso 1-3.

\[\begin{array}{c} \text {Quedando v(x,y) = k_{x} u_{x} (x) + k_{y} u_{y} (y)y k_{x} + k_{y} = 1; usando v(x,y_{*}) = v(x_{*},y^{*}), queda} \\k _ {x} u _ {x} (x) + k _ {y} u _ {y} (y _ {*}) = k _ {x} u _ {x} (x _ {*}) + k _ {y} u _ {y} (y ^ {*}) y k _ {x} + k _ {y} = 1, \end{array}\]

que es un sistema de dos ecuaciones con dos incógnitas. Al resolver este sistema está determinada la función valor para dos medidas de efectividad.

Si sólo quedan dos medidas de efectividad se evalúan, usando la función valor, todas las alternativas y la que tenga el valor mayor es la mejor opción y termina el algoritmo. Si no es así, habrá que continuar con el Paso 3-4.

Paso 3-4 El analista con la función valor de dos medidas de efectividad, puede hacer que todas las alternativas tengan igual una de esas medidas. En cuanto logra ello se elimina la medida de efectividad que es igual para todas las alternativas, y se regresa al Paso 3-1.

Se presenta ahora un ejemplo donde se aplican los pasos del algoritmo.

IV. EJEMPLO. ELECCIÓN DEL SITIO DONDE SE VA A CONSTRUIR UN AEROPUERTO

Parte uno: Definición de objetivos, medidas de efectividad y función utilidad

Paso 1-1 El analista le pide al decisor que defina el horizonte de planeación, los objetivos, las medidas de efectividad y sus rangos, estableciendo si las preferencias de las medidas son monotónicas crecientes o decrecientes.

El decisor definió como horizonte de planeación 30 años, y cuatro objetivos:

Minimizar el costo de construcción y de operación del aeropuerto. Su medida de efectividad, , es dicho costo en miles de millones de pesos en valor presente; su rango [a, b] = [50, 200] mil millones de pesos; su preferencia es monotónica decreciente.

Minimizar el tiempo de acceso al aeropuerto. Su medida de efectividad, , es el tiempo promedio para llegar al aeropuerto desde diferentes zonas en minutos; su rango [a, b] = [15, 90] minutos; su preferencia es monotónica decreciente.

Maximizar la capacidad del aeropuerto. Su medida, , es el número de despegues y aterrizajes por hora; su rango [a, b] = [50, 250] operaciones por hora; su preferencia es monotónica creciente.

Minimizar el efecto de la contaminación del ruido causado por los aviones. Su medida, , es el número de miles de personas sujetas a un nivel de ruido de 80 decibeles o más; su rango [a, b] = [10, 200] miles de personas; su preferencia es monotónica decreciente.

Paso 1-2 El analista para cada medida de efectividad forma la lotería L cuyas consecuencias son los extremos del rango proporcionado en el paso anterior, le asigna probabilidades de 0.5 a cada extremo y le pregunta al decisor el equivalente bajo certeza (EBC) de cada lotería.

En la Tabla 3 están las respuestas del decisor.

Tabla 3: EBC de la lotería L = [a, 0.5; b, 0.5]

Medida de efectividadabEBC
$x_{1}$ : costo50200140
$x_{2}$ : tiempo de acceso159060
$x_{3}$ : capacidad50250120
$x_{4}$ : personas afectadas por el ruido10200120

Paso 1-3 El analista calcula el valor esperado (VE) de cada lotería. Con los valores de EBC y VE, así como el conocimiento si la medida de efectividad es monotónica creciente o decreciente, utilizando las Tablas 1 y 2, determina la función utilidad, u(x). El resultado está en la Tabla 4.

Tabla 4: u(x) para cada medida de efectividad

Medida de efectividadLoteríaVEFunción monotónicaVE vs EBCu(x)
$x_1$ L(50, 0.5; 200, 0.5)125decrecienteVE<EBC $-e^{cx}$
$x_2$ L(15, 0.5; 90, 0.5)52.5decrecienteVE<EBC $-e^{cx}$
$x_3$ L(50, 0.5; 250, 0.5)150crecienteVE>EBC $-e^{-cx}$
$x_4$ L(10, 0.5; 200, 0.5)105decrecienteVE<EBC $-e^{cx}$

Paso 1-4 El analista, para cada medida de efectividad donde el decisor tiene aversión o propensión al riesgo, formula la ecuación u(EBC) = u(L), que es una ecuación con una sola incógnita, el parámetro c, y resuelve dicha ecuación. (Puede utilizar el algoritmo Newton-Raphson)

Para la primera medida de efectividad,

De manera que

Su derivada

\[c _ {n + 1} = c _ {n} - \frac {f (c)}{f ^ {\prime} (c)}\]

Proporciona c = 0.00548, luego

Para la segunda medida de efectividad,

De manera que

Su derivada

Utilizando el algoritmo Newton-Raphson

\[c _ {n + 1} = c _ {n} - \frac {f (c)}{f ^ {\prime} (c)}\]

Proporciona c = 0.01096, luego

Para la tercera medida de efectividad,

De manera que

Su derivada

Utilizando el algoritmo Newton-Raphson

\[c _ {n + 1} = c _ {n} - \frac {f (c)}{f ^ {\prime} (c)}\]

Proporciona c = 0.00639, luego

Para la cuarta medida de efectividad,

De manera que

Su derivada

Utilizando el algoritmo Newton-Raphson

\[c _ {n + 1} = c _ {n} - \frac {f (c)}{f ^ {\prime} (c)}\]

Proporciona , luego

Parte dos: Estipulación de alternativas

Paso 2-1 El analista le pide al decisor que estipule las alternativas, precisando para cada una de ellas la función densidad de probabilidad de cada medida de efectividad, su rango, su media y su desviación estándar. Información que se presenta en las Tablas 5 a 8.

Tabla 5: Medida de efectividad 1. Costo (miles de millones de pesos)

AlternativaFunción densidadRangoMedia $m = \frac{a+b}{2}$ Desviación estándar $\sigma = \frac{b-a}{\sqrt{12}}$
ab
Sitio 1Uniforme5070605.773503
Sitio 2Uniforme601008011.54701
Sitio 3Uniforme8012010011.54701
Sitio 4Uniforme10020015028.86751
Sitio 5Uniforme6090758.660254
Sitio 6Uniforme9015012017.32051
Sitio 7Uniforme13019016017.32051
Sitio 8Uniforme1201401305.773503
Sitio 9Uniforme13017015011.54701
Sitio 10Uniforme1501801658.660254

Tabla 6: Medida de efectividad 2. Tiempo de acceso (minutos)

AlternativaFunción densidadRangoMedia $m = \frac{a+b}{2}$ Desviación estándar $\sigma = \frac{b-m}{3.6}$
ab
Sitio 1Normal1535252.777778
Sitio 2Normal3555452.777778
Sitio 3Normal5575652.777778
Sitio 4Normal759082.52.083333
Sitio 5Normal2040302.777778
Sitio 6Normal4060502.777778
Sitio 7Normal6080702.777778
Sitio 8Normal1575458.333333
Sitio 9Normal359062.57.638889
Sitio 10Normal2030251.388889

Tabla 7: Medida de efectividad 3. Capacidad (operaciones por hora)

AlternativaFunción densidadRangoMedia $m = \frac{a+b}{2}$ Desviación estándar $\sigma = \frac{b-m}{3.6}$
ab
Sitio 1Normal50100756.944444
Sitio 2Normal1001501256.944444
Sitio 3Normal1502001756.944444
Sitio 4Normal2002502256.944444
Sitio 5Normal70120956.944444
Sitio 6Normal1401601502.777778
Sitio 7Normal1802202005.555556
Sitio 8Normal2402502451.388889
Sitio 9Normal2202302251.388889
Sitio 10Normal2002102051.388889

Tabla 8: Medida de efectividad 4. Personas afectadas por el ruido (miles de personas) Paso 2-2 El analista calcula el EBC de cada alternativa y medida de efectividad, utilizando las tablas 3 y 4. Y presenta los EBC calculados en una Tabla cuyos renglones son las alternativas y las columnas las medidas de efectividad. El resultado está en la Tabla 9.

AlternativaFunción densidadRangoMedia $m = \frac{a+b}{2}$ Desviación estándar $\sigma = \frac{b-m}{3.6}$
ab
Sitio 1Normal1050305.555556
Sitio 2Normal50100756.944444
Sitio 3Normal1001501256.944444
Sitio 4Normal1502001756.944444
Sitio 5Normal4070554.166667
Sitio 6Normal60110856.944444
Sitio 7Normal1301701505.555556
Sitio 8Normal3070505.555556
Sitio 9Normal70120956.944444
Sitio 10Normal1401601502.777778
AlternativasCostoTiempo de accesoCapacidadAfectación por ruido
Sitio 160.0925.0274.8530.05
Sitio 280.3745.03124.8575.08
Sitio 3100.3765.03174.85125.08
Sitio 4152.2882.51224.85175.08
Sitio 575.2130.0394.8555.03
Sitio 6120.8250.03149.9885.08
Sitio 7160.8270.02199.90150.05
Sitio 8130.0945.22244.9950.05
Sitio 9150.3762.69224.9995.08
Sitio 10165.2125.01204.99150.01
AlternativasCostoTiempo de accesoCapacidadAfectación por ruido
Sitio 160.0925.0274.8530.05
Sitio 280.3745.03124.8575.08
Sitio 3100.3765.03174.85125.08
Sitio 4152.2882.51224.85175.08
Sitio 575.2130.0394.8555.03
Sitio 6120.8250.03149.9885.08
Sitio 7160.8270.02199.90150.05
Sitio 8130.0945.22244.9950.05
Sitio 9150.3762.69224.9995.08
Sitio 10165.2125.01204.99150.01

Paso 3-1 El analista encuentra y elimina las alternativas dominadas. La alternativa A domina a la B cuando A es mejor que B en algunas medidas de efectividad y en todas los demás es igual que B.

El sitio 8 domina a los sitios 4 y 7, por lo cual quedan eliminados. Quedando así la Tabla 10.

Tabla 10: Equivalentes bajo certeza

AlternativasCostoTiempo de accesoCapacidadAfectación por ruido
Sitio 160.0925.0274.8530.05
Sitio 280.3745.03124.8575.08
Sitio 3100.3765.03174.85125.08
Sitio 575.2130.0394.8555.03
Sitio 6120.8250.03149.9885.08
Sitio 8130.0945.22244.9950.05
Sitio 9150.3762.69224.9995.08
Sitio 10165.2125.01204.99150.01
AlternativasCostoTiempo de accesoCapacidadAfectación por ruido
Sitio 160.0925.0274.8530.05
Sitio 280.3745.03124.8575.08
Sitio 3100.3765.03174.85125.08
Sitio 575.2130.0394.8555.03
Sitio 6120.8250.03149.9885.08
Sitio 8130.0945.22244.9950.05
Sitio 9150.3762.69224.9995.08
Sitio 10165.2125.01204.99150.01

Paso 3-2 El analista ve las alternativas que no han sido eliminadas. Si sólo queda una, esa es la mejor opción y termina el algoritmo. Si no es así, habrá que continuar con el Paso 3-3.

Como aún hay 8 sitios no dominados, se continúa con el Paso 3-3

Paso 3-3 El analista le hará una pregunta al decisor para determinar la función valor para dos medidas de efectividad.

El analista escoge dos medidas de efectividad. En este ejemplo seleccionó : costo con rango [50, 200] y : tiempo de acceso con rango [15,90].

El analista le pregunta al decisor el valor de que le hace indiferente (200, 15) con . El decisor contestó .

Por lo anterior,

Como entonces donde que fueron determinadas en el Paso 1-4.

Esta ecuación junto con constituye un sistema de dos ecuaciones con dos incógnitas, que al resolverlo se tiene ,

Como el número de medidas de efectividad que permanecen es diferente de 2, habrá que continuar con el paso siguiente.

Paso 3-4 El analista con la función valor de dos medidas de efectividad, puede hacer que todas las alternativas tengan igual una de esas medidas. En cuanto se logra ello se elimina la medida de efectividad que es igual para todas las alternativas, y el analista regresa al Paso 3-1.

El analista fijó el tiempo de acceso minutos y calculó el costo que deben tener las alternativas. Para hacer este cálculo, empleó

Que proporciona

Los costos nuevos, calculados con la fórmula anterior, se muestran en la Tabla 11.

Tabla 11: Nuevos equivalentes bajo certeza

AlternativasCostoTiempo de accesoCapacidadAfectación por ruido
Sitio 160.0925.0274.8530.05
Sitio 2109.7025.02124.8575.08
Sitio 3155.2425.02174.85125.08
Sitio 582.5825.0294.8555.03
Sitio 6150.9825.02149.9885.08
Sitio 8153.0825.02244.9950.05
Sitio 9190.7525.02224.9995.08
Sitio 10165.1925.02204.99150.01
AlternativasCostoTiempo de accesoCapacidadAfectación por ruido
Sitio 160.0925.0274.8530.05
Sitio 2109.7025.02124.8575.08
Sitio 3155.2425.02174.85125.08
Sitio 582.5825.0294.8555.03
Sitio 6150.9825.02149.9885.08
Sitio 8153.0825.02244.9950.05
Sitio 9190.7525.02224.9995.08
Sitio 10165.1925.02204.99150.01

Como todas las alternativas tienen el mismo tiempo de acceso se elimina esa medida de efectividad quedando la Tabla 12.

Tabla 12: Nuevos equivalentes bajo certeza Se regresa al paso 3-1

AlternativasCostoCapacidadAfectación por ruido
Sitio 160.0974.8530.05
Sitio 2109.70124.8575.08
Sitio 3155.24174.85125.08
Sitio 582.5894.8555.03
Sitio 6150.98149.9885.08
Sitio 8153.08244.9950.05
Sitio 9190.75224.9995.08
Sitio 10165.19204.99150.01
AlternativasCostoCapacidadAfectación por ruido
Sitio 160.0974.8530.05
Sitio 2109.70124.8575.08
Sitio 3155.24174.85125.08
Sitio 582.5894.8555.03
Sitio 6150.98149.9885.08
Sitio 8153.08244.9950.05
Sitio 9190.75224.9995.08
Sitio 10165.19204.99150.01

Paso 3-1 El analista encuentra y elimina las alternativas dominadas.

El sitio 8 domina a los sitios 9 y 10, por lo que se eliminan los sitios dominados, quedando la Tabla 13.

Tabla 13: Nuevos equivalentes bajo certeza

AlternativasCostoCapacidadAfectación por ruido
Sitio 160.0974.8530.05
Sitio 2109.70124.8575.08
Sitio 3155.24174.85125.08
Sitio 582.5894.8555.03
Sitio 6150.98149.9885.08
Sitio 8153.08244.9950.05
AlternativasCostoCapacidadAfectación por ruido
Sitio 160.0974.8530.05
Sitio 2109.70124.8575.08
Sitio 3155.24174.85125.08
Sitio 582.5894.8555.03
Sitio 6150.98149.9885.08
Sitio 8153.08244.9950.05

Paso 3-2 El analista ve las alternativas que no han sido eliminadas. Si sólo queda una, esa es la mejor opción y termina el algoritmo. Si no es así, habrá que continuar con el Paso 3-3.

Como aún hay 6 sitios no dominados, se continúa con el Paso 3-3

Paso 3-3 El analista le hará una pregunta al decisor para determinar la función valor para dos medidas de efectividad.

El analista escoge dos medidas de efectividad. En este ejemplo seleccionó : costo con rango [50, 200] y : capacidad con rango [50, 250].

El analista le pregunta al decisor el valor de que le hace indiferente (200, 250) con . El decisor contestó .

Por lo anterior,

Como entonces donde que fueron determinadas en el Paso 1-4.

Esta ecuación junto con constituye un sistema de dos ecuaciones con dos incógnitas, que al resolverlo se tiene ,

\[v\left(x_{1},x_{3}\right) = 0.246411\left(-e^{0.00548x_{1}}\right) + 0.753589\left(-e^{-0.00639x_{3}}\right)\]

Como el número de medidas de efectividad que permanecen es diferente de 2, habrá que continuar con el paso siguiente.

Paso 3-4 El analista con la función valor de dos medidas de efectividad, puede hacer que todas las alternativas tengan igual una de esas medidas. En cuanto se logra ello se elimina la medida de efectividad que es igual para todas las alternativas, y el analista regresa al Paso 3-1.

El analista fijó la capacidad operaciones por hora y calculó el costo que deben tener las alternativas. Para hacer este cálculo, empleó

Que proporciona

Los costos nuevos, calculados con la fórmula anterior, se muestran en la Tabla 14.

Tabla 14: Nuevos equivalentes bajo certeza Como todas las alternativas tienen la misma capacidad se elimina esa medida de efectividad quedando la Tabla 15.

AlternativasCostoCapacidadAfectación por ruido
Sitio 1177.60244.9930.05
Sitio 2171.70244.9975.08
Sitio 3181.44244.99125.08
Sitio 5174.46244.9955.03
Sitio 6189.25244.9985.08
Sitio 8153.08244.9950.05

Tabla 15: Nuevos equivalentes bajo certeza

AlternativasCostoAfectación por ruido
Sitio 1177.6030.05
Sitio 2171.7075.08
Sitio 3181.44125.08
Sitio 5174.4655.03
Sitio 6189.2585.08
Sitio 8153.0850.05

Paso 3-1 El analista encuentra y elimina las alternativas dominadas.

El sitio 8 domina a los sitios 2, 3, 5 y 6, por lo que se eliminan los sitios dominados, quedando la Tabla 16.

Tabla 16: Nuevos equivalentes bajo certeza

AlternativasCostoAfectación por ruido
Sitio 1177.6030.05
Sitio 8153.0850.05

Paso 3-2 El analista ve las alternativas que no han sido eliminadas. Si sólo queda una, esa es la mejor opción y termina el algoritmo. Si no es así, habrá que continuar con el Paso 3-3.

Como aún hay 2 sitios no dominados, se continúa con el Paso 3-3

Paso 3-3 El analista le hará una pregunta al decisor para determinar la función valor para dos medidas de efectividad.

Sólo quedan dos medidas de efectividad : costo con rango [50, 200] y : afectación por ruido con rango [10,200].

El analista le pregunta al decisor el valor de que le hace indiferente (200, 10) con . El decisor contestó .

Por lo anterior,

Como entonces donde que fueron determinadas en el Paso 1-4.

Esta ecuación junto con constituye un sistema de dos ecuaciones con dos incógnitas, que al resolverlo se tiene ,

\[v\left(x_{1},x_{4}\right) = 0.75\left(- e^{0.00548x_{1}}\right) + 0.25\left(- e^{0.00338x_{4}}\right)\]

Como sólo quedan dos medidas de efectividad, usando la ecuación anterior, se calcula el valor de los sitios 1 y 8 y se elige el que tenga mayor valor. Así, , por lo que el sitio 8 es el mejor.

V. CONCLUSIONES

En este artículo se presenta un algoritmo donde interactúan un decisor y un analista. El analista le hace preguntas al decisor y con sus respuestas efectúa los cálculos requeridos.

El algoritmo siempre converge a una solución óptima.

Los análisis deterministas son más sencillos que los aleatorios, con lo que se logran simplificaciones importantes que permiten llegar a la solución con un esfuerzo mínimo tanto del decisor como del analista.

Este algoritmo es útil cuando existe un solo decisor, una sola etapa y las distribuciones de probabilidad son: uniforme, normal, exponencial, Cauchy, Chi-cuadrada, Erlang, Gamma y Laplace. Si en el problema que se estuviese analizando las distribuciones fuesen diferentes, este artículo presenta el procedimiento para que el analista pueda determinar los valores requeridos de los equivalentes bajo certeza.

Conflict of Interest

The authors declare no conflict of interest.

Ethical Approval

Not applicable

Data Availability

The datasets used in this study are openly available at [repository link] and the source code is available on GitHub at [GitHub link].

Funding

This work did not receive any external funding.

References

7 Cites in Article

Cite this article

Generating citation...

Related Research

  • Code: 003
  • Version of record

    v1.0

  • Issue date

    05 August 2024

  • Language

    en

Article Placeholder
Open Access
Research Article
CC-BY-NC 4.0
Views 658
Downloads 29
Special Issue

Launch a focused special issue to highlight research, emerging trends, and expert insights in your academic field.

Support