title: Graphical output of density for the function gammamixEM (package mixtools)
body: <p>I'm using the function <code>gammamixEM</code> from the package <code>mixtools</code>. How can I return the graphical output of density as in the function <code>normalmixEM</code> (i.e., the second plot in <code>plot(...,which=2)</code>) ? </p>
  
  <p><strong>Update</strong>:</p>
  
  <p>Here is a <a href="http://www.rdocumentation.org/packages/mixtools/versions/1.0.4/topics/gammamixEM" rel="nofollow noreferrer">reproducible example</a> for the function <code>gammamixEM</code>:</p>
  
  <pre><code>x &lt;- c(rgamma(200, shape = 0.2, scale = 14), rgamma(200, 
       shape = 32, scale = 10), rgamma(200, shape = 5, scale = 6))
  out &lt;- gammamixEM(x, lambda = c(1, 1, 1)/3, verb = TRUE)
  </code></pre>
  
  <p>Here is a <a href="http://www.rdocumentation.org/packages/mixtools/versions/1.0.4/topics/normalmixEM" rel="nofollow noreferrer">reproducible example </a>for the function <code>normalmixEM</code>:</p>
  
  <pre><code>data(faithful)
  attach(faithful)
  out &lt;- normalmixEM(waiting, arbvar = FALSE, epsilon = 1e-03)
  plot(out, which=2)
  </code></pre>
  
  <p><a href="https://i.stack.imgur.com/t1V54.jpg" rel="nofollow noreferrer"><img src="https://i.stack.imgur.com/t1V54.jpg" alt="enter image description here"></a></p>
  
  <p>I would like to obtain this graphical output of density from the function <code>gammamixEM</code>.</p>