• R语言

  • R怎么画一个跨期迭代含有时间t的图？

λt=γλt-1+μt 其中γ=0.95，μ1=0.01
怎么用R画出这个图，横坐标为时间t，纵坐标为λ。?

λt=γλt−1+μt\lambda_{t} = \gamma \lambda_{t-1} + \mu_{t}

还缺一个λ0\lambda_0 跟μt\mu_{t}的函数，总之我先把μt\mu_{t}拿掉，画一个λt=γλt−1\lambda_{t} = \gamma \lambda_{t-1}

library(purrr)

x=1:1000
dummy <- function(acc,nxt){
  return (0.95 * acc )
}
y= x %>% accumulate(dummy)

plot(x,y)

    tctcab 是这样只有t=1的时候，有μ=0.01，之后就没有μ了

      tctcab 如果加上初始值λ=0，t=1时的μ=0.01，怎么画啊？我思考了好久，还是没有出来

          ttdds

          library(purrr)
library(ggplot2)
library(dplyr)
x=0:1000

dummy <- function(acc,nxt){
  if (nxt ==1) mu =0.01
  else mu = 0

  return (0.95 * acc + mu)
}
y= x %>% accumulate(dummy)

df <- data.frame(x=x,y=y)

gg <- ggplot(data=df %>% filter(x > 0),aes(x=x,y=y))+
  geom_point()+
  xlab("t") +
  ylab(expression(paste(lambda[t])))
print (gg)

              tctcab 谢谢！如果再嵌套一个新的函数K，请问怎么将K Y C画在一个图里？
              λt=0.95λt-1+μt μ的初始值第一期还是0.01，和前面一样

              Kt+1=0.9Kt+0.1λt （K Y C的初始值都是0）
              Yt=0.2Kt+1.4λt
              Ct+1=0.6Kt+0.4λt