平岡 和幸, 堀 玄. プログラミングのための線形代数. オーム社 2004. 自分が読んだ線型代数の本でピカイチ。 これを読んで、行列の変換としての意味が完璧に理解できた。 学部1年では東京大学出版会の線型代数入門 (緑のやつ, 解析学は青) が使われていた。 …

パッケージをもりもり作っていたら、 * checking Rd cross-references ...が表示されたまま止まってしまう。 自分の環境だと何故か Sys.setenv( ) 関数を .Rprofile に書いておいても働いてくれない、これも不思議。 しょうがないので check.pl の240行目ぐ…

ご要望があったので、機能追加版を作ってみました。 プロットに用いたデータを、グラフを描画したシートに移動するように変更 箱の幅を指定できるように変更 箱の幅指定に関しては、Excel 2007の致命的なバグで、書式設定画面から幅を指定できないため、VBA…

Rcmdr package自体のインストールは、car packageを追加した上で、「ローカルから...」を選べばよいので簡単でした。 しかし、Rcmdrの初回起動時に、依存関係があるパッケージをディレクトリを指定して一括でインストールする画面がでます。 ここで、 must i…

せっかくなので触ってみました。 data logistic; call streaminit(88745623); a = 0; b = 1; do i = 1 to 200; do x = 0 to 1; logit = a +b * x; p = exp(logit) / (1 + exp(logit)); y = rand("Bernoulli", p); output; end; end; keep x y; run; ods list…

意外に色々あるんですね。 ods listing close; ods pdf file = "c:\out.pdf" style = Journal; ods graphics on; proc xxxx; run; ods graphics off; ods pdf close; ods listing; その他使いそうな形式。 ods csvall file = "c:\out.csv"; ods rtf file = "…

09/11/29版「W32TeX」のdvipdfmxの挙動について (TeX Q&A) 私もほぼ同じ症状が出てしまいました。 数式の一部と、文章中の斜体 (要は斜体の一部?) が歯抜けになって表示できない状態でした。 斜体のアルファベットとギリシャ文字は歯抜けして、subscriptの大…

関数を作ってみました。 HS.tapply <- function(x, index, type) { # create index indicate matrix index <- as.numeric(factor(index)) mindex <- NULL for(i in 1:length(unique(index))) { mindex <- cbind(mindex, index-(i-1)) } if (type == 5 || typ…

うーん、こりゃ参ったねｗ n <- 20 rep <- 100000 d <- matrix(rnorm(n*rep), nrow=n, ncol=rep ) # apply Rprof() varA <- apply(d, 2, var) Rprof(NULL) print(a <- summaryRprof()$sampling.time) # matrix operation (1) Rprof() varB <- NULL for(i in …

前回書いた結果をもとにして、ちゃんと近似できているのか試してみることにしました。 今回は Welch (1938) の (2) 式の統計量 \[ v=\frac{\bar{X}_1-\bar{X}_2}{\sqrt{\frac{V_1}{n_1}+\frac{V_2}{n_2}}} \] を導出する過程で出てくる分布の近似を考えてみ…

[1] の説明が分からなかったので、考えたり調べたりしていたら [2] を見つけることが出来た。 ある尤度関数 $p(y \mid \theta)$ に対して、無情報事前分布を構成する事を考える。 どんなパラメータ $h(\theta)$ に対して一様な事前分布を考えるかによって、$…

Ntzoufras I. Bayesian Modeling Using WinBUGS (Wiley Series in Computational Statistics). New York: John Wiley & Sons, Inc. 2009 ぱらぱらっと読ませてもらいました。 うーん、欲しいけど高い。

Biometrics Bulletin 1946 [2] を最初に読んだけど、詳しい話が書いてなかったので Psychometrika 1941 [1] を読んだ。概要は、不偏分散の推定量 \[ V_i=\frac{1}{\nu_i}\sum_{j=1}^{n_i}(Y_{ij}-\bar{Y}_i)^2 \] の和の分布および、カイ二乗分布に従う確率…

シミュレーション中に進捗状況を監視。 path <- "C:/"; set.seed(091021) r <- 20 m <- c(0, 0, 0.3) s <- c(.1, .1, .1) n <- c(10, 10, 10) d <- do.call("rbind", mapply(function(n, m, s) matrix(rnorm(r * n, mean = m, sd = s), nrow=n, ncol=r), n, …

SkyDrive を導入してみました。 ついでに昔公開していたファイルを移動しました。 依頼があったらExcelマクロとか作ってアップロードすると思います。

松原望. 入門ベイズ統計―意思決定の理論と発展. 東京図書 2008. 中古で買ってみたら結構新しい本でした。 難易度としてはさほど高くないと思います。珍しくタイトルに入門がついていて、ホントに入門っぽい本ですね。 カルマン・フィルターの話は結構面白か…

はじめてインターネットにアクセスしてから10年ぐらい。 もう記憶は定かではないけど、一番古い記録がこれ？

線型モデル $\mathbf{Y}=\mathbf{X}\mathbf{\beta}+\boldsymbol{\epsilon}$ において、不等分散下 \[ \mathrm{E}\left(\mathbf{Y}\right)= \mathbf{X}\mathbf{\beta}= \begin{pmatrix} x_{11} & x_{12} & \cdots & x_{1m} \\ x_{21} & x_{22} & \cdots & x_{…

Windows Update「一部の更新をインストールできませんでした」 タスクスケジューラを無効にしていると失敗するらしい。 Info: Office Genuine Advantage Notifications (KB949810) failure when Task Scheduler Not Running

SAS/Graphのannotate機能のBAR関数を使います。 proc datasets lib = work kill; run; option linesize = 130 pagesize = 500 mprint; dm 'log; clear; output; clear'; %let execpath = " "; %let Path = " "; %macro setexecpath; %let execpath = %sysfun…

野田 一雄, 宮岡 悦良. 数理統計学の基礎. 共立出版. 一番好き、非常に丁寧に書かれている上に、カバーしてくれる範囲が広い。 確率の話はあっさりめに書いてあって、収束の話、確率モデル、検定、推定などを中心に勉強するのに適していると思います。 内容…

Rでシミュレーションを行おうと思って、「群変数が共通で、データを乱数として r セット生成する」スクリプトをループを使わずに書こうと思ったら、なかなかうまくいかない。 ↑を逆(「データが共通で、群変数を乱数として r セット生成する」)にするとリサン…

Kotz S, Nadarajah S. Multivariate T-Distributions and Their Applications. Cambridge University Press, 2004. 未知パラメータを消すために、最近買いましたｗ

インストーラが9.1.3と全く違いました。 Windowsで統計解析ソフトウェアとして使うだけなら9.1.3の方が使いやすかった。 サーバ型ソフトウェアとか、他の使い方なら使いやすいんだろうか？ SAS Deployment Wizardの使い方がよく分からず(インストーラみたい…

デフォルトではリスト記号に画像が埋め込まれるので、番号付きリストで入力したい。 "\beamertemplatetextbibitems"というオプションを入れると、番号付きリストで出力してくれる。 \begin{frame}[t,allowframebreaks, allowdisplaybreaks] \frametitle{Refe…

私事になるが、SAS 9.2が導入された(まだインストールしていないけど)。 グラフ関係とベイズ関係の procedure は早めに触っておこうと思う。

OSを入れ直したPCにThunderbirdを入れたところ、メール画面の右側に表示できるカレンダー画面にToDoが表示されなくなってしまった。 何でだろう、と真剣に悩んだ。 むちゃくちゃ簡単だった。 念のためですが、ToDoを表示するためには、Lightningが必要です。…

検索してきた方がいたのでメモ。 proc datasets lib = work kill; run; option linesize = 130 pagesize = 9999 mprint; dm 'log; clear; output; clear'; %let execpath = " "; %let Path = " "; %macro setexecpath; %let execpath = %sysfunc(getoption(s…

自民に入れました、増税嫌だな。 今回の選挙で個人的に楽しみだったことは、朝日RDD[1]の精度を見たかった所ぐらいかな。 選挙結果と朝日RDD調査による予測 (全体) 政党名 選挙結果 朝日RDD (全体) 推定値 信頼区間？* 民主308321[307, 330] 自民119103[ 89,…

便利かなと思って作成してみました。 data palet; input color depth code $ @@; cards; 1 1 cxA60000 1 2 cxBF3030 1 3 cxFF0000 1 4 cxFF4040 1 5 cxFF7373 2 1 cx85004B 2 2 cx992667 2 3 cxCD0074 2 4 cxE6399B 2 5 cxE667AF 3 1 cx48036F 3 2 cx602580 …