junctionScore.Rd
Function to calculate the JCC score for a given gene
junctionScore(uniqreads, mmreads, predcovs, g, beta = 1, ...)
Vector with the number of uniquely mapping junction-spanning reads for each junction in the gene.
Vector with the number of multimapping junction-spanning reads for each junction in the gene.
Vector with the predicted junction coverages for each junction in the gene.
Junction weight function that takes as input the fraction of uniquely mapping reads for the junction, and returns a weight for the junction contribution to the score.
Scaling exponent. If beta=1, predicted coverages are scaled to have the same (weighted) sum as the observed coverages. If beta=0, no scaling is done.
Additional parameters passed to g
.
The value of the JCC score for the gene