The University of Sheffield
Neil Lawrence ML@SITraN
Probabilistic Dynamical Model for Inference of the Regulatory Mechanism of Transcription - MATLAB Code

Dynamical Model for Transcription

This page describes examples of how to use the Dynamical Model for Transcription detailed in this paper. The code is available for download here.

Release Information

Current release is 0.11.

Release 0.11 is the second release associated with the paper, containing scripts for recreating the results given.

The toolboxes required to run the code are listed below.

The CHIPDYNO software can be downloaded here.

Release Information

Current release is 0.121.

As well as downloading the CHIPDYNO software you need to obtain the toolboxes specified below. These can be downloaded using the same password you get from registering for the CHIPDYNO software.
Toolbox Version
NDLUTIL 0.159
NETLAB 3.3

Version 0.12

Release version with some files missing (blank files for several demSpellman scripts).

Finally you will also need the NETLAB toolbox in your path.

Results in the Paper

The main scripts are demTu and demSpellman which run the code on the metabolic and cell cycle data sets. The function chipdynoTransFact computes the posterior, gene specific, transcription factor activities for a given transcription factor. The function chipdynoTransFactNoise does the same thing for expression levels with associated noise (such as computed by mmgMOS). The function chipdynoGeneAct sorts the transcription factor activities for a particular gene. The function chipdynoTransFactCorrel normalises the correlations and presents the most correlated TFs. Some scripts have been added to recreate the further validatory experiments detailed in the Supplementary material. Specifically,demSpellmanp1, demSpellmanp2 and demSpellmanp5 run the code on the cell cycle data set with connectivity matrix obtained using different cut-offs at p=0.001, p=0.0002 and p=0.005 respectively. demSpellmanp1p2 and demSpellmanp5p2 do the same thing but only use the genes and transcription factors active at p=0.0002, so that the systems have the same dimensions. The script chipDynoCompareX can then be used to compare the matrices of effective regulations in the different cases.

Page updated on Fri Mar 14 17:11:24 2008