将图表 (ggplot2) 中的标题与另一个文件的标签作为下标

Gre*_*Rov 3 linux command loops r ggplot2

在我的程序中,我有两个主要文件,第一个包含数据,第二个包含标签(或图形标题):

\n\n

文件总数1(数据)

\n\n
3   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\n
Run 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* \n
Run Code Online (Sandbox Code Playgroud)\n\n

我在同一张 bash 卡上同时使用 linux 命令和 R 来绘制图表。我正在这个脚本中循环:

\n\n
\n

!/bin/bash

\n\n

对于 {1..9} 中的 l;做

\n\n

R -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
\n\n

程序没问题,我只需要在我的图形标题中以下标形式写入“pH24 Final”的24和“pH0 Inicial”的0,以从我使用命令的标签文件中写入标题ggtitle(leg[$l,1])。有谁知道我如何通过使用另一个文件的标签来做到这一点?

\n

Bin*_*ing 7

看看这是否可以给你一些想法。在基本 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)

在此输入图像描述