最早是在 小丫畫圖 公眾號發(fā)現(xiàn)的這個推文 這套代碼，完全復(fù)現(xiàn)這篇Nature Communications，當(dāng)時點了收藏，然后就讓他去吃灰色了。最近又給想起來了。大體看了一下論文和原始代碼，決定花時間重復(fù)一下。論文里一共有11幅圖，4幅圖放在正文里，7幅圖放在了補充材料里。數(shù)據(jù)和代碼全部公開在了github，對應(yīng)的鏈接是https://github.com/karkman/crassphage_project 。接下來的時間希望每天能夠重復(fù)其中的一幅圖。今天的推文算是立個小旗。后面如果沒有做到歡迎大家催更哈。
下面是來自論文中的圖
image.png
圖一是用基礎(chǔ)繪圖函數(shù)畫的，重復(fù)完原始代碼爭取使用ggplot2來畫一下。
image.png
圖二散點圖然后分組添加擬合曲線
image.png
圖三和圖四和圖二差不多，但是更改了圖例的位置
image.png
補充材料圖一 散點圖，分組添加擬合曲線，然后組圖，共享圖例。
image.png
補充材料圖二
image.png
補充材料圖三
image.png
補充材料圖四和圖五 散點圖和箱線圖
image.png
補充材料圖六合圖七
