步骤一:打开链接并选择物种
http://www.genome.jp/kegg-bin/get_htext?hsa00001+3101
步骤二:对文件进行解析
步骤三:统计信息
一级结构(6大类):
grep ^A dre00001.keg |les
二级结构(下面三种方法都行,54小类):
grep ^B dre00001.keg |awk '$0 != "B" {print}' |les grep ^B dre00001.keg |grep -v 'B$' |les grep ^B dre00001.keg |sort -u |sed '1d' |les
三级结构(515 信号通路):
grep ^C dre00001.keg |les
四级结构(27101个基因):
grep ^D dre00001.keg |wc -l