2019-01-01から1年間の記事一覧

乱れた列を揃えたい(主にSQL用途)

create table hoge ( hoge string hageeeeeeeeeeee int moge string ) を create table hoge ( hoge string hageeeeeeeeeeee int moge string ) のように揃えたいときは、↑を一旦適当なファイルに突っ込んでおいて、column するといい $ column -t -s $' ' <filename></filename>

ggplot2の軸とかに数式を使いたい

expression 使えばよろし。 library("mmetrics") df <- gsummarise(iris, Species, metrics = define(y=sum(Petal.Width))) ggplot(df, aes(x = Species, y = y)) + geom_bar(position = "dodge", color="black", stat="identity") + labs(x = expression(al…

特異値分解の検算

R

すぐ忘れる。 svd関数の戻り値をxとして x$u %*% diag(x$d) %*% t(x$v) で元通り。 > b <- array(rnorm(10), dim=c(5,2)) > b [,1] [,2] [1,] -0.5036082 -1.218417405 [2,] -0.7905285 1.331412278 [3,] 0.8311744 -0.147847237 [4,] -0.8628976 0.68170753…

ある関数のデフォルト引数がある引数のみを、そのデフォルト引数とともに抜く

例えばここにこんな関数があったとする。 f <- function(x, y = 100, args = list()){ y } そしてこの関数のデフォルト引数がある引数のみを、そのデフォルト引数とともに抜きたいと思うことがある、俺にはある。 こういう処理をしてくれる関数 get_default_…

各大グループにおいて、各小グループの値がどの程度の割合を占めているのかを出す

R

昔、 r-wakalangで教えてもらったがメモってなかった話なのでメモっておく。 こういうデータ df があった時に > df <- data.frame( + big_group=rep(letters[1:2], 30), + small_group=rep(1:3, 10), + value=runif(30) + ) > head(df, 10) big_group small_…

データ増やした時の標準誤差の減り具合

掲題の件、これはふつー「1/√データ数」ですが、適当にデータとして0だけ突っ込んでいったらどうなるんじゃいと思ってやってみたら大体似たようなもんだった。 理論的な話はあとで考える。 > size_unit <- 10^3 > x <- c() > y <- c() > for (i in 0:9){ + …

各Bin・Cellに同じ数だけデータが入っているHistgramを作る

R

調べた感じ引数ではやれなそうだったので、やむなく quantile 関数で対応することにした。 > x <- rnorm(1000) > tmp <- hist(x, breaks=quantile(x, c(0, 0.2, 0.4, 0.6, 0.8, 1))) > tmp $breaks 0% 20% 40% 60% 80% 100% -3.4892322 -0.9159784 -0.265114…