R：在輪廓中填充封閉區域

Question

在R中使用繪圖繪製knn分類圖以繪製樣本和輪廓以繪製將平面分類的線。

這裏是我的代碼：

k<-1 
datax<-rbind(matrix(rnorm(30,-1,5.25),15,2),matrix(rnorm(36,1,5.25),18,2)) 
datay<-rbind(matrix(1,15,1),matrix(0,18,1)) 
plot(datax[,1], datax[,2],pch = datay+1,axes=FALSE,ann=FALSE) 
box() 
n <- 1000 
xp <- seq(length=n, from = min(datax[,1]), to = max(datax[,1])) 
yp <- seq(length=n,from = min(datax[,2]) ,to = max(datax[,2])) 
gr <- expand.grid(xp, yp) 
library(class) 
z <- as.numeric(knn(datax, gr, datay,k))-1 
zM <- matrix(z, n, n, byrow = FALSE) 
contour(xp, yp, zM, xlab="x",ylab="",nlevels = 1 ,lwd=2, add=TRUE, drawlabels =FALSE) 

我的問題是：我怎樣才能顏色在劇情的封閉區域？我試過filled.contour，但沒有添加參數。我只想要分類器= 0的區域爲白色，並且其分類爲藍色的區域爲1。我應該怎麼做？ 感謝

Answer 1

相反的contour，您可以使用contourLines保持輪廓線的邊緣的座標，並與polygon繪製出來。

plot(datax[,1], datax[,2],axes=FALSE,ann=FALSE, type="n") 
box() 
cL <- contourLines(xp, yp, zM,nlevels = 1) 
lapply(cL,function(x)polygon(x$x,x$y,col="red")) 
points(datax[,1], datax[,2],pch = datay+1) 

但是它並不完美與達到情節的邊緣輪廓線（見第二個圖的左下角），所以它需要一些手工製作的調整：

編輯：在嵌套的輪廓線，我不認爲有一個簡單的方法來對付它，但這裏的情況是單向的：

library(splancs) 
ord <- sapply(lapply(cL,function(x)datay[inout(datax,cbind(x$x,x$y))]), 
       median) #Check what values are present in the polygon and 
         #take the most common one 
plot(datax[,1], datax[,2],axes=FALSE,ann=FALSE, type="n") 
box() 
lapply(cL[ord==1],function(x)polygon(x$x,x$y,col="blue")) 
lapply(cL[ord==0],function(x)polygon(x$x,x$y,col="white")) 
points(datax[,1], datax[,2],pch = datay+1) 

2日編輯：當然還有也是在你的情況使用功能image的可能性：

image(xp, yp, zM, col=c("transparent","blue")) 
points(datax[,1], datax[,2],pch = datay+1)