昔、マルコフ連鎖インベンションみたいな事をやったり、ルールベースの対位法もどきみたいな事をやって放置していたけどようやく再開。 やり方は,楽譜ソフトに音を入力して music21: a Toolkit for Computer-Aided Musicology というPythonのライブラリで読…

状態遷移の式 α(t+1)=T*αt+ε (記号はKFASから） でTの中に2種類のQが入るだけなので、何も問題がないはずです。 nn=3 a=rep(c(1,2,3),5*nn) b=rep(c(5,5,0,0,0),3*nn) d=10+a+b ts.plot(d[1:30]) scode=" data{ int N; real y[N]; } parameters{ real sigma;…

UCI Machine Learning Repository: Poker Hand Data Set ポーカーの手の判定。プログラミングの課題でよくあるけど、統計的な分類で判定します。 > head(df) C 1 C 2 C 3 C 4 C 5 C 6 C 7 C 8 C 9 C 10 C 11 1 3 12 3 11 3 13 3 10 3 1 9 2 4 10 4 11 4 1 4 …

Pythonとjavascriptで http://biones.dip.jp/space/space.py （鯖が止まってると動かない可能性あり） 的なモノを昔作って何か分析しようとおもいつつ放置していたのに再着手。 青色の点が目的地。 歩行軌跡間の「近さ」はデータの長さが違っても類似度が計…

Kawazoe (@riverplus) n 個のリール（数字が描かれている部分です）を持つスロットマシンを回します。 各リールには、0 から 9 のいずれかの数字がランダムに出ます。 このとき、「最も多く出現した数字の出現回数」に等しいドルが賞⾦として得られます。 例…

練習。 UCIのデータセットより。 UCI Machine Learning Repository: Wine Quality Data Set hoge=sample(1:1599,900) #900件をトレーニング用に train=wine[hoge,] test=wine[-hoge,] # 0.66くらい相関がある特徴があったけど、とりあえずただの線形回帰 tra…

下の例で計算すると0.96。 線形トレンド 観測方程式 の状態空間で、x,y（変数名の方）の残差の相関をstanで計算すると-0.12程度になった。 N.t=60 mu.x = numeric(N.t) x = numeric(N.t) mu.y = numeric(N.t) x=y=c() mu.x[1] = mu.y[1] = 10 s1 = 2 s2 = 1 …