R 函数行()不起作用

R 函数行()不起作用,r,plot,lines,R,Plot,Lines,我对功能行有问题 这是我到目前为止写的: 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]<

我对功能行有问题

这是我到目前为止写的:

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

model.ew您确实需要阅读一些关于R的介绍性手册。转到本页,选择一个说明如何使用R进行线性回归的手册:

首先,我们需要制作一些数据:

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

工资从哪里来?请包括示例数据。
predict.lm
已矢量化。无需为
循环设置
。研究
?predict.lm
中的示例。此外,您希望将x和y值传递给
行,而不仅仅是x。
model.ew <- lm(Empl~Wage)
summary(model.ew)
plot(Empl~Wage) # Note. You had the axes flipped here
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"))