Gre*_*Rov 3 linux command loops r ggplot2
在我的程序中,我有两个主要文件,第一个包含数据,第二个包含标签(或图形标题):
\n\n文件总数1(数据)
\n\n3 10000 3 32039232 1 0.0017290351 2 0.0002781092\n3 10001 3 32101193 1 0.0045398899 2 0.0032875689\n3 1000 1 60233253 1 0.0022057964 2 6.747e-06\n3 10002 3 32108182 1 0.0219913914 2 0.0102120679\n3 10003 3 32133994 1 0.0007025013 2 0.0010197563\n3 10004 3 32192498 1 0.0029210855 2 0.0036980008\n3 10005 3 32230041 1 0.0005408603 2 0.0015782048\n3 10006 3 32271305 1 1.099e-07 2 0.0033466856\n3 10007 3 32289336 1 0.0185812303 2 0.0027349589\n3 10008 3 32453784 1 0.0080117379 2 0.0003596759\nRun Code Online (Sandbox Code Playgroud)\n\n锉腿(标签)
\n\n\xc3\x81re de olho de lombo\nEspessura de gordura subcut\xc3\xa2nea\npH0 inicial \npH24 final\nPerda por coc\xc3\xa7\xc3\xa3o\nFor\xc3\xa7a de cisalhamento\nCor L*\nCor a*\nCor b* \nRun Code Online (Sandbox Code Playgroud)\n\n我在同一张 bash 卡上同时使用 linux 命令和 R 来绘制图表。我正在这个脚本中循环:
\n\n\n\n\n!/bin/bash
\n\n对于 {1..9} 中的 l;做
\n\nR -q -e "leg<-read.table(\'leg\',\n header=F,sep=\'\\t\');write.table(leg,\'lef1.txt\', sep=\'\\t\', row.names =\n F,col.names=F);uni$l<-read.table(\'./var/chrsnpvar_uni$l\',header=F) ;bi$l<-read.table(\'./var/chrsnpvar_bi$l\',\n header=F);map<-read.table(\'snp_map_clean\',sep=\'\\t\ ',header =\n F);uni1$l=uni$l[ 其中(uni$l[,1]==1 & uni$l[,2]==3), ]; bi1$l=bi$ l [\n 其中(bi$l[,1]==2 & bi$l[,2]==3), ];total$l <- 合并(uni1$l,bi1$l,\n by= c(2,4,5,6), all=T);write.table(total$l,\'total$l.txt\', sep=\'\\t\',\n row.names = \n F,col.names=F);library(ggplot2);tiff(\'./solution/cor$l.tiff\', 宽度\n =10 , 高度 = 6, 单位 = \'in\', res = 75 ); ggplot(总计$l,aes(x=总计$l[,6], y=总计$l[,8])) + geom_point() +\n scale_x_continuous(breaks = round(seq(min) (0), max(0.08), by =\n 0.01),2),limits=c(0,0.08))+ scale_y_continuous(breaks = round(seq(min(0), max(0.08), by = 0.01) ),2),limits=c(0,0.08)) +\n geom_smooth(method=lm , color=\'grey35\', se=FALSE)+\n geom_vline(xintercept =\n 分位数(总计$l[ ,6],0.95),lty=2,size=1,颜色=\'灰色\')+ \n geom_hline(yintercept = 分位数(总计$l[,8], .95),\n lty=2,大小=1,颜色=\'灰色\') + ylab(\'Bicaacter\xc3\xadstica\') +\n xlab(\'Unicaracter\xc3\xadstica\')+ theme_classic()+\n theme(axis .line.x=element_line(颜色=\'黑色\',size=1),axis.line.y=element_line(颜色=\'黑色\',size=1),\n axis.ticks.length=单位(4,\'mm\'),legend.title =\n element_blank(),axis.title.y = element_text(size = 22), axis.title.x =\n element_text(size = 22),axis. text.x=element_text(size =\n 22),axis.text.y=element_text(size = 22), legend.text=element_text(size\n = 22), title=element_text(size =22)) + 注释(\'文本\', x = 0.005, y=0.07, 标签 = \'A\',size=8) + 注释(\'文本\', x = 0.06, y=0.07, 标签\n = \' B\',size=8)+ 注释(\'文本\', x = 0.005, y=0.005, 标签=\'C\', 颜色=\'白色\', size=8)+ 注释(\'文本\', x = 0.06, y=0.005, 标签=\n \'D\',size=8)+ggtitle(leg[$l,1]); dev.off()"
\n\n完毕
\n
程序没问题,我只需要在我的图形标题中以下标形式写入“pH24 Final”的24和“pH0 Inicial”的0,以从我使用命令的标签文件中写入标题ggtitle(leg[$l,1])。有谁知道我如何通过使用另一个文件的标签来做到这一点?
看看这是否可以给你一些想法。在基本 R 中使用数学表达式似乎适用于 ggplot 标题。
library(ggplot2)
df=data.frame(x=1:10, y=rnorm(10))
vn=expression(paste(pH[10]," inicial"))
ggplot(df, aes(x,y))+geom_point()+ ggtitle(vn)
Run Code Online (Sandbox Code Playgroud)
| 归档时间: |
|
| 查看次数: |
5728 次 |
| 最近记录: |