Добавил:
rn Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Учебное пособие Немирко Манило.pdf
Скачиваний:
106
Добавлен:
14.06.2020
Размер:
897.06 Кб
Скачать

1

 

 

1

0,95

 

 

0,95

0,9

N

 

0,85

 

 

 

 

 

0,6

 

 

0,6

A1

S*A

A2

 

 

Рис. 3.4

 

 

По табл. 3.14 построим из исходной матрицы M c

матрицу M f . Она бу-

дет иметь вид табл. 3.20. Можно легко убедиться, что данная матрица не имеет седловой точки. Воспользуемся для решения геометрической интерпретацией игры. Для этого выполним необходимые построения. Они представлены на рис. 3.4, откуда видно, что в точке решения N пересекаются

только две прямые, соответствующие стратегиям S1 и S3* , поэтому для точки N игру можно представить в вице игры 2 ×2 с матрицей M f в виде

табл. 3.21.

Пользуясь формулой (3.6), получаем решение в виде смешанной стратегии хирурга

p1

=

f22

f21

=

10,85

= 0,4 .

p

f

f

0,95 0,6

 

 

 

2

 

11

12

 

 

 

Таким образом, F (A1 )F (A2 )= 0,4 , и хирургу надо почти в два раза чаще отказываться от операции, чем оперировать больного.

3.5. Минимизация риска хирургического вмешательства в онкологии

Злокачественная опухоль – это неуклонно прогрессирующее заболевание с безусловно плохим прогнозом. Будем рассматривать только те из них, при которых нет конкурирующих методов лечения, а рекомендуемое хирургическое вмешательство сопряжено с непосредственным хирургическим риском. Последний может выражаться в виде послеоперационной летальности q , которая в случае онкологических заболеваний зависит от локализации

86

опухоли и характера заболевания и нередко достигает 20…40 % [8]. В этом случае клиническая операбельность – величина, равная вероятности выживания больного в случае успешного выполнения радикального вмешательства при резектабельной опухоли, равна 1q . Она оценивает возможность больного перенести в данном лечебном учреждении показанную ему тяжелую радикальную операцию по поводу рака, выжить и быть выписанным.

Рассматривая тактику хирурга при неосложненных опухолях, будем решать вопрос о том, предлагать или не предлагать больному радикальную операцию при имеющемся риске, полагая, что она целесообразна по онкологическим соображениям и может быть выполнена технически. При этом в качестве цели радикальной операции при раке рассмотрим максимизацию продолжительности жизни онкологического больного.

В этом случае мы имеем два состояния больного («природы»): S1 – больной операбелен, S2 – больной неоперабелен; вероятность состояния S2 равна q , а вероятность состояния S1 (1q). В распоряжении хирурга две стратегии: A1 – предложить больному радикальную операцию и A2 – отказаться от вмешательства. Выигрыш хирурга aij обозначим следующим обра-

зом: D – математическое ожидание продолжительности жизни данного больного при отказе от радикальной операции, G – в случае успешного исхода радикальной операции и O – при летальном ее исходе. Очевидно, что G > D . Эти величины могут быть найдены приближенно на основании статистических данных для каждого учреждения либо получены экспертным пу-

тем для каждого больного. Матрицы Ma , Mr и

M f для данного случая

имеют вид табл. 3.22, 3.23 и 3.24 соответственно.

 

 

 

 

 

 

 

Таблица 3.22

 

Таблица 3.23

 

 

Таблица 3.24

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ai

 

S j

 

Ai

 

S j

 

 

Ai

 

S j

 

 

 

 

 

 

 

 

 

 

 

 

S1

 

S2

 

S1

 

S2

 

S1

 

S2

 

 

 

 

 

 

 

 

 

A1

G

 

0

 

A1

0

 

0

 

 

A1

G

 

D

A2

D

 

D

 

A2

G D

 

D

 

 

A2

2D G

 

D

87

 

 

Для этих данных по Ma при

известном q

можно вычислить

 

 

= G (1q),

 

= D . В случае

 

 

>

 

 

 

получаем G(1q)> D, q <1D G и

 

a1

a2

a1

a2

принимаем решение A1. Если

 

<

 

, что соответствует

q >1D G , выби-

a1

a2

раем A2 . При равенстве a1 = a2 или q =1DG можно применять обе тактики.

Таким образом, при величине операционного риска q <1DG будет оптимальным предложить больному оперативное вмешательство, при q >1DG оптимальным решением является отказ от радикальной операции. При q =1DG оптимально любое решение, и показания к операции или отказ от нее могут быть продиктованы чисто клиническими соображениями.

Можно легко убедиться, что при известном q , используя матрицу M f вместо Ma , мы придем к тому же алгоритму принятия решения. Относительно эквивалентности матриц Ma и Mr в этом же смысле мы уже упоминали.

Пример 3.6. Определить стратегию хирурга, если ожидаемая средняя продолжительность жизни больного в случае успешной радикальной операции равна 18 месяцам. Без операции он может прожить в среднем 12 месяцев. Вероятность летального исхода радикальной операции для данного больного равна 0,4.

По условию задачи, согласно принятым ранее обозначениям q = 0,4, G =18, D =12. Так как 1DG =11218 = 0,33, то q >1DG , и оптимальной стратегией хирурга будет отказ от операции (стратегия A2 ).

Пример 3.7. Решить задачу, сформулированную в предыдущем примере, при условии: q = 0,6, Gгода= 3, 1годD = .

Так как 1DG =113 = 0,66 и q <1DG , то оптимальной стратегией хирурга будет проведение операции (стратегия A1).

Рассмотрим теперь случай, когда величина операционного риска q врачу неизвестна. В матрице Ma имеется седловая точка α =β = a22 = D (в табл. 3.22 она обведена), и поэтому оптимальной стратегией хирурга, согласно критерию Вальда, является чистая стратегия A2 – отказ от операции. Это

очень осторожная стратегия, приводящая к тому, что никто из больных не проживет более D лет, хотя среди них есть и такие, которым хирург мог бы

88

продлить жизнь до G лет. Для этой части больных неоправданные потери в виде нереализованной продолжительности жизни равны G D .

При использовании матрицы M f , учитывающей и необходимость ми-

нимизации неоправданных потерь, оказывается, что у нее седловой точки нет и решение ищется в смешанных стратегиях. В соответствии с (3.6) оптималь-

ное соотношение частот применения стратегий A1 и A2 этом случае

 

 

F (A1)

= G D .

(3.11)

 

F (A )

 

G + D

 

2

 

 

Таким образом, при незнании величины операционного риска q

опти-

мальным является применение хирургом смешанной стратегии с соотношением частот, даваемым выражением (З.11). Это правило менее точно, чем выбор оптимальной чистой стратегии при известной величине q , но оно все же лучше, чем отсутствие какого бы то ни было решения в аналогичной ситуации.

Если в силу разных причин хирург не может уверенно назвать ни величину D , ни величину G , то можно утверждать, что в такой ситуации вследствие априорно известного условия G + D > G D оптимальным является применение смешанной стратегии, причем A2 должна использоваться чаще, чем A1.

Для практической реализации приведенных выше типовых алгоритмов решения существует специальная карта [8] (см. прил. 2), которая вкладывается в историю болезни.

3.6. Решение игр m × n

Для игр m ×n геометрическая интерпретация неприменима. Здесь применяются чисто расчетные методы. Можно показать, что решение любой конечной игры m ×n может быть сведено к задаче линейного программирования.

Рассмотрим игру m ×n . У игрока A имеется m стратегий: A1, A2, ..., Am ; у игрока B есть n стратегий: B1, B2, ..., Bn . Такая игра задается матрицей иг-

ры

m ×n

a . Нужно

найти

две

оптимальные

смешанные стратегии

 

 

 

ij

 

 

 

 

 

 

 

 

S*

= (p ,

p , ..., p

) и S*

= (q ,

q , ..., q ), где

p , p

, ..., p

и q , q , ..., q

A

1

2

m

B

1 2

n

1 2

m

1 2

n

вероятности применения соответствующих чистых стратегий

A1, A2, ..., Am и

 

 

 

m

n

 

 

 

 

 

 

 

B1, B2, ..., Bn

и pi

=1, q j =1.

 

 

 

 

 

 

 

 

i=1

j=1

 

 

 

 

 

 

 

 

 

 

 

 

 

89

 

 

 

 

Нахождение S*A . Положим, что цена игры γ положительна, γ ≥ 0 . Это

всегда можно сделать, добавив ко всем членам матрицы игры достаточно большое положительное число М. При этом решение игры не изменится, а найденную величину γ нужно будет в конце также увеличить на М.

Если мы применяем S*A , а противник – чистую стратегию B j , то наш средний выигрыш будет равен

a j = p1a1 j + p2a2 j +... + pmamj

( j =1, 2, ..., n) .

Так как мы применяем S*A ,

то наш средний выигрыш не может быть

меньше цены игры γ, т. е. a j ≥ γ,

j =1, 2, ...,

n , поэтому

p1a11 + p2a21 +... + pmam1 p1a12 + p2a22 +... + pmam2

≥ γ;

 

≥ γ;

(3.12)

.................................................

p1a1n + p2a2n +... + pmamn ≥ γ.

Внимание! Строки вышеприведенной системы пишутся по столбцам матрицы игры!

Разделим все вышеприведенные неравенства на положительную величину γ и введем обозначения

x1 = p1γ, x2 = p2 γ, …, xm = pm γ .

Тогда система (3.12) превращается в следующую:

a11x1 +a21x2 +...

+am1xm 1;

 

a12x1 +a22x2 +...

+am2xm 1;

(3.13)

.............................................

 

a1n x1 +a2n x2 +... +amn xm 1.

Так как p1 + p2 +... + pm =1, то

x1 + x2 +... + xm = 1γ .

Мы хотим сделать наш гарантированный выигрыш максимально возможным. При этом величина 1γ принимает минимальное значение.

Таким образом мы получаем следующую задачу линейного программирования: найти такие неотрицательные значения переменных x1, x2, ..., xm , которые удовлетворяли бы линейным ограничениям (3.13) и обращали бы в минимум линейную функцию

90

L = x1 + x2 +... + xm

Решив эту задачу линейного программирования, мы можем найти оптимальную стратегию S*A игрока A.

Нахождение SB* . Оптимальная стратегия SB* находится аналогично. Разница заключается в том, что игрок B стремится не максимизировать, а минимизировать выигрыш, а значит максимизировать величину 1γ. Следовательно, вместо условий (3.13) должны соблюдаться условия

 

 

a11y1 +a12 y2

+... +a1n yn 1;

 

 

 

a21y1 +a22 y2

+... +a2n yn 1;

(3.14)

 

 

...............................................

 

 

 

am1y1 +am2 y2 +... +amn yn 1,

 

где y j = q j γ,

j =1,

2, ..., n .

 

 

Требуется

так

выбрать неотрицательные значения

переменных

y1, y2, ..., yn , чтобы они удовлетворяли условиям (3.14) и обращали в максимум линейную функцию

L = y1 + y2 +... + yn =1γ

или, что то же самое, обращали в минимум линейную функцию L′= −L:

L′= −y1 y2 ... yn = −1γ.

Таким образом, любая конечная игра m ×n сводится к паре задач линейного программирования.

Пример 3.8. Рассмотрим игру «Три пальца», которая формулируется следующим образом. Два игрока одновременно и независимо показывают 1, 2 или 3 пальца. Выигрыш равен сумме показанных пальцев. Если это число четное, то выигрывает игрок А, а если нечетное, то выигрыш достается игроку В. Платежная матрица этой игры имеет вид, как в табл. 3.25.

91

 

 

 

Таблица 3.25

 

 

 

 

 

Ai

 

Bj

 

 

 

 

 

 

B1

B2

 

B3

 

 

 

 

 

 

 

A1

2

–3

 

4

A2

–3

4

 

–5

A3

4

–5

 

6

 

 

 

Таблица 3.26

 

 

 

 

 

Ai

 

Bj

 

 

 

 

 

 

B1

B2

 

B3

 

 

 

 

 

 

 

A1

7

2

 

9

A2

2

9

 

0

A3

9

0

 

11

Прибавив ко всем членам этой матрицы число М = 5, сделаем их неотрицательными. Тогда матрица примет вид, как в табл. 3.26.

Найдем стратегию S*A . Для этого по полученной матрице составим с и- стему линейных неравенств (3.13):

7x1 +2x2 +9x3 1;

2x1 +9x2 1;

9x1 +11x3 1,

а функцию L = x1 + x2 + x3 устремим к минимуму.

Внимание! Неравенства составляются по столбцам этой матрицы. Введя новые переменные y1, y2, y3, которые, как и x1, x2, x3 , должны

быть неотрицательными, и перейдем от условий неравенств к условиям равенств, т. е. к основной задаче ЛП. Тогда

y1 = 7x1 +2x2 +9x3 1;

 

y2 = 2x1 +9x2 1;

 

y3 = 9x1 +11x3 1.

(n = 6) – и 3 уравнения

Получаем 6 переменных: x1, x2, x3, y1, y2, y3

(m = 3), следовательно, n m = 3. Как видно из пр

иведенных выражений,

нужны 3 свободные переменные. Выберем в качестве свободных x1, x2, x3 . Тогда y1, y2, y3 – базисные переменные. Пробуем первое опорное решение: x1 = x2 = x3 = 0 . Это решение недопустимо, так как все переменные y1, y2, y3 в этом случае становятся отрицательными.

Для нахождения опорного решения нужно шаг за шагом обменивать базисные переменные со свободными (так, чтобы число отрицательных сво-

92

бодных членов уменьшалось или при этом убывали их абсолютные величины). Осуществим замену переменных:

x3 y3

и переразрешим систему уравнений относительно новых базисных переменных y1, y2, x3. Тогда получим систему уравнений

y1 = −114 x1 +2x2 +119 y3 112 ; y2 = 2x1 +9x2 1;

x3 = −119 x1 +111 y3 +111 .

Пробуем следующее опорное решение: x1 = x2 = y3 = 0 . Оно также не подходит, поскольку y1 = −112 , y2 = −1. Но это решение лучше предыдуще-

го, так как число отрицательных переменных уменьшилось. Снова осуществим замену переменных

x2 y2 ;

 

 

 

y

= −

80 x

+ 2 y

2

+

 

9

 

 

y

3

+

 

4

 

;

 

 

 

 

 

 

 

 

 

 

 

11

99

 

 

 

 

 

 

 

 

 

 

 

 

1

 

99

1

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

= −

2 x +

1 y

2

 

+

1

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

9

 

1

 

9

 

 

 

 

9

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

= −

9

 

x

+

1

 

y

3

 

+

 

1

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

11

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

1

 

11

 

 

 

 

 

11

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

Пробуем опорное решение

x

 

= y

2

 

= y

3

= 0 ,

при этом x

=

,

x

=

,

 

 

 

9

11

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

3

 

 

y =

 

4

. Все переменные положительны,

поэтому

это решение

является

99

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

опорным. Вычислим целевую функцию

L = x1 + x2 + x3 = − 994 x1 + 19 y2 +111 y3 + 9920 .

В нашей опорной точке L = 9920 . Так как коэффициент при x1 имеет от-

рицательный знак, то, сделав x1 больше нуля, можно еще значительнее уменьшить L. Однако при этом надо сделать равной нулю какую-то переменную из базисных ( y1, x2, x3) . Из последней системы уравнений видно, что

93

при увеличении x1 быстрее всего станет равным нулю переменная y1. Поэтому нужно снова осуществить замену переменных:

x1 y1;

x1 = − 8099 y1 + 1140 y2 + 8081 y3 + 201 ;

 

 

 

x =

11 y +

 

1

 

 

y

2

 

9

 

y +

 

1

;

 

 

 

 

 

 

20

40

10

 

 

 

 

 

 

2

 

40

1

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

x =

81 y +

9

 

 

y

2

59 y

3

+

1

;

 

 

 

 

 

 

40

 

20

 

 

 

 

 

 

3

 

80

1

 

 

 

 

80

 

 

 

 

 

 

 

 

 

 

L =

 

1

y +

 

1

 

 

y

2

+

 

1

 

y

3

+

1 .

 

 

 

 

 

 

 

20

10

 

20

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

5

 

 

1

 

Полагая y

= y

2

= y

3

= 0 ,

получаем

 

опорное

решение x

= x =

,

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

3

20

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2 = 101 , L = 15 . Это – оптимальное решение, так как в L все коэффициенты

положительны и никакое увеличение y1, y2, y3 сверх нуля не приведет к

уменьшению значения целевой функции.

Итак, оптимальное решение основной задачи ЛП найдено. Теперь найдем решение нашей игры 3×3:

 

γ′=

1

= 5; γ = γ′−5 = 0 ,

 

 

L

 

 

1

 

min

 

 

 

1

 

 

 

1 .

а вероятности p1 = x1γ′=

; p2 = x2γ′=

;

p3

= x3γ′=

 

4

 

 

 

 

 

2

 

 

 

4

 

 

*

 

1

,

1

,

1

 

 

Оптимальная стратегия SA =

 

4

2

4

, а цена игры γ = 0 .

 

 

 

 

 

 

 

 

 

Аналогично решается задача поиска оптимальной стратегии SB* игрока

*

 

1

,

1

,

1

 

B . Можно показать, что SB =

4

2

4

.

 

 

 

 

 

94