Этот пример показывает, как выполнить скан параметра путем симуляции модели несколько раз, каждый раз изменяя значение параметра.
В модели, описанной в модели дрожжевого гетеротримерного цикла G-белка, скорость инактивации G-белка (kGd
) намного ниже по мутантному штамму по сравнению со штаммом дикого типа (kGd = 0.004
от kGd = 0.11
), что объясняет более высокие уровни активированного G-белка (Ga
) в мутантном штамме. Для получения подробной информации о том, как изменяется уровень kGd
влияет на уровень Ga
, выполните скан параметров над различными значениями kGd
.
Загрузите gprotein.sbproj
проект, который включает переменную m1
, объект модели.
sbioloadproject gprotein
Создайте вектор с пятью равномерно расположенными значениями для kGd
в диапазоне от 0.001
на 0.15
.
kGdValues = linspace(1e-3,0.15,5)';
Создайте SimFunction
объект, где kGd
- параметр входа для сканирования и Ga
- наблюдаемый вид. Передайте в пустом массиве [] как последний входной параметр для обозначения отсутствия дозированных видов.
simfunc = createSimFunction(m1,{'kGd'},{'Ga'},[]);
Симулируйте модель несколько раз с различными значениями кГд. Установите время остановки равным 1000.
sd = simfunc(kGdValues,1000);
Постройте график результатов симуляции, чтобы увидеть, как изменяется уровень kGd
влияет на уровень Ga
.
sbioplot(sd);
createSimFunction
| SimFunction object