CRISPR Screen QC, Normalization, and Advanced Methods
Screen Quality Metrics
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
from scipy import stats
def gini_index(counts: np.ndarray) -> float:
"""Gini index of guide count distribution. Perfect evenness = 0, one guide all reads = 1.
Good screens: Gini < 0.1. Problematic: Gini > 0.2.
"""
sorted_counts = np.sort(counts)
n = len(sorted_counts)
cumulative = np.cumsum(sorted_counts)
[Description truncada. Veja o README completo no GitHub.]