dir="/data/zhangyong/yyp/yypold/hic/chenhebing/cs" ##定义工作的文件夹
sedwd=(dir)
options(scipen = 20)
Args <- commandArgs(T)
library(GENOVA)
Ha40 <- load_contacts(signal_path = '/data/zhangyong/yyp/yypold/hic/chenhebing/cs/hic_results/matrix/GSE125129_WT_HiC1/iced/40000/GSE125129_WT_HiC1_40000_iced.matrix',
indices_path = '/data/zhangyong/yyp/yypold/hic/chenhebing/cs/hic_results/matrix/GSE125129_WT_HiC1/raw/40000/GSE125129_WT_HiC1_40000_abs.bed',
sample_name = "WT",
colour = "black")
Hb40 <- load_contacts(signal_path = '/data/zhangyong/yyp/yypold/hic/chenhebing/cs/hic_results/matrix/GSE125129_AdnpKO1/iced/40000/GSE125129_AdnpKO1_40000_iced.matrix',
indices_path = '/data/zhangyong/yyp/yypold/hic/chenhebing/cs/hic_results/matrix/GSE125129_AdnpKO1/raw/40000/GSE125129_AdnpKO1_40000_abs.bed',
sample_name = "WAPL",
colour = "red")
jpeg(file = "style.jpg") ##plot函数输出的文件名字
cisChrom_out <- cis_trans( list(Ha40, Hb40) )
barplot(cisChrom_out$cis, names.arg = cisChrom_out$sample, ylim = c(0, 100) )
abline(h = 90, col = 'red', lty = 3)
abline(h = 93, col = 'red', lty = 3)
dev.off()
今天,我想用R画图并保存成jpg文件。应用的方法如下:
setwd("c://")
plot(1:10)
rect(1, 5, 3, 7, col="white")
savePlot("CTplot", type=c("jpg"),device=dev.cur(),restoreConsole=TRUE)
在R编辑器下运行成功。
可是我在cmd中,使用Rscript 去运行就无法生成图片文件。
问题应该是,命令行无法在窗口中画图,所以也就无法保存图片。
后来,在网上找到一种直接生成图片的方法。
setwd("c://")
jpeg(file="myplot.jpeg")
plot(1:10)
rect(1, 5, 3, 7, col="white")
dev.off()
命令行下运行成功。同时查到另外的几种格式生成函数
第一种png格式
png(file="myplot.png", bg="transparent")
...
dev.off()
第二种jpeg格式
jpeg(file="myplot.jpeg")
...
dev.off()
第三种pdf格式
pdf(file="myplot.pdf")