But we need to have some code to examine if we want to talk about algorithms. So, if you would like to follow the discussion that's coming later, please paste the code below into a textfile on your desktop (my copy is named gibbs.py but you can do whatever you want). You will also need the code from here, saved as motif.py, and the code from here, saved as scoring.py. You can run it in the usual way by doing "python gibbs.py."
Here is the code, and output for the particular settings. Discussion coming soon.
# Lawrence et al. Gibbs Sampler algorithm
numSeq = 10