This function implements the mathematical core of FracFixR: fitting a non-negative least squares (NNLS) regression to estimate fraction weights and correct individual transcript abundances.
ProcessReplicate(RepMat, transcriptlist)List containing:
Propestimates: Proportion estimates for each transcript
Coefficients: NNLS regression coefficients (fraction weights)
Fractions: Normalized fraction proportions
plot: Diagnostic plot of fitted vs residuals
Data frame with transcripts as rows, samples as columns. Must include a "Total" column representing the whole cell lysate
Character vector of transcript IDs to use for regression. These should be informative transcripts in the appropriate abundance range
Mathematical basis: $$Total = \alpha_0 + \alpha_1 \times Fraction1 + \alpha_2 \times Fraction2 + ... + \epsilon$$ Where \(\alpha_0\) represents the "lost" fraction and other \(\alpha_i\) are fraction weights