Hi,
I would like to calculate the Gini coefficient for my TCR repertoire data. I have seen that this can be calculated using the DescTools and ineq packages in R, however from the package documentation, I am unclear what input is required in the context of TCR repertoire data.
Is anyone able to help with this?
Best wishes,
Lucy
1 answer
Hi,
The Gini coefficient is used in many diverse ways, it seems.
In the context of TCR-seq, you could use it to check for 'inequality' among clonotypes on a per sample or group (condition) basis.
Which data have you got, currently?
Kevin
I see. The function in DescTools appears to just literally accept, at minimal, a vector of values, like c(45,33,67, ..., 45). I suppose that, technically, a Gini coefficient can be calculated from this, where each number represents usage (normalised read count?) of a particular clonotype.
Based on this vignette (https://www.statology.org/gini-coefficient-in-r/), I think the required input would be the clonotype assignment for each cell rather than the abundances.
Log in to answer this question.
-
ESTIMATE -tidyestimate package- tumour purity
written by Nona •
-
Specify the order of the sets in VennDiagram R?
written by Lucy •
-
How does --primary work with -M and --fractional options featureCounts?
written by Lucy •
-
featureCounts parameters
written by Lucy •
-
Seurat sctransform - does anyone know of a good explanation for biologists?
written by Lucy •
-
Goodness-of-fit test - observed vs. expected (some expected values < 5)
written by Lucy •
-
Identification of modules of genes that vary over pseudotime (scRNA-seq pseudotime analysis)
written by Lucy •
-
Negbinom vs. logistic regression for Seurat FindMarkers()
written by Lucy •
-
Seurat sctransform return.only.var.genes = FALSE
written by Lucy •
-
Error with fitCentipede
written by Lucy •
