1

データファイル: http://ubuntuone.com/3x5z4kFEUcVUB8KXYjICpD

上記のデータを次のようにプロットすると、その一部が切り取られます。

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)

これはしません：

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()

プロットされていないポイント (ざっと見ただけ) は、岩相列の下で「ドロマイト」として分類される 2 つのポイントです。それらはスケール外であるため、ggplot2はそれらをプロットしていませんか? それとも2点しかないから？

4

1 に答える 1

2

このような何かがあなたを始めるかもしれません

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()

なぜそれがどういうわけかトリミングされているように見えるのかわからない、私はよく知らないdensity()

于 2012-04-12T08:11:02.430 に答える