#delimit ; use ~/work/optmat/revisited/plugin/birthseq; keep if _n<=100; matrix subsmat = (0,1,2,3 \ 1,0,1,2 \ 2,1,0,1 \ 3,2,1,0 ); scalar indelcost = 2; set matsize 700; matrix resmat0 = J(_N,_N,0); gen len=73; timer on 2; oma len state1-state73, subs(subsmat) indel(2) pwdist(resmat0) length(73); timer off 2; timer list 2; mata:; B1=st_matrix("resmat0") mata matsave bspure B1, replace end clustermat wards resmat0,add; cluster generate g8=groups(8); gen str73 stripe = ""; foreach x of varlist state1-state73 {; replace stripe = stripe + "." if `x'==4; replace stripe = stripe + "=" if `x'==3; replace stripe = stripe + "x" if `x'==2; replace stripe = stripe + "M" if `x'==1; };