data(Seq1093)
Seq1093
nchar(Seq1093)
data(HydroScore)
seqScore=CharSequence2ScoreSequence(Seq1093,HydroScore)
seqScore[1:50]
localScoreC(seqScore)$localScore
LS=localScoreC(seqScore)$localScore[1]
prob1 = scoreSequences2probabilityVector(list(seqScore))
daudin(localScore = LS, sequence_length = nchar(Seq1093),
score_probabilities = prob1,
sequence_min = min(seqScore),
sequence_max = max(seqScore))
karlin(localScore = LS, sequence_length = nchar(Seq1093),
score_probabilities = prob1,
sequence_min = min(seqScore),
sequence_max = max(seqScore))
Run the code above in your browser using DataLab