Documentation

This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English verison of the page.

Note: This page has been translated by MathWorks. Please click here
To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.

Perform a Parameter Scan

This example shows how to perform a parameter scan by simulating a model multiple times, each time varying the value of a parameter.

In the model described in Model of the Yeast Heterotrimeric G Protein Cycle, the rate of G protein inactivation (kGd) is much lower in the mutant strain versus the wild-type strain (kGd = 0.004 versus kGd = 0.11), which explains higher levels of activated G protein (Ga) in the mutant strain. For a detailed look at how varying the level of kGd affects the level of Ga, perform a parameter scan over different values of kGd.

Load the gprotein.sbproj project, which includes the variable m1, a model object.

sbioloadproject gprotein

Create a vector of five evenly spaced values for kGd ranging from 0.001 to 0.15.

kGdValues = linspace(1e-3,0.15,5)';

Create a SimFunction object, where kGd is the input parameter to scan, and Ga is the observed species. Pass in an empty array [] as the last input argument to denote there are no dosed species.

simfunc = createSimFunction(m1,{'kGd'},{'Ga'},[]);

Simulate the model multiple times with different kGd values. Set the stop time to 1000.

sd = simfunc(kGdValues,1000);

Plot the simulation results to see how varying the level of kGd affects the level of Ga.

sbioplot(sd);

See Also

|

Related Topics

Was this topic helpful?