Error(s) in re-building vignettes:
  ...
--- re-building ‘MLdist.Rmd’ using rmarkdown
Warning in xy.coords(x, y, xlabel, ylabel, log) :
  182 y values <= 0 omitted from logarithmic plot
Warning in xy.coords(x, y, xlabel, ylabel, log) :
  165 y values <= 0 omitted from logarithmic plot
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(function(th) exp(-g(th)), lower = l.th, upper = u.th,  :
  the integral is probably divergent
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in uniroot(function(th) log(g(th)), lower = l.th, upper = u.th,  :
  -Inf replaced by maximally negative value
Warning in .integrate2(g1, lower = a, upper = b, subdivisions = subdivisions,  :
  roundoff error is detected in the extrapolation table
--- finished re-building ‘MLdist.Rmd’

--- re-building ‘parametrisation.Rmd’ using rmarkdown
pandoc-citeproc: Error in $: Incompatible API versions: encoded with [1,20] but attempted to decode with [1,17,5,4].
CallStack (from HasCallStack):
  error, called at ./Text/Pandoc/JSON.hs:111:48 in pandoc-types-1.17.5.4-5tHZ3B61A58JaKOMxwGQR4:Text.Pandoc.JSON
Error running filter /usr/bin/pandoc-citeproc:
Filter returned error status 1
Error: processing vignette 'parametrisation.Rmd' failed with diagnostics:
pandoc document conversion failed with error 83
--- failed re-building ‘parametrisation.Rmd’

--- re-building ‘probsNquantiles.Rmd’ using rmarkdown
pandoc-citeproc: Error in $: Incompatible API versions: encoded with [1,20] but attempted to decode with [1,17,5,4].
CallStack (from HasCallStack):
  error, called at ./Text/Pandoc/JSON.hs:111:48 in pandoc-types-1.17.5.4-5tHZ3B61A58JaKOMxwGQR4:Text.Pandoc.JSON
Error running filter /usr/bin/pandoc-citeproc:
Filter returned error status 1
Error: processing vignette 'probsNquantiles.Rmd' failed with diagnostics:
pandoc document conversion failed with error 83
--- failed re-building ‘probsNquantiles.Rmd’

SUMMARY: processing the following files failed:
  ‘parametrisation.Rmd’ ‘probsNquantiles.Rmd’

Error: Vignette re-building failed.
Execution halted