title: Multiple power law plots with package poweRlaw
body: <p>I tried using the package <a href="https://github.com/csgillespie/poweRlaw" rel="nofollow noreferrer">poweRlaw</a> to plot a few powerlaw fits. It seems working for a single curve. But I could not plot multiple plots on the same graph.
  <a href="https://stackoverflow.com/questions/14476426/how-to-plot-degree-distribution-in-r">Ref:</a> Is there a way within this package?
  [P.S. I am a novice]</p>
  
  <pre><code>set.seed(1)
  x1 &lt;- ceiling(rlnorm(1000, 4))
  
  x2 &lt;- ceiling(rlnorm(1000, 2))
  
  library(poweRlaw)
  
  pl_d = pl_data$new(x1)
  plot(pl_d)
  
  #Now fit the powerlaw
  m = displ$new(pl_d)
  #Estimate the cut-off
  #estimate_xmin(m)
  aaa &lt;- estimate_xmin(m)
  aaa &lt;- as.data.frame(aaa)
  aaa &lt;- aaa[2,1]
  x_min &lt;- min(table(x))
  m$setXmin(aaa); m$mle()
  #Plot the data and the PL line
  #plot(m)
  lines(m, col=2)
  
  # next POWER LAW graph
  
  #Plot the data
  pl_d = pl_data$new(x2)
  points(pl_d)
  
  #Now fit the powerlaw
  m = displ$new(pl_d)
  #Estimate the cut-off
  #estimate_xmin(m)
  aaa &lt;- estimate_xmin(m)
  aaa &lt;- as.data.frame(aaa)
  aaa &lt;- aaa[2,1]
  x_min &lt;- min(table(x))
  m$setXmin(aaa); m$mle()
  #Plot the data and the PL line
  #points(m)
  lines(m, col=3)
  </code></pre>