2016-10-29 9 views
1

私は、a〜kの文字で表される11変数を持つデータフレームを2点散布図にプロットしています。私はすべての点を識別するために符号Rをプログラムしたい私は点(D、E、FおよびG、i)が30未満ユークリッド距離と軸の値に基づいて点を検索するR

d<-dist(cor, method = 'euclidean');d 

距離を有することが観察ユークリッド距離を用い

cor<-data.frame(X=c(0.36187115, -0.54755904, -0.82417308, -0.70806545, -0.77422866, -0.70003404, 
        -0.70043884, 0.73602124,-0.89909694, -0.05937341, 0.93496883), 
      Y=c(-0.54354070,-0.81211142, -0.52775892, 0.40191296, 0.36820779, 0.28163131, 
        -0.26161395, -0.26386668,-0.31894766, -0.91541962, -0.04548996), 
      row.names = letters[1:11]);cor 

a<-seq(0,2*pi, length=100) 
plot(cos(a),sin(a), type="l", lty=2, xlab = "X", ylab = 'Y') 
points(cor[cor$X<0 & cor$Y<0,-3], pch=20, col='blue') 
points(cor[cor$X<0 & cor$Y>0,-3], pch=20, col='forestgreen') 
points(cor[cor$X>0 & cor$Y<0,-3], pch=20, col='red') 
abline(v = 0, h = 0) 
text(cor, rownames(cor), pos = 3, cex = 0.8) 

距離は< 30で、X軸とY軸の値に基づいてこれらの点の1つを取得します。

例:点iとgの距離値は、軸点の基準に基づいて0.206であり、X軸とY軸の値が小さいため除外する必要があります。

しかし、どこから始めるべきなのかわかりません。

誰かがコードを開始するのを助けることができますか?

ありがとう

答えて

1

この機能はうまくいくはずですが、おそらくこれを行うより良い方法があります。

nearby <- function(data, d){ 
    dist <- as.matrix(dist(data)) 
    dist[upper.tri(dist, diag = TRUE)] <- NA 
    pairs <- which(dist < d ,arr.ind = TRUE) 
    for (i in 1:nrow(pairs)){ 
     for (j in 1:2){ 
      pairs[i,j] <- letters[as.numeric(pairs[i,j])] 
     } 
    } 
    rownames(pairs) <- NULL 
    colnames(pairs) <- NULL 
    pairs[,2:1] 
} 

ので機能は2点のだけの座標(平面上の点)とのポイントのために働くことを

> nearby(data = cor, d = 0.3) 
    [,1] [,2] 
[1,] "c" "g" 
[2,] "c" "i" 
[3,] "d" "e" 
[4,] "d" "f" 
[5,] "e" "f" 
[6,] "g" "i" 
[7,] "h" "k" 

注意を行い、それらの間の距離が0.3未満であるためにペアを取得します。

+0

ループを 'cbind(rownames(d)[pairs [、1]]、colnames(d)[pairs [、2]]' – rawr

関連する問題