r - R: 不等式関数 (3 変数) のグラフのプロット方法

Question

私の不等式:

x^2 + y^2 + z^2 < 1 AND x^2 + y^2 < z^2

この不等式の論理的な組み合わせを R でプロットする方法は? (3 つの軸すべての範囲 = -1,1)

Answer 1

これは、3Dプロットを作成し、条件によって色付けされたポイントを追加するためのrglソリューションです。

# some data on a grid
x = seq(-1, 1, len = 10); 
df <- expand.grid(x=x,y=x,z=x)
# indicator for color
df$ind <- with(df, x^2 + y^2 + z^2 < 1 & x^2 + y^2 < z^2)

require(rgl)
# empty plot
plot3d(df$x, df$y, df$z, type = 'n')
# add points
with(df[df$ind, ], points3d(x, y, z, color = 'red', size = 10))
with(df[!df$ind, ], points3d(x, y, z, color = 'blue', size = 10))

Answer 2

使用rgl:

x <- y <- z <- seq(-1, 1, by=0.01)
df <- setNames(expand.grid(x, y, z), c("x", "y", "z"))
df <- transform( df, ueq = (x^2 + y^2 + z^2 < 1) & (x^2 + y^2 < z^2))
df$color <- ifelse(df$ueq == TRUE, "green" , "red")
require(rgl)
with(df[df$ueq == TRUE, ], plot3d(x=x, y=y, z=z, col=color, type="p", size=5))
grid3d(c("x", "y+", "z"))