/>
小さな工夫と発見の蓄積
[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))