Error(s) in re-building vignettes:
--- re-building ‘MLdist.Rmd’ using rmarkdown
Warning in xy.coords(x, y, xlabel, ylabel, log) :
  20 y values <= 0 omitted from logarithmic plot
Warning in xy.coords(x, y, xlabel, ylabel, log) :
  173 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
! LaTeX Error: Environment cslreferences undefined.

Error: processing vignette 'parametrisation.Rmd' failed with diagnostics:
Failed to compile parametrisation.tex. See https://yihui.name/tinytex/r/#debugging for debugging tips. See parametrisation.log for more info.
--- failed re-building ‘parametrisation.Rmd’

--- re-building ‘probsNquantiles.Rmd’ using rmarkdown
! LaTeX Error: Environment cslreferences undefined.

Error: processing vignette 'probsNquantiles.Rmd' failed with diagnostics:
Failed to compile probsNquantiles.tex. See https://yihui.name/tinytex/r/#debugging for debugging tips. See probsNquantiles.log for more info.
--- failed re-building ‘probsNquantiles.Rmd’

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

Error: Vignette re-building failed.
Execution halted