localScore (version 1.0.6)

karlin: Karlin [p-value] [iid]

Description

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

Usage

karlin(
  localScore,
  sequence_length,
  score_probabilities,
  sequence_min,
  sequence_max
)

Arguments

localScore

the observed local score

sequence_length

length of the sequence (at least several hundreds)

score_probabilities

the probabilities for each unique score from lowest to greatest

sequence_min

minimum score

sequence_max

maximum score

Value

A double representing the probability of a localScore as high as the one given as argument

Details

This method works the better the longer the sequence is.

Examples

# NOT RUN {
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)
# }