Learn R Programming

localScore (version 1.0.8)

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
)

Value

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

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

Details

This method works the better the longer the sequence is.

Examples

Run this code
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