ジョンとヨーコのイマジン日記

キョウとアンナのラヴラブダイアイリー改め、ジョンとヨーコのイマジン日記です。

2020-01-01から1年間の記事一覧

事後分布を正規分布で近似する変分推論のアルゴリズムを Julia で書く

アルゴリズム 機械学習スタートアップシリーズ ベイズ推論による機械学習入門作者:須山敦志講談社Amazon 須山『ベイズ推論による機械学習入門』ではロジスティック回帰とニューラルネットのところで近似事後分布として1次元の正規分布を仮定して、変分パラメ…

Albert (2008) 打者の調子の波のモデル化

モチベーション 以前にAlbert (2008)、 https://www.stat.berkeley.edu/~aldous/157/Papers/albert_streaky.pdf を読んでやってみたことのJulia版です。次の0と1の羅列はカルロス・ギーエン(カルロス・ギーエン - Wikipedia)という選手の2005年の打撃成績…

Julia によるかんたんなパーティクルフィルタの例

ローカルレベルモデルです。(システムモデル) (観測モデル)パーティクルフィルタで固定パラメータの とか を推定する方法は色々提案されてたりするんだけど、シミュレーテッドアニーリング法に丸投げです。「こう書くともっと速くなるよ」みたいなご指摘…

Julia でかんたんラプラス近似

データを生成した分布を「真の分布」と呼ぶことにする。確率モデルを として、確率モデルにパラメータ の最尤推定量 を代入し、 「真の分布はおよそ だろう」 と推測することを最尤推定による予測分布(略して、MLE)と呼ぶことにする。最尤推定量(またはMA…

区間打ち切りデータの最尤推定のシミュレーション:JuliaとR

Julia です。 @time using Optim @time using Distributions @time using StatsFuns function make_dat_gamma0(n,a,b) t = rand(Gamma(a,b),n) E_len = rand(Exponential(),n) u = rand(n) E_R = t - E_len .* (1.0.-u) E_L = t + E_len .* u E_L, E_R end f…

変分ベイズで RFM 指標から顧客生涯価値を計算してみる

R

モデルと尤度 RFM 指標から将来の購買回数を予測する Pareto / NBD モデルから派生したモデルに BG / NBD モデルがあります。Pareto / NBD モデルより計算がかんたんです。モデル: 顧客の購買はレート の定常ポアソン過程に従う 顧客の購買のたびに確率 で…

ドミノでピノークル

Domino Pinochle - domino game rules より、ドミノを使ったトリックテイキングゲーム、ドミノピノークルを紹介します。 人数 3人 必要なもの W6のドミノ 配り方 各人に9枚ずつ。ディーラーは余った一つの牌を表にします。この牌をウィドウと呼びます。 プレ…

花札の遊び方(新潟県のある地域のルール)

花札にはローカルルールが大量にあります。 ここでは私が子供の頃遊んでいたルールを紹介します。 任天堂 花札 丸福天狗 赤 メディア: おもちゃ&ホビー 任天堂 花札 丸福天狗 黒 メディア: おもちゃ&ホビー ゲームの流れ 花合わせ - Wikipedia ゲームの流…

ドミノゲーム『テキサスムーン』ルール紹介

はじめに ムーン(テキサスムーン、シュート・ザ・ムーン)はドミノを使ったトリックテイキングゲームです。 ムーンのルールは ムーン - ゲームばっかりやってないで で解説されていますが、アプリ版のルールとは違うようなので、あらためて紹介します。 Tex…

ツーペン(Toepen)というゲームのアプリを買った

ToepenHD - leukste kaartspel! GraafICT ゲーム ¥150 Toepen - leukste kaartspel! GraafICT ゲーム ¥250 ゲームの遊びかたは以下に書きます。 大まかな流れ speel でゲームスタート。 ツーペン(Toepen;トゥーペン?)は 2 から 6 を抜いた 32 枚のトラン…

7-8(2人用トリックテイキングゲーム)

How To Play 3 2 5 を参考に 7-8 (Sath-Aath) というトランプゲームを紹介します。インドやパキスタンでメジャーなゲームだそうです。 ブラックバイスクル メディア: スポーツ用品 人数 2人 使用するカード スペード: 7, 8, 9, 10, J, Q, K, Aダイヤ: 8, 9…

シディバラニ・ヤス(4人用トリックテイキングゲーム)

Rules of Card Games: Sidi Barrani Jass を参考にシディバラニ・ヤスというゲームを紹介します。 バイスクル ブラックブルー 発売日: 2018/11/09 メディア: おもちゃ&ホビー 準備 使用するカード:ふつうの52枚のトランプから2〜5をのぞいた36枚プレイヤー…

将校シャーフコップ(Officers' Schafkopf):手札公開の二人用トリックテイキングゲーム

Officers' Schafkopf - Wikipedia を参考にして紹介します。 将校スカートと同じようなゲームです。 gokurakism.com 将校スカートはスマホのアプリもあります。 play.google.com Officer's Skat Victor Olaru Games Free apps.apple.com 準備 使用するカード…

ペニーチ(Penneech):二人用トリックテイキングゲーム

古いイギリスのゲームだそうです。 Pennech: historical card game を参考にして紹介します。 任天堂 トランプ ナップ 818 赤 メディア: おもちゃ&ホビー 準備 使用するカード:ふつうの52枚のトランプ(ジョーカー除く)プレイヤー:2人カードの強さ:強い…

カラーブラインド(3人用トリックテイキングゲーム)

デヴィッド・パーレット氏の創作ゲームからカラーブラインド(Colour-blind: an original card game)を紹介します。 バウンサー(http://abrahamcow.hatenablog.com/entry/2017/09/14/023529)に似た感じのルールです。 用意するもの 人数: 3カード: 52枚の…

バウンサー(3人用トリックテイキングゲーム)

バウンサー(BOUNCERS)というトランプゲームを紹介します。 デヴィッド・パーレット氏の創作ゲームで、ハーツ(http://abrahamcow.hatenablog.com/entry/2015/01/08/014215)のバリエーションの一つと見ることもできます。 Bouncers: an original card game…

フランク(3人用トリックテイキングゲーム)

フランク(Flunk)というトランプゲームを紹介します。 Flunk: an original card game by David Parlett デヴィッド・パーレット氏の創作ゲームです。 ブラックバイスクル メディア: スポーツ用品 人数 3人 カードの強さ 強い方から順に、A>K>Q>J>10>9>8>7>6…

クロスパルパス(2人用トリックテイキングゲーム)

デヴィッド・パーレット氏の創作ゲームからクロスパルパス(Cross Purposes: an original card game)を紹介します。 4人用のルールもありますが、ここでは2人用のルールについてのみ書きます。 用意するもの 普通の52枚のトランプ. ルール ディールごとにプ…

2人用スペード(トランプゲーム)の遊び方

スペードはもともと4人で遊ぶトランプゲームですが、2人用のルールを紹介します。 How to Play Spades With 2 People の説明を参照しています。 ブラックバイスクル メディア: スポーツ用品 ジャンル トリックテイキングゲーム:切り札あり、マストフォロー…

マジョレ(Marjolet):2人用トリックテイキングゲーム

べジークやピノークルに近いフランスのゲームだそうです。ブラックバイスクルメディア: スポーツ用品 カード 2から6を除く32枚のトランプを使います。 ランク カードの強さは強い方から順にA、10、K、Q、J、9、8、7です。 カードのポイントはAと10が10点、他…

ggplot2でヒストグラムを箱ひげ図風に並べてプロットするパッケージggbrickを書いている

R

これです:GitHub - abikoushi/ggbrick: An R package for brick plotggtetrisって名前にしようかと思ったんだけど、もうある(GitHub - EmilHvitfeldt/ggtetris: Create Tetris Chart Visualizations in R)みたいなので、ggbrickにしました。brickはレンガ…

藤原香織、渡辺澄夫(2006)ベイズ尤度比検定による変化点検出のシミュレーション

R

これです:http://watanabe-www.math.dis.titech.ac.jp/users/fujiwara/doc/fujiwara_ibis2006.ppt.pdfこの研究成果は論文化されているようですが, オープンアクセスではないみたいです:https://search.ieice.org/bin/summary.php?id=j91-d_4_889&category=…

今日の川柳

ぼくの以前のブログには「今日の川柳」というモジュールが設置されていた。こんな風だった: function OnButtonClick() { var a=["多分A明日BならやつはC","歯ブラシを眺めて泣いた夜もある","本当はラッパーになりたかった","ていねいなしらない人の顔写真"…

ggplot2 で floating bar chart を書く R のパッケージ ggfloatbar を作っている

R

GitHub - abikoushi/ggfloatbar: Floating bar chart on ggplot2です。インストールは devtools::install_github("abikoushi/ggfloatbar") でたぶんいけます。たぶんまだ不具合とかあると思うし、説明とかぜんぜん書いてないです。なにか要望とか意見、感想…