水産資源学におけるベイズ統計の応用ワークショップのスライドの22,23枚目
となっているのだが、23枚目のRのコードは
n2 <-sample(length(x1), 1000, replace=T, prob=dlnorm(x1,1.1,0.6))
とwを求めずにいきなりサンプリングしているので？となったが、一様分布だから重みを計算する必要ないのね。
わざわざ書くならこうか。
x1 <-runif(100000,0,60)
w <- dlnorm(x1,1.1,0.6) / (dunif(x1,0,60))
n2 <-sample(length(x1), 1000, replace=T, prob=w)
x2 <-x1[n2]
plot(density(x2))
18082009 life
あーなるほど。
小選挙区制度には、数々の欠陥が存在する。このような欠陥は高校の授業でだいたい教える。一番の欠陥は死票の多さであろう。しかし、それでも小選挙区制度を採用しているのは、選挙の目的が「政権交代」であり、「合法的な革命」に他ならないからだ。
イジングモデルとかエルゴード性を思い浮かべた。
フルーティーで香りがよい。
地元でしか出回らないらしいんだけど、今回は多めにつくられたようで、こっちのほうでも売ってるらしい。
あと、波瀬正吉はしょうきちじゃなくてまさきちだという小ネタも仕入れたが使い道がないなぁ。
18082009 life
分かりやすかった。良書。
偏微分方程式からテイラー展開、ヤコビアン。 線積分、面積分をやってストークスの定理へ。 最後がラグランジュの未定乗数法。
特に最後のラグランジュの未定乗数法の章が良かった。
18082009 chemoinformatics bioinformatics q-chem
ドッキングシミュレーション関係の本だろうか。
あとGaussian
Gaussianは使ってないのでこっちはあまり興味をそそられないが。
17082009 MCMC
夏休み前にBUGSの具体的な事例が載っていたので買って積んでおいた本を読みはじめた。さっきまでWinBUGSをmacから使えるように頑張っていたけど、JAGSもpymcもあるし、無理して使わんでもまぁええかと諦めた。
統計科学のフロンティアシリーズのやつも薦められたので読んでみたい。
伊庭先生と言えば「ベイズ統計と統計物理」にはお世話になった。
で、この本買ったのが、GIWとなんかの共同開催の時で、基調講演が赤池先生だったのは覚えている。
U隊長がオーブンとドライヤーを一緒に使ったら局所的にブレーカーが落ちて、とばっちりでサーバーが落っこちた。
以前にちょっと色々やっててサーバーの電源をUPSのバッテリ側に差してなかったっていう僕の落ち度もあるんだけど、Estraierのindexが壊れたっぽい。
今日はもうやる気が起きないので、来週なんとかする。
15082009 work
出るまでもないけど、ちょっと把握しておきたい会議ってのは少なからず存在していて、そういうのはイントラのustreamみたいなもので流れてくれると嬉しい。あと、議事録で、言った言わないとか面倒だし。実際、議事録に記載する数行の文章をどうするかでメールが5往復とかあったりするしな。
ニコニコ動画みたいな仕組みでオンライン議事録とか会議録取れるようなものが欲しいなぁと思う。コメントのとこってどのタイミングでどういうこと言ったかメモれるし便利ですよなぁ。さらに時系列でまとまっててみやすいし。したのほうに、wikiっぽい仕組みのものでもくっつけておけばさらによいだろう。
Kamaitachiとか使えばやれるのかなぁ。
パターン認識と機械学習 下 - ベイズ理論による統計的予測
キーポイント多変数の微分積分 (理工系数学のキーポイント (7))
タンパク質計算科学 ―基礎と創薬への応用― [CD-ROM付]
Gaussianプログラムによる量子化学計算マニュアル 計算入力法から実験値との比較まで
マルコフ連鎖モンテカルロ法 (統計ライブラリー)
計算統計 2 マルコフ連鎖モンテカルロ法とその周辺 (統計科学のフロンティア)
岩波講座物理の世界 物理と情報 （3） ベイズ統計と統計物理
OMRON BZ50T UPS(無停電電源装置)