Question

我想在下面的代码中对两个底线之间的区域进行着色。但我所能达到的只是连接abline的两个截距。有一个简单的方法来处理这个？更复杂的解决方案也非常受欢迎！

x1s <- c(.5,1,1,2,3,3.5,     1,3.5,4,5,5.5,6)
x2s <- c(3.5,1,2.5,2,1,1.2,  5.8,3,4,5,4,1)
ys <- c(rep(+1,6),          rep(-1,6))
my.data <- data.frame(x1=x1s, x2=x2s, type=as.factor(ys))
my.data

library('e1071')
svm.model <- svm(type ~ ., data=my.data, type='C-classification', kernel='linear',scale=FALSE)

plot(my.data[,-3],col=(ys+3)/2, pch=19, xlim=c(-1,6), ylim=c(-1,6))
points(my.data[svm.model$index,c(1,2)],col="blue",cex=2) # show the support vectors

w <- t(svm.model$coefs) %*% svm.model$SV
b <- -svm.model$rho
p <- svm.model$SV

abline(a=-b/w[1,2], b=-w[1,1]/w[1,2], col="black", lty=1)
abline(a=(-b-1)/w[1,2], b=-w[1,1]/w[1,2], col="orange", lty=3)
abline(a=(-b+1)/w[1,2], b=-w[1,1]/w[1,2], col="orange", lty=3)
polygon(c(0,0), c((-b-1)/w[1,2],(-b+1)/w[1,2]), col = 'grey80', border = NA)

Answer 1

我这样做的方法是使用函数curve而不是abline，以便保持两条线的路径坐标并将它们提供给`polygon：

c1 <- curve((-b-1)/w[1,2] - x*w[1,1]/w[1,2], from=-10, to=10, 
                              col="orange",add=TRUE,lty=3)
c2 <- curve((-b+1)/w[1,2] - x*w[1,1]/w[1,2], from=-10, to=10, 
                              col="orange",add=TRUE,lty=3)
polygon(c(c1$x,rev(c2$x)), c(c1$y, rev(c2$y)),col="grey80", border=NA)

修改：正如Julian Urbano正确建议的那样，您可以使用par('usr')自动定义参数from和to的值：

u <- par('usr')
c1 <- curve((-b-1)/w[1,2] - x*w[1,1]/w[1,2], from=u[1], to=u[2], 
                              col="orange",add=TRUE,lty=3)
c2 <- curve((-b+1)/w[1,2] - x*w[1,1]/w[1,2], from=u[1], to=u[2], 
                              col="orange",add=TRUE,lty=3)
polygon(c(c1$x,rev(c2$x)), c(c1$y, rev(c2$y)),col="grey80", border=NA)

R中两个'ablines'之间的阴影区域

1 个答案: