|
Size: 2772
Comment:
|
← Revision 14 as of 2020-06-01 18:40:49 ⇥
Size: 2724
Comment:
|
| Deletions are marked like this. | Additions are marked like this. |
| Line 11: | Line 11: |
| {{{ | {{{#!sagecell |
| Line 13: | Line 13: |
| vv = []; nn = 0 | vv = [] nn = 0 |
| Line 15: | Line 16: |
| def foo(pts = checkbox(True, "Show points"), refresh = checkbox(False, "New random walk every time"), steps = (50,(10..500))): |
def foo(pts=checkbox(True, "Show points"), refresh=checkbox(False, "New random walk every time"), steps=slider([10..500],default=50)): |
| Line 19: | Line 20: |
| # checking or unchecking the points checkbox doesn't change # the random walk. |
# checking or unchecking the points checkbox does not change # the random walk. |
| Line 23: | Line 24: |
| if refresh or len(vv) == 0: s = 0; v = [(0,0)] |
if refresh or not vv: s = 0 v = [(0, 0)] |
| Line 31: | Line 33: |
| s = vv[-1][1]; j = len(vv) | s = vv[-1][1] j = len(vv) |
| Line 34: | Line 37: |
| vv.append((i+j,s)) | vv.append((i + j, s)) |
| Line 39: | Line 42: |
| if pts: L += points(v, pointsize=10, rgbcolor='red') show(L, xmin=0, figsize=[8,3]) |
if pts: L += points(v, pointsize=10, rgbcolor='red') show(L, xmin=0, figsize=[8, 3]) |
| Line 45: | Line 49: |
| {{{ | {{{#!sagecell |
| Line 47: | Line 51: |
| def rwalk3d(n=(50,1000), frame=True): pnt = [0,0,0] |
def rwalk3d(n=slider([50..1000]), frame=True): pnt = [0, 0, 0] |
| Line 51: | Line 55: |
| pnt[0] += random()-0.5 pnt[1] += random()-0.5 pnt[2] += random()-0.5 |
pnt[0] += random() - 0.5 pnt[1] += random() - 0.5 pnt[2] += random() - 0.5 |
| Line 55: | Line 59: |
| show(line3d(v,color='black'),aspect_ratio=[1,1,1],frame=frame) | show(line3d(v, color='black'), aspect_ratio=[1, 1, 1], frame=frame) |
| Line 63: | Line 67: |
| based on the classic example in "Biological Sequence Analysis" by Durbin et al. To get the dice picture you need to download this {{worksheet copy|dishonest_casino.sws}} {{{ |
Based on the classic example in "Biological Sequence Analysis" by Durbin et al. {{{#!sagecell |
| Line 70: | Line 75: |
| a = os.system('cp '+DATA+'dice_sm.png ./dice.png') | |
| Line 72: | Line 76: |
| html('<h3>The Occasionally Dishonest Casino</h3>') html('<img src="'+DATA+'dice_sm.png">') html('Emissions:'+str(test[0]).replace(',','').replace(' ','')[1:-1]) |
pretty_print(html('<h3>The Occasionally Dishonest Casino</h3>')) pretty_print(html('Emissions:'+str(test[0]).replace(',','').replace(' ','')[1:-1])) |
| Line 77: | Line 80: |
| html('Viterbi: '+vit_str) | pretty_print(html('Viterbi: '+vit_str)) |
| Line 80: | Line 83: |
| html('Actual: '+ actual_str) | pretty_print(html('Actual: '+ actual_str)) |
Sage Interactions - Statistics and Probability
goto interact main page
Contents
A Random Walk
by William Stein
3D Random Walk
Hidden Markov Model: The Occasionally Dishonest Casino
by Marshall Hampton
Based on the classic example in "Biological Sequence Analysis" by Durbin et al.
