EKMCMC (version 0.1.0)

fore_p: Simulation plot of enzyme kinetics model

Description

The function depicts the overlayed two plots; one is observed data, the other is simulation result using fitted MM constant and catalytic constant.

Usage

fore_p(method = TRUE, CL = TRUE, time, species, enz, subs, MM, catal,
  nrepeat = 100, ti = 1)

Arguments

method

method selection: T=TQ model, F=SQ model(default = T)

CL

Adding empircal 95% confidence interval (default = T)

time

observed time interval

species

observed trajectory of product

enz

enzyme concentration

subs

substrate concentration

MM

true value of MM constant

catal

initial value of catalytic constant

nrepeat

total number of simulation (default=100)

ti

tme interval for descreted simulatin result (default =1)

Value

This functio has no returned object.

Details

Basically this function draws overlayed picture: The one is trajectory of given data of products for enzyme kinetics model. The other is trajectory of products from simulation result of Gillespie algorithm with estimated two constants, Michaelis-Menten constant and catalytic constant. CL option controls the plot that depicts the observed data and mean of simulated series only or adding 95 confidence interval with 10 samples of simulated trajectory.

Examples

# NOT RUN {
data("Chymo_low")
time1=Chymo_low[,1]
species1=Chymo_low[,2]
fore_p(method=TRUE, CL=TRUE, time=time1,species=species1,enz=4.4e+7, subs=4.4e+7
       ,MM=4.4e+8, catal=.051)
# }