Clipboard

生き物 / R / 研究 / 美味しいもの

QGISのOpenLayersプラグインにSTAMEN MAPタイルを追加する

QGISのベースマップにWebからstamen mapを利用するのは便利だけれど、地名なども表示されるので少し邪魔。
maps.stamen.com

stamen mapのToner backgroundやTerrain backgroundを使いたい。
下記の方法でできた。

cecilialiao.com

Qgis2.18だと、gdalTMS='stamen_toner.xml'というxmlファイルもhtmlと同じように編集する必要があった。
phyton全然分からないけど、書いてある通りに編集したら使えるようになった。

R: ggplotでのエラーバープロット

エラーバーはgeom_errorbarで描画できる。x軸方向のエラーバーはgeom_errorbarhを使用する。

data(iris)
d <- summaryBy(Sepal.Length + Sepal.Width ~ Species, data=iris, FUN=c(mean,sd))

g <- ggplot(d, aes(x=Sepal.Width.mean, y=Sepal.Length.mean, col=Species)) + 
geom_point(size=5) +
geom_errorbar(aes(ymax = Sepal.Length.mean + Sepal.Length.sd, ymin= Sepal.Length.mean - Sepal.Length.sd), data=d, width=0.05) +
geom_errorbarh(aes(xmax = Sepal.Width.mean + Sepal.Width.sd, xmin= Sepal.Width.mean - Sepal.Width.sd), data=d, height=0.05)

g

f:id:yumico:20140528121327j:plain

R : データの集計

library(doBy)のsummaryBy()関数が便利。
library(plyr)でも集計できるけれど、こちらの方が直観的に書きやすい感じ。

data(dietox)
dietox12 <- subset(dietox,Time==12)

summaryBy(Weight+Feed~Evit+Cu, data=dietox12,FUN=c(mean,var,length))