Skip to contents

Identify optimal reads per cell range for power analysis grid (C++)

identify_library_size_range_cpp.Rd

C++ implementation that determines the minimum and maximum reads per cell values for power analysis grid generation using binary search on the S-M curve. Uses saturation-based thresholds (10% and 80%) instead of platform-specific minimums.

Usage

identify_library_size_range_cpp(experimental_platform, rSAC_fn_wrapper)

Arguments

experimental_platform

String. Experimental platform identifier (kept for compatibility, not used).

rSAC_fn_wrapper

List. Parameters from library_estimation containing:

  • method_used: "ZTNB" or "RFA"

  • reads_norm: Normalization constant

  • n_cells: Number of cells

  • UMI_per_cell_at_saturation: Maximum UMI per cell

  • For ZTNB: L, size, mu

  • For RFA: valid_estimator, coefs_real, coefs_imag, poles_real, poles_imag, or constant_value

Value

List with min_reads_per_cell and max_reads_per_cell elements.

Details

This C++ implementation uses efficient binary search to find the reads per cell range for power analysis. Uses saturation-based thresholds:

  • Minimum reads: 10% UMI saturation (dynamic based on UMI_per_cell_at_saturation)

  • Maximum reads: 80% UMI saturation (diminishing returns beyond this point)

See also

identify_library_size_range for R version