RでBase64のエンコード/デコード
金曜日, 3 月 13th, 2009ログを確認していたら、えらく長いリファラが付いているアクセスを発見しました。文字列から察するに、Base64でエンコードされたもののようだったので、デコードしてみようと思いました。普段ならPerlでワンライナーでも書こ [...]
FizzBuzz問題をRでやってみる
火曜日, 2 月 10th, 2009Slashdotを読んでいて、ふと目に止まったFizzBuzz問題。”R fizzbuzz”等でググれば、当然Rでやっている先人たちのサイトがちらほら見つかります(一例)。以下の図は、7種のコードで100,000回ループをまわした際の実行時間を10回計測し、箱ひげ図にしたものです。ちなみにy軸は対数。
単純に速度を追求するなら、Cなりアセンブラで書け、という話ですが…。
さて、Rのソースは以下の通りです。
# fizzbuzz.r
# 一番シンプルな形
fb0 <- function(len=100) {
a <- b <- 1:len
a[b %% 3 == 0] <- "fizz"
a[b %% 5 == 0] <- "buzz"
a[b %% 15 == 0] <- "fizzbuzz"
a
}
# ワンライナー風、x %% 5を2回評価するのが若干無駄
fb1 <- function(len=100) {
x <- 1:len
z <- ifelse(!x %% 3, ifelse(!x %% 5, "fizzbuzz", "fizz"), ifelse(!x %% 5, "buzz", x))
z
}
# さらに短く、Rらしさのある怪しいコード
fb2 <- function(len=100) {
z [...]
線形SVMの識別関数を求める (ksvmパッケージ)
土曜日, 2 月 7th, 2009サポートベクターマシン(SVM)は、変数に重みを付けるのではなく、サンプルに重みを付けることで識別超平面を求めます。RBFカーネル等を用いた非線形モデルでは識別関数が複雑な形(RBFならいくつかのガウス関数の和)となる [...]
Rでテキストマイニング
木曜日, 2 月 5th, 2009PubMedのアブストラクトで遊ぼうと思い、調査中。
CRAN Task View: Natural Language Processing
徳島大学 石田基広先生のスライド [pdf]
最近は「Rによる [...]
Affymetrix GeneChip 薀蓄
木曜日, 2 月 5th, 2009ググれば他でも出てくるような話ですが、一応書いておきます。
[とりあえずRMAでノーマライズ]
setwd()なりメニューなりから、CELファイル(gz圧縮可)のあるディレクトリに移動して、
> e <- [...]
Rまめ知識
月曜日, 2 月 2nd, 2009随時更新予定です。
“a <-1″は”1 -> a”とも書ける (コマンドラインで重宝)
でかいファイルはgzで圧縮しておき、read.delim(gzfile(”xxx.gz”))などとして読む
で [...]
« Previous Entries