旧问题新问？

huang6894

一个旧问题：
原文：
=======================================================================
好吧。。。我不知道怎么描述问题，大神们见谅。。。

两个文本：
文本1：

1. 
   
2. CYP21A2 =>      *29
   
3. [chr6::32008338::32008339::C/T] =>      *99
   
4. [chr6::32008695::32008696::G/A]&&[chr6::32006860::32006861::C/G]        =>      *56 
   
5. [chr6::32008338::32008339::C/T]&&[chr6::32006860::32006861::C/G]&&[chr6::320088::320089::C/T]        =>      *86
   
6. [chr6::32008338::32008339::C/T]&&[chr6::320088::320089::C/T]  =>  *96
   
7. [chr6::320088::320089::C/T]     =>      *66
   
8.  

复制代码

文本2：

1. 
   
2. CYP21A2   NM_000016       c.134C>T        chr6    32008338        32008339        1p31.1  1       .      C       snv    C/T   .   HET  
   
3. CYP21A2  NM_000016       c.145C>T        chr6    32008695        32008696        1p31.1  1       .       C       snv     G/A   .   HOM
   
4. CYP21A2  NM_000016       c.145C>G        chr1    32006860        32006861        1p31.1  1       .       C       snv     C/G   .   HET
   
5. CYP21A2   NM_000016       c.14C>T        chr6    320088        320089        1p31.1  1       .      C       snv    C/T   .   HOM
   
6. CYP21A2   NM_000016       c.1234C>T        chr16    32011088        32011089        1p31.1  1       .      C       snv    C/T   .   HET

复制代码

怎么做能够当
=====================================================================
引自怎么严格限定条件后得到对应的结果?
http://bbs.chinaunix.net/forum.p ... mp;fromuid=29097174
感谢@yestreenstars@rulebook两位大神的回答，可是我回过头研究的时候，想问问
--假如我现在需要引入一个新的变量：文本二的最后一列（由两种情况HET、HOM）作为条件需求匹配，应该怎么做呢？
新条件如下：
==========================================================

1. 
   
2.   if[文本2中4、5、6、12列满足于文本1 => 前面的所有条件（条件a)]&&[文本2的最后一列是HOM];then
   
3.          输出文本1 => 后面的字符串;
   
4.   elsif[文本2中1列满足于文本1 => 前面的条件（条件b)]&&[文本2的最后一列是HOM];then
   
5.           输出文本1 => 后面的字符串;
   
6.   elsif[文本2的最后一列是HET];then
   
7.         if[该行满足条件a]&&[条件a只有一个];then
   
8.            输出一个条件b的文本1 => 后面的字符串外加一个满足条件a的文本1 => 后面的字符串  #满足条件a必然满足条件b
   
9.         elsif[该行满足条件a]&&[条件a不止一个];then
   
10.             输出同时满足多个条件a时的文本1 => 后面的字符串
    
11.          fi
    
12.    elsif[该行满足条件a]&&[条件a不止一个]&&[文本2的最后一列既有HOM也有HET];then
    
13.          首先输出最后一列为HOM的文本1 => 后面的字符串 ，HET的再组合考虑
    
14.    else
    
15.          输出".";
    
16.     fi
    
17. 
    

复制代码

举例：
----------------------------------------------------------------------------------------------------------------------
当文本2中有一行第4列是chr6、第5列是32008338、第六列是32008339、第12列是C/T，第一列是CYP21A2，且最后一列是HET时，输出*29*99
当文本2中有一行第4列是chr6、第5列是32008695、第六列是32008696、第12列是G/A的同时还存在一行第4列是chr6、第5列是32006860、第六列是32006861、第12列是C/G，两行最后一列都是HET时，输出*56
当文本2中有一行第4列是chr6、第5列是32008338、第六列是32008339、第12列是C/T的同时还存在一行第4列是chr6、第5列是32006860、第六列是32006861、第12列是C/G的同时还存在一行第4列是chr6、第5列是320088、第六列是320089、第12列是C/T，最后一列都是HET时，输出*86
--------------------------------------------------------------------------------------------------------------------
当文本2中有一行第4列是chr6、第5列是32008338、第六列是32008339、第12列是C/T，且最后一列是HOM时，输出*99
当文本2中有一行第4列是chr6、第5列是32008695、第六列是32008696、第12列是G/A的（第一行）同时还存在一行第4列是chr6、第5列是32006860、第六列是32006861、第12列是C/G（第二行），即使两行最后一列都是HOM。输出*56，因为在文本1中没有第一行、第二行能单独决定的字符串；
当文本2中有一行第4列是chr6、第5列是32008338、第六列是32008339、第12列是C/T（第一行）的同时还存在一行第4列是chr6、第5列是32006860、第六列是32006861、第12列是C/G（第二行）的同时还存在一行第4列是chr6、第5列是320088、第六列是320089、第12列（第三行）是C/T，最后一列都是HOM时，第一行输出*99、第二行输出"."（因为切掉第一第三行的记录，单独的第二行或者第二第三行组合不决定字符串）、第三行输出 *66
----------------------------------------------------------------------------------------------------------------------
当文本2中有一行第4列是chr6、第5列是32008695、第六列是32008696、第12列是G/A的（第一行）同时还存在一行第4列是chr6、第5列是32006860、第六列是32006861、第12列是C/G（第二行），第一行是HOM、第二行是HET。输出*56，因为在文本1中没有第一行、第二行能单独决定的字符串；
当文本2中有一行第4列是chr6、第5列是32008338、第六列是32008339、第12列是C/T（第一行）的同时还存在一行第4列是chr6、第5列是32006860、第六列是32006861、第12列是C/G（第二行）的同时还存在一行第4列是chr6、第5列是320088、第六列是320089、第12列（第三行）是C/T，第一行是HET、第二行是HOM、第三行是HET时，第一行输出*96、第二行输出"."（因为切掉第一第三行的记录，单独的第二行或者第二第三行组合不决定字符串）、第三行输出 *96
-----------------------------------------------------------------------------------------------------------------------
当文本2中有一行第4列是chr6、第5列是32008338、第六列是32008339、第12列是C/T（第一行）的同时还存在一行第4列是chr6、第5列是32006860、第六列是32006861、第12列是C/G（第二行）的同时还存在一行第4列是chr6、第5列是320088、第六列是320089、第12列（第三行）是C/T，第一行是HOM、第二行是HOM、第三行是HET时，第一行输出*99、第二行输出"."（因为切掉第一第三行的记录，单独的第二行或者第二第三行组合不决定字符串）、第三行输出 *29*26
&&表示前后都要满足，文本1中可能同时存在多个&&

也就是最后得到：

1. 
   
2. CYP21A2   NM_000016       c.134C>T        chr6    32008338        32008339        1p31.1  1       .      C       snv    C/T   .   HET  *29*99
   
3. CYP21A2  NM_000016       c.145C>T        chr6    32008695        32008696        1p31.1  1       .       C       snv     G/A   .   HET   *56 
   
4. CYP21A2  NM_000016       c.145C>G        chr1    32006860        32006861        1p31.1  1       .       C       snv     C/G   .   HET     *56 
   
5. CYP21A2   NM_000016       c.14C>T        chr6    320088        320089        1p31.1  1       .      C       snv    C/T   .   HOM   *66
   
6. CYP21A2   NM_000016       c.1234C>T        chr16    32011088        32011089        1p31.1  1       .      C       snv    C/T   .   HET   .
   

复制代码

sorry~跪求大神们原谅我的描述能力~

Herowinter

个人理解能力有限，呼唤熟悉问题的两位大神
@yestreenstars@rulebook。

huang6894

回复 2# Herowinter
Herowinter大神啊~为什么我看了很多遍论坛上面说的提问艺术我还不会提问呢？

   

Herowinter

回复 3# huang6894

和你的问题描述无关，是我想偷懒，问题描述实在太多。
更何况以前有人帮你写过代码了，让他们改一下更快吧。
   

huang6894

回复 4# Herowinter


    谢谢大神给我留了面子。。。哈哈，我也觉得描述好多。。。。

yestreenstars

信息量太大~

Herowinter

回复 1# huang6894

看了好久，文本1 文本2得到这个结果的算法没看懂

1. CYP21A2   NM_000016       c.134C>T        chr6    32008338        32008339        1p31.1  1       .      C       snv    C/T   .   HET  *29*99
   
2. CYP21A2  NM_000016       c.145C>T        chr6    32008695        32008696        1p31.1  1       .       C       snv     G/A   .   HET   *56
   
3. CYP21A2  NM_000016       c.145C>G        chr1    32006860        32006861        1p31.1  1       .       C       snv     C/G   .   HET     *56
   
4. CYP21A2   NM_000016       c.14C>T        chr6    320088        320089        1p31.1  1       .      C       snv    C/T   .   HOM   *66
   
5. CYP21A2   NM_000016       c.1234C>T        chr16    32011088        32011089        1p31.1  1       .      C       snv    C/T   .   HET   

复制代码

从这结果看，好象是文本2的每一行加点查文本1得到的字符串，但你的算法描述里还要结合其他行，
2-3行综合考虑，实在看不懂。
   

expert1

这描述，看的蛋疼

关阴月飞

进来恶补一下.....

huang6894

回复 6# yestreenstars


    是呀，，哈哈~大大如果有时间帮忙看看？