根据B站小明的数据分析笔记本的rscu堆积图代码改的。侵删
以下为原视频链接
R语言ggplot2画堆积柱形图展示密码子偏向性的RSCU值_哔哩哔哩_bilibili
以下为改动之后的图片效果。颜色和密码子末尾碱基对应。这样每个氨基酸可以看到偏好哪个碱基结尾的密码子。配色至今没搞明白。就找了色轮来看,数还是对不上。
目前的,A=15,T=135,C=195,G=315
以下为代码部分,在原基础上略有改动。x5是y轴高度。x6是数字颜色。
library(xlsx)
ray = read.xlsx('/Users/shmilyy/Desktop/1234567.xlsx',1,header=F,stringsAsFactors = F)
library(ggplot2)
library(aplot)
上半部分的p1
p1<-ggplot(ray,aes(x=X3,y=X4,fill=as.character(X6)))+geom_bar(stat="identity",position = "stack")+
theme_bw()+theme(legend.position = "")+
scale_y_continuous(expand = c(0,0),limits = c(0,6.2))+
labs(x=NULL,y="RSCU")
下半部分的P2
p2<-ggplot(ray,aes(x=X3,y=X5,fill=as.character(X6)))+
geom_label(aes(label=X2),size=2)+
theme_minimal()+
ylim(1.8,6)+
theme(legend.position = "none",
axis.title = element_blank(),
axis.text = element_blank(),
panel.grid = element_blank())
insert_bottom(p1,p2,height = 0.8)
写在末尾。 装rstudio是2022.9.09装的。图是今天转完愈伤板子改的。卡壳的点也有很多,莫名其妙报错最多。昨天java崩了。重新装的java。r还是很友好的, 出问题告诉你哪里出的,改起来目的性还是比较强的。比试验友好一些。试验出问题从来都是试试哪里出了问题。。。
另外x11模块报错的问题,可以用最上面两行代码。读取文件数据来就好了。
纯小白一天把图粘出来了。我觉得我已经很好了。哈哈哈哈哈。