fore_p: Simulation plot of enzyme kinetics model

In EKMCMC: MCMC Procedures for Estimating Enzyme Kinetics Constants

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)

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.

Value

This functio has no returned object.

Examples

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)

Example output

There were 50 or more warnings (use warnings() to see the first 50)

EKMCMC documentation built on May 2, 2019, 3:47 p.m.