小さな工夫と発見の蓄積
[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。
sudo apt-get install myspell-en-us使用するスペルチェッカーを変更するには、Option -> Configure Texmakerから行う。
> head(x)
game.id picher.id type speed batter.name runner picher.name
1 2005061405 11715 カーブ 111 赤星 憲広 0 松坂 大輔
2 2005082405 11715 カーブ 111 石本 努 0 松坂 大輔
3 2005032601 11715 カーブ 112 北川 博敏 0 松坂 大輔
4 2005050504 11715 カーブ 112 セギノール 0 松坂 大輔
5 2005053102 11715 カーブ 112 ウッズ 10 松坂 大輔
6 2005053102 11715 カーブ 112 アレックス 0 松坂 大輔
visitor.id home.id hv
1 5 7 home
2 8 7 home
3 11 7 home
4 8 7 home
5 7 4 visitor
6 7 4 visitor
aes(x)に、分布を調べる変数を指定する。colourで外枠、fillで塗りつぶす色を指定する。
library(ggplot2)
p <- ggplot(x, aes(x=speed)) +
ggtitle('松坂大輔投手の球速 2005') + xlab('球速 (km/h)')
p + geom_histogram(binwidth=1, colour='grey5', fill='grey35')
..density..というのは定数か何か。こう書くものと覚える。 sizeで線の太さを変える。 p + geom_histogram(binwidth=1, aes(y=..density..),
colour='grey5', fill='grey35') +
geom_density(colour='azure', size=2) +
geom_density(colour='black', size=1.5)
alphaを1より小さい値に指定する（薄く塗りつぶす）と重なった時に綺麗。 facet_grid()を用いて、球種ごとに描画する。p + geom_density(aes(colour=type, fill=type), alpha=.3, size=1)
p + geom_histogram(binwidth=2, aes(y=..density..),
colour='grey5', fill='grey35') +
geom_density() + facet_grid(type ~. )
p + geom_histogram(aes(fill=type),
binwidth=1, stackgroup=T, colour='grey')
p <- ggplot(subset(x, !is.na(hv)), aes(x=speed)) +
ggtitle('松坂大輔投手の球速 2005') + xlab('球速 (km/h)')
p + geom_histogram(binwidth=2, aes(y=..density..),
colour='grey5', fill='grey35') +
geom_density() + facet_grid(type ~ hv)
p <- ggplot(subset(x, !is.na(hv)), aes(type, speed)) +
ggtitle('松坂大輔投手の球速 2005') + xlab('球種') + ylab('球速 (km/h)')
p + geom_boxplot(aes(fill=hv), alpha=0.5) +
theme(legend.position='bottom', legend.title=element_blank())
x <- read.table('http://blog.cnobi.jp/v1/blog/user/ca2e456143c0d20195537cc5daa5fd14/1396483677',
as.is=T, header=T)
names(x) <- c('game.id', 'picher.id', 'type', 'speed', 'batter.name',
'runner', 'picher.name', 'visitor.id', 'home.id')
x$hv <- NA
x$hv[x$home.id %in% 7] <- 'home'
x$hv[x$visitor.id %in% 7] <- 'visitor'
library(ggplot2)
p <- ggplot(x, aes(x=speed)) +
ggtitle('松坂大輔投手の球速 2005') + xlab('球速 (km/h)')
p + geom_histogram(binwidth=1, colour='grey5', fill='grey35') # ...(1)
p + geom_histogram(binwidth=1, aes(y=..density..),
colour='grey5', fill='grey35') +
geom_density(colour='azure', size=2) +
geom_density(colour='black', size=1.5) # ...(2)
p + geom_density(aes(colour=type, fill=type), alpha=.3, size=1)
# ...(3a)
p + geom_histogram(binwidth=2, aes(y=..density..),
colour='grey5', fill='grey35') +
geom_density() + facet_grid(type ~. ) # ...(3b)
p + geom_histogram(aes(fill=type),
binwidth=1, stackgroup=T, colour='grey') # ...(3c)
p <- ggplot(subset(x, !is.na(hv)), aes(x=speed)) +
ggtitle('松坂大輔投手の球速 2005') + xlab('球速 (km/h)')
p + geom_histogram(binwidth=2, aes(y=..density..),
colour='grey5', fill='grey35') +
geom_density() + facet_grid(type ~ hv) # ...(4)
p <- ggplot(subset(x, !is.na(hv)), aes(type, speed)) +
ggtitle('松坂大輔投手の球速 2005') + xlab('球種') + ylab('球速 (km/h)')
p + geom_boxplot(aes(fill=hv), alpha=0.5) +
theme(legend.position='bottom', legend.title=element_blank())
# ...(5)
pdftk ps2pdf pdf2pspdftk input1.pdf input2.pdf input3.pdf cat output merged.pdf
pdftk input.pdf cat 1-endeast output rotated.pdf
pdftk input.pdf cat 1east 2-3west 4 5north 6south output mess.pdf
pdftk input.pdf cat 5-6 1-4 output reordered.pdf pdftk input.pdf cat 2-5 output extracted.pdf
pdf2ps input.pdf tmp.ps ps2pdf tmp.ps small.pdf
for file in *.[pP][dD][fF];
do
echo "$file";
pdftk "$file" cat 1-endeast output "${file%.*}-rotated.pdf";
done
for file in *.[pP][dD][fF];
do
echo "$file";
pdf2ps "$file" __tmp.ps;
ps2pdf __tmp.ps "${file%.*}-small.pdf";
done
ggplot2はRのグラッフィク・ライブラリの１つだ。あまり詳しい設定をいじらなくても、結構きれいな画像を描画してくれる。しかも、しっかりと仕組みを理解すればかなり自由にカスタマイズすることも可能らしい。一念発起して、ggplot2の習得を目指す。p <- ggplot(data, aes(...)) # オブジェクトの作成 p + geom_bar() ... # 描画
* ggplot関数は、必要なデータが１つのデータフレームに入っているようなときに便利。よりシンプルなqplot関数も用意されているが、より汎用性の高いggplot2を覚えれば、qplotでできることは全てできる。
* ggplotオブジェクトをまず作り、そこにどんどん描画（レイヤー）を足していくイメージ。
* aes (エステティック)は、グラフの構成要素と考えればよい。xは横軸変数、fillは塗りつぶしに使う色など。
* geom_bar関数は、カテゴリごとの個体数を表示するようにデザインされた。このカテゴリに対応する変数を、aes(x)として指定する。
* 棒の長さに対応する値はデータから個体数を数えて求めるものであって、今回のように値がすでに与えられているようなケースでは少し工夫が必要：横軸変数（順位）をaes(x)に指定し、値変数（賞金額）をaes(weight)に指定する。
* aes変数は、あとから追加できるし、指定を変更することもできる。
* 上から描画を足すことができる。aes変数は原則としてそのまま使われるが、別のものを指定してもいい。
x <- read.table('http://blog.cnobi.jp/v1/blog/user/ca2e456143c0d20195537cc5daa5fd14/1396338447',
header=F, as.is=T)
names(x) <- c('sei', 'mei', 'shokin', 'shozoku')
x$shokin <- x$shokin/10000
x$rank <- as.factor(1:nrow(x)) # use as category
x$sei <- factor(x$sei, levels=x$sei) # use as category
barwidth <- .3
library(ggplot2)
library(scales) # comma() function
p <- ggplot(x, aes(rank, weight=shokin))
p + geom_bar(width=barwidth) + xlab('') + ylab('賞金額（万円）') +
ggtitle('日本棋院賞金ランキング 2013') +
scale_y_continuous(labels=comma) # ...(1)
p + geom_bar(aes(fill=shozoku), width=barwidth) +
xlab('') + ylab('賞金額（万円）') +
ggtitle('日本棋院賞金ランキング 2013') +
scale_y_continuous(labels=comma) +
theme(legend.position='bottom', legend.title=element_blank())
# ...(2)
p + geom_bar(aes(x=sei), width=barwidth) +
xlab('') + ylab('賞金額（万円）') +
ggtitle('日本棋院賞金ランキング 2013') +
scale_y_continuous(labels=comma) # ...(3)
p <- p + geom_bar(width=barwidth) +
xlab('') + ylab('賞金額（万円）') +
ggtitle('日本棋院賞金ランキング 2013') +
scale_y_continuous(labels=comma) # 中間地点
p + geom_text(aes(y=shokin, label=sei),
col='white', size=5.2, fontface=2) +
geom_text(aes(y=shokin, label=sei)) # ...(4)
p + geom_bar(width=barwidth) +
xlab('') + ylab('賞金額（万円）') +
ggtitle('日本棋院賞金ランキング 2013') +
scale_y_continuous(labels=comma) +
geom_text(aes(y=shokin, label=sei), vjust=-0.1) # ...(5)
p + geom_bar(width=barwidth) + xlab('') + ylab('賞金額（万円）') +
ggtitle('日本棋院賞金ランキング 2013') +
scale_y_continuous(labels=comma) +
geom_text(aes(y=shokin, label=sei), hjust=-0.3) # ...(6)