Re: [问题] 捞资料

楼主: yaowei2010 (yaowei)   2014-08-20 18:28:10
不好意思 让回复的大大看不懂我的问题 只好重新回一篇
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
楼主: yaowei2010 (yaowei)   2014-08-20 18:31:00
补充 Probe和gene的数值全都是numeric位子范围就是指start~end
作者: koai (可笑的自以为...)   2014-08-20 18:31:00
请确认 A gene 有match到 a3 吗??
楼主: yaowei2010 (yaowei)   2014-08-20 18:33:00
阿 没有@@" 眼残惹那大大得应该是对的 谢谢指正与指教!!文已修正

Links booklink

Contact Us: admin [ a t ] ucptt.com