關(guān)于GBD數(shù)據(jù)庫的實(shí)操視頻已經(jīng)陸續(xù)更新，大家可以關(guān)注B站的“小明學(xué)習(xí)室”。

如需要加入GBD數(shù)據(jù)交流群，請(qǐng)加小編微信“Endoscopy_1991”,小編拉你進(jìn)群

接著上一期我們重現(xiàn)了文中的表格1后，我們繼續(xù)復(fù)現(xiàn)文中的圖1。

復(fù)現(xiàn)這個(gè)數(shù)據(jù)，需要我們下載每個(gè)國家的數(shù)據(jù)。同樣的我們登錄GBD網(wǎng)站，下載好國家層面的數(shù)據(jù)，按照前面限制的條件，將location限定在select only countries and territories上。

按照上節(jié)課的內(nèi)容，我們已經(jīng)得到地區(qū)層面的Table 1，現(xiàn)在我們按照同樣的方法需要得到國家層面的數(shù)據(jù)，再將數(shù)據(jù)映射到地圖上。

我們首先復(fù)現(xiàn)圖1中的A圖，也就是ASR標(biāo)準(zhǔn)發(fā)病率的地圖展示，我們讀取需要的數(shù)據(jù)

我們首先設(shè)置工作路徑，安裝以下幾個(gè)包后，讀取需要的包

我們讀取2019年的ASR數(shù)據(jù)

我們讀取地圖數(shù)據(jù)，

查看下worldData數(shù)據(jù)集，可以看到各個(gè)經(jīng)緯度對(duì)應(yīng)的國家名稱（region）

由于region里的部分國家名稱可能會(huì)和ASR_2019國家數(shù)據(jù)集的國家名稱（location）不一致，因此我們需要對(duì)不一致的內(nèi)容進(jìn)行調(diào)整。以下代碼是我核對(duì)兩者數(shù)據(jù)集后不一樣的國家名稱后對(duì)ASR_2019國家數(shù)據(jù)集的國家名稱（location）進(jìn)行調(diào)整后的數(shù)據(jù)結(jié)果，這里我已經(jīng)調(diào)整完畢，大家可以直接取用。

接著我們來畫圖，采用以下代碼，大家調(diào)整自己數(shù)據(jù)的時(shí)候需要對(duì)break進(jìn)行調(diào)整，不然圖出來會(huì)比較難看。

運(yùn)行后就可以得到以下圖1A，然后大家導(dǎo)出為矢量圖格式，然后用AI調(diào)整后圖片后就可以得到文章中的圖1A。

同樣的道理，我們來繪制文章中的圖1B,根據(jù)圖1B的描述，我們需要得到2019年以及1990的發(fā)病數(shù)，從而得到圖B。

首先，我們?cè)诶L圖之前需要得到目標(biāo)的數(shù)值，change in cancer cases.

運(yùn)行上述代碼即可獲得我們的目標(biāo)結(jié)果，接著我們需要匹配好兩個(gè)數(shù)據(jù)集的國家名字要一致

然后調(diào)整好數(shù)據(jù)集后，我們進(jìn)行繪圖

運(yùn)行后就可以得到圖1B的類似結(jié)果，然后大家導(dǎo)出為矢量圖格式，然后用AI調(diào)整后圖片后就可以得到文章中的圖1B。

同樣的思路來繪制EAPC的地圖圖1C

首先獲取EAPC結(jié)果

校正好國家名稱，繪制圖像

最后得到如下結(jié)果

好了，大家接著在AI里調(diào)整3張圖的格式后合并成1張圖后就可以得到想要的圖1了，大家按照上述思路可以實(shí)操起來了

大家如果需要上述代碼，關(guān)注公眾號(hào)“小明學(xué)習(xí)室”，回復(fù)關(guān)鍵詞“GBD”即可下載全部代碼。