2012-06-18 4 views
-1

Я использовал следующий код. Но вместо начала с -5000 ось x начинается с 0. Кроме того, я хотел бы указать интервал 100 для диапазонов x и y-осей. Я пытался использовать xaxt='n' и axis(side=1), но это не помогло. Следующий код является тем, который я сейчас работаю наНастройка интервалов оси X 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(), поэтому вы можете заменить их на ось() в следующей строке.

+0

Привет, я попытался следующие 'ч <-hist (х, оси = FALSE, то xaxt = 'N', yaxt = 'N', Col = "зеленый", основная = "GATA3-TPMDistribution", xlab = «TPMValues», вероятность = TRUE) ось (сторона = 1, at = seq (-162000,162000, 1)) ось (сторона = 2, at = seq (0,0.008,0.000001)) ' I Я получаю гистограмму с отрицательными координатами. Но начальная точка кривой и пик усекаются. Любые мысли о том, почему это происходит? –

Смежные вопросы