2013-05-11 43 views
1

我有大量的X和Y變量以及它們的集羣成員列。羣集成員資格(羣組)可能並不總是正確的(限制完善羣集算法),我想交互式地顯示羣集並操縱羣集成員身份以識別點。與xy點圖集羣交互工作 - r中的組操作

我試過rggobi,以下是我能夠得到的觀點(我不是說我需要使用rggobi/ggobi,如果有更好的選項可用,歡迎您提出建議)。

# data 
set.seed (1234) 
c1 <- rnorm (40, 0.1, 0.02); c2 <- rnorm (40, 0.3, 0.01) 
c3 <- rnorm (40, 0.5, 0.01); c4 <- rnorm (40, 0.7, 0.01) 
c5 <- rnorm (40, 0.9, 0.03) 
Yv <- 0.3 + rnorm (200, 0.05, 0.05) 
myd <- data.frame (Xv = round (c(c1, c2, c3, c4, c5), 2), Yv = round (Yv, 2), 
cltr = factor (rep(1:5, each = 40))) 

require(rggobi) 
g <- ggobi(myd) 
display(g[1], vars=list(X="Xv", Y="Yv")) 

enter image description here

可以看到五組,與CLTR可變顏色不同。我手動確定了異常點,並且我想在cltr變量中將它們的值設爲NA。他們有沒有簡單的方法來解除這種會員資格並寫入文件。

回答

1

你可以嘗試identify拿到手動之分指數:

## use base::plot 
plot(myd$Xv, myd$Yv, col=myd$cltr) 

exclude <- identify(myd$Xv, myd$Yv) ## left click on the points you want to exclude (right click to stop/finish) 

myd$cltr[exclude] <- NA 
+0

涼,謝謝。可能是我只是將cltr值替換爲NA(不從繪圖中移除Xv,Yv,但是myd $ cltr [exclude] < - NA不工作,有什麼想法? – jon 2013-05-11 15:50:14

+0

請不要說「...不工作」 ,而是提供錯誤信息。 – 2013-05-11 18:06:37

+0

現在適合我,謝謝 – jon 2013-05-12 14:02:16