Bernardo JM, Smith AFM. Bayesian Theory (Wiley Series in Probability and Statistics). John Wiley & Sons Inc, 2000. 図書館で読んだけど、理論的な説明が非常に詳細で、これ全部ちゃんと読んだら相当効きそう。 応用(MCMCとか計算)に関してはあまり書…

交換可能性を満たす場合、$n$ 個の確率変数 $X_i$ ($i=1, 2, \ldots, n$) の同時確率 (密度) 関数において \[ p(x_1, x_2, \ldots, x_n)= p(x_{\s(1)}, x_{\s(2)}, \ldots, x_{\s(n)}), \] が成り立つ。 ただし、$\s$ は集合 $\{1, 2, \ldots, n\}$ に対する…

衝動買いした、24,000円だった。 竹内 啓 (編集). 統計学辞典. 東洋経済新報社, 1989. ここで検索すると多分安く売ってるよ (2010/02/04)。

○○支持率：メディアによる洗脳効果を計る指標の一つ 電話調査：固定電話を日中に受けられる集団の結果を得るための調査方法 街頭アンケート：昼間に商店街をぶらついているような、非常に偏った集団の結果を得るための調査方法 日本のジャーナリズム：情報を…

Word 2007の新しい数式の入力方法を調べた。 WordやPowerPointは環境でズレまくり&文字化けしまくりなので、Excelしか使わない事にしてるけど、Wordで書けといわれちゃーしょうがない。 MathTypeとか数式エディタだとどうにも行列が綺麗に入力できず、2007な…

調べものをするついでに、図書館で借りてきて読んだ。 高橋 一. 経済学とファイナンスのための数学 (新経済学ライブラリ). 新世社 1999. 経済学もファイナンスもサッパリ分からないけど、自分が使いそうな領域を含んでいたので買うことにした。 中身は微積分…

Rでもバッチ実行ディレクトリにファイルを保存したかったので調べてみた。 Rscriptを使うと出来るようだ。 Rscript [--options] [-e expr] inputfile [args]のように、argsを引数として渡すことができる。 Rscript -e "commandArgs(); proc.time()" arg1 arg…

Windowsにインストールされているフォントを適当に指定してfamilyとして定義すると使えるようになるみたい。 windowsFonts(Times="Times New Roman") par(family="Times") 定義済みのフォントは windowsFonts() で表示される。 n <- 30 y <- 2 Bin <- functi…

不完全ベータ関数の近似はまとめるほど理解できなかったので、また時間があるときにやろう。 今日はラプラス近似が何をしてるかを少しグラフィカルに検討した。 \[ \int_{\mathbf b\in R^q}\exp\{h(\mathbf b)\}d\mathbf b \simeq \exp\{h(\hat{\mathbf b})\…

\[ p(\theta \mid y)=\int_{\theta_1}^{\theta_2}\begin{pmatrix}n \\ y\end{pmatrix}\theta^y(1-\theta)^{n-y}d\theta \] Laplace が使った近似を調べてるんだけど、$\theta^y(1-\theta)^{n-y}$ を $y/n$ の周りで展開して・・・、なんて言う話と、正規分布…

Bayesian inference 勉強中。 prior distribution を一様分布で尤度が二項分布の時の posterior distribution は \[ p(\theta \mid y)=\frac{\displaystyle \int_{\theta_1}^{\theta_2} \begin{pmatrix}n \\ y\end{pmatrix} \theta^y(1-\theta)^{n-y} \mathr…

090217 の記事で挙げた Asymptotic methods in analysis を読んできた、1st ed. だった。 4.2 A general case, 4.6 Multiple integrals にきれいに書いてあった、前提条件も書いてある。 積分近似自体は Laplace のアイデアだから、The Laplace method for i…

積分のラプラス近似についてです。 だいぶ昔に [1] を引用して非線形混合効果モデルについて書いたことがあったのですが、積み残しを消化 (詳細を理解) するためにラプラス近似について調べてみました。ラプラス近似を見てなんじゃこりゃと思ったのが [1] p.…

今まで使ってたmaglogからお引っ越し。 編集中にサーバ負荷で書きかけの記事消滅とかはもう勘弁して欲しいので・・・。 おや、早速重くなってきたな・・・。 インポートできなかったから手作業でコツコツと、はずかしいところを修正しつつ・・・ TeXが使える…

ネタがなかったんで、ログとアウトプットの制御について。 printto procedureを使うと色々制御できます。 /* logとoutputを消す */ filename fdummy dummy; proc printto log = fdummy print = fdummy; run; /* logとoutputをファイル出力 */ proc printto l…

メモが目的なので要点だけ。 dev.copy2eps()だと何故かうまくいかなかったので、postscript()を使った。 postscript( file="c:/test.eps", width = 9, height = 9, horizontal = FALSE, onefile = FALSE, paper = "special", family="AvantGarde" ) family="…

いつの間にかASAから出てる論文がWebでフリーで読めるようになってる(2001年以降のものだけ)。 JASA, The American Statistician辺りはそこそこ読むので取り寄せないで済むから助かる。 後はProject Euclidも外せない。 Annals系は余り読む機会もないが、Ber…

Excel 2007 を導入しました。 マウスで何カ所も設定をクリックするの面倒くさいのは相変わらずでしたが、グラフは意外ときれいになっていました。 論文の査読コメントについて相談を受けたついでにざっと全体を見せてもらったら、毛付き棒グラフを描いていた…

やろうと思っていること 研究 線型モデル 特殊なパラメータ変換をしたときの conditional exact test (積多項分布モデル) リサンプリング・ブートストラップを使った方法 ベイズ流の解析方法、MCMC EMアルゴリズム ネットワークアルゴリズム Rパッケージを作…