tracciare oggetto elenco utilizzando ggplot

Question

Sto eseguendo alcune simulazioni che mi stavo chiedendo di tracciare i risultati in un bellissimo ggplot, ma sembra che ggplot non possa occuparsi degli oggetti lista. Qualcuno sa come incollare i risultati in ggplot chart?

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 
    } 

soluzione Ugly:

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

Answer 1

ggplot2 bisogno di un data.frame come dati di origine. Quindi è necessario:

1. trasformare i dati con reshape2 (o plyr o molti altri strumenti)

2. trama usando qplot o ggplot

   esempio

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

Il meglio che potevo fare il tuo indizio è:

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

ma ancora non riflette le probabilità; i punti dovrebbero essere impilati come una sorta di istogramma per riflettere le probabilità.

Un approccio diverso utilizza la funzione di densità, ma può creare problemi se ho molte categorie di campioni da tracciare.

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