powered by
Calculates an approximated p-value of a given local score value and a long sequence length in the identically and independantly distributed model for the sequence. See also mcc() function for another approximated method in the i.i.d. model
karlin( localScore, sequence_length, score_probabilities, sequence_min, sequence_max )
A double representing the probability of a localScore as high as the one given as argument
the observed local score
length of the sequence (at least several hundreds)
the probabilities for each unique score from lowest to greatest
minimum score
maximum score
This method works the better the longer the sequence is.
karlin(150, 10000, c(0.08, 0.32, 0.08, 0.00, 0.08, 0.00, 0.00, 0.08, 0.02, 0.32, 0.02), -5, 5)
Run the code above in your browser using DataLab