Resolved: Set Axis Limits of mixfitEM plot

rm(list=ls())<br />
library(mixR)<br />
library("tidyverse")<br />
set.seed(07062022)</p>
<p>N <- 1000
lbda=.5
mu1=1
Del=3
mu2=mu1+Del
components <- sample(1:2,prob=c((1-lbda), lbda),size=N,replace=TRUE)
mus <- c(mu1,mu2)
sds <- sqrt(c(1,.5))

Y <- rlnorm(n=N,meanlog=mus[components],sdlog=sds[components])
[/code]

</p>



Then I try to fit two log normals into this data using <code>mixfit()</code> from <code>mixR</code>. In the first attempt I tried to use the <code>plot()</code> function from the <code>base</code> package; but <code>plot</code> ignores the <code>xlim</code> argument.<br> 


<p>[code]mod4 <- mixfit(Y, ncomp = 2, family = 'lnorm')

plot(mod4, title = 'Log-Normal Mixture (2 components)', xlim=c(0,200))
[/code]

</p>



Then I tried to plot with ggplot, which, in theory should be possible according to the <a href="https://search.r-project.org/CRAN/refmans/mixR/html/plot.mixfitEM.html" rel="nofollow noreferrer">mixR manual</a>. But ggplot does not understand the <code>mixfitEM</code> class.<br> 


<p>[code]ggplot(mod4)+<br />
  + coord_cartesian(xlim = c(0, 200))

> ggplot(data.frame(mod4))+
+   + coord_cartesian(xlim = c(0, 200))
Error in as.data.frame.default(x[[i]], optional = TRUE, stringsAsFactors = stringsAsFactors) : 
  cannot coerce class ‘"mixfitEM"’ to a data.frame
> ggplot(mod4)+
+   + coord_cartesian(xlim = c(0, 200))
Error in `fortify()`:
! `data` must be a data frame, or other object coercible by `fortify()`, not an S3 object with class mixfitEM.
Run `rlang::last_error()` to see where the error occurred.

plot(mod4, title = 'Log-Normal Mixture (2 components)') + xlim(c(0, 200))