我在库(ks)中使用plot.kde来提取核密度图的轮廓水平。我想叠加多个图,这样我就可以使轮廓填充半透明。然而,有一个边界/轮廓线,其颜色我似乎无法控制。
我尝试更改plot.kde函数中所有不同的col,cont.col,color(etc)选项,但似乎无法磨练轮廓本身的颜色。我可以通过从kde对象中提取轮廓坐标然后使用polygon()函数绘制它来使用一些解决方法,但我真的想从plot.kde中控制它。这是我将要经营多次的事情。
这可能非常简单,但我只是错过了它!在下图中,它是我试图控制的较粗的红线。
感谢!!!
library(ks)
data(unicef)
H.scv <- Hscv(x=unicef)
fhat <- kde(x=unicef, H=H.scv, compute.cont=TRUE)
plot(fhat, display="filled.contour2", cont=c(10),col=c(NA,rgb(1,0,0,0.5)))