R：ggplot距離式

Question

ggplotを使用して散布図の下に等濃度線をプロットしようとしていますが、stat_functionを正しく使用する方法がわかりません。

アイソクラインは、距離の式に基づいている：

sqrt((x1-x2)^2 + (y1-y2)^2) 

中心プロットの原点であろう除き、これら 同心円ようになります。

私が今まで試したことは、ggplot内のdistance関数をそう呼んでいることです（注：実際の問題で固定値があるので、x1 = 1とy1 = 1を使用します）

distance <- function(x, y) {sqrt((x - 1)^2 + (y - 1)^2)} 
ggplot(my_data, aes(x, y))+ 
    geom_point()+ 
    stat_function(fun=distance) 

が、Rがエラーを返した：

Computation failed in 'stat_function()': argument "y" is missing, with no default

それが原点に中心に、距離の式の一般的なプロットをプロットするようにどのように私は正しくstat_functionにxとyの値を養うのですか？

Answer 1

少し複雑なものについては、私はstat機能の使用を避けます。彼らは主に迅速な計算を目指しています。通常、xに基づいてyと計算されます。私はちょうど代わりstat_contourとデータおよびプロットを事前に計算します：

distance <- function(x, y) {sqrt((x - 1)^2 + (y - 1)^2)} 
d <- expand.grid(x = seq(0, 2, 0.02), y = seq(0, 2, 0.02)) 
d$dist <- mapply(distance, x = d$x, y = d$y) 

ggplot(d, aes(x, y)) + 
    geom_raster(aes(fill = dist), interpolate = T) + 
    stat_contour(aes(z = dist), col = 'white') + 
    coord_fixed() + 
    viridis::scale_fill_viridis(direction = -1)