McCullagh and Nelder の GLM が 2 日前にとどいた！ 送料込みで $79.10 だった、結構安いと思う。 入門としては Dobson の An introduction ... を読むといいらしいが、僕は読んでない。多分 SAS の Help で勉強したと思う。 2. outline, 7. conditional li…

set.seed(314159265) height <- rnorm(30, 170, 7) factor(cut(height, breaks = quantile(height, c(0, 0.25, 0.50, 0.75, 1)), include.lowest = TRUE)) [1] (175,181] (167,172] (167,172] (167,172] (167,172] [158,167] (175,181] (175,181] (167,172] …

ggplot2-0.8.7 で、テーマ theme_gray() を日本語向けに直してみよう。 数ヶ月前に ggplot2 でフォントファミリー・サイズを変更できるようなテーマを作りました。 ggplot2 は当然英字出力で調整しているはずなので、日本語を含む場合はファミリーやサイズに…

吉田朋広．数理統計学．朝倉書店，2006． 少し図書館を散歩していたらこの本を発見。 残念ながら目次の詳細が Web 上で見付からなかったのですが、漸近展開の辺りの記述が参考になると思いました。 購入予定です。

m3 の海外ジャーナル日本語サマリー NEJM, Lancet, JAMA, Ann Intern Med, BMJ の新着論文リスト。 日本語に要約した文章が 2 行付いてきてとても便利。 登録 (医療従事者のみ) するともう少し詳細な要約が読める。試しにやってみたところ、条件には合致して…

分散のオンラインアルゴリズム (逐次更新計算)

モンテカルロ法で積分値を推定する問題では、いくつかの分散減少法がしられています。 そのうちの一つである、antithetic variates method (負相関変量法・対称変量法) についてまとめてみようと思います。では、積分 \[ \int_a^b f(x) \mathrm{d}x= \int_0^…

今回も \[ \int_a^b f(x) \mathrm{d}x= \theta \] の様な積分値を推定する問題を考えます。このままでは $[0,\, 1]$ の一様乱数を利用できないので、 \[ x=a+(b-a)r \] \[ \frac{\mathrm{d}x}{\mathrm{d}r}=b-a \] \[ 0\leq r \leq 1 \] の置換積分を考えま…

単純モンテカルロ積分なんていう基礎的な事を復習してみました。自由度 $\nu=1$、分散共分散行列が \[ \mathbf{R}=\begin{pmatrix}1.0&0.7\\0.7&1.0\end{pmatrix} \] の 2 変量 $t$ 分布について、 \[ \int_{a_1}^{b_1}\int_{a_2}^{b_2}\mathrm{MVT}_2(t_1,\…

巷では噂になっている、Amazon.com を装ったフィッシングメールが届いた。 HTML メールで、リンク先のドメインが caveidea.ru / superkub.com だった。 Thunderbird 3.1 は感知してくれました。

十進数を二進数に変換する演算が必要になったので調べてみました。 wle package に binary() 関数というものがあり、そいつを少し弄ればよいことが分かった。 それだけではちょっとつまらなかったので、再帰を使って自分でも書いてみた。 # 再帰で書いた関数…

以前は単なる対応データの割合が一致するかどうかの検定ぐらいの認識だったので、ちゃんと勉強してみようと思いました。 でも McNemar の原典 [3] を読むとさっぱり分からなかったので Kbauth 1985 [2] なども読みました。McNemar's test は、 \[ H_0:\, \pi…

お役立ちリンク情報。NCI The Biometric Research BranchReprints & Presentations から様々な文献がダウンロードできる。本もダウンロードできる、という情報を後輩からもらったので紹介。 Statistics SurvaysASA, IMS, Bernoulli, Statistical Society of …

正規線型モデルとは書いたのですが、結局の所 $\mu$ なのか $X\beta$ なのかの違いだけなので、確率モデル \[ Y_{i} \overset{\mathrm{i.i.d.}}{\sim} \mathrm{N}(\mu, \s^2); i=1, \ldots, n, \] を考えます。この確率モデルのもとで、$\mathbf{Y}$ の同時…

TeX の PDF viewer を Sumatra PDF viewer に変えてみました。 PDF ファイルをロックしない仕様で、ファイルが更新されたら読み直してくれる、dviout ライクに使えてとっても便利になりました。

Beamer & pdfpagelabels warning (LaTeX Community) \documentclass[hyperref={pdfpagelabels=false}]{beamer} ドキュメントクラスのオプションに入れておくらしい。 気になっちゃう人向け。 関連記事 Beamer 3.10 にアップグレードした Beamer 3.10 の場合…

最近勉強していると、必ずこの問題にぶつかるので、今のうちに整理しておこうと思う。 ここで、統計モデル \[ f(\mathbf{Y} \mid \boldsymbol{\theta}, \boldsymbol{\phi}), \] を考えよう ($\mathbf{Y}, \boldsymbol{\theta}, \boldsymbol{\phi}$ は全てベ…

印刷したら薄くなってしまった。 bibtex などをちゃんと使っている時は、エントリーの種類毎に \setbeamerXXX できる。 bibitem べた打ちの場合は {bibliography entry author} が \setbeamerXXX の対象になるみたい。 \setbeamercolor{bibliography entry a…

なんか悔しかったので、ループを使わない方法を考えてみた。 @echo off call set CPATH=%CD% cd ../ call set UPATH=%CD% call set CDN=%%CPATH:%UPATH%=%% call set CDN=%CDN:~1% cd %CPATH% echo %CDN% PAUSE ルートディレクトリだった場合とかは考えてい…

プロットの作製 基本プロットを作る Geoms Aesthetics 違う種類のグラフを重ねる 参照線の追加 グループ分け 層別プロット スケールと軸 Scales 軸ラベルやタイトルの変更 軸の表示範囲を変更する 軸の左右の余白を削除する 軸表示の修飾 日時の軸スケール …

古い記事を移動しました。 基本プロットを作る 基本的な手順は、 ggplot() にデータフレームと各軸や層に対応する変数名を指定する 書きたいグラフに対応する geom_XXX() を足す の二つが分かっていれば OK です。 これに加えて、統計処理をしてくれる stat_…

Usage Note 32464: Registering file types via the SAS Filetype Registration Utility for SAS 9.2 Phase 1 and SAS Enterprise Guide を使う。 上の記事内の [1. Navigate to Index of /techsup/download/pc.] のリンクからユーティリティーをダウンロー…

gdipp ってなに 最近話題なので入れてみました。 gdipp - Project Hosting on Google Code 窓の杜 - インストールするだけであらゆるソフトのテキストを美しく滑らかに描画「gdipp」 gdipp_x86_0.7.6.msi を program loader でインストール。 MS Wordだと こ…

%macro DELSHEET(rc); %local vlen nlen k var; data _null_; rc = "&rc."; s = find(rc, "R"); e = find(rc, "C"); r1 = substr(rc, s+1, e-s-1); s = e; e = find(rc, ":"); c1 = substr(rc, s+1, e-s-1); rc = substr(rc, e+1, length(rc)-e); s = find(r…

Hide/Show Google left-hand menu 1.0.4 右なら良かったのにねぇ。 ↓ 慣れてきた。 旧検索にも戻せるらしい。 http://www.google.co.jp/webhp?hl=all 消えちゃった？

IARC Publications - PDFs online 個人で非商用の研究に利用する場合は OK だそうです。 Volume I はペーパーバック版を持っています。内容は最新のものではないですが、case-control study に関連した統計手法について、恐ろしいほど多彩な情報を網羅してい…

[1] の文献をかなり参考にさせていただきました。 少し修正して、配列を演算に使えるようにしています。 でも、かなり小標本じゃないと使えないですよねぇ。 shuffer <- function(n, m) { m <- as.integer(m) n <- as.integer(n) if (m > n) { tmp <- m; m <…

R の survival package には、条件付きロジスティック回帰分析を行うための clogit() 関数なるものがあります。 最近この関数を使っていて、ちょっと気づいた点がありましたので、メモを残しておこうと思います。 method="approximate" は危ない clogit() 関…

scale_colour_manual を使う事で、ggplot2 で用いられる色のパターンを自分で設定出来るらしいです。 library(ggplot2) colours <- c("#D92121", "#9999FF", "#D92121", "#21D921", "#FFFF4D", "#FF9326") d <- data.frame( ID = floor(seq(1, 6 + 7/8, by =…

高速 tapply(X, INDEX, var) もどき関数 を作ったときに、簡単な実装を考えていたのでまとめてみます。 xxx.test <- function(d, g) { n <- as.numeric(table(g)) index <- as.numeric(factor(g)) mindex <- NULL for(i in 1:length(unique(index))) { minde…