出るみたい

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

普通の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]))

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

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

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

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

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

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

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

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

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

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

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

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

満足。

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

初日

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

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

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

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

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

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

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

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

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

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

二日目

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

大きなゾイドがある。

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

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

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

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

結構楽しんだ。

三日目

高瀬観光やなへ。

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

鮎こくもなかなか。

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

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

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

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

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

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

弱い紐帯というか。

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

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

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

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

PRMLも再読してるが、手動かしながら読むほうがしっくりくる。 p.70のベータ分布の逐次学習ってのが面白そうなのでやってみた。

data <- rbinom(30,1,0.5)

plot_dbeta <- function(data){
a <- 0
b <- 0
x <- seq(0.01, 1.0, len = 500)
for (i in 1:length(data)) {
    if(data[i] == 1) { a = a + 1 } 
    else             { b = b + 1 }
    if(i>3){
    if(i==4){ plot(x, dbeta(x,a,b),type = "l",col=i,ylim=c(0,5)) }
    else    { plot(x, dbeta(x,a,b),type = "l",col=i,xlab="",ylab="",axes=F,ylim=c(0,5)) }
    par(new=T)
    }
}
}

二項分布で30個のリストを発生させて読みながらa,bのデータを更新していきつつその時のベータ分布の密度を求めつつプロットしていくという操作。

> data
 [1] 1 0 1 1 0 0 0 1 0 0 0 1 0 0 1 1 1 0 1 0 0 0 1 0 0 1 0 1 0
[30] 1

データはこんな感じで、4番目のデータからプロットしていくとだんだん二項分布に近づいていく様がみてとれる。