gillespied(smfsb)
gillespied()所属R语言包:smfsb
Simulate a sample path from a stochastic kinetic model described by a stochastic Petri net
从一个随机Petri网的随机动力学模型所描述的,模拟的样本路径
译者:生物统计家园网 机器人LoveR
描述----------Description----------
This function simulates a single realisation from a discrete stochastic kinetic model described by a stochastic Petri net and discretises the output onto a regular time grid.
此功能实现从离散的随机Petri网的随机动力学模型和模拟的是单discretises输出到一个固定的时间格。
用法----------Usage----------
gillespied(N, T=100, dt=1, ...)
参数----------Arguments----------
参数:N
An R list with named components representing a stochastic Petri net (SPN). Should contain N$M, a vector representing the initial marking of the net, N$Pre, a matrix representing the LHS stoichiometries, N$Post, a matrix representing the RHS stoichiometries, and N$h, a function representing the rates of the reaction processes. N$h should have first argument x, a vector representing the current state of the system, and second argument t, a scalar representing the current simulation time (in the typical time-homogeneous case, N$h will ignore this argument). N$h may posses additional arguments, representing reaction rates, for example.
R列表代表一个随机Petri网(SPN)的命名的组件。应包含N$M,一个向量,表示净值的初始标识,N$Pre,代表了LHS化学计量矩阵,N$Post,矩阵的RHS化学计量,和N$h的,代表的反应过程的速率的函数。 N$h应该有第一个参数x,一个向量,表示该系统的当前状态,和第二个参数t,一个标量,表示当前的仿真时间(在典型的时间均匀的情况下, N$h将忽略此参数)。 N$h可能拥有额外的参数,代表反应速率,例如。
参数:T
The required length of simulation time. Defaults to 100 time units.
所需长度的模拟时间。默认为100个时间单位。
参数:dt
The grid size for the output. Note that this parameter simply determines the volume of output. It has no bearing on the correctness of the simulation algorithm. Defaults to one time unit.
输出的网格大小。请注意,这个参数决定的产出量。它有没有影响的模拟算法的正确性。默认为一个单位时间。
参数:...
Additional arguments will be passed into the function N$h.
额外的参数将被传递到函数N$h。
值----------Value----------
参见----------See Also----------
simpleEuler, rdiff,
simpleEuler,rdiff,
实例----------Examples----------
# load LV model[加载LV模型]
data(spnModels)
# simulate and plot a realisation[模拟和绘制实现]
plot(gillespied(LV,T=100,dt=0.01))
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|