kde2dとperspを使用して、2 つの変数 x と y の結合密度をプロットしています。
どうすれば次のことができますか:
1) x と y の値に応じて、プロットの色を変えたい。つまり、値 x>.5 および y>.5 の場合、色は青、x<=.5 および y<=.5 の場合、色は赤になります。現在、グラフ全体が青色になっています (以下のコードを参照)。
2) 表面下の領域 (つまり、確率質量) に色を付けるにはどうすればよいですか? 青と赤でそれぞれ色付けしたいと思います。
library(MASS)
i<-1
x<-NULL
y<-NULL
while (i<=30) {
x1 <- rnorm(1, 0, 1)
x2 <- rnorm(1, 0, 1)
x <- c(x, x1)
y<-c(y,x2)
i<-i+1
}
den3d <- kde2d(x, y)
persp(den3d,xlab="Var1", ylab="Var2", zlab="Density", shade = 0.1, col = "blue", box=TRUE)
助けていただければ幸いです。