2012-06-18 376 views
-1

我一直在使用下面的代碼。但是,不是從-5000開始,x軸從0開始。另外,我想爲x和y軸範圍指定100的間隔。我試圖使用xaxt='n'axis(side=1),但這沒有幫助。以下代碼是我目前正在使用的代碼R直方圖中X和Y軸的間隔設置

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

你要指定斷裂,可能是在一個陣列,在這裏看到: 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

您可能需要使用軸= FALSE在你的hist()中,所以你可以在後續行中用axis()替換它們。

+0

嗨,我試過以下 'h <-hist(x,axes = FALSE,xaxt ='n',yaxt ='n',col =「green」,main =「GATA3-TPMDistribution」,xlab = 「TPMValues」,概率= TRUE) 軸(側= 1,在= SEQ(-162000,162000,1)) 軸(側= 2,在= SEQ(0,0.008,0.000001))' 我我得到的直方圖與負面的座標。但是,曲線的起點和峯值正在被截斷。任何想法爲什麼發生這種情況? –