# NOT RUN {
data(LongSeq)
LongSeq
nchar(LongSeq)
data(dico)
LongSeqScore=CharSequence2ScoreSequence(LongSeq,dico)
LongSeqScore[1:50]
localScoreC(LongSeqScore)$localScore
LS=localScoreC(LongSeqScore)$localScore[1]
prob1 = scoreSequences2probabilityVector(list(LongSeqScore))
daudin(localScore = LS, sequence_length = nchar(LongSeq),
score_probabilities = prob1,
sequence_min = min(LongSeqScore),
sequence_max = max(LongSeqScore))
karlin(localScore = LS, sequence_length = nchar(LongSeq),
score_probabilities = prob1,
sequence_min = min(LongSeqScore),
sequence_max = max(LongSeqScore))
# }
Run the code above in your browser using DataCamp Workspace