GNU R で、image() を使って3次元データをカラーマップとして可視化するとき、その横に色と値の対応を示したいですよね。R はそのような機能を用意してくれているのかな？ 少し探したけど発見できなかったので、以下のようにしています。もっと良い方法があ…

las オプションで指定する。las=1 で常に水平、las=2 で常に垂直など。 詳しくは help(par) を参照すること。

R に翻訳してみた．そういう気分だった．でも，以下のやりかたじゃない別のもっと良い方法があるような気がする．R の神様降ってこい！ critics <- data.frame(matrix( # LW, SP, JML, SR, YMD, TNL c( 2.5, 3.5, 3.0, 3.5, 2.5, 3.0, # Lisa Rose 3.0, 3.5,…

明日28日から emobile の新しい料金プラン「スーパーライトデータプラン」がスタートするようですね．私は現在，ライトデータプランに契約していて，今月末まで無料キャンペーンの対象になっています．emobile は，プランの変更は自由にできるのでユーザフレ…

R には標準でベキ分布を生成する関数がありませんね．以下の関数を定義すると，Zipf-Mandelbrot 分布 に従う乱数を runif と同じ感覚で生成できます． rzipf <- function (n, s=1, q=0, supp=c(1,n)) { x <- supp[1]:supp[2] p <- 1/(x + q)**s y <- sample(…

lsoda ODE (常微分方程式) の求解系

http://www.tctv.ne.jp/members/mkim/geography/city_population_2001.htm に関するはてブのエントリに次のようなコメントがありました： # 2007年03月05日 caramelly 地理 最新のは→http://uub.jp/rnk/rnk.cgi?T=c&S=jとかhttp://uub.jp/rnk/rnk.cgi?T=cktv…