This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Are logr ratios supposed to be normally distributed?

In a CNV analysis workflow, do you expect the logr ratios along the genome to be from a normal (or approximately normal) distribution? Why?

statistics cnv

3 answers

The read depths are technically a mixture of Poisson or log-Poisson distributions (depending on whether PCR was done), I think, but usually modelled well enough with the inverse binomial distribution. In log scale it is a bit more overdispersed than Poisson, and the second parameter in the density function for the inverse binomial distribution can be used to account for that.

It may be "normal enough" for your needs, with the understanding that there will be more "outliers" or extreme values than you'd expect to see under the normal distribution. Robust statistics (e.g. MAD or Tukey's biweight midvariance instead of standard deviation) can be used to insulate your results from these extreme values if you'd like to use statistical methods that assume normality.

you can easily know by running normality tests like shapiro..

I know. I also did a qqplot. Both methods show that my data considerably deviates from normality. My question was whether that is what you would expect for logr ratios (i.e. should I be concerned). Sorry if I failed to make that clear in the original question.

Or in R draw a histogram with "hist" on the underlying variable to take a look.

Log in to answer this question.