ggplot2 Clips Daten in stat_density (Dichteplot)

Question

Datendatei: http://ubuntuone.com/3x5z4kFEUcVUB8KXYjICpD

die obigen Daten ähnliche Clips einige davon Plotten:

tc = read.table('tertiary-tc.csv', header=T, sep=',') 
library(ggplot2) 
old <- theme_set(theme_bw()) 
pg <- ggplot(tertiary, aes(Conductivity)) + stat_density(geom = 'path' 
    , position = 'identity') + facet_wrap(~Lithology) 
print(pg) 

während dies nicht:

library(lattice) 
densityplot(~Conductivity | Lithology 
    , data = tc 
    , groups=Lithology 
    , plot.points=T 
    , ref = T 
    , main="Density plot" 
    , auto.key=F 
    , scales=list(tck=-1)) # set to a negative to plot inside ticks? Works! 
dev.off() 

Die Punkte, die nicht geplottet wurden (nur schnell betrachtet), sind zwei Punkte, die in der Spalte "Lithologie" als "Dolomit" klassifiziert sind. Zeichnet ggplot2 sie nicht auf, weil sie nicht skaliert sind? Oder weil es nur zwei Punkte gibt?

Answer 1

So etwas könnte man

density2 <- function(x, select, ...){ 
    n <- nrow(x) 
    range <- range(x[[select]]) 
    xgrid <- seq(range[1], range[2], length=200) 

    weight <- rep(1, n)/n 
    d <- density(x[[select]], ..., 
        weight=weight, from=range[1], to=range[2]) 

    as.data.frame(d[c("x","y")]) 

} 


summaries <- ddply(tc, "Lithology", density2, select="Conductivity") 

ggplot(summaries, aes(x, y)) + facet_wrap(~Lithology) + 
    geom_path() 

nicht sicher erhalten begonnen, warum es irgendwie beschnitten sieht, ich bin nicht vertraut mit density()