DDEを使ってExcelに出力するときに、特殊文字を出力する方法が無いか調べてみました。 & (&)等はコードを使わなくても普通に出力出来ますが、– (–)等はうまく出力出来ません。 文字をSASのエディタにコピーすると "?" に化けてしまいます(画像の16…

Madrid themeでオプションを比較してみよう。 Madridデフォルト ナビゲーションバーがない、フッター情報が入っている。 \useoutertheme[subsection=false]{smoothbars} \setbeamertemplate{footline}[page number] ナビゲーションバーが欲しいので追加する…

大和総研発表の「法則」ドリカム人気→株価上昇？@フジサンケイ ビジネスアイ(魚拓) ネタなのか本気なのか判断が難しいところだ。 データ元がビデオリサーチという時点で怪しいもんだが・・・ サザエさんだけを切り出している所が理解不能だし、ドリカムは偶…

Adobe Acrobatです、Adobe Readerでも多分同じだと思います。 MS PowerPointで言うスライドショーの開始。→Ctrl + L MS PowerPointで言うプレゼン中のページ番号指定移動。→Ctrl + Shift + N

はてなダイアリー市民になってました。

書式 PERCENTILE(配列, 率) 配列: 相対的な位置を決定するため、数値データを含む配列またはセル範囲を指定します。 率: 0 〜 1 の範囲で、目的の百分位の値を指定します。 解説 ・配列にデータが含まれていないとき、または配列のデータ数が 8,191 個を超え…

minipage等で段組みした部分にblock環境を使うと、幅が調整されずに横に突き抜けてしまいます。 LaTeX-Community.org - Beamer: Block width? ではvarblock環境を新しく定義することで、問題を解消しています。便利です。 \begin{frame}[t] \frametitle{varb…

この発想はアートだと思う。 確率変数 $X_1, X_2, \ldots, X_i, \ldots, X_n$ が互いに独立に $\mathrm{N}(\mu, \s^2)$ に従うとき、統計量 \[ T=\frac{\bar{X}-\mu}{\sqrt{\frac{U^2}{n}}} \] は自由度 $n-1$ の $t$ 分布に従う。 ただし、 \[ \bar{X}=\fra…

セクション等をPDFのしおりには出力したいが、(ヘッダーの)ナビゲーションバーには出力したくない場合の書き方を調べました。 \documentclass[dvipdfm]{beamer} \usetheme{Madrid} \useoutertheme[subsection=false]{smoothbars} \setbeamertemplate{navigat…

SASスクリプトをバッチ実行したときや、実行元のSASスクリプトファイルがあるディレクトリにファイル出力したり、グラフ出力したりするために作った物です。 多分2006年頃に作ったものですが、書く事もなかったので挙げなおしました。 マクロ本体 %macro set…

通常フォントファミリーで指定してしまうので、フォントを混在させることはあまりありません。 mtext関数を駆使する事で、複数のフォントを使ってグラフを作成することができます。 postscript出力では未チェックですが、Windows環境でemf出力するか、ビット…

expression()関数中で書体の切り替えができたらしい、何で今まで気づかなかったんだろう・・・。 plain(x) draw x in normal font bold(x) draw x in bold font italic(x) draw x in italic font bolditalic(x) draw x in bolditalic font 他の基本的な数学…

あるパッケージを入れてフォントを変更してみたところ、本文のフォントは変更されたものの、数式のフォントは変更されませんでした。 数式全体のフォントごと変更してしまってもいいですが、記号がずれたり、ちょっとデザインが凝りすぎなのもマッチしません…

Windwosバッチファイルで相対パスって使えないんですね。 ディレクトリ構成 .\hoge.bat .\dir1\ .\dir2\ バッチファイル cd dir1 set D1=%CD% cd ..\dir2 set D2=%CD% for /F "delims= " %%i in ('dir /s /b %D1%\*.txt') do ( echo %D2%\%%i ) 'dir /s /b %…

理論的な説明はとりあえずすっ飛ばして、 今回は、果たしてどういう計算をしているのか、という視点でまとめます。今回使うのは [1] に書かれている有名な方法で、一部では min P とも呼ばれるらしいです。 リサンプリングによって全ての比較が帰無仮説の時 …

この本図書館になかったし、近々購入したい。 と思ってamazon.co.jpで検索したら17000円か！高いな！ amazon.comは$122.98、ちょっと悩む。→注文した(09/06/22)→届いた！(09/06/27)→面白い！役に立つ！(その後) Westfall PH, Young SS. Resampling–Based Mul…

正しいのかどうか分かりませんが、一応できたのでメモっておきます。 URW Garamondを使ってみよう。 まず、Rをインストールしたディレクトリの中の /R/R-X.X.X/library/grDevices/afm に使いたいフォントのafmファイルを一つずつgzファイルにしてからコピー…

ネットワークアルゴリズム。 日本語の解説ってあったっけ？ 参考文献 [1] Mehta CR, Patel NR. A network algorithm for performing Fisher's Exact Test in r × c contingency tables. Journal of the American Statistical Association 1983; 78(382): 427…

今新しくページを作成中。 共通メニュー+コンテンツ(ブログ的な)形式にしようとし、手でhtmlとcssを書いた。 最初はまとめて作っていたが、メニューを更新すると全部書き換えるのがものすごく面倒な事に気づいて、分離して作って適当にバッチファイルで結合…

スコア信頼区間の求め方を復習。積二項分布モデルの尤度から、リスク比 $\phi=p_1/p_2$ のスコア信頼区間を求める。 スコア関数と Fisher 情報行列から $\mathbf{S}^t(\phi)\mathbf{I}^{-1}\mathbf{S}(\phi)=z^2$ をみたす $\phi_{L}$ と $\phi_{U}$ を求め…

TeX Wikiを参考にしてバッチファイルを書いてみました。 Ghost Script/binとTeX/binにパスを通しておきます。 for %%i in ("*.eps") do ( gswin32c -q -dSAFER -dNOPAUSE -dBATCH -sDEVICE=pngalpha -dEPSCrop -r600 -dTextAlphaBits=4 -dGraphicsAlphaBits=…

共通オッズ比の検定の帰無仮説は、層の数を $k$ として第 $i$ 層のオッズ比を $\theta_i$ とすると、 \[ H_0: \psi_1=\psi_2=\ldots=\psi_k=\psi=1 \] \[ H_1: \mbox{not } H_0 \] 以下である。 これは、オッズ比均質性の検定よりは厳しい帰無仮説になってま…

数学セミナー(81年10月号)を古本屋で買いました。 大学でニセ数学しか勉強していなかった私にとっては非常に分かりやすかった。 何で必要になったのかということと、何が出来るようになるのかも厳密にではないが理解できた。 後は確率にどう繋がっていくのか…

共通オッズ比の検定 (test for common odds ratio) は Mantel-Haenszel test を中心に、入門書でも結構詳細に解説してあるものが多いと思う。 Mantel-Haenszel test の原典[1]は30ページ近くあり、かなり読むのが大変。純粋に統計の論文ではなく、後ろ向き研…

有用そうなマクロが公開されています。 ライセンスもGPLだそうです。 survival analysisはあると助かるかも？Mayo Clinic. Locally Written SAS Macros. (cited Apr. 6, 2009) http://mayoresearch.mayo.edu/mayo/research/biostat/sasmacros.cfm

barplotでやると盛大にはみ出ます・・・ Ety <- function(y) { n <- 20 return( (y+1)/(n+2) ) } par(mar=c(3.2, 3.2, 1.5, 0.5), mgp=c(2.0, 0.7, 0)) par(bty="l") names <- replace(seq(0, 20), which(seq(0, 20) %% 5 != 0), NA) barplot( Ety(seq(0,20)…

一楽 重雄. 集合と位相そのまま使える答えの書き方. 講談社 2001. Amazonマーケットプレイスで購入、400+340=740円。 値段も内容も非常に良い買い物でした。

plot()→abline()などで参照線を引くと、参照線のはずなのにグラフの線より前に表示されてしまう。 なので、panel.firstオプションを使うと、参照線→グラフ線の順に描画してくれる。 xx <- runif(100) yy <- runif(100) par(mar=c(3.2, 3.2, 1.5, 0.5), mgp=c…

listings.styでRのスクリプトで関数を取り込んでいたら、表示がプロポーショナルになってくれない。 RjpWikiの利用例3からコピーして使っていたが、 \lstset{% ...略... columns=[l]{fullflexible}, ...略... } このオプションで、コードの単語間のスペース…

Rのデフォルトだとちょっと軸ラベルが離れすぎているし、余白も多めなので調べてみた。 グラフィックスパラメータの mpg で、軸ラベル・軸目盛・軸線のバランスを整える事ができる様だ。 ヘルプには以下のような感じに記述してあった。 軸からラベル、軸から…