myPar <- bsdPar(10, 2)
M <- genSeq(myPar, 2)
getProb(M)
# all Sequences
par <- pbrPar(bc=c(2,2))
refSet <- getAllSeq(myPar)
probs <- getProb(refSet)
# sequences with probabilities
cbind(probs, refSet$M)
Run the code above in your browser using DataLab