R言語の条件式について

Question

R言語を用いて，プログラミングを行っているのですが条件式の書き方が分からず迷っています．
二次元カーネル密度推定を行った結果に閾値を設定したいため条件式を用いようと思ったのですが調べても出ず迷っています．

プログラム内容は以下の通りです．
x=scan("clipboard")
y=scan("clipboard")
plot(x,y,xlab="latitude",ylab="longitude")
install.packages("MASS", repos="http://cran.ism.ac.jp/")
library(MASS)
d=kde2d(x,y,c(bandwidth.nrd(x),bandwidth.nrd(y)),n=80)
a<-mean(d$z)

aよりも値が大きいものだけを抽出したいのですが方法がわからないため教えていただきたいです．

Answer 1

d[d > a] で良いかと思います。

R では
ベクトルx > ベクトルy
で「xの各成分がyの対応する成分よりも大きいか」という条件式になっており
TRUE, FALSE, TRUE, TRUE, ...
のようなものが格納されます。

つまり、「d[d > a]」で
「dの成分のうち、aの対応する成分よりも大きいもの」
を抽出することができます。

Comments

[qtakusora1014]

> d[d>a]
Error: (list) object cannot be coerced to type 'double'
となりエラーが起きてしまうんです．．．・

[qtakusora1014]

お忙しい中申し訳ございませんが，アドバイスをいただけたら幸いです．．

[Victorian-Vat]

では
　z <- unlist(d$z)
とした上で
　z[z > a]
としてみてください。

[qtakusora1014]

抽出することはその方法でできたのですが，それに対応するｘ，ｙも保管したいのですが，できますか？
d$xに緯度，d$yに経度が入っています．
persp(x1, y1, z1, theta = 0, phi = 0, col = "red")
のような形で，最終的にグラフを絵画したいです，，