使用ggplot绘制列表对象

Question

我正在运行一些模拟，我想知道在一个漂亮的ggplot中绘制结果，但似乎ggplot无法处理列表对象。有谁知道如何将结果粘贴到ggplot图表中？

   N <- 8619170         
   nn <- c(1000, 1200, 3000)
   p <- .27     
   nsim <- 100

    phat <- list()
    for (i in 1:length(nn)) {
    n <- nn[i]
    x <- rhyper(nsim, N * p, N * (1 - p), n)
    phat[[i]] <- x / n
    }

丑陋的解决方案：

    names(phat) <- paste("n=", nn)
    stripchart(phat, method="stack")
    abline(v=p, lty=2, col="red")

Answer 1

ggplot2需要data.frame作为源数据。所以你需要：

1. 使用reshape2（或plyr或许多其他工具）转换数据

2. 使用qplot或ggplot

   绘图

   e.g

    ## transform data
    require(reshape2)
    h <- do.call(cbind, phat)
    h.melt <- melt(h)
   
    ## rename variables so they look nicer on plots
    names(h.melt) <- c("test","N","value")     
   
    ## stripchart (not shown)
    qplot(data = h.melt, x = value,y = N,color=N)+geom_point()
   
    ## histogram (not shown)    
    ggplot(h.melt,aes(x=value,fill=N))+geom_histogram()+facet_grid(N~.)
   
    ## dotplot with rug (not shown)
    ggplot(h.melt,aes(x=value,fill=N))+geom_dotplot()+facet_grid(N~.)+geom_rug()  
   
    ##density plot with rug (shown below)
    ggplot(h.melt,aes(x=value,fill=N))+geom_density()+facet_grid(N~.)+geom_rug() 
   

   

Answer 2

我能做的最好的事情就是：

qplot(data = h.melt, x = value,y = Var2)+ geom_point(shape=1, size=5)

但它仍然没有反映概率;这些点应该作为一种直方图堆叠以反映概率。

另一种方法是使用密度函数，但如果我有许多样本类别可以绘制，那么它可能会弄乱。

ggplot(h.melt, aes(x=value, fill=Var2)) + geom_density(alpha=.5, position="identity")