TBtools進(jìn)行GO富集分析

1、下載GOA文件,下載網(wǎng)址如下:
http://current.geneontology.org/annotations/index.html
2、下載人的文件,如圖:

image.png

3、然后利用腳本選擇你需要的列,腳本如下:

def main():
    """
    main function
    """
    f = open('goa_human.gaf', 'r')
    result = open('goa_human_result.txt', 'w')

    f_lines = f.readlines()
    for i in f_lines:
        i = i.strip()

        if 'ensembl:ENSMUSP' in i:
            i_list = i.split('\t')
            i_list_ensmbl = i_list[4:9]
            str_ensmbl = '\t'.join(i_list_ensmbl)
            result.write(str_ensmbl + '\n')


    f.close()
    result.close()

if __name__ == '__main__':
    main()

4、生成背景文件Human_GO_background.txt。將需要進(jìn)行GO富集分析的基因整理成列表,下載go-basic.obo文件,直接點(diǎn)擊保存即可,如圖:


image.png

5、輸入的背景數(shù)據(jù)格式Human_GO_background.txt,如圖:


image.png

6、輸入的基因數(shù)據(jù)格式ITGB1_GO.txt,如圖:
image.png
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

友情鏈接更多精彩內(nèi)容