r - 功能行（）不工作

Question

我的功能线有问题。

这是我到目前为止所写的：

model.ew<-lm(Empl~Wage)

summary(model.ew)

plot(Empl,Wage)

mean<-1:500 

lw<-1:500

up<-1:500

for(i in 1:500){

  mean[i]<-predict(model.ew,data.frame(Wage=i*100),interval="confidence",level=0.90)[1]

  lw[i]<-predict(model.ew,data.frame(Wage=i*100),interval="confidence",level=0.90)[2]

  up[i]<-predict(model.ew,data.frame(Wage=i*100),interval="confidence",level=0.90)[3]

}

plot(Wage,Empl)

lines(mean,type="l",col="red")

lines(up,type="l",col="blue")

lines(lw,type="l",col="blue")

我的问题是我的情节上没有出现任何线条，我不知道为什么。

有人可以帮助我吗？

Answer 1

您确实需要阅读 R 的一些介绍性手册。转到此页面，然后选择一个说明使用 R 进行线性回归的手册： http: //cran.r-project.org/other-docs.html

首先我们需要做一些数据：

set.seed(42)
Wage <- rnorm(100, 50)
Empl <- Wage + rnorm(100, 0)

现在我们运行您的回归并绘制线条：

model.ew <- lm(Empl~Wage)
summary(model.ew)
plot(Empl~Wage) # Note. You had the axes flipped here

你的第一个问题是你翻转了轴。因变量 (Empl) 在垂直轴上。这就是你在情节上没有得到任何台词的主要原因。要获得预测线根本不需要循环，只需使用 matlines() 进行一次绘图调用：

xval <- seq(min(Wage), max(Wage), length.out=101)
conf <- predict(model.ew, data.frame(Wage=xval), 
    interval="confidence", level=.90)
matlines(xval, conf, col=c("red", "blue", "blue"))

这里的所有都是它的。