Hi,

I have a file of the following style

chr1 122333 1223344 1.87887877778
chr1 2324343 433433 1.98794334343
chr2 343545 54545454 9.0009123232
chr3 843749 89449379 10.004343722
................................................
................................................

setwd("/path/")
data<-read.table("input.txt",sep="\t",header=F)
x<-data$V4
h<-hist(x, axes=FALSE, xaxt='n', yaxt='n',col="green",main="GATA3-TPMDistribution",xlab="TPMValues",probability=TRUE)
axis (side=1, at=seq(0,162450, 1000))
axis (side=2, at=seq(0,0.008,0.0001))
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()

In what way is it being truncated? What are the results you expect, and what are the results you get?

How large is your input data in bytes, not records?

Hi Corona,

I could only see the end of the normal curve. But, I would like to see the whole normal curve, I mean its start, peak and stop.

I would like to look at the distribution of the values in the 4th column using an histogram. The results I see are fine enough, but without properly looking into the curve, I can't infer anything.

My input file is around 30MB.

Please let me know if you would like to know anything in detail.

I would like to have my histogram with a curve like this -




But, I am getting it this way -