私的メモ VOLUMEコマンド後に当該ディレクトリに作成したファイルは作成されない http://qiita.com/namutaka/items/f6a574f75f0997a1bb1d 環境変数はコンテナ起動時にその時のユーザに割り当てられる（多分） http://stackoverflow.com/questions/36978880/w…
pacman使えよって話だが。 # ほしいパッケージ全部書いとく packages <- c("ggmap", "dplyr", "tidyr", "data.table", "rgl", "stringr", "R.utils", "scales") # ないやつはインストール new_packages <- packages[!(packages %in% installed.packages()[,"…
要するに、”あるグループごとのデータフレームに分解してデータを持つ”ということをしたいんだ、そしてこんなんでた。 > x <- iris %>% + group_by(Species) %>% + do(val=data.frame(.)) > x Source: local data frame [3 x 2] Groups: <by row> Species val (fctr)</by>…
使うたびにはまる系なのでめもる＆答えはgeom_rect()でfill=NAとする。 library("ggplot2") ggplot(mtcars, aes(qsec, mpg)) + geom_point()+ geom_rect(aes(xmin=18, xmax=20, ymin=21, ymax=23), fill=NA, color="red")
dplyr でグループごとに複数カラムを追加したい #rstatsj に対する１回答。 > library("dplyr") > library("tidyr") > probs <- seq(0, 1, by=0.25) > iris %>% group_by(Species) %>% + summarize(tmp=list(data.frame(prob=probs, quantile=quantile(Sepal…
そういうことです。今のところの私の答えは以下で breaksでちゃんとX軸のどれを書くかを指定 trans_formatでべき乗表示に変換 です。 ggplot(data=data.frame(x=10^(1:4), y=1:4), aes(x=x, y=y)) + scale_x_log10(breaks=10^(1:4), labels=trans_format('lo…
はじめに 主に俺の健康を管理するために、ヘルスプラネットというタニタの運営するサービスのAPIを叩いてデータを取得するパッケージを作った。 とりあえずinnerscan系の情報 体重 (kg) 体脂肪率 (%) 筋肉量 (kg) 筋肉スコア 内臓脂肪レベル 基礎代謝量 (kca…
そういう時もある。書き方結構あるなぁ。 書き方1 import scala.io.StdIn.readLine object Main { def main(args: Array[String]){ val Array(a, b) = readLine().split(" ").map(_.toInt) println((a,b)) val s = readLine() println(s) } } 実行例 $ scala…
掲題の件、そういうことです。 この上場投信(ETF)のサイトは NEXT FUNDS NOMURA 原油インデックス連動型上場投信 《愛称》原油先物ETF（1699） ここではついでに、 この上場投信価格が「WTI(West Texas Intermediate)原油先物 × USDJPY」でちゃんと計算でき…
Linuxだとwc使えばいいんだけど、Windowsだとどうしたらいいんだろうかと思って見つけたメーリングリスト。 比較的大きいファイルを扱うので、効率的に処理したいところ。 http://r.789695.n4.nabble.com/Fast-way-to-determine-number-of-lines-in-a-file-t…
これは君らにとって小さな一歩だが、私にとって偉大な飛躍なんだ。 library("ggplot2") library("scales") #元データ(適当) df <- data.frame(x=rep(1:10, 2), y=rnorm(20), type=rep(LETTERS[1:2], each=10)) ggplot(df, aes(x=x, y=y)) + #フォントを大き…
これは君らにとって小さな一歩だが、私にとって偉大な飛躍なんだ。 library("ggplot2") library("scales") #元データ(適当) df <- data.frame(x=rep(1:10, 2), y=rnorm(20), type=rep(LETTERS[1:2], each=10)) ggplot(df, aes(x=x, y=y)) + #%+replace%でテ…
なぜWindows？とい無粋な話はどうでもいい。 最新版をビルドしているありがたい方がいらっしゃるので、これを拝借する。 https://www.barik.net/archive/2015/01/19/172716/ のリンクにあるtar.gzファイルをダウンロードして、そこをHADOOP_HOMEと環境変数に…
scala> for(x <- 1 to 3; y <- 1 until 5){println(x + " " + y)} 1 1 1 2 1 3 1 4 2 1 2 2 2 3 2 4 3 1 3 2 3 3 3 4 これはイケてる。
モンテカルロ法でコール・オプション価格 - My Life as a Mock Quantを流行りの分散処理環境Apache Sparkでやってみたという話。 コード書くのは一瞬なんだけど、依存関係周りの解消が辛かったので、積極的にコード残しておきたい。下記コード+sbtファイルを…
とりあえず慣れてるのからということで・・・Wolfram Alphaで出した結果と大体あってるので、これでよさげだ。 import scala.util.Random import scala.math object App { def main(args: Array[String]) = { val random = new Random() // Monte Carlo sett…
そういうことっぽい。 scala> Seq(1,2,3) res0: Seq[Int] = List(1, 2, 3) scala> Seq.apply[Int](1,2,3) res1: Seq[Int] = List(1, 2, 3) 参考 第22章：Scalaの抽出子
TVにつないでるRaspberry piの画面を引っ張りたい。そんなとき。 sudo apt-get install x11vnc sudo apt-get install iceweasel browser-plugin-gnashからのGUI起動時設定 sudo vim /etc/xdg/lxsession/LXDE-pi/autostart で中身を以下のようにして、Raspber…
R言語の中でfactor(因子型）は数値として扱われており、外っつらが同じでも中が異なる場合があるので注意が必要だ。 例えば >|r| > rlove <- c("R", "L", "O", "V", "E") > x <- factor(rlove) > x [1] R L O V E Levels: E L O R V > as.numeric(x) [1] 4 2…
Scalaおじさんにもなってきたので、マメにハマりをめも。Typesafe Activator使ってScalaしてて、Hello PlayFramework的なTutorialが動かなくてこまった。 同じ症状の人はこちら。 https://github.com/typesafehub/activator/issues/1060 要するに、build.sbt…
塗りつぶした長方形じゃなくて、淵にだけ色つけたい＆そもそもgeom_tileの指定方法知らなかった。 df <- data.frame( x = c(1,2), y = c(3,4), z = c(1,2) ) ggplot(df, aes(x, y)) + geom_tile(aes(fill = z)) ggplot(df, aes(x, y)) + geom_tile(aes(fill …
関数名覚えられないマン。 > today <- Sys.Date() > class(today) [1] "Date" > today [1] "2015-12-26" > strftime(today, format="%Y%m%d") [1] "20151226" > class(strftime(today, format="%Y%m%d")) [1] "character"
この記事はR Advent Calendar 2015の9日目の記事です。ggvisパッケージを使ってちょっとだけ凝った可視化をしてみたので、その備忘録＆情報共有です。 内容としては (適当な)予測モデルを作って 予測値からの実際のずれが上に(下に)大きいほど、描画される点…
連続する数値をたくさん作ろうとコードを書いたら以下のように怒られた。 > seq(0, 1:15) Error in seq.default(0, 1:15) : 'to' は長さ 1 でなければなりません 世間はなんて辛いんだ。 そこで、↑が動くor等価な出力をする結果がほしいとr-wakalangで質問し…
自分で作ったパッケージについて、パッケージをCRANにあげるのは大変な苦労をともなう。 それを少しでも軽くしたい、手元や継続的CI環境ではテストに成功しているのにCRANでもテストするだなんて面倒くさい。そんな時にはtestthatパッケージで単体テストを書…
ロジックメモ ある確率変数とその写像で定義される２つの確率変数、特に確率変数の確率分布について考えたい。 の確率分布関数は と確率変数の確率分布関数を用いてあらわすことができる。ここで、は単調増加な関数であると仮定している。単調減少の場合は不…
あたまだし 掲題の件、最近出たわかりやすいデータ分析の入門書、 において 売上げ v.s. 気温 として回帰分析をしていたケースがあった。これが時系列データに対するやってはいけない回帰、すなわち、みせかけの回帰に相当するケースではないか？という話を…
"グループ毎にXXXしたい"という時のXXXの内容が複雑な場合、関数として処理を切り出したくて、以下のように書いてみた。 処理の内容は"Speciesのグループ毎にその最終行のPetal.Lengthの値を取得する"だ。 > process <- function(df){ + df[nrow(df),]$Petal…
したい。 関数によってミリ秒だったり秒だったり、OSによって実時間だったりCPU時間だったりする点に注意。 #include <iostream> #include <vector> #include <chrono> #include <time.h> int main() { const int size = 1000 * 1000; //C++11版, ミリ秒単位、実時間ベース std::vector<int> x; auto</int></time.h></chrono></vector></iostream>…
gccのバージョンを4.9系にしないとC++の正規表現ライブラリが使えないので、バージョンアップ。そのためにはScientific Linux Cern 6 repositoryを使う。 cd /etc/yum.repos.d wget http://linuxsoft.cern.ch/cern/scl/slc6-scl.repo 本家のレポジトリだと重…
