小ネタ Windows 版 SAS の Enhanced Editor では、文字列を選択し "Ctrl + ?" を押すと 選択範囲をコメント化できる。 "Ctrl + Shift + ?" を押すと 選択範囲のコメントを解除できる。

ggplot2 で geom_bar(stat = "identity") を利用して棒グラフを作成するときに、横軸の順序を変える方法をまとめてみます。 特に要約データなどを持っていて、ヒストグラムをプロットしたい時などに利用できるかと思います。 方法としては、データフレームに…

ある本を読んでいたら、箱ひげ図に jitter を付加した散布図を重ねて描いていたので、再現してみようと思い適当に作りました。 散布図のプロットの際に、 近くに点があり重なりそうな場合には、jitter を強めに付加する 特に tie になる場合はちゃんと重なら…

おもしろそうなコードを 教えてもらった。 require(ggplot2) DS <- data.frame(speed=rnorm(10), dist=rnorm(10)) f <- function(DS, x, y) { aes <- eval(substitute(aes(x, y), list(x = substitute(x), y = substitute(y)))) p <- ggplot(DS, aes) + geom…

最近 FREQ procedure をよく使ってます、便利ですよね。 FREQ procedure は便利なのですが、table statement 中の out オプションで分割表の集計結果をデータセットに出力するときに、デフォルトでは度数が 0 のセルを出力してくれません。 0 セルを出力した…

追記 最新バージョン (R commander 用プラグイン) が CRAN から利用できます！ RcmdrPlugin.KMggplot2_0.0-3 is on CRAN now 旧バージョン RcmdrPlugin.KMggplot2_0.0-2.zip: Windows 用バイナリ (Dropbox) RcmdrPlugin.KMggplot2_0.0-2.tar.gz: Linux & Mac…

論文風な Kaplan-Meier plot を書いてみよう で試してみましたが、grid package の viewport を使うと、柔軟な感じにレイアウトを決めることができました。 普通に png や pdf でファイルを保存することもできますが、ggplot2 pacakge の ggsave が便利なの…

R 2.13.0 から OpenMP 用のマクロ SHLIB_OPENMP_CFLAGS SHLIB_OPENMP_CXXFLAGS SHLIB_OPENMP_FCFLAGS SHLIB_OPENMP_FFLAGS が Makevars ファイル中で使えるようになっていたらしい (Writing R Extensions)。パッケージの src ディレクトリ中に OpenMP に対応…

追記 マウス操作のみで作図が可能な R commander 用プラグインが CRAN から利用できます！ RcmdrPlugin.KMggplot2_0.1-0 is on CRAN now - Triad sou. テスト版 require(ggplot2) require(survival) `kmg2.theme_gray` <- function (base_size = 12, base_fa…

ggplot2-0.9.x では修正済です require(ggplot2) t <- c("a", "b", "c") x <- c(1, 2, 3) y <- c(2, 2, 2) d <- data.frame(x, y, t) windowsFonts(serif = "Times New Roman") p <- ggplot(data = d, aes(x = x, y = y, label = t)) + geom_text(size = 16 …

グループ変数が複数セットあり、測定値が 1 セットある様なデータに対して、標本平均、総和、標本分散、偏差平方和をたぶん高速に計算できる関数を作りました。 以前作った 高速 tapply(X, INDEX, var) もどき関数 (グループ変数が1セットで、測定値が複数セ…

R で OpenMP を使って並列演算できるかどうか試してみたのでメモを書いておこうと思います。 まず、Windows 環境では、RTools (MinGW) で OpenMP が使えなさそう (Re: Missing libpthread in RTools) でした。 しかし、R devel - gomp library with Rtools21…

Firefox Add-on の論文検索リンク 1.0を復活させよう で修正したものに加えて xpi の定義ファイルのバージョンの部分を修正しただけですが、一応おいておきます。 ダウンロード (Dropbox) 地震でメイン PC の HDD が死んで、復旧作業中。 私はこの程度で済ん…

SGCライブラリ 18 臨時別冊・数理科学2002年11月「独立成分分析」〜多変量データ解析の新しい方法〜サイエンス社 2002.甘利俊一, 村田昇.詳細を見る 独立成分分析の本だが、推定関数の記述が非常に良い。 中身は、一夫と次郎と花子の3人の会話で内容を解説す…

最近お勉強中のものをまとめてみよう。 推定関数、推定方程式 セミパラメトリックモデル Model miss-specification KLD Missing data 因果推論 情報幾何 数値解析 今までやってきたことが、いかに貧弱で、低レベルだったかということを再認識した。

ods pdf のデフォルト出力は解像度が残念すぎて、とても使う気にはなれません。 という事で、解像度の設定について調べてみたら、dpi というオプションがありました。 ods results off; ods listing close; ods pdf file = "&Path.filename.pdf" dpi = 600 s…

mixi にログインして、html ファイルから一部を抜き出すソフトを作ってみた。 ほとんど書いたこと無かったけど Java ってラクだねｗ import java.awt.BorderLayout; import java.awt.Container; import java.awt.Frame; import java.awt.event.ActionEvent; …

Firefox 4 用の更新版 Firefox 4 で論文検索リンク 1.0 Firefox 3 用 グーグルの検索画面の上に、論文検索 (Google Scholar) へのリンクを追加するアドオン 論文検索リンク 1.0 が、Google ショッピングが追加されたときに使えなくなってしまいました。 この…

数式オートコレクトのオプションが外れていると、linear format で入力した後にスペースで変換できない。

SAS で ods output するときに解析結果が大量あると、Results のツリー出力に時間がかかることがあります。 Output 画面への出力制御は ods exclude statement で設定出来ますが、ツリー出力の場合は ods results statement で設定できます。 data sim; do r…

フリーで使える読み上げソフト ReadPlease 2003 を教えていただきました。 国際学会発表の準備に使っておられるそうなので、僕も使ってみようと思います。

togetter Rユーザー会2010 初日 Professor Hadley Wickham の発表資料 Rユーザー会2010 二日目 Professor Uwe Ligges の発表資料 Rユーザー会2010 三日目 鈴木 了太 (株式会社ef-prime) 「Rによるデータ解析ツールの新提案」 R Analytic Flow。 SAS EG みた…

いくつかの PC に SAS 9.2 をインストールしてあるのですが、そのうち 1 台だけ起動すると、 Enhanced Editor Unable to load managed editor component. Error: 0x80040154. The editor may not be installed properly.OLE: Out of memory.このようなエラー…

SAS Enterprise Guide 4.2 を Windows 7 に入れてみたところ SAS Enterprise Guide にエラーが発生しました。 作業中のファイルを保存して SAS Enterprise Guide を再起動してください。という、まるでだめなエラーメッセージを吐いて終了してしまった。 詳…

以下を入手した。 LaTeX Beamer class 3.10 A (CTAN から) pgf 2.00 xcolor hyperref oberdiek /beamer/base/art で ebb する。 hyperref, xcolor, oberdiek を platex *.ins する。 とりあえず動くようになった。 公式サイトのファイルでエラーが出ることが…

追記 最新バージョン(R commander 用プラグイン) が CRAN から利用できます！ RcmdrPlugin.KMggplot2_0.1-0 is on CRAN now - Triad sou. 旧バージョン RcmdrPlugin.KMggplot2_0.0-2.zip: Windows 用バイナリ (Dropbox) RcmdrPlugin.KMggplot2_0.0-2.tar.gz:…

Windows7 での動作確認のついでに、日本語切り替えを追加してみた。

リリース板 最新バージョン(R commander 用プラグイン) が CRAN から利用できます！ RcmdrPlugin.KMggplot2_0.1-0 is on CRAN now - Triad sou. RcmdrPlugin.KMggplot2_0.0-3 is on CRAN now - Triad sou. 試作品 ggplot2 を使った RcmdrPlugin を試作してみ…

2010年度 統計数理研究所共同研究集会 「データ解析環境Rの整備と利用」開催のお知らせ 2010年度の統計数理研究所共同研究集会「データ解析環境Rの整備と利用」 を、下記の日程で開催いたします。興味のある方は、ぜひご参加ください。 日程：2010年11月25日…

基礎的すぎるから知らなかったのか confint(ModelName, level = 1-alpha) という関数があったらしい、プロファイル尤度に基づく信頼区間を計算できる (2012/6/19 記述を修正)。 d <- c( 1,1,1,1,0,0,0,0,0,0, 1,1,1,1,1,1,1,1,0,0, 1,0,0,0,0,0,0,0,0,0, 1,1…