この記事はマルチポストなアドベントカレンダー記事なんだ、すまない諸君 この記事は R Advent Calendar 2013 : ATNDの27日目の記事で、26日目は@berobero11氏で「RStanで『予測にいかす統計モデリングの基本』の売上データの分析をトレースしてみた」でした…

plotした後に、contour関数をadd=TRUEとして使えばOK。contour関数に食わせるデータはouter関数で作るヨロシ。 x <- seq(-4, 4, 0.1) y <- x z <- outer(x, y, function(x,y){exp(-0.5*(x^2+y^2))}) plot(rnorm(200), rnorm(200), xlim=c(-4, 4), ylim=c(-4,…

とりあえず試行錯誤メモって感じで。やりたい事は関数h内のxの値をf内から変更するって話。 x <- 1:3 g <- function() { f <- function(){x <<- 10} h(f) } h <- function(f) { x <- 1 print(x) f() print(x) } これを実行すると > g() [1] 1 [1] 1 > x [1] …

symnum関数について、http://www.okada.jp.org/RWiki/?%C3%CE%A4%C3%A4%C6%A4%A4%A4%EB%A4%C8%A4%A4%A4%C4%A4%AB%CC%F2%A4%CB%CE%A9%A4%C4%28%3F%29%B4%D8%BF%F4%C3%A3#i1187283にちょいとだけ日本語情報があるだけで、その他の日本語情報が全くなかったので…

掲題の件、そういうことです。 組み込まれてるアルファベット文字一覧(letters変数)でテスト。 > letters [1] "a" "b" "c" "d" "e" "f" "g" "h" "i" "j" "k" "l" "m" "n" "o" "p" "q" "r" "s" "t" "u" "v" "w" "x" "y" "z" > noquote(letters) [1] a b c d e…

setNames関数というのがありまして、"names"属性を持つオブジェクトに対して名前付初期化をするのに便利です。 > setNames( 1:3, c("foo", "bar", "baz") ) foo bar baz 1 2 3 > # 以下のコードと同じ > tmp <- 1:3 > names(tmp) <- c("foo", "bar", "baz") …

F#上でRを使って回帰するにはどうしたらええんじゃい？ってのを簡易メモ。 実行前に#rつかってRProvider & RDotNetにリファレンス張る必要がある。流れとしては リファレンス張って〜 回帰(R.lm)して〜 結果(R.summary)表示 って感じ。最後のsummaryの結果は…