Question

I run a bunch of simulations to evaluate type I error, so the result is a vector such as

pdata = c(0,0,0,0,0,0,0,0,0,0.07,0,0.02,0.03)

The mean of the simulated vector should be 0.05. Now I am thinking of a way to display the results via boxplots. The default function in R

boxplot(pdata)

gives a boxplot that is rather hard to see the typical value as there are many 0’s. In addition, it shows the median, but what I really want is the mean to be displayed on the plot. Are there any graphical display that is effective in such situation? I know that I can simply report the numerical values, but because my simulation involves other factors which I hope to compare, a boxplot-like graph will be ideal. Thanks!

Answer 1

Something like this maybe :

 plot(table(pdata))

Here a ggplot2 version :

   ggplot(as.data.frame(table(pdata)),aes(x=pdata,y=Freq))+geom_bar()