Fit probability law on genes length distribution · Biostar

This is a test version of Biostars. For the public version, visit https://www.biostars.org.

Fit probability law on genes length distribution

Hi, I am trying to fit a probability law on this distribution representing the genes length of E. rectale. I have tried to fit a skew normal distribution on it (in red) Do you think this is the most appropriate law? Thanks in advance. Genes length distribution of E. rectale

genes length distribution

• 1,640 views • link

updated 9.7 years ago by Steven Lakin

• written 9.7 years ago by Florian Plaza Oñate •

Try plotting the log(length).

• 0 views • link

written 9.7 years ago by Devon Ryan

1 answer

You may want to try an overdispersed Poisson distribution or the negative binomial. For a variety of reasons, these distributions tend to fit biological data well in cases of discrete counts (length can likely be modeled this way with success). Distribution choice largely depends on what you intend to model.

• 0 views • link

written 9.7 years ago by Steven Lakin

Log in to answer this question.

More posts like this

biological vs. technical definition of "insert size"
written 3.5 years ago by sovrappensiero •

I'm trying to simulate metagenomic reads using Grinder. I want to simulate paired-end reads, so I need to select arguments for insert size and distribution. …
Forum: Is there ever a good reason to use a log transformation instead of Box-Cox?
written 4.8 years ago by Mensur Dlakic

I often deal with non-normal data distributions, and for some downstream application it is needed to transform the data to at least resemble a normal …
WGCNA: scale-free topology fit indexes are negative
written 5.1 years ago by Erica •

Hi all, I am trying to perform WGCNA on my RNA-sequencing data (29 samples). First, I obtained the TPM and did a log2 transformation. After …
Using R package to implement Bayesian phase I/II dose-finding design for three outcomes
written 5.2 years ago by Ran •

Dear all, I am trying to implement the work done by Suyu Liu, "A Bayesian Phase I/II Trial Design for Immunotherapy", using R, since the …
ChromoMap Error in 1:nrow(temp.df) : argument of length 0
written 6.1 years ago by kousi31 •

I am trying to plot the distribution of genes on chromosome using `chromoMap`. I prepared a chromosome_file. txt as per [manual][1] NC_037545.1 1 202105980 NC_037546.1 …
Question regarding deFuse and fusion calling
written 7.5 years ago by CY

I recently tried [deFuse][1] and raised several questions after reading the paper. If I understand correctly, two necessary conditions set on discordant reads are hard …
remove batch effect for multiple tumor data
written 9.4 years ago by hellocita •

Dear all, I am dealing with a data contains samples with 7 tissues(normal control vs. tumor), and 2 experimental batches. These control and tumor of …
MA plots DESeq: strange MA plots
written 11.1 years ago by VHahaut

Hello! I am currently analyzing with DESeq2 a group 32 ovine tumors against 3 healthy controls. 80% of the tumor samples are technical replicates. I …
estimation of size factors in DESeq2 analysis
written 12.1 years ago by Assa Yeroslaviz

Hi all, in relation to a mail from January this yearin the r-help community, I followed Simon's advice to do my analyses in DESeq2 instead …
Fitting Distribution Of The Gene Expression Levels Of Rna Seq Data
written 12.5 years ago by k.nirmalraman

<p>I am trying to fit a distribution of the gene expression (DESeq Normalized counts) of each sample from a RNA Seq data. And I came …