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