论文

Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic

image.png

本地文件 s41564-020-0771-4.pdf

代码和数据下载链接

https://github.com/plemey/SARSCoV2origins

今天的推文我们来重复一下论文中的 Figure5 a

image.png

加载R包

library(pheatmap)
library(tidyverse)
library(viridis)

读取数据

all_rscu <- read_csv('all_rscu_codonBiasReanalysis.csv')

选取数据集的2到26列

all_df <- as.data.frame(all_rscu[,2:26])

给数据集赋予行名

row.names(all_df) <- all_rscu$codon

画热图

pheatmap(all_df,
         cluster_rows = F,
         scale = 'none',
         fontsize_row = 5,