在R中绘制函数

Question

我正在对我的数据进行QDA，需要在R或MATLAB中绘制决策边界，这看起来像下面的函数。

0.651 - 0.728(x_1) - 0.552(x_2) - 0.006(x_1 * x_2) - 0.071(x_1)^2 +.170 (x_2)^2 = 0

任何人都可以帮助我吗？我一直在搜索互联网，似乎无法找到解决方案。

Answer 1

R没有绘制这样的隐式函数的函数，但你可以用等高线图很好地伪造它。假设您想绘制区域[0,1] ^ 2（可以很容易地更改），您需要一些像这样的代码：

#f is your function
f <- function (x_1, x_2) 0.651 - 0.728 * x_1 - 0.552 * x_2 - 0.006 * x_1 * x_2 - 0.071 * x_1^2 +.170 * x_2^2

#length=1000 to be safe, could be set lower
x_1 <- seq(0, 1, length=1000)
x_2 <- seq(0, 1, length=1000)
z <- outer(x_1, x_2, FUN=f)
contour(x_1, x_2, z, levels=0)

Answer 2

使用ggplot同样的事情：

library(ggplot2)
library(reshape2)   # for melt(...)
## same as first response
x <- seq(0,10,length=1000)
y <- seq(-10,10, length=1000)
z <- outer(x,y,FUN=f)
# need this extra step 
gg <- melt(z, value.name="z",varnames=c("X","Y"))
# creates the plot
ggplot(gg) + stat_contour(aes(x=X, y=Y, z=z), breaks=0)