1 | p.path.init(s, x, theta, dimen, Jst.t = 0.64, Jst.rat = 0.5776972)
|
s |
|
x |
|
theta |
|
dimen |
|
Jst.t |
|
Jst.rat |
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | ##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (s, x, theta, dimen, Jst.t = 0.64, Jst.rat = 0.5776972)
{
passage <- matrix(0, dimen, 6)
colnames(passage) <- c("dimen", "tau", "y", "minI", "l",
"u")
for (i in 1:dimen) {
ith.passage <- bm.pass(s, x[i], theta[i], Jst.t, Jst.rat)
passage[i, ] <- c(i, ith.passage$tau, ith.passage$y,
ith.passage$minI, x[i] - theta[i], x[i] + theta[i])
}
list(curr.time = s, pass.mat = passage, cyc.mat = mat.sort.c(passage,
dimen, "tau"))
}
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.