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

   例えば

   ## 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")