代码之家  ›  专栏  ›  技术社区  ›  Karl

如何绘制观测时间序列上的拟合模型

  •  1
  • Karl  · 技术社区  · 16 年前

    model <- lm(x~y+z)
    plot(x)
    lines(fitted(model))
    

    3 回复  |  直到 9 年前
        1
  •  8
  •   Dirk is no longer here    16 年前

    abline(model)

     z <- lm(dist ~ speed, data = cars)
     plot(cars)
     abline(z) # equivalent to abline(reg = z) or
     abline(coef = coef(z))
    
        2
  •  5
  •   Mikko    8 年前
    x <- rnorm(100)
    y <- rnorm(100)
    z <- rnorm(100)
    
    model <- lm(x~y+z)
    plot(x,type="l",col="green")
    lines(fitted(model),col="blue")
    

        3
  •  0
  •   adamsss6    8 年前

    n = 100; mi = 0; sigma = 2
    x = rnorm(n,mi,sigma)
    e = rnorm(n,0,1)
    b0 = 1; b1 = 2
    y = b1*x + b0 + e
    #plot observations
    plot(x,y)
    #model
    lm_res= lm(y~x)
    summary(lm_res)
    arg= c(min(x),max(x))
    out = coef(lm_res)[2]*arg+ coef(lm_res)[1]
    #plot model line
    lines(arg, out, col = 'red')