-5

関数行に問題があります。

これは私がこれまでに書いたものです：

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

私の問題は、プロットに線が表示されず、その理由がわからないことです。

誰か助けてくれませんか？

4

1 に答える 1

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() を使用してプロットを 1 回呼び出すだけです。

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

それだけです。 ここに画像の説明を入力

于 2013-01-13T23:32:09.697 に答える