我一直在使用以下代码。但是,不是从-5000开始,x轴从0开始。另外,我想为x轴和y轴范围指定100的间隔。我尝试使用xaxt='n'
和axis(side=1)
,但这没有帮助。以下代码是我目前正在处理的代码
setwd("/path/")
data<-read.table("input.txt",sep="\t",header=F)
x<-data$V4
h<-hist(x, col="green",xlim=c(-5000, 162000),ylim=c(0,0.0003),main="TPMDistribution",xlab="TPMValues",probability=TRUE)
s=sd(x)
m=mean(x)
curve(dnorm(x,mean=m,sd=s),add=TRUE,lwd=3,col="red")
lines(density(x),col="blue")
abline(v=mean(x),col="blue")
mtext(paste("mean",round(mean(x),1),";sd",round(sd(x),1),";N",length(x),sep=""),side=1,cex=.75)
dev.off()
所有人都很感激。提前谢谢。
答案 0 :(得分:0)
您可能希望指定中断,可能在数组中,请参阅此处: http://stat.ethz.ch/R-manual/R-patched/library/graphics/html/hist.html
使用以下方法调整轴: http://stat.ethz.ch/R-manual/R-patched/library/graphics/html/axis.html
您可能需要在hist()中使用axes = FALSE,因此可以在后续行中用axis()替换它们。