This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Gini Coefficient and Lorenz curve to measure uniformity of coverage

Hello, would you suggest me a nice tool to compute Gini index / plot Lorenz curve to estimate uniformity of coverage over the genome other than the R Package ‘ineq’ or the old R Package 'htSeqTools' ?

Thank you in advance!

gini lorenz curve coverage uniformity
• 1,987 views • link
updated by abascalfederico abascalfederico • written by am • am

2 answers

Hi Federico, I kept using ineq or python code: https://github.com/auroramaurizio/my_DNA_seq_pipelines/blob/main/GINI_index/Snakefile. I also suggest you to check the bioconductor R package SCOPE https://www.bioconductor.org/packages/devel/bioc/vignettes/SCOPE/inst/doc/SCOPE_vignette.html#gini-coefficient.

HTH!

• 0 views • link
written by am • am

Thank you! I finally used some code from my supervisor for the Lorenz curves, it's relatively easy, and DescTools::Gini to compute the gini coefficient. The input data was just 10Mbp of samtools depth

• 0 views • link
written by abascalfederico abascalfederico

Hi, did you find a program to do this?

• 0 views • link
written by abascalfederico abascalfederico

Log in to answer this question.

More posts like this