Re: [問題] 撈資料

[yaowei2010 (yaowei)]

不好意思 讓回復的大大看不懂我的問題 只好重新回一篇


 gene_a  Chromosome  gene_start gene_end |  Probe_b Chromosome2 Chr_s  Chr_e
 A           1          25000    50000   |  a1         2        175   200
 B           X           1000     2000   |  a2         4        600   625
 C           2           0         800   |  a3         1       23575  23600
                                         |  a4         X        1010   1035
                                         |  a5         1       30000  30025



我想做的是

先對Chromosome是不是在一個→再對位子→假設兩者都相符，收data(不符丟掉

例如上面兩data

probe a3和a5都為chromosome 1 兩位子分別23575~23600及30000~30025

A gene為chromosome 1 ，其位子為25000~50000，上述a5 probe都符合其chromosome

以及gene範圍，歸類到gene A裡的match中


我的想法是兩關，第一關篩選出和gene一樣的chromosome的probe

第二關比較probe的數值要落入gene的位子範圍中(gene_start~gene_end)

大概是這樣> <


以此類推


最後才變成

 gene_a   match_probe
  A        a5
  B        a4
  C        a1

--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.112.129.5
※ 文章網址: http://www.ptt.cc/bbs/R_Language/M.1408530493.A.512.html
→ yaowei2010: 補充 Probe和gene的數值全都是numeric
→ yaowei2010: 位子範圍就是指start~end
→ koai: 請確認 A gene 有match到 a3 嗎??
→ yaowei2010: 阿 沒有@@" 眼殘惹
→ yaowei2010: 那大大得應該是對的 謝謝指正與指教!!
※ 編輯: yaowei2010 (140.112.129.5), 08/20/2014 18:35:22
→ yaowei2010: 文已修正