Imperial College/Courses/Fall2008/Synthetic Biology (MRes class)/'R' Tutorial/Practical 2: Difference between revisions

From OpenWetWare
Jump to navigationJump to search
No edit summary
Line 19: Line 19:


==Exercise 1 (Data frames and inputs/ouputs)==
==Exercise 1 (Data frames and inputs/ouputs)==
In R, files can be used for different purposes such as storing your command script, storing your results, or load data previously generated.
* '''Storing your command script:'''
** create a .txt file
** write R commands (1 per line), like creating a random vector of uniform data between [-10,10], and create a sequence of integers between [1,100], and finally plot the resulting scatterplot.
** save the file
** in R, use the source("path_to_file") command to execute the commands stored into the file.
* '''Storing your results (with a CSV format for example):'''
* '''Loading data previously generated:'''


==Exercise 2 (ODE Solver)==
==Exercise 2 (ODE Solver)==
Line 37: Line 46:


==Exercise 3 (Non linear regression)==
==Exercise 3 (Non linear regression)==
Some synthetic data has been generated to mimic the characterisation of a repressed gene construct. The experiment consisted in measuring the reporter concentration at steady-state, expressed by the repressed construct, as a function of different repressor concentration.
The measurements are as follow:
Repressor, Reporter
'''Questions''':
# Plot the measurements: Reporter = f(Repressor)
# In order to analyse the data, you need to build an hypothetical underlying model of the repressed gene construct. Using ODEs, write down the equations governing the dynamics of the system. In order to model the repressor effect on the transcription, you will use a Hill function.
# From the previous ODEs, define the steady-state equilibrium of the system. What are the parameters that characterise the repressed gene construct ?
# Using a non linear regression analysis, fit your model to the measurements. Plot on top of the measurements, the resulting fit.
# Convert your system of ODEs so that it becomes representative of a PoPS device. Plot its transfer function.

Revision as of 02:56, 5 December 2008

Fall 2008 - Synthetic Biology (MRes class)

Home        'R' Tutorial        Resources        Literature

<html> <body> <!-- Start of StatCounter Code --> <script type="text/javascript"> var sc_project=3315864; var sc_invisible=0; var sc_partition=36; var sc_security="8bb2efcd"; </script>

<script type="text/javascript" src="http://www.statcounter.com/counter/counter_xhtml.js"></script><noscript><div class="statcounter"><a class="statcounter" href="http://www.statcounter.com/"><img class="statcounter" src="http://c37.statcounter.com/3315864/0/8bb2efcd/0/" alt="blog stats" /></a></div></noscript> <!-- End of StatCounter Code -->

</body> </html>

Introduction to 'R'



Exercise 1 (Data frames and inputs/ouputs)

In R, files can be used for different purposes such as storing your command script, storing your results, or load data previously generated.

  • Storing your command script:
    • create a .txt file
    • write R commands (1 per line), like creating a random vector of uniform data between [-10,10], and create a sequence of integers between [1,100], and finally plot the resulting scatterplot.
    • save the file
    • in R, use the source("path_to_file") command to execute the commands stored into the file.
  • Storing your results (with a CSV format for example):
  • Loading data previously generated:

Exercise 2 (ODE Solver)

  • Model: Constitutive gene expression
    • ODE System:
      • d[mRNA]/dt = k_1 - d_1*[mRNA]
      • d[protein]/dt = k_2[mRNA] -d_2*[protein]
  • Questions:
    • 1.Considering that the half-life for the mRNA and the protein is respectively 3min and 1h. Workout the values of (k_1, d_1, k_2, d_2) so that, at steady state, [mRNA] = 3 and [protein] = 500.
    • 2.Using the previously found parameter values, plot on the same graph [mRNA](t) and [protein](t) for t=[0, 5h].
    • 3. Parameter scanning: Consider parameter d_2 to vary [-20%, 20%] of its nominal value (10 uniformly spread values). Generate on the same graph the 10 [protein](t) trajectories (t=[0, 5h]).

Exercise 3 (Non linear regression)

Some synthetic data has been generated to mimic the characterisation of a repressed gene construct. The experiment consisted in measuring the reporter concentration at steady-state, expressed by the repressed construct, as a function of different repressor concentration. The measurements are as follow:

Repressor, Reporter



Questions:

  1. Plot the measurements: Reporter = f(Repressor)
  2. In order to analyse the data, you need to build an hypothetical underlying model of the repressed gene construct. Using ODEs, write down the equations governing the dynamics of the system. In order to model the repressor effect on the transcription, you will use a Hill function.
  3. From the previous ODEs, define the steady-state equilibrium of the system. What are the parameters that characterise the repressed gene construct ?
  4. Using a non linear regression analysis, fit your model to the measurements. Plot on top of the measurements, the resulting fit.
  5. Convert your system of ODEs so that it becomes representative of a PoPS device. Plot its transfer function.
Retrieved from "https://openwetware.org/mediawiki/index.php?title=Imperial_College/Courses/Fall2008/Synthetic_Biology_(MRes_class)/%27R%27_Tutorial/Practical_2&oldid=267505"