User:Vincent Rouilly/Computational Biology With R: Difference between revisions
From OpenWetWare
Jump to navigationJump to search
No edit summary |
|||
Line 10: | Line 10: | ||
# Load library 'SBMLR' in R by typing: | # Load library 'SBMLR' in R by typing: | ||
## > library("SBMLR") | ## > library("SBMLR") | ||
# | # Download sample SBML file describing a A-->B reaction: [[media:here]] | ||
# Read SBML file into R | |||
## > myModel <- readSBML("AtoB.xml") | |||
# Simulate model between [0,100] with 100 points | |||
## > results=simulate(myModel,seq(0,100,1)) | |||
# Plot results | |||
## > attach(results) | |||
## > par(mfrow=c(2,1)) | |||
## > plot(time,s1,type="l") | |||
## > plot(time,s2,type="l") | |||
## > detach(results) | |||
===Perturbation analysis on compound concentration=== | |||
===Perturbation analysis on kinetic rate=== |
Revision as of 10:39, 14 December 2009
Computational Biology with R
Tutorials
Running SBML models in R
- Install SBMLR package
- Start R, and enter:
- > source("http://bioconductor.org/biocLite.R")
- > biocLite("SBMLR")
- Load library 'SBMLR' in R by typing:
- > library("SBMLR")
- Download sample SBML file describing a A-->B reaction: media:here
- Read SBML file into R
- > myModel <- readSBML("AtoB.xml")
- Simulate model between [0,100] with 100 points
- > results=simulate(myModel,seq(0,100,1))
- Plot results
- > attach(results)
- > par(mfrow=c(2,1))
- > plot(time,s1,type="l")
- > plot(time,s2,type="l")
- > detach(results)