Тема: ПОСТРОЕНИЕ ОПТИМАЛЬНЫХ И ЛОКАЛЬНО ОПТИМАЛЬНЫХ ПЛАНОВ ДЛЯ РЕГРЕССИОННЫХ МОДЕЛЕЙ

Теория планирования — раздел математической статистики, занимающийся оптимальным планированием условий эксперимента. Основание этой теории было заложено Р. Фишером в 1935 году в его работе [1]. Результат эксперимента зависит от некоторых параметров и часто стоит задача нахождения оценки этих параметров или проверки некоторой гипотезы относительно них. Когда проведение опыта связано со значительными временными или материальными затратами, требуется осуществлять рациональный выбор плана эксперимента, чтобы за наименьшее число опытов получить наиболее точную оценку.
Существуют множество критериев оптимальности, наиболее исследованным из которых является D-критерий, минимизирующий объем доверительного эллипсоида. Также популярным является критерий L-оптимальности, который позволяет минимизировать среднюю дисперсию оценок параметров (определения соответствующих критериев будут даны в пункте1.2).
Целью данной работы является построение L-оптимального плана для квадратической и кубической модели, а также D-оптимального плана для обощенной модели Михаэлиса-Ментен и дробно-рациональной модели с четырьмя параметрами (в случае большого и малого промежутков). Работа осуществлена в рамках темы СПбГУ 6.38.435.2015.
Сначала будут рассмотрены L-оптимальные планы на отрезке [-1,1] для некоторых полиномиальных моделей (D-оптимальные планы для них уже были найдены, см. [2, с. 64]).
Далее в работе будут рассмотрены локально D-оптимальные планы для некоторых нелинейных по параметрам моделей. Например, для обобщенной модели Михаэлиса- Ментен, которая находит широкое применение во многих областях науки: например, сельское хозяйство [3], биология сохранения живой природы [4], безопасность окружающей среды [5] — лишь некоторые из них. В биохимии с помощью этой модели описываются уравнения ферментативной кинетики — зависимости скорости протекания химической реакции от концентрации субстрата.
Наконец, мы рассмотрим дробно-рациональную модель с четырьмя параметрами, для которой будет более подробно описан метод нахождения локально D-оптимального плана, представленный в работе [6], а также найдены некоторые новые результаты, позволяющие численно находить локально D-оптимальный план в том случае, когда выразить его явно не представляется возможным. Данная модель также представляет практический интерес, поскольку применяется в химической кинетике [7], [8], а также в биологии и сельском хозяйстве

Купить

Таким образом, мы построили L-оптимальный план для полиномиальной модели,
т.е. для /(ж) = (1, х,..., хт1)Tдля т = 2 :
и для т = 3 :
-1 -0.439545
<з* = I
0.145348 0.354652
Также мы рассмотрели нелинейную обобщенную модель Михаэлиса-Ментен с тремя параметрами ц(х,4) = 41+
х +43
и построили локально D-оптимальный план на отрезке X = [0,4] :
Далее, перейдя к рассмотрению дробно-рациональной модели с четырьмя параметрами:
41 4з
4(х,4)= —-г- + —-г-,
X + 42 X + 44
мы подробно описали способ получения явной формулы для опорных точек в случае достаточно большого промежутка (когда d — достаточно большое число):
х2,4 = -1(1 + у Д(1 + (1 )2- 4). х* = 1.
Наиболее интересным представляется новый результат относительно поиска локально D-оптимального плана для дробно-рациональной модели в случае малого промежутка (когда 4 меньше некоторого критического значения 4*). Применив доказанную в работе [6, с. 41] теорему, мы получили метод нахождения численного решения для задачи. Был написан алгоритм на языке R,реализующий данный метод.
Частью алгоритма является минимизация специальной функции, описанной в утверждении теоремы2. Для поиска минимума этой функции были получены необходимые для минимизации ограничения (2.39) и схема построения начального приближения.

Купить