Box-MullerをPythonで

PRML11章。一様乱数から正規乱数を発生させる。

from math import sin,cos,pi,log,sqrt
from random import random

def box_muller():
    z1 = random()
    z2 = random()
    return sqrt(-2*log(z2)) * sin(2*pi*z1)

if __name__ == "__main__":
    from pylab import *
    x = [box_muller() for i in range(10000)]
    hist(x,20)
    savefig("box_muller.png")

box_muller

ProductName パターン認識と機械学習 下 - ベイズ理論による統計的予測
C. M. ビショップ
シュプリンガー・ジャパン株式会社 / ¥ 8,190 ()
在庫あり。

Blu Mar Tenの新譜

出るみたい

これは今聴いても良いですな。

カーネル主成分分析

普通のPCAとちがうのは非線形主成分を扱えること。でも、カーネル関数をうまく選ばないといけないことを意味する。

Rとカーネル法を参考に。

> library(kernlab)
> x <- as.matrix(iris[,1:4])
> iris.kpc1 <- kpca(x,kernel="rbfdot",feature=2,kpar=list(sigma=0.1))
> plot(pcv(iris.kpc1),col=as.integer(iris[,5]))

kpca

あータニモトカーネルを使ったkernel PCAとかなんかいいかもしれない。

ProductName パターン認識と機械学習 下 - ベイズ理論による統計的予測
C. M. ビショップ
シュプリンガー・ジャパン株式会社 / ¥ 8,190 ()
在庫あり。

「読む数学」を読んだ

はじめての現代数学の著者による、数学の本。

小難しい数式はほとんど出てこないので寝転がって読めた。

抽象的な感覚というか、もっと化学とか生物よりの言い方をすれば本質を見極めるってことなんだろうけど、そういう感覚が改めて得られるのは気持ちいい。

となると、自分が子供の時にこういう本に出会っておけばなぁなんてことを考えてしまうんだけど、一方でそのころの自分は抽象の美しさとか感じることが出来たのだろうか?なんてことも思ってしまうんだよなぁ。

子供にはそういう子供のころ読んでいればなぁ、、、的な本は用意しておいてやりたいなぁとは思うんだけど。

PRML

PRMLは下巻になるといきなり難易度が上がるのは気のせいだろうか?

MCMCとSVMをやっつけたい今日この頃。

ProductName パターン認識と機械学習 下 - ベイズ理論による統計的予測
C. M. ビショップ
シュプリンガー・ジャパン株式会社 / ¥ 8,190 ()
在庫あり。

王様のカレー

富士宮のビオデリの近くにカレー屋があるらしいのでいってみた。店内はゆったりしていて居心地がよい。さらに店員さんの接客態度が凄く良かった。

1250074766 1250074769

サービスセットのカレー二種とチーズナンのセット。

ナンが美味い。あとほうれん草のカレーが気に入った。辛さは4くらいにしてもらったんだけどこれだったら6くらいまであげてもよいかも。カレーのもう一方のはトマトとナスのカレーでこっちはわりと甘め。

1250074762

U隊長の頼んだレディースセットはこの他にマンゴーソースがけのバニラアイスがついてた。カレーもあまり辛くなかったので娘がばくついてた。

1250074773

満足。

栃木帰省

ひいばあさんがひ孫の顔をみたいと子供みたいなわがままを言うので、急遽帰省することにした。

初日

佐野ラーメンでもと、佐野に立ち寄り、山銀本店へ。

オーソドックスに手打ちラーメンを注文。

佐野ラーメンってこしっていうよりは平べったい縮れ麺なのね。ニュルッとしててワンタンみたい。これはこれでありかも。

1249988035 1249987993

ラーメン食べたら、車をちょっと走らせて足利学校でも見学。多分初めて。儒教らしい。

1249988014 1249988049

帰りに娘のばあさんと宇都宮の東武で待ち合わせて、総菜を買ってかえるが、7時閉店なのにはびびった。

B1Fの和食の総菜屋で買ったトマトの梅酒煮が美味。あとメークイーンのそぼろがけもなかなか

1249987985 1249988017

同じく、東武で買った東力士のどすこいなんとかという銘柄。

1249988021

親には「東力士は甘くてうまくねーっていっておいてなんでまた買うんだ?」と散々いわれたが、仙禽売ってないんだもんしょうがないじゃんみたいな。

澤姫とか鳳凰美田でも良かったんだけどそういう銘柄すらないのね。あそこのデパート。

結局仙禽は飲めずじまいだった (大那は買ってかえったけど)。

二日目

壬生町にあるおもちゃの博物館とかいうところに行った。

大きなゾイドがある。

1249988000

中に入って尾っぽから出られる。

1249988052 1249988045

昼は併設されているユウガオ館という蕎麦屋さんへ(JAがやってる感じ)

1249988038

天ざる。蕎麦は悪くはないけどちょっとやわかったかも。

1249987996 1249988007

トラクターの路面電車も乗せてくれる。

1249988042

結構楽しんだ。

三日目

高瀬観光やなへ。

鮎は刺身よりも塩焼きのほうが美味い。

1249988010 1249988003

鮎こくもなかなか。

1249988028

最終日は十一正宗の純米酒を。地元の材料で醸したから純矢板らしい。もうちょっとコクがあってもよいかなと思った。なんちゅーか薄く感じた。

1249988024

実家に帰ったら飾ってあったフォトフレーム

1249987989

これはなかなか好評だった。今回撮った写真もSDカードに詰めて送ってやればよいが、フォトフレームもなんかオンラインのストレージに対応して自動で同期するようにならんかなと思った。

JavaでPBPK

PKQuestというものがあるらしい

「仕事するのにオフィスはいらない」を読んだ

クラウドが重要というよりは、クラウドのように抽象化された状態が重要なのではなかろうか。

ProductName 仕事するのにオフィスはいらない (光文社新書)
佐々木俊尚
光文社 / ¥ 798 ()
在庫あり。

会社はジョインしてる感が強くなっていくのかな?成果主義とかそういう指向だよなぁと。

弱い紐帯というか。

ASA-CHANG&巡礼「花」

白隠正宗 蔵元日記で知った。

うお、pearls/sadeじゃねーか。

この曲、内省的で超好きなんだけど、「花」はさらに内省度があがってる。リリックが非常に素晴らしい、というか堕ちまくる。

ProductName
ASA-CHANG&巡礼
ミュージックマインID / ¥ 1,890 (2001-03-28)
在庫あり。

ProductName Love Deluxe
Sade
Sony Mid-Price / ¥ 817 (2000-11-23)
通常7~12日以内に発送

あー、lifeをどう描くかだよなぁ。