R

よい検定、わるい検定

R

この二つの解説は非常に分かりやすい。 http://aoki2.si.gunma-u.ac.jp/lecture/BF/index.html http://aoki2.si.gunma-u.ac.jp/lecture/BF/sankouzu.html 2群はwelchのt検定、3群はoneway.test(ANOVA(aov)はNG)が良いとのこと。なるほど。 一つ気になるのは…

R

http://www.kyoritsu-pub.co.jp/series/arudemanabu.html 仕事の半分をRでデータ解析している身としては、completeすべきか。 (社会科学系はいらない様な気もするが。。。) しかし、このURL、"あるでまなぶ"。センスに脱帽。 機械学習が3冊でオーバーラップ…

Locale

R

今日のはまり所 > lc <- Sys.setlocale("LC_ALL", "C") ; sort(c(" ", "10", " 9", "9")) [1] " " " 9" "10" "9" > lc <- Sys.setlocale("LC_ALL", "en_US.UTF-8") ; sort(c(" ", "10", " 9", "9")) [1] " " "10" "9" " 9"USでは自分の常識が通用しないと思…

PDF2

R

久しぶりにアッブデート。 https://r-forge.r-project.org/projects/pdf2/ といっても、R-2.11.1用に少し修正しただけ。

PDF2

R

R-2.10.0 -> R-2.10.1 grDeviceライブラリの変更は2行のみ、久しぶりに平穏なアップデート。

Rユーザー会

R

普段、交流のない方々と話ができました。 hyperlinkやpopup textだけでなく、javascriptが入れられるはず。これは面白いかも。 cairoもannotationを検討中だけど、実装はまだみたい。http://lists.cairographics.org/archives/cairo/2009-February/016620.ht…

Rユーザー会

R

Rユーザー会でPDF2の紹介をする予定。

ふたつのR

最近こんな記事がありました。 Rubyはpseudo codeを実行出来る言語だと書かれてあるけど、まさにその通りで、 実際このコードは、わずかな修正でRubyで動いてる。 (もちろんRGLのおかげなんだけれども。。。) Rも負けてない。集合やベクトル、行列演算のアル…

pdf2

R

しばらく放っておいたPDF2のアップデート。 http://r-forge.r-project.org/projects/pdf2 (ビルドエラーがなければ、明日ぐらいにはパッケージがダウンロード出来るはず) R-2.8.1への対応のみ、機能追加なし。 自分が使っていないとモチベーションが上げられ…

RSRuby

どうやら、パッチが取り込まれたらしい。 http://github.com/alexgutteridge/rsruby/tree/master ruby 1.9.1への対応だけでなく、shutdownメソッドのfixも入ったので良かった。 rubyfogeのgemも更新されているみたい。 アップデートの際は、下記のようにオプ…

RSRubyのruby 1.9.1への対応

RSRubyがruby1.9.1で動かないのでバッチを作った。 http://github.com/tadakado/rsruby/tree/master 以下はgithubを使うための手順 ssh-keygen -t rsa cat ~/.ssh/id_rsa.pub githubにユーザ登録(id_rsa.pubを登録する) rsrubyのサイトでforkをクリック mkdi…

pdf2

R

一部のコアディベロッパには、パッケージングの方法に関して印象が悪いが、、、 なんとかWindows版もリリースのめどが立った。手元のPCでのビルドは成功。 あとはR-Forgeのマシンでビルド出来るかどうかだけ。 T君にはいろいろ手伝ってもらったのだけど、結…

pdf2

R

R教授に指摘を受けたので、修正版をアップロードした。 眠い。。。

pdf2

R

β版をリリースした。 http://r-forge.r-project.org/projects/pdf2/

A patch for extending pdf device to embed popup text and web links

R

Sample script & patch to R-2.6.2 pdf("sample.pdf") par(oma=c(2,2,2,2)) plot(1:10) pdf.link.on.box(1, 5, 3, 8, url="http://www.r-project.org") for(i in 1:10) pdf.text.on.box(i-0.1, i-0.1, i+0.1, i+0.1, LETTERS[i], border=c(0,0,0)) text(5, 9…

A patch for extending pdf device to embed popup text and web links

R

patch to Revision: 44977 Last Changed Date: 2008-03-29 diff -r -c R.org/src/library/grDevices/NAMESPACE R/src/library/grDevices/NAMESPACE *** R.org/src/library/grDevices/NAMESPACE 2008-03-29 23:36:12.000000000 +0900 --- R/src/library/grDev…

PDF link patch

R

R-2.6.2で、PDFにアノテーションとして、メモやURLのリンクを埋め込むことができるパッチです。 今のところマニュアルはありません、パッチを適用してコンパイル後、下記のスクリプトを走らせてみて、 作成されたPDFから、どんなことが出来るか確認して下さ…

いいものめっけ(PNG出力)

R

png()を使うとX11の通信が生じる。 数十万点のプロットとかさせると、リモートでの作業はつらい。 これがあれば、画期的に改善されそう。 早速、インストールしよう。http://d.hatena.ne.jp/dancing_infobio/20071026/p2 install.packages(c("GDD")) library…

小技

R

capture.outputでコンソールの出力を保存出来る。 > data(iris) > capture.output(summary(iris), file="output.txt") sink()なんてのもある http://www.okada.jp.org/RWiki/?cmd=read&page=R%BD%D0%CE%CF%A4%CE%B5%AD%CF%BF&word=%BD%D0%CE%CF%C0%E8

1000倍早くする方法

R

同じデータで縦横を入れ替えたファイルでは、1000倍読み込み時間が異なる。 縦長なデータの方が有利。 > system.time(a1 [1] 859.060 4.620 866.269 0.000 0.000 > system.time(a2 [1] 0.610 0.010 0.682 0.000 0.000 > dim(a1) [1] 10 22283 > dim(a2) [1] …

3倍早いファイル読み込み

R

Rのread.table(というかread.delimが好みだけど)は結構遅い。 大量のファイルを処理する必要があって、改良してみた。 単に数値データであることを指定してるだけだが、効果てきめん。 でも、3倍ってのは言い過ぎか。。。 > system.time(m1 system.time(m2 d…