Pharmacovigilance Analysis — phv

Pharmacovigilance, also known as drug safety. In the context of pharmacovigilance studies, disproportionality analysis primarily served as a tool to evaluate possible association between a specific adverse event and a particular drug which can then be investigated through clinical assessment of individual case reports.

phv_signal(
  a,
  b,
  c,
  d,
  methods = NULL,
  alpha = 0.05,
  correct = TRUE,
  n_mcmc = 100000L,
  alpha1 = 0.5,
  alpha2 = 0.5,
  theta_init = NULL,
  squashing = TRUE,
  BPPARAM = SerialParam()
)

phv_ror(a, b, c, d, alpha = 0.05)

phv_prr(a, b, c, d, alpha = 0.05)

phv_chisq(a, b, c, d, correct = TRUE, BPPARAM = SerialParam())

phv_fisher(a, b, c, d, alpha = 0.05, BPPARAM = SerialParam())

phv_bcpnn_norm(a, b, c, d, alpha = 0.05)

phv_bcpnn_mcmc(
  a,
  b,
  c,
  d,
  alpha = 0.05,
  n_mcmc = 100000L,
  BPPARAM = SerialParam()
)

phv_obsexp_shrink(
  a,
  b,
  c,
  d,
  alpha = 0.05,
  alpha1 = 0.5,
  alpha2 = 0.5,
  n_mcmc = 100000L,
  BPPARAM = SerialParam()
)

phv_ebgm(a, b, c, d, alpha = 0.05, theta_init = NULL, squashing = TRUE)

Arguments

a: also referred to as n11 as this is the count of event of interest under exposure of interest.
b: also referred to as n10 as this is the count of not event of interest under exposure of interest.
c: also referred to as n01 as this is the count of event of interest under not exposure of interest.
d: also referred to as n00 as this is the count of not event of interest under not exposure of interest.
methods: An atomic character, specifies the method used to signal mining. Currently, only "ror", "prr", "chisq", "bcpnn_norm", "bcpnn_mcmc", "obsexp_shrink", "fisher", and "ebgm" are supported. If NULL, all supported methods will be used.
alpha: Level of significance, for construction of the confidence intervals.
correct: A bool indicating whether to apply Yates's continuity correction when computing the chi-squared statistic.
n_mcmc: Number of MCMC simulations per (a,b,c,d)-tuple to calculate confidence intervals.
alpha1: Numerator shrinkage parameter >=0, default 0.5.
alpha2: Denominator shrinkage parameter >=0, default 0.5.
theta_init: A data frame of initial hyperparameter guesses with columns ordered as: alpha1, beta1, alpha2, beta2, P. See openEBGM::autoHyper
squashing: A bool, whether do automated data squashing. If any zeros found in a, will always be TRUE.
BPPARAM: An optional BiocParallelParam instance defining the parallel back-end to be used during evaluation.

Value

A data.table with columns of estimated value and it's confidence interval (ci_low and ci_high). Estimated column are as follows:

phv_ror: reporting odds ratio (ror).
phv_prr: proportional reporting ratio (prr). Signal defined as a prr of at least 2, chi-squared with Yates's correction of at least 4 and a >= 3. An equivalent alternative to chi-squared is to calculate a confidence interval around the prr.
phv_bcpnn_norm: information component (ic).
phv_bcpnn_mcmc: information component (ic).
phv_obsexp_shrink: observed to expected ratio (oe_ratio).
phv_ebgm: Empirical Bayes Geometric Mean (ebgm).

Details

Note that the a, b, c, d inputs can be an atomic vectors of equal length, for which the function will perform the calculations for each individual (a,b,c,d)-tuple moving across the vectors.

It is assumed that the contingency table under consideration has drugs/exposures in the rows and outcomes/events in the columns. See contingency table section.

We use the distinct patient count method to obtain the frequency counts of patients exposed to each interested drug, those reporting interested event. As illustrated in the Contingency table, n equals the total number of patients in the database, n11 is the number of patients with exposure to the interested drug during the model period and reporting interested events, n10 is the number of patients that have used the interested drug but did not experience interested event during any of the model periods associated with the drug, n01 is the number of patients that did not use the interested drug but experienced interested event, and n00 is the number of patients that were not exposed to the interested drug and did not report interested condition.

Contingency table

	ADR of interest	Other ADRs	Total
Drug of interest	a=n11	b=n10	a+b=n1.
Other drugs	c=n01	d=n00	c+d=n0.
Total	a+c=n.1	b+d=n.0	a+b+c+d=n

phv_obsexp_shrink

The observed to expected (OE) ratio with approximate confidence intervals are constructed on the log2 scale as outlined in Norén et al. (2013).

Expected value was estimated by (a + b) / (a + b + c + d) * (a + c).

The OE ratio with shrinkage estimates is calculated as (O + alpha1) / (E + alpha2).

If (O + alpha1) < 1, then the exact uncertainty limits should be used. That is the confidence intervals as implemented in phv_bcpnn_mcmc (Norén et al., 2013).

log2(OE) approximates the Bayesian confidence propagation neural network information component (IC) with reasonable accuracy when alpha1 = alpha2 = 0.5 (Norén et al., 2013).

phv_ebgm

An implementation of the Gamma-Poisson Shrinker (GPS) model for identifying unexpected counts in large contingency tables using an empirical Bayes approach. The Empirical Bayes Geometric Mean (EBGM) and quantile scores are obtained from the GPS model estimates. The GPS was proposed by DuMouchel as a signal detection tool for large frequency tables with both observed (O) and expected (E) counts for each drug-outcome pair. It assumes the observed count of any drug-outcome pair follows the Poisson distribution.

For each drug-outcome pair, the primary parameter of interest was the risk ratio. Rather than using the observed over expected (O/E), GPS uses the empirical Bayesian geometric mean (EBGM) posterior distribution of the risk ratio and the surrounding confidence interval for each drug-outcome pair to identify statistical signals of excess risk. To prevent spurious false positives due to implausibly high risk ratios, GPS implements a Bayesian framework that “shrinks” O/E estimates towards a value which is close to the average O/E values for all drug-event pairs at each level of granularity.

References

David Olaleye, SAS Institute Inc. (2019), Real-World Evidence and Population Health Analytics: Intersection and Application, https://support.sas.com/resources/papers/proceedings19/3361-2019.pdf

Evans, S.J.W., Waller, P.C. and Davis, S. (2001), Use of proportional reporting ratios (PRRs) for signal generation from spontaneous adverse drug reaction reports. Pharmacoepidem. Drug Safe., 10: 483-486. https://doi.org/10.1002/pds.677

Norén GN, Hopstadius J, Bate A. Shrinkage observed-to-expected ratios for robust and transparent large-scale pattern discovery. Statistical methods in medical research. 2013 Feb;22(1):57-69.

https://journal.r-project.org/archive/2017/RJ-2017-063/RJ-2017-063.pdf

Examples

phv_signal(122, 1320, 381, 31341)
#> ℹ Running `phv_ror()`
#> ℹ Running `phv_prr()`
#> ℹ Running `phv_chisq()`
#> ℹ Running `phv_bcpnn_norm()`
#> ℹ Running `phv_bcpnn_mcmc()`
#> ℹ Running `phv_obsexp_shrink()`
#> ℹ Running `phv_fisher()`
#> ℹ Running `phv_ebgm()`
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: ✖ Cannot estimate hyperparameter using openEBGM
#> ℹ Using `nlminb()` directly to estimate hyperparameter
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#>    expected      ror ror_ci_low ror_ci_high      prr prr_ci_low prr_ci_high
#>       <num>    <num>      <num>       <num>    <num>      <num>       <num>
#> 1: 21.87088 7.602804   6.155362    9.390613 7.044175   5.784965    8.577477
#>       chisq  chisq_pvalue bcpnn_norm_ic bcpnn_norm_ic_ci_low
#>       <num>         <num>         <num>                <num>
#> 1: 481.7814 8.752383e-107       2.41958             2.126359
#>    bcpnn_norm_ic_ci_high bcpnn_mcmc_ic bcpnn_mcmc_ic_ci_low
#>                    <num>         <num>                <num>
#> 1:              2.712802      2.453088             2.224909
#>    bcpnn_mcmc_ic_ci_high oe_ratio oe_ratio_ci_low oe_ratio_ci_high odds_ratio
#>                    <num>    <num>           <num>            <num>      <num>
#> 1:              2.660812 2.453088        2.153455         2.669561   7.603072
#>    odds_ratio_ci_low odds_ratio_ci_high fisher_pvalue     ebgm ebgm_ci_low
#>                <num>              <num>         <num>    <num>       <num>
#> 1:          6.102962           9.418401  3.925584e-56 4.416748        3.69
#>    ebgm_ci_high
#>           <num>
#> 1:         5.25
phv_signal(122, 1320, 381, 31341, "ror")
#> ℹ Running `phv_ror()`
#>    expected      ror ror_ci_low ror_ci_high
#>       <num>    <num>      <num>       <num>
#> 1: 21.87088 7.602804   6.155362    9.390613
phv_ror(122, 1320, 381, 31341)
#>         ror   ci_low  ci_high
#>       <num>    <num>    <num>
#> 1: 7.602804 6.155362 9.390613
phv_signal(122, 1320, 381, 31341, "prr")
#> ℹ Running `phv_prr()`
#>    expected      prr prr_ci_low prr_ci_high
#>       <num>    <num>      <num>       <num>
#> 1: 21.87088 7.044175   5.784965    8.577477
phv_prr(122, 1320, 381, 31341)
#>         prr   ci_low  ci_high
#>       <num>    <num>    <num>
#> 1: 7.044175 5.784965 8.577477
phv_signal(122, 1320, 381, 31341, "chisq")
#> ℹ Running `phv_chisq()`
#>    expected    chisq  chisq_pvalue
#>       <num>    <num>         <num>
#> 1: 21.87088 481.7814 8.752383e-107
phv_chisq(122, 1320, 381, 31341)
#>       chisq        pvalue
#>       <num>         <num>
#> 1: 481.7814 8.752383e-107
phv_signal(122, 1320, 381, 31341, "bcpnn_norm")
#> ℹ Running `phv_bcpnn_norm()`
#>    expected bcpnn_norm_ic bcpnn_norm_ic_ci_low bcpnn_norm_ic_ci_high
#>       <num>         <num>                <num>                 <num>
#> 1: 21.87088       2.41958             2.126359              2.712802
phv_bcpnn_norm(122, 1320, 381, 31341)
#>         ic   ci_low  ci_high
#>      <num>    <num>    <num>
#> 1: 2.41958 2.126359 2.712802
phv_signal(122, 1320, 381, 31341, "bcpnn_mcmc")
#> ℹ Running `phv_bcpnn_mcmc()`
#>    expected bcpnn_mcmc_ic bcpnn_mcmc_ic_ci_low bcpnn_mcmc_ic_ci_high
#>       <num>         <num>                <num>                 <num>
#> 1: 21.87088      2.453088             2.224909              2.660832
phv_bcpnn_mcmc(122, 1320, 381, 31341)
#>          ic   ci_low  ci_high
#>       <num>    <num>    <num>
#> 1: 2.453088 2.224876 2.660825
phv_signal(122, 1320, 381, 31341, "obsexp_shrink")
#> ℹ Running `phv_obsexp_shrink()`
#>    expected oe_ratio oe_ratio_ci_low oe_ratio_ci_high
#>       <num>    <num>           <num>            <num>
#> 1: 21.87088 2.453088        2.153455         2.669561
phv_obsexp_shrink(122, 1320, 381, 31341)
#>    oe_ratio   ci_low  ci_high
#>       <num>    <num>    <num>
#> 1: 2.453088 2.153455 2.669561
phv_signal(122, 1320, 381, 31341, "fisher")
#> ℹ Running `phv_fisher()`
#>    expected odds_ratio odds_ratio_ci_low odds_ratio_ci_high fisher_pvalue
#>       <num>      <num>             <num>              <num>         <num>
#> 1: 21.87088   7.603072          6.102962           9.418401  3.925584e-56
phv_fisher(122, 1320, 381, 31341)
#>    odds_ratio   ci_low  ci_high       pvalue
#>         <num>    <num>    <num>        <num>
#> 1:   7.603072 6.102962 9.418401 3.925584e-56
phv_signal(122, 1320, 381, 31341, "ebgm")
#> ℹ Running `phv_ebgm()`
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: ✖ Cannot estimate hyperparameter using openEBGM
#> ℹ Using `nlminb()` directly to estimate hyperparameter
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#>    expected     ebgm ebgm_ci_low ebgm_ci_high
#>       <num>    <num>       <num>        <num>
#> 1: 21.87088 4.416748        3.69         5.25
phv_ebgm(122, 1320, 381, 31341)
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN replaced by maximum positive value
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: ✖ Cannot estimate hyperparameter using openEBGM
#> ℹ Using `nlminb()` directly to estimate hyperparameter
#> Warning: NA/NaN function evaluation
#> Warning: NaNs produced
#> Warning: NaNs produced
#> Warning: NA/NaN function evaluation
#>        ebgm ci_low ci_high
#>       <num>  <num>   <num>
#> 1: 4.416748   3.69    5.25