Linear Regression

Figure2.1 in The Elements of Statistical Learning を書いてみました。 前回同様Input SpaceはX = {\bf R}^{2}であることに注意!

library(ElemStatLearn)

x = mixture.example$x
g = mixture.example$y
xnew = mixture.example$xnew

X = cbind(matrix(1, dim(x)[1], 1), x)
Xnew = cbind(matrix(1, dim(xnew)[1], 1), xnew)
beta = solve(t(X) %*% X) %*% t(X) %*% g

yhat = Xnew %*% beta 
px1 = mixture.example$px1
px2 = mixture.example$px2
contour(px1, px2, matrix(yhat, length(px1), length(px2)), levels=0.5, axes=FALSE, xlab="", ylab="", labels="", main="Linear Regression of 0/1 Response")
points(x, col=ifelse(g==1, "coral", "cornflowerblue"))
points(expand.grid(px1, px2), pch=".", cex=1.2, col=ifelse(yhat>0.5, "coral", "cornflowerblue"))
box()

file:///Users/hiroyukitanaka/Dropbox/myR/plots/esl2.1.png

コメントを残す

以下に詳細を記入するか、アイコンをクリックしてログインしてください。

WordPress.com ロゴ

WordPress.com アカウントを使ってコメントしています。 ログアウト / 変更 )

Twitter 画像

Twitter アカウントを使ってコメントしています。 ログアウト / 変更 )

Facebook の写真

Facebook アカウントを使ってコメントしています。 ログアウト / 変更 )

Google+ フォト

Google+ アカウントを使ってコメントしています。 ログアウト / 変更 )

%s と連携中