lunix的昵称 发表于 2015-11-10 13:12 
很聪明的方法,谢谢!
不过在I2=findc(d,"ACJ");这个命令中ACJ是你根据我提供的三行数据总结出来的三种情 ...
简单一点:
data a;
input;
chb = scan(scan(tranwrd(_infile_,"CHB=","|"),2,"|"),1,";");
datalines4;
ASN=2,3;CEU=2;CHB=4,5,6;CHS=6,7
ASN=3;CEU=2,3;AFR=2,3;CHB=4,5;CHS=7,8
ASN=2;CEU=2,5;JPN=4,5,6;AFR=4,5,6;CHB=5,6,7
;;;;
run;