# NOT RUN {
data(MidSeq)
MidSeq
nchar(MidSeq)
data(dico)
MidSeqScore=CharSequence2ScoreSequence(MidSeq,dico)
MidSeqScore[1:30]
localScoreC(MidSeqScore)$localScore
prob1 = scoreSequences2probabilityVector(list(MidSeqScore))
daudin(localScore = 52, sequence_length = nchar(MidSeq),
score_probabilities = prob1,
sequence_min = min(MidSeqScore),
sequence_max = max(MidSeqScore))
score=-5:5
prob2=c(0.15,0.15,0.1,0.1,0.0,0.05,0.15,0.05,0.2,0.0,0.05)
daudin(localScore = 52, sequence_length = nchar(MidSeq),
score_probabilities = prob2,
sequence_min = min(MidSeqScore),
sequence_max = max(MidSeqScore))
# }
Run the code above in your browser using DataCamp Workspace