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

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

plot r

1

Karl · 技术社区 · 16 年前

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

3 回复 | 直到 10 年前

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 9 年前

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')