あんまよくわかってなくて、酷い目に会いましたよと。まず、以下のようなシミュレーション結果「ある経路を表すx, y列を持つデータ」が手もとにあったとする。 > head(df) x y 1 119.2424 14.39025 2 118.9490 14.94898 3 118.6767 16.27511 4 117.9835 17.2…

東京で尻を洗うカレンダープロット お、おお、こんなのあるんだ。昔、似たような感じの話で 2011年・TOPIXのヒートマップを作ってみた - My Life as a Mock Quant で試したコードよりも簡素化されていて大変よい。大変良いので真似てみた。ここでは、今目下…

やってみよう分析！おまけ 2 - 1： Excel VBAで強化学習(Q-learning, ε-greedy / softmax 行動選択) のR実装版。強化学習自体の解説は上の記事読んどいたらいい。めんどいのでとりあえずε-greedyのみやった。計算結果は > Qlearning() [,1] [,2] [1,] 5 20.0…

なんか殺伐とした空気を吹き飛ばせるんじゃないか、そう強く感じさせるパッケージを２つみつけた。 まず、Githubからのインストール。 library(devtools) install_github("wrathematics/Rfiglet") install_github("sckott/cowsay") Rfigletから使ってみる。 …

@teramonagi 1:length(x) の代わりに seq_along(x) を使うと良いってごみ箱が言ってた。— kos59125 (@kos59125) 2014, 8月 18seq_alongの方が、空っぽのベクターに対しても安全に動作しますよっと。 > x <- 1:3 > 1:length(x) [1] 1 2 3 > seq_along(x) [1] …

掲題の件、そういうことです。動的にインタラクティブに地図いじりたい！って人はgoogleVisパッケージでいいけど、静的な地図ならggmapなな印象だね！・・・と以下に書いてあった。 R help - googleVis plot and knitr/sweave というわけで、今、静的な地図…

他の部分に対する突っ込みはよくわからないが、とりあえずここはこうだろうと。 データサイエンスの本を買ったが一部アツすぎる内容で萎えた話 | JUMPERZ.NET Blog ↑この本を買ったのだが、95ページ周辺にむちゃくちゃな事が書いてあって萎えた。統計学では…

以下のプラグイン plugins=c("cpp11")をcppFunctionに指定するだけ。Winでもいけた。C++11なんで戻り値をdecltypeにしてやろうと思ったら、できなかった。windowsのgccが古いせいかな？ library(Rcpp) cppFunction(' std::vector<double> twoTimes(std::vector<double> xs) {</double></double>…

Rでの例外処理(tryCatchの使い方メモ) を参考にもっと簡略化したもの。 > #文字列型⇒日時型へ変換 > as.POSIXct("2013-3-2 13:45:32") [1] "2013-03-02 13:45:32 JST" > #文字列型⇒日時型へ変換？？？ > as.POSIXct("2013年3月2日") Error in as.POSIXlt.cha…

ちょいとRcppをがっつり使ってみようと、その試行錯誤記録が続く予定。 速度検証のテストコードとして、ランダムウォークする系列データを生成するコードを書いた。 library(Rcpp) sourceCpp(code=' #include<vector> #include <functional> #include <numeric> #include <Rcpp.h> using namespace</rcpp.h></numeric></functional></vector>…

継続渡しスタイル - Wikipedia 末尾再帰 - Wikipedia ここらに書いてるのはこういうことだろうという。以下はリストの中の最大値を求めるコード。 let maximum xs = let rec maximum proc xs = match xs with | [] -> failwith "List is empty!!!" | x::[] -…

俺たちの最強Rクオンツブロガー・Timely Portfolioことタイポさんの記事をみると、ファーマ・フレンチの３ファクターをサクッとゲットできるなんて方法を解説されていたので、パク真似てみた。上の記事で、彼はzooパッケージによる可視化をしていたので、こ…

やりたいことへの試行錯誤 例えばdplyrとmagrittrを活用して、mtcarsデータセットのmpg, cylの両列に対してround関数を適用するには、以下のように書くわけです。 > library(dplyr) > mtcars %>% select(mpg, cyl) %>% sapply(round) mpg cyl [1,] 21 6 [2,]…

どうやら as.Date(character()) として空の日付を作るのが大事っぽい。 > x <- data.frame(as.Date(character()), as.double()) > str(x) 'data.frame': 0 obs. of 2 variables: $ as.Date.character...:Class 'Date' num(0) $ as.double.. : num > rbind(x,…

これがオーソドックスな正しいやり方なのかわからないけど、とりあえずやってみた。 optim関数で数値的に逆関数の計算したい(一次元版) とりあえず一次元の関数をと書いた時の、が与えられた下でのの推定を考えたい。これはもちろん解析解があるとハッピーな…