首页 > 解决方案 > R图中的Hmisc minor.tick范围

问题描述

我想使用RNA contigs vs Time 的 R plot plot 为随着时间(timept)映射到 mtDNA 基因组(mtDNA_rds)的数量 DNA contigs 的简单散点图添加小刻度线。我在 Hmisc 包中使用了 minor.tick ,它有效,除了它沿着 xlim 和 ylim 范围两侧的轴放置小刻度 - 我绝对不想要 - 特别是对于 sub 0 区域。我尝试将 xaxs 选项从“r”(常规)更改为“i”(内部),这解决了问题,但在第 0 小时和“Contigs Mapped ..”时,我的数据点符号减半。值是 0。我知道我还有一些额外的绘图格式要做 - 但首先,我如何摆脱那些无关紧要的小刻度?这是我的数据和 R 脚本。

install.packages("Hmisc")
library(Hmisc)
timepts = c(0,2,4,6,8,12,24,36,48,72,96,120)
mtDNA_rds = c(1876,282,82,327,26,11,3,2,1,0,1,0)
Zf.0 = data.frame(timepts,mtDNA_rds)  
plot(Zf.0$timept,Zf.0$mtDNA_rds,bty="l",xlim=c(0,120),ylim=c(0,2000),pch=19,xaxs="r",yaxs="r",xaxt="n",yaxt="n",xlab="Hours Post Removal of Zf",ylab="RNA Contigs Mapped to Zf mtDNA Genome",cex.lab = 2,cex = 1.75) 
axis(1,at=c(0,5,10,15,20,25,30,35,40,45,50,55,60,65,70,75,80,85,90,95,100,105,110,115,120))
axis(2,at=c(0,250,500,750,1000,1250,1500,1750,2000))
minor.tick(nx=20,ny=10,tick.ratio=0.65)

标签: rplothmisc

解决方案


推荐阅读