u011919863的博客

1978年，Margaret Belle (Oakley) Dayhoff教授给出了一个PAM矩阵，这个矩阵实际上是一个概率表，就是每一个氨基酸突变成另一个氨基酸的概率是多少，由于没有足够多的序列和理想的比对算法，她选择了相似程度非常高（≥85%）的蛋白质序列集。在统计了1572个突变点之后，她总结出了每一个氨基酸突变成另一个氨基酸的概率表，就是PAM1矩阵。      

data<-read.csv("cdr15.csv",header=T)library(ggplot2)library(reshape2)data_m=melt(data)# variable和value为矩阵melt后的两列的名字，内部变量, variable代表了点线的属性，value代表对应的值。#箱线图#p1<-ggplot(data_m,aes(x=variable,y=va