亚洲乱码中文字幕综合,中国熟女仑乱hd,亚洲精品乱拍国产一区二区三区,一本大道卡一卡二卡三乱码全集资源,又粗又黄又硬又爽的免费视频

Python 利用Entrez庫篩選下載PubMed文獻(xiàn)摘要的示例

 更新時(shí)間:2020年11月24日 16:24:11   作者:v林三歲  
這篇文章主要介紹了Python 利用Entrez庫篩選下載PubMed文獻(xiàn)摘要的示例,幫助大家更好的理解和使用python,感興趣的朋友可以了解下

作者:xiaolanLin

聲明 :本文版權(quán)歸作者和博客園共有,來源網(wǎng)址:https://www.cnblogs.com/xiaolan-Lin

一個(gè)不是學(xué)生物的孩子來搞生物,當(dāng)真是變成了一塊廢鐵啊,但也是讓我體會(huì)到了一把生物信息的力量。

廢話不多說,開整!

任務(wù):快速高效從PubMed上下載滿足條件的文獻(xiàn)PMID、標(biāo)題(TI)、摘要(AB)。

PubMed官網(wǎng) https://pubmed.ncbi.nlm.nih.gov

此處有幾種選擇可以達(dá)到目的:

(1)官網(wǎng)上匹配篩選條件(注:匹配快速,但是下載下來的數(shù)量受到限制,每次只能下載10000條數(shù)據(jù),甚至更少。)

 可以看到,我需要的數(shù)據(jù)是有三十多萬條,但是每次只能下載10000條,那我豈不是要手動(dòng)n次。。很明顯,在大批量下載文獻(xiàn)的情況下,官網(wǎng)不是很友好。

(2)R語言有個(gè)R包,叫做easyPubMed,這里我也給大家貼上學(xué)習(xí)指南(https://cran.r-project.org/web/packages/easyPubMed/vignettes/getting_started_with_easyPubMed.html

由于我不喜歡用R寫代碼,所以我寫一半還是換了Python,熟練R的小伙伴可以自行根據(jù)指南走通需求。

(3)重量級(jí)庫來了,Python自帶的Bio包中的Entrez檢索庫,簡(jiǎn)直就是我的救星,以下是我的代碼:

  注:Entrez在Bio包中,Bio的安裝請(qǐng)移步 https://www.cnblogs.com/xiaolan-Lin/p/14023147.html

 import numpy as np
 from Bio import Medline, Entrez # 一般是通過BioPython的Bio.Entrez模塊訪問Entrez
 from collections import Counter
 
 Entrez.email = "(此處寫你自己在官網(wǎng)注冊(cè)的郵箱賬號(hào))" # 應(yīng)用自己的賬號(hào)訪問NCBI數(shù)據(jù)庫
 
 # 此處需將服務(wù)器協(xié)議指定為1.0,否則會(huì)出現(xiàn)報(bào)錯(cuò)。http.client.IncompleteRead: IncompleteRead(0 bytes read)
 # 服務(wù)器http協(xié)議1.0,而python的是1.1,解決辦法就是指定客戶端http協(xié)議版本
 import http.client
 http.client.HTTPConnection._http_vsn = 10
 http.client.HTTPConnection._http_vsn_str = 'HTTP/1.0'
 
 """
 Entrez 是一個(gè)檢索系統(tǒng),可以用其訪問NCBI數(shù)據(jù)庫,比如說PubMed,GenBank,GEO等。
 獲得有關(guān) global PBDE 的所有文獻(xiàn)的PubMed IDs
 """
 # handle_0 = Entrez.esearch(db="pubmed", term="drug therapy[Subheading] AND adverse effects[Subheading] AND humans[MeSH Terms]", retmax=306431)
 handle_0 = Entrez.esearch(db="pubmed", term="drug therapy[MeSH Subheading] AND adverse effects[MeSH Subheading] AND humans[MeSH Terms] AND (2000/01/01[Date - Publication] : 2021/12/31[Date - Publication])",
              ptyp="Review", usehistory="y", retmax=306431)
 record = Entrez.read(handle_0) # 獲取檢索條件的所有文獻(xiàn)
 idlist = record["IdList"] # 提取出文獻(xiàn)id
 print ("Total: ", record["Count"])
 No_Papers = len(idlist)  # 共306431篇文獻(xiàn) 2000-01-01:2021-12-31
 webenv = record['WebEnv']
 query_key = record['QueryKey']
 
 total = No_Papers
 step = 1300
 print("Result items:", total)
 with open("./Data_PubMed/PBDE1.txt", 'w') as f:
   for start in range(0, total, step):
     print("Download record %i to %i" % (start + 1, int(start + step)))
     handle_1 = Entrez.efetch(db="pubmed", retstart=start, rettype="medline", retmode="text",
                 retmax=step, webenv=webenv, query_key=query_key) # 獲取上述所有文獻(xiàn)的PubMed IDs
     records = Medline.parse(handle_1)
     records = list(records) # 將迭代器轉(zhuǎn)換至列表(list)
     for index in np.arange(len(records)):
       id = records[index].get("PMID", "?")
       title = records[index].get("TI", "?")
       title = title.replace('[', '').replace('].', '') # 若提取的標(biāo)題出現(xiàn)[].符號(hào),則去除
       abstract = records[index].get("AB", "?")
       f.write(id)
       f.write("\n")
       f.write(title)
       f.write("\n")
       f.write(abstract)
       f.write("\n")

話不多說,結(jié)果跑出來了我真的很快樂~

最后的結(jié)果是存放在txt文件中,大伙兒根據(jù)自己的需求改變代碼所需字段啊。

現(xiàn)在我來解釋一下,我貼上的這串代碼的實(shí)現(xiàn)原理,首先是通過Entrez檢索到符合我篩選條件的文獻(xiàn),里邊的限制條件包括了幾個(gè)詞匯匹配以及時(shí)間限制,時(shí)間我限制在了2000年1月1日到2021年的12月31日(這里的時(shí)間我選用的是Date - Publication,時(shí)間選取Date - Completion、Date - Modification還是Date - Publication其實(shí)還是有爭(zhēng)議的,大家自行考慮選取)。

Entrez.esearch的作用就是用來檢索的,里邊的參數(shù)db指向你要檢索的數(shù)據(jù)庫,代碼中的注釋也寫了,Entrez作為一個(gè)接口檢索,除了能夠檢索PubMed中的文獻(xiàn),也能去到別的數(shù)據(jù)庫檢索文獻(xiàn);term是寫你的篩選語句,注意你寫的檢索語句不能帶有引號(hào),單引號(hào)也不行,否則會(huì)檢索不到,如果不知道檢索語句怎么寫,或者是不知道字段是否被定義,可以在官網(wǎng)的檢索那里https://pubmed.ncbi.nlm.nih.gov/advanced/選擇字段輸入內(nèi)容自動(dòng)生成query,但是生成的語句是不太智能的,會(huì)有很多括號(hào)是你不需要的,自己寫代碼的時(shí)候要適當(dāng)去掉;ptyp我這里用的是Review,usehistory是y,意思是后邊我的檢索要記住這個(gè)語句,根據(jù)歷史查詢來檢索;retmax如果不進(jìn)行設(shè)置的話,默認(rèn)給你的最大數(shù)據(jù)量好像是只有1000,我要的檢索內(nèi)容是超過這個(gè)值的,因此我需要自定義檢索的數(shù)量。

Entrez.read是對(duì)Entrez.esearch檢索到的內(nèi)容進(jìn)行讀取,里邊包含了9種元素,我們主要是想從中得到文獻(xiàn)的id號(hào),只有拿到了文獻(xiàn)的id號(hào),我們后面進(jìn)行摘要的提取才能準(zhǔn)確定位。

最后是循環(huán)當(dāng)中步長(zhǎng)的設(shè)置,這里就要根據(jù)自己的需求來定義了,包括內(nèi)容的提取,因?yàn)槲抑恍枰狿MID、標(biāo)題(TI)、摘要(AB),所以我就沒有加載別的內(nèi)容進(jìn)來,這里也有一點(diǎn)要注意,標(biāo)題下載下來是大部分帶有[ ].的,方便操作我直接就在下載的時(shí)候給去除了,這也是上面replace代碼的由來。

附上我參考的鏈接,如果我這篇文章解決不了你的問題,那么希望下面的渠道能夠幫助到你

https://zhuanlan.zhihu.com/p/54611852

https://zhuanlan.zhihu.com/p/262957260

以上就是Python 利用Entrez庫篩選下載PubMed文獻(xiàn)摘要的示例的詳細(xì)內(nèi)容,更多關(guān)于Python Entrez庫下載PubMed文獻(xiàn)的資料請(qǐng)關(guān)注腳本之家其它相關(guān)文章!

相關(guān)文章

  • python 七種郵件內(nèi)容發(fā)送方法實(shí)例

    python 七種郵件內(nèi)容發(fā)送方法實(shí)例

    這篇文章主要介紹了python 七種郵件內(nèi)容發(fā)送方法實(shí)例,需要的朋友可以參考下
    2014-04-04
  • python獲取當(dāng)前計(jì)算機(jī)cpu數(shù)量的方法

    python獲取當(dāng)前計(jì)算機(jī)cpu數(shù)量的方法

    這篇文章主要介紹了python獲取當(dāng)前計(jì)算機(jī)cpu數(shù)量的方法,涉及Python操作計(jì)算機(jī)硬件的技巧,代碼簡(jiǎn)單易懂,非常具有實(shí)用價(jià)值,需要的朋友可以參考下
    2015-04-04
  • PyTorch中可視化工具的使用

    PyTorch中可視化工具的使用

    本文主要介紹了PyTorch中可視化工具的使用,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧
    2023-04-04
  • python列表推導(dǎo)式 經(jīng)典代碼

    python列表推導(dǎo)式 經(jīng)典代碼

    這篇文章主要介紹了python列表推導(dǎo)式,文章內(nèi)容附有經(jīng)典的詳細(xì)代碼,感興趣的小伙按可以參考下面文章具體內(nèi)容
    2021-10-10
  • 解決阿里云郵件發(fā)送不能使用25端口問題

    解決阿里云郵件發(fā)送不能使用25端口問題

    這篇文章主要介紹了解決阿里云郵件發(fā)送不能使用25端口問題,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧
    2020-08-08
  • Python3 replace()函數(shù)使用方法

    Python3 replace()函數(shù)使用方法

    這篇文章主要介紹了Python3 replace()函數(shù)使用方法,需要的朋友可以參考下
    2018-03-03
  • python讀取Windows注冊(cè)表的示例代碼

    python讀取Windows注冊(cè)表的示例代碼

    本文主要介紹了python讀取Windows注冊(cè)表的示例代碼,文中通過示例代碼介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友們下面隨著小編來一起學(xué)習(xí)學(xué)習(xí)吧
    2022-08-08
  • 實(shí)例講解Python中浮點(diǎn)型的基本內(nèi)容

    實(shí)例講解Python中浮點(diǎn)型的基本內(nèi)容

    在本文里小編給大家整理了關(guān)于Python中浮點(diǎn)型的基本知識(shí)點(diǎn)內(nèi)容,有興趣的朋友們學(xué)習(xí)下。
    2019-02-02
  • python+selenium實(shí)現(xiàn)自動(dòng)搶票功能實(shí)例代碼

    python+selenium實(shí)現(xiàn)自動(dòng)搶票功能實(shí)例代碼

    Selenium是ThoughtWorks公司的一個(gè)強(qiáng)大的開源Web功能測(cè)試工具系列,采用Javascript來管理整個(gè)測(cè)試過程,包括讀入測(cè)試套件、執(zhí)行測(cè)試和記錄測(cè)試結(jié)果。這篇文章主要介紹了python+selenium實(shí)現(xiàn)自動(dòng)搶票,需要的朋友可以參考下
    2018-11-11
  • 使用django的objects.filter()方法匹配多個(gè)關(guān)鍵字的方法

    使用django的objects.filter()方法匹配多個(gè)關(guān)鍵字的方法

    今天小編就為大家分享一篇使用django的objects.filter()方法匹配多個(gè)關(guān)鍵字的方法,具有很好的參考價(jià)值,希望對(duì)大家有所幫助。一起跟隨小編過來看看吧
    2019-07-07

最新評(píng)論