Документация

Поиск и просмотр папки с помощью браузера «Текущая папка» matlab\toolbox\curvefit\curvefit.

Щелкните файл правой кнопкой мыши OpioidHypnoticSynergy.txtи выберите «Импорт данных». Появится мастер импорта.

Кроме того, данные можно импортировать программным путем. Введите следующий код для считывания данных «доза-реакция» из файла в рабочую область MATLAB.

data = importdata( 'OpioidHypnoticSynergy.txt' );
Propofol      = data.data(:,1);
Remifentanil  = data.data(:,2);
Algometry     = data.data(:,3);
Tetany        = data.data(:,4);
Sedation      = data.data(:,5);
Laryingoscopy = data.data(:,6);

Чтобы создать поверхности отклика, необходимо выбрать два препарата для входов X и Y и один из четырех эффектов для выхода Z. После загрузки переменных в рабочее пространство можно либо открыть инструмент и выбрать переменные в интерактивном режиме, либо указать переменные начального вписывания с помощью cftool команда.

Введите следующее, чтобы открыть приложение «Фитинг кривой» (при необходимости) и создать новую поверхность ответа для Algometry:

cftool(Propofol, Remifentanil, Algometry)

Просмотр приложения «Фитинг кривой» X, Y, и Z входные и выходные органы управления. Инструмент отображает выбранные переменные Propofol, Remifentanil и Algometry, с посадкой поверхности. Аппроксимация по умолчанию - это интерполированная поверхность, которая проходит через точки данных.

Создайте копию текущей посадки поверхности одним из следующих способов:

Выберите тип подгонки «Пользовательское уравнение» в раскрывающемся списке, чтобы определить собственное уравнение для подгонки данных.

Выберите и удалите пример пользовательского текста уравнения в поле редактирования.

Можно использовать пользовательское поле редактирования уравнений для ввода кода MATLAB для определения модели. Уравнение, определяющее модель, должно зависеть от входных переменных x и y и список фиксированных параметров, оцениваемых параметров или и того, и другого.

Модель из документа:

где _CA и _CB - концентрации лекарственного средства, а IC50A, IC50B, альфа и n - коэффициенты, подлежащие оценке.

Это можно определить в коде MATLAB как

Effect = Emax*( CA/IC50A + CB/IC50B + alpha*( CA/IC50A )...
         .* ( CB/IC50B ) ).^n ./(( CA/IC50A + CB/IC50B + ...
         alpha*( CA/IC50A ) .* ( CB/IC50B ) ).^n  + 1);

Указание инструменту, какие переменные следует уместить и какие параметры оценивать, требует перезаписи имен переменных CA и CB кому x, и y. Необходимо включить x и y при вводе пользовательского уравнения в поле редактирования. Принять Emax = 1 поскольку выходные данные эффекта нормализованы.

Введите следующий текст в поле редактирования пользовательской формулы.

( x/IC50A + y/IC50B + alpha*( x/IC50A ) .* ( y/IC50B ) ).^n
  ./(( x/IC50A + y/IC50B + alpha*( x/IC50A ) .* 
  ( y/IC50B ) ).^n  + 1);

Приложение «Фитинг кривой» подгоняет поверхность к данным с помощью пользовательской модели уравнений.

Задайте некоторые опции подгонки, щелкнув Опции подгонки (Fit Options) в пользовательской формуле.

В диалоговом окне «Параметры посадки» выполните следующие действия.

Просмотрите область «Результаты». Просмотрите (и, при необходимости, скопируйте) любой из следующих результатов:

Отобразите график остатков для проверки распределения точек относительно поверхности, нажав кнопку панели инструментов или выбрав Вид (View) > График остатков (Residuals Plot).

Чтобы создать код для всех посадок и графиков в сессии приложения «Фитинг кривой», выберите «Файл» > «Создать код».

Приложение «Фитинг кривой» генерирует код из сессии и отображает файл в редакторе MATLAB. Файл включает все посадки и графики в текущей сессии.

Сохраните файл с именем по умолчанию, createFits.m.

Можно воссоздать посадки и графики, вызвав файл из командной строки (с исходными данными или новыми данными в качестве входных аргументов). В этом случае исходные данные по-прежнему отображаются в рабочей области.

Выделите первую строку файла (исключая слово function) и вычислите его, щелкнув правой кнопкой мыши и выбрав «Вычислить выбор» в окне команд, нажав клавишу F9 или скопировав и вставив в командную строку следующее:

[fitresult, gof] = createFits(Propofol,...
 Remifentanil, Algometry)

Функция создает окно фигуры для каждого вписывания в сессии. На пользовательском рисунке подгонки показаны графики поверхности и остатков, созданные в интерактивном режиме в приложении «Фитинг кривой».

Создайте новое вписывание в ответ тетаны вместо алгометрии, введя:

[fitresult, gof] = createFits(Propofol,...
 Remifentanil, Tetany)

Файл необходимо отредактировать, если требуется новая метка ответа на графиках. Созданный код можно использовать в качестве начальной точки для изменения посадок поверхности и графиков в соответствии с требованиями пользователя. Список методов, которые можно использовать, см. в разделе sfit.