从Gompertz方程得到中点值

Question

time <- 1:12
y <- c(0,0,0,0,0,0,0,12,34,69,100,100)
mdl <- gcFitModel(time,y,control = grofit.control(fit.opt = "m", model.type = "gompertz"))

我从上面的mdl得到我的参数以适应gompertz方程式

y <- A* exp(-exp(mu * e * (lambda - time)/mu + 1))

mu <- 36.162016
lambda <- 7.9800164
A <- 100
time <- 1:12

time这是15天的一步。例如，time = 1表示15天的中午，time 2表示接下来15天的中间日期，时间3表示暗示接下来15天的中午，依此类推。

我安装了以下功能：

  e <- exp(1)
  y <- 100 * exp(-exp(mu * e * (lambda - time)/mu + 1))
  plot(time,y)

lambda控制沿x轴的运动。 我希望修改这条曲线，以便通过将ID转换为数周而不是每15天的中点来获得更多数据点，我希望每7天获得一次y。我怎样才能做到这一点？

Answer 1

您的代码并未包含任何适合例程，因此我假设这是一个关于绘图的问题。

以下是time = 1和time = 12之间100点的情节示例。

time <- seq(1, 12, length.out = 100);

mu <- 34.55844
y <- 100 * exp(-exp(mu * e * (3 - time)/mu + 1))
plot(time,y)