廿TT

譬如水怙牛過窓櫺 頭角四蹄都過了 因甚麼尾巴過不得

2017-03-29から1日間の記事一覧

pymc3 で状態空間(ローカルレベル)モデル

Python ハローワールドです。Rodeo(Yhat End-to-End Data Science Platform: Rodeo)っていう RStudio みたいなやつをいれました。なかなか快適です。 import numpy as np import pymc3 as pm import matplotlib.pyplot as plt import statistics as stat i…

R によるすごくかんたんなパーティクルフィルタの実装例

これであってるのかあんまり自信ない。主に http://daweb.ism.ac.jp/koza/koza2008/PF_Nakano20081030.pdf を参考にしました。パーティクルフィルタのアルゴリズムは、 一期先の予測を乱数でばらまく 得られた乱数に尤度の重みをつけてまき直す だと理解しま…

Rcpp で独立メトロポリスヘイスティングス

独立メトロポリス・ヘイスティングス法を用いたベイズ推測の簡単な例題 - 廿TT でやったのと同じことを Rcpp で書いてみた。 ハローワールド。C++ のコードはこう。 #include <Rcpp.h> using namespace Rcpp; double lik(double lambda, NumericVector dat) { double</rcpp.h>…