В этом примере показано, как идентифицировать важные параметры модели для модели [1] роста опухоли. В этом примере вы вычисляете локальную зависящую от времени чувствительность роста опухоли относительно параметров модели для лекарства от рака.
Модель, используемая в этом примере, является реализацией SimBiology фармакокинетической/фармакодинамической модели (PK/PD) Simeoni и др. Это определяет количество эффекта лекарств от рака на кинетике роста опухоли от в естественных условиях исследований на животных. Фармакокинетика препарата описана моделью 2D отсека с болюсным введением IV и линейным устранением (ke) из отсека Central. Рост опухоли описан как двухфазный процесс с начальным экспоненциальным ростом, сопровождаемым линейным ростом. Темп роста распространяющихся опухолевых клеток описан
.
Здесь, L0, L1 и Ψ являются параметрами роста опухоли, x1 является весом распространяющихся опухолевых клеток, и w является общим весом опухоли. В отсутствие любых наркотиков опухоль состоится только распространяющихся ячеек, то есть, w = x1. В присутствии противоракового агента это принято, что часть распространяющихся ячеек преобразовывается в нераспространяющиеся ячейки. Уровень этого преобразования принят, чтобы быть функцией концентрации препарата в плазме и факторе эффективности k2. Нераспространяющиеся ячейки, x2, проходят серию транзитных этапов (x3 и x4) и в конечном счете очищены от системы. Поток через транзитных отсеков моделируется как процесс первого порядка с константой скорости k1.
Некоторые корректировки были внесены в фармакодинамику роста опухоли.
Вместо того, чтобы задать вес опухоли как сумму x1, x2, x3 и x4, вес опухоли определяется реакцией под названием Increase, null → tumor_weight, со скоростью реакции .
tumor_weight является общим весом опухоли, x1 является весом распространяющихся опухолевых клеток и L0, и L1 является параметрами роста опухоли.
Точно так же уменьшение в весе опухоли задано реакцией под названием Decay, tumor_weight → null, со скоростью реакции k1*x4. k1 является параметром форвардного курса, и x4 является последними разновидностями в ряду транзитных сокращений веса опухоли.
ke задан как функция разрешения и объем центрального отсека: ke = Cl_Central/Central.
Анализ чувствительности позволяет вам определить, какие константы скорости и концентрации в модели имеют значительное влияние на полное поведение модели. В этом примере показано, как заняться расследованиями, какие параметры модели чувствительны к росту опухоли с помощью встроенной программы анализа чувствительности. А именно, этот пример показывает вам, как вычислить чувствительность разновидностей tumor_weight относительно параметров модели можно следующим образом:
Предположим, что вы уже калибровали параметры модели роста опухоли к экспериментальным данным PK/PD и сохраненным оценкам параметра и начальным условиям как вариант под названием parameterEstimates в вашей модели. Как последующий анализ, вы хотите узнать, какие параметры модели чувствительны к росту опухоли. Следующие шаги показывают, как сделать настолько использующий встроенную программу анализа чувствительности.
Откройте приложение SimBiology Model Analyzer путем ввода simBiologyModelAnalyzer в командной строке MATLAB® или путем нажатия на SimBiology Model Analyzer на вкладке Apps.
На вкладке Home нажмите Open. Перейдите к папке matlabroot\examples\simbio, где matlabroot является папкой, где MATLAB установлен. Откройте файл с именем проекта SimBiology tumor_growth_vpop_sa.sbproj.
Если вы используете Mac, нажмите Command +Shift+G в диалоговом окне File Browser и введите полный путь в папку.
На вкладке Home выберите Program> Calculate Sensitivities. Новая программа открывается.
В разделе Variants на шаге Model программы кликните по опции, чтобы просмотреть варианты, которые будут применены. Затем выберите parameterEstimates. Этот вариант содержит ранее оцененные значения параметров и начальные условия для препарата.
В разделе Doses выберите interval_dose. Эта повторная доза применяет 30 мг препарата каждые 4 дня, начиная в день 7 в течение в общей сложности пяти раз.
В разделе States To Log выберите [Tumor Growth Model].tumor_weight только.
Контекстное меню таблицы содержит опции, чтобы добавить разновидности, параметры или отсеки, и позволяет вам сделать выборы. Доступ к контекстному меню путем щелчка правой кнопкой по таблице.
На шаге Simulation (шаг выполнения), по умолчанию, метод нормализации установлен в Full (full dedimensionalization), означая, что приложение полностью нормирует чувствительность так, чтобы они могли сравниться друг с другом. Для получения дополнительной информации смотрите Normalization.
Каждый шаг выполнения программы может быть запущен отдельно. Шаг выполнения отмечен значком запуска наверху шага. Выполнение отдельного шага полезно особенно, если программа содержит несколько шагов, и вы хотите видеть, что промежуточное звено следует из конкретного шага. Путем выполнения так, можно внести изменения по мере необходимости прежде, чем запустить следующий шаг программы или целую программу. Чтобы запустить целую программу, нажмите кнопку Run на вкладке Home.
В разделе Sensitivities to Compute щелкните правой кнопкой где угодно по таблице и выберите Add All Constant Parameters. Очистите ke, потому что он задан как функция Central объем и Cl_Central, который уже является входом. Кликните по последней пустой строке и введите tumor. Затем выберите [Tumor Growth Model].tumor_weight из списка. Разновидность tumor_weight затем добавляется как выход.
На вкладке Home нажмите Run.
После анализа чувствительности программа автоматически генерирует два графика. Кликните по вкладке Plot1, чтобы просмотреть график курса времени веса опухоли и других значений чувствительности d[tumor_weight]/d[parameter]. Графики временной зависимости для d[tumor_weight]/d[k21] и d[tumor_weight]/d[k12] останьтесь относительно плоскими по сравнению с другой чувствительностью, указав, что рост опухоли не чувствителен к параметрам k12 и k21. Property Editor позволяет вам выбрать который графики отобразиться.
Графики поддерживаются данными, которые в настоящее время присутствуют в рабочей области приложения. Графики не являются снимками состояния. Когда данные (или экспериментальные данные или результаты симуляции) удалены или изменены, графики также обновляются согласно изменениям в базовых данных.
Можно также исследовать те же данные о чувствительности путем проверки величины вычисленной чувствительности, интегрируемой в зависимости от времени. Второй сгенерированный график (Plot2) показывает столбиковую диаграмму, содержащую такую информацию. График показывает, что вес опухоли чувствителен к L0, L1, w0, k2, k1 и Cl_Central, но не к k12 и k21.
Сохраните результаты анализа чувствительности в отдельной папке.
Щелкните правой кнопкой по папке LastRun программы в панели Browser.
Выберите Save Data.
В диалоговом окне Save Data введите sa_data как имя данных.
Результаты анализа чувствительности показывают, что вес опухоли не чувствителен к значениям модели k12 и k21. Исследуйте, приводят ли, варьируясь эти значения параметров к различным результатам чувствительности.
Варьируйтесь k12. Во-первых, варьируйтесь значения для параметра k12.
Кликните по вкладке Program1. Кликните (+) плюс наверху программы и выберите Generate Samples.
Шаг Generate Samples появляется. В разделе Parameter Set шага дважды кликните пустую ячейку в Component Name и введите k12. Установите следующие опции:
Ввод Область значений значений
Spacing — линейный
Min 1e-3
Max 1
# Of Steps — 25
Отключите генерацию графика по умолчанию путем нажатия кнопки графика наверху шага.
На шаге Simulation, под Sensitivities to Compute, очищают все входные параметры кроме k12. Сохраните tumor_weight как выход. Кроме того, отключите генерацию графика по умолчанию путем нажатия кнопки графика наверху шага.
На вкладке Home нажмите кнопку Run, чтобы запустить целую программу.
Если симуляция заканчивается, чувствительность, results хранится в папке LastRun. Расширьте папку LastRun. Нажмите results. Затем выберите time из раздела Plot по вкладке Home.
Responses в графике соответствует весу опухоли и чувствительности веса опухоли относительно k12. Каждый ответ построен с помощью различного стиля линии. Scenarios соответствует каждому сканированию параметра (сценарий симуляции). Каждый сценарий построен с помощью различного цвета.
Настройте график так, чтобы график использовал различный цвет в каждом ответе вместо этого. В таблице Slice Data очистите Style Scenarios путем выбора пустой. Установите Responses на Color.
В таблице Responses очистите tumor_weight и сохраните только результаты чувствительности.
График предполагает, что вес опухоли чувствителен к некоторым значениям параметров k12, когда варианты модели, отличающиеся, чем предполагаемые значения параметров, исследуются.
Постройте чувствительность k12 вместе с предыдущими результатами чувствительности, насколько чувствительный рост опухоли к k12 относительно других параметров модели. В панели Browser расширьте папку sa_data что вы сохраненный ранее. Затем расширьте results.
Мультивыберите (Ctrl + Click) элементы данных чувствительности для L0, L1, w0, k2, k1 и Cl_Central. Затем перетащите их на графике чувствительности k12 Plot3.
График предполагает, что рост опухоли еще более чувствителен к большинству других параметров, чем k12.
Сохраните результаты в отдельной папке. Щелкните правой кнопкой по LastRun и выберите Save Data. Введите k12_data как имя данных.
Варьируйтесь k21. Можно выполнить подобный анализ различным k21 и сравнив его результаты с другой чувствительностью параметра.
Возвратитесь к Program1. В нижней части шага Generate Samples нажмите Add parameter set to scan. Появляется новый раздел Parameter Set.
Дважды кликните пустую ячейку в Component Name и введите k21. Установите те же опции как k12:
Ввод Область значений значений
Spacing — линейный
Min 1e-3
Max 1
# Of Steps — 25
Отключите первый набор параметра (PS1). Укажите на верхний угол раздела Parameter Set PS1 и кликните по кнопке действий. Затем выберите Disable Parameter Set.
Parameter Set затем grayed, чтобы указать, что вы успешно отключили его.
На шаге Simulation выберите k21 как единственный вход и tumor_weight как выход.
Поразите Run от вкладки Home, чтобы запустить программу снова.
Если симуляция заканчивается, чувствительность, results хранится в папке LastRun. Расширьте папку LastRun. Нажмите results. Затем выберите time из раздела Plot по вкладке Home.
В таблице Slice Data очистите Style Scenarios путем выбора пустой. Установите Responses на Color.
В таблице Responses очистите tumor_weight и сохраните только результаты чувствительности.
График предполагает, что вес опухоли чувствителен к некоторым значениям параметров k21, когда варианты модели, отличающиеся, чем предполагаемые значения параметров, исследуются.
Постройте чувствительность k21 вместе с предыдущими результатами чувствительности выдержать сравнение. Нажмите Plot3. В панели Browser расширьте папку LastRun. Затем расширьте results.
Перетащите данные о чувствительности k21 d[tumor_weight]/d[k21] на Plot3.
График предполагает, что рост опухоли относительно более чувствителен к большинству других параметров, чем k21 и k12.
В этом примере показано, как выполнить анализ чувствительности, чтобы найти важные параметры модели, к которым рост опухоли чувствителен. Начальный анализ приходит к заключению, что вес опухоли не чувствителен к некоторым параметрам. Пример затем исследует пространство параметров этих менее чувствительных параметров и сравнивает результаты чувствительности с теми из других параметров модели. Можно использовать такую информацию о чувствительности в других исследованиях, например, выполнить симуляции Монте-Карло путем варьирования чувствительных параметров, чтобы исследовать динамику модели и биологическую изменчивость.
[1] Simeoni, M., П. Магни, К. Кэммия, Г. Де Николао, V. Шафраны, Э. Пезенти, М. Джермани, я. Poggesi и М. Роккетти. 2004. Прогнозирующее фармакокинетическо-фармакодинамическое моделирование кинетики роста опухоли в моделях ксенотрансплантата после администрирования противораковых агентов. Исследования рака. 64:1094-1101.