> ggscatter(data,x="NuCy1_FPKM",y="NuCy2_FPKM",add="reg.line",cor.coeff.args = list(method="pearson",label.x=3,label.sep="\n"))+stat_cor(method = "pearson",label.x = 1500,label.y = 5000)+xlab("NuCy1_FPKM")+ylab("NuCy2_FPKM")> ggscatter(data,x="NuCy1_FPKM",y="NuCy2_FPKM",add="reg.line",cor.coeff.args = list(method="pearson",label.x=3,label.sep="\n"))+stat_cor(method = "pearson",label.x = 1500,label.y = 5000)+xlab("NuCy1_FPKM")+ylab("NuCy2_FPKM")+xlim(0,6000)+ylim(0,6000)
詳細(xì)看收藏文章
相關(guān)系數(shù)：
>a<-read.table("FPKM.txt")
> corr <- cor(a)
> corrplot(corr=corr,addCoef.col = "white")
