1

可能であれば信頼区間を使用して、回帰の非パラメトリック部分を散布図に描画しようとしています。やってみても失敗した、ここで質問する理由

library(np)
library(foreign)  

wage1 <- read.dta(file="PATH HERE\\wage1.dta")
reg.np <- npplreg(lwage ~ female + married + educ + tenure | exper, data = wage1)

plot(reg.np)

plot(wage1$exper, wage1$lwage, xlab = "age", ylab = "log(wage)")
lines(wage1$exper, fitted(reg.np), lty = 1, col = "blue")

これは、私が試みた多くの試みの 1 つです。誰かが助けてくれれば本当にありがたいです。

賃金1 データベース: http://www.sendspace.com/file/txksv3

4

1 に答える 1

1

wage1データにアクセスできませんが、適切なx軸値を含む新しいデータフレームを作成し、これを引数としてpredict（）を使用して適切なy軸値を取得するために何をしたいかを考えます。これらのx軸とy軸の値を使用して、モデルをプロットに追加します。

fited（）を使用する場合は、x変量の昇順でデータを並べ替える必要がある場合があります。

于 2012-11-21T20:11:34.473 に答える