/>
小さな工夫と発見の蓄積
[PR]上記の広告は3ヶ月以上新規記事投稿のないブログに表示されています。新しい記事を書く事で広告が消えます。
ただいまコメントを受けつけておりません。
library(ggplot2) library(reshape2) x <- read.csv('http://blog.cnobi.jp/v1/blog/user/ca2e456143c0d20195537cc5daa5fd14/1399760770', sep='\t', as.is=T) z <- read.csv('http://blog.cnobi.jp/v1/blog/user/ca2e456143c0d20195537cc5daa5fd14/1399761108', sep='\t', as.is=T, header=F) fit <- loess(rent ~ lon + lat, data=x, span=0.5) gd <- 100 lon <- seq(min(x$lon), max(x$lon), length=gd) lat <- seq(min(x$lat), max(x$lat), length=gd) D <- expand.grid(lon=lon, lat=lat) v <- melt(predict(fit, D), value.name='rent') v$lon <- as.numeric(gsub('lon=', '', v$lon)) v$lat <- as.numeric(gsub('lat=', '', v$lat)) p <- ggplot(v, aes(lon, lat, z=rent)) p + geom_point(aes(col=rent), linetype='blank', size=5, shape=15) + stat_contour(bins=10, col='white', size=1.2) + stat_contour(bins=10) + geom_line(data=z, aes(V4, V3, z=V3), size=2, col='lightgreen') + geom_line(data=z, aes(V4, V3, z=V3), size=1.2, color='coral4') + geom_text(data=z, aes(V4, V3, z=V3, label=V1), size=5.12, col='white', fontface=2) + geom_text(data=z, aes(V4, V3, z=V3, label=V1), fontface=2, size=5) + scale_colour_gradientn(colours=rich.colors(10)) p + geom_tile(aes(fill=rent)) + stat_contour(bins=10, col='white', size=1.2) + stat_contour(bins=10) + geom_line(data=z, aes(V4, V3, z=V3), size=2, col='lightgreen') + geom_line(data=z, aes(V4, V3, z=V3), size=1.2, color='coral4') + geom_text(data=z, aes(V4, V3, z=V3, label=V1), size=5.12, col='white', fontface=2) + geom_text(data=z, aes(V4, V3, z=V3, label=V1), fontface=2, size=5) + scale_fill_gradientn(colours=rich.colors(10))