1

2 つの fasta シーケンスの dotplot (seqinr パッケージの dotPlot() を使用) があり、プロットからいくつかの値 (x、y) を抽出する必要があります。

「通常の」プロットで機能し、プロットからクリックされた値を含むベクトルを返す identify() を使用しようとしましたが、dotPlot では機能しません。

シェルまたはPythonでドットを使用するなど、可能な解決策は大歓迎です。ありがとう

4

2 に答える 2

2

dotPlotを使用imageしてポイントを描画します。identifyしたがって、いつ呼び出されるかを認識する「ポイント」はありません。2 つのシーケンスをプロットし、同じロジックを使用すると、dotPlotより多くの成功が得られます。

seq2 <- seq1 <- paste(sample(c("A","G","T","C"), 100, repl=TRUE), collapse="")
for ( i in sample(1:100, 20) ) { seq2 <- 
  paste0(substr(seq2, 1, i-1), "Z", substr(seq2, i+1, nchar(seq2) ) )}
plot(NA, xlim=c(0,100), ylim=c(0,100)); 
for ( i in 1:100 ) { 
           if(substr(seq1,i,i)==substr(seq2, i,i) ){ points(i,i)}
pts <- identify(1:100, 1:100, n=5)
于 2013-10-25T05:33:31.910 に答える