Question

我想问的问题是如何在经过一段时间之后更改我的变量的一个值（在我的列表中记为LO $ M）。

我想要达到的目标是，经过20,000秒后，我想将Lac的值改为Lac的值，时间为20,0000 +10,000

因此在t = 20,000时，Lac = Lac + 10,000

我对我的代码所遇到的问题是，如果tt＆gt; = 20000，我的if命令就会出现问题，但这会导致Lac的每个值在20,000之后增加10,000，我想要的是FIRST 20,000之后的价值增加10,000。

基本上，在我的实验经过20,000次后，我试图在实验中注入10,000多个Lac。

我的代码如下：

LO = list()
LO$M = c(i = 1, ri = 0, I = 50, Lac = 20, ILac = 0, o = 1, Io = 0, RNAP = 100, RNAPo = 0, r = 0, z = 0)
LO$Pre = matrix(c(1,0,0,0,0,0,0,0,0,0,0,
                  0,1,0,0,0,0,0,0,0,0,0,
                  0,0,1,1,0,0,0,0,0,0,0,
                  0,0,0,0,1,0,0,0,0,0,0,
                  0,0,1,0,0,1,0,0,0,0,0,
                  0,0,0,0,0,0,1,0,0,0,0,
                  0,0,0,0,0,1,0,1,0,0,0,
                  0,0,0,0,0,0,0,0,1,0,0,
                  0,0,0,0,0,0,0,0,1,0,0,
                  0,0,0,0,0,0,0,0,0,1,0,
                  0,0,0,1,0,0,0,0,0,0,1,
                  0,1,0,0,0,0,0,0,0,0,0,
                  0,0,1,0,0,0,0,0,0,0,0,
                  0,0,0,0,1,0,0,0,0,0,0,
                  0,0,0,0,0,0,0,0,0,1,0,
                  0,0,0,0,0,0,0,0,0,0,1), ncol=11, byrow=TRUE)
LO$Post = matrix(c(1,1,0,0,0,0,0,0,0,0,0,
                   0,1,1,0,0,0,0,0,0,0,0,
                   0,0,0,0,1,0,0,0,0,0,0,
                   0,0,1,1,0,0,0,0,0,0,0,
                   0,0,0,0,0,0,1,0,0,0,0,
                   0,0,1,0,0,1,0,0,0,0,0,
                   0,0,0,0,0,0,0,0,1,0,0,
                   0,0,0,0,0,1,0,1,0,0,0,
                   0,0,0,0,0,1,0,1,0,1,0,
                   0,0,0,0,0,0,0,0,0,1,1,
                   0,0,0,0,0,0,0,0,0,0,1,
                   0,0,0,0,0,0,0,0,0,0,0,
                   0,0,0,0,0,0,0,0,0,0,0,
                   0,0,0,1,0,0,0,0,0,0,0,
                   0,0,0,0,0,0,0,0,0,0,0,
                   0,0,0,0,0,0,0,0,0,0,0), ncol=11, byrow=TRUE)
LO$h = function(x,t,th=c(0.02,0.1,0.005,0.1,1,0.01,0.1,0.01,0.03,0.1,1e-05,0.01,0.002,0.01,0.001))
{
  with(as.list(c(x, th)), {
    return(c(th[1]*i, th[2]*ri, th[3]*I*Lac, th[4]*ILac, th[5]*I*o, th[6]*Io, th[7]*o*RNAP,
             th[8]*RNAPo, th[9]*RNAPo, th[10]*r, th[11]*Lac*z, th[12]*ri, th[13]*I, 
             th[13]*ILac, th[14]*r, th[15]*z))
  })
}

gillespie1 = function (N, n, ...) 
{
  tt = 0
  x = N$M
  S = t(N$Post - N$Pre)
  u = nrow(S)
  v = ncol(S)
  tvec = vector("numeric", n)
  xmat = matrix(ncol = u, nrow = n + 1)
  xmat[1, ] = x
  for (i in 1:n) {
    h = N$h(x, tt, ...)
    tt = tt + rexp(1, sum(h))
    j = sample(v, 1, prob = h)
    x = x + S[, j]
    tvec[i] = tt
    xmat[i + 1, ] = x
    if( tt >=20000){
      x[4] = x[4] +10000
    }
  }
  return(list(t = tvec, x = xmat))
}

newout = gillespie1(LO,200000)
matplot(newout$x[,4], type="l", lwd=0.25, col="grey")

我没有足够的声誉来附加图片，但它看起来应该是这样的： https://gyazo.com/0ffd940a22df23b2ccfdf4a17e85dca8

对不起，如果不清楚的话。感谢

Answer 1

在此示例中，您具有函数myTask()。当您致电execMyTask()时，您将执行myTask()一次，之后，您将以1到max_wait毫秒之间的随机间隔执行它。当您感到疲倦时，可以使用tclTaskDelete()终止任务。

library(tcltk2)

myTask <- function() cat("some task!\n")

id = "execMyTask"

execMyTask <- function(max_wait = 3000) {
  id <- toString(match.call()[[1]])
  myTask()
  wait = sample(1:max_wait, 1)
  cat("Waiting", wait, "miliseconds\n")  # replace with your function

  if (is.null(tclTaskGet(id))) {
    tclTaskSchedule(wait=wait, execMyTask(), id=id, redo = TRUE)
  } else {
    tclTaskChange(wait=wait, execMyTask(), id=id, redo = TRUE)
  }

}

execMyTask()
tclTaskDelete(id)

到目前为止，这种方法存在一些问题，因为我们无法在tclTaskChange()中为函数 fun 提供参数。

R：在一定时间后改变变量值

1 个答案: