Any suggestions for how I can plot mixEM type data using ggplot2

by

Look at the structure of the returned object (this should be documented in the help):

> # simple mixture of normals:
> x=c(rnorm(10000,8,2),rnorm(10000,17,4))
> xMix = normalmixEM(x, lambda=NULL, mu=NULL, sigma=NULL)

Now what:

> str(xMix)
List of 9
 $ x         : num [1:20000] 6.18 9.92 9.07 8.84 9.93 ...
 $ lambda    : num [1:2] 0.502 0.498
 $ mu        : num [1:2] 7.99 17.05
 $ sigma     : num [1:2] 2.03 4.02
 $ loglik    : num -59877

The lambda, mu, and sigma components define the returned normal densities. You can plot these in ggplot using qplot and stat_function. But first make a function that returns scaled normal densities:

sdnorm =
function(x, mean=0, sd=1, lambda=1){lambda*dnorm(x, mean=mean, sd=sd)}

Then:

qplot(x,geom="density") + stat_function(fun=sdnorm,args=list(mean=xMix$mu[1],sd=xMix$sigma[1], lambda=xMix$lambda[1]),fill="blue",geom="polygon")  + stat_function(fun=sdnorm,args=list(mean=xMix$mu[2],sd=xMix$sigma[2], lambda=xMix$lambda[2]),fill="#FF0000",geom="polygon")

enter image description here

Or whatever ggplot skills you have. Transparent colours on the densities might be nice.

ggplot(data.frame(x=x)) + 
 geom_histogram(aes(x=x,y=..density..),fill="white",color="black") +
 stat_function(fun=sdnorm,
    args=list(mean=xMix$mu[2],
             sd=xMix$sigma[2],
             lambda=xMix$lambda[2]),
             fill="#FF000080",geom="polygon") +
 stat_function(fun=sdnorm,
    args=list(mean=xMix$mu[1],
             sd=xMix$sigma[1],
             lambda=xMix$lambda[1]),
             fill="#00FF0080",geom="polygon")

producing:

enter image description here

More Related Contents:

  1. Using ggplot2, can I insert a break in the axis?
  2. Fixing the order of facets in ggplot
  3. Grouped bar plot in ggplot
  4. Plot with conditional colors based on values in R [duplicate]
  5. Putting mathematical symbols and subscripts mixed with regular letters [duplicate]
  6. Create discrete color bar with varying interval widths and no spacing between legend levels
  7. Overlaying histograms with ggplot2 in R
  8. Increase distance between text and title on the y-axis
  9. Align ggplot2 plots vertically
  10. It is possible to create inset graphs?
  11. R + ggplot2 => add labels on facet pie chart [duplicate]
  12. How to add a number of observations per group and use group mean in ggplot2 boxplot?
  13. Setting Defaults for geoms and scales ggplot2
  14. Normalizing y-axis in histograms in R ggplot to proportion
  15. R Plotting confidence bands with ggplot
  16. R – emulate the default behavior of hist() with ggplot2 for bin width
  17. Add percentage labels to a stacked barplot
  18. How to jitter both geom_line and geom_point by the same magnitude?
  19. How to show only part of the plot area of polar ggplot with facet?
  20. Different legend-keys inside same legend in ggplot2
  21. Plot two graphs in a same plot
  22. ggplot multiple grouping bar
  23. Showing equation of nls model with ggpmisc
  24. pie chart with ggplot2 with specific order and percentage annotations
  25. R stat_smooth all points
  26. Color one point and add an annotation in ggplot2?
  27. Preserve proportion of graphs using grid.arrange
  28. Remove unused factor levels from a ggplot bar plot
  29. How to display emojis in ggplot2 using emo package in R?
  30. adding x and y axis labels in ggplot2

Leave a Comment