2014-09-11 3 views
8

Ich bin neu bei ggplot2 und versuche ein Diagramm zu replizieren, das ich mit filled.contour mit ggplot2 erstellt habe.filled.contour vs. ggplot2 + stat_contour

unter meinem Code:

require(ggplot2) 
require(reshape2) 

#data prep 
scale <- 10 

xs <- scale * c(0, 0.5, 0.8, 0.9, 0.95, 0.99, 1) 
ys <- scale * c(0, 0.01, 0.05, 0.1, 0.2, 0.5, 1) 

df <- data.frame(expand.grid(xs,ys)) 
colnames(df) <- c('x','y') 
df$z <- ((scale-df$x) * df$y)/((scale-df$x) * df$y + 1) 

#filled contour looks good 
filled.contour(xs, ys, acast(df, x~y, value.var='z')) 

#ggplot contour looks bad 
p <- ggplot(df, aes(x=x, y=y, z=z)) 

p + stat_contour(geom='polygon', aes(fill=..level..)) 

ich nicht herausfinden kann, wie ggplot Kontur bekommt die Polygone des ganzen Weg bis zu der oberen linken Seite (es gibt einen Punkt, an (0,10 zu füllen mit) z = 0,99) ... alles, was ich diese seltsamen Dreiecke sind erhalten

Antwort

3

um eine ggplot Version des filled.contour Plot erstellen Sie eine größere data.frame als die df Objekt in Ihrem Beispiel und mit geom_tile haben müssen produzieren die Handlung, die du suchst. Betrachten Sie das folgende:

# a larger data set 
scl <- 10 
dat <- expand.grid(x = scl * seq(0, 1, by = 0.01), 
        y = scl * seq(0, 1, by = 0.01)) 
dat$z <- ((scl - dat$x) * dat$y)/((scl - dat$x) * dat$y + 1) 

# create the plot, the geom_contour may not be needed, but I find it helpful 
ggplot(dat) + 
aes(x = x, y = y, z = z, fill = z) + 
geom_tile() + 
geom_contour(color = "white", alpha = 0.5) + 
scale_fill_gradient(low = "lightblue", high = "magenta") + 
theme_bw()