Running ‘test_asis.R’ [4s/4s]
  Running ‘test_bayesqgcomp.R’ [3s/4s]
  Running ‘test_boot_ints.R’ [4s/5s]
  Running ‘test_boot_vs_noboot.R’ [5s/5s]
  Running ‘test_bootchooser.R’ [47s/49s]
  Running ‘test_cox_msmtest.R’ [6s/6s]
  Running ‘test_zi.R’ [20s/21s]
Running the tests in ‘tests/test_zi.R’ failed.
Complete output:
  > cat("# zi test\n")
  # zi test
  > library(qgcomp)
  > n=100
  > dat <- data.frame(y=rbinom(n, 1, 0.5)*rpois(n, 1.2), x1=runif(n), x2=runif(n), z=runif(n))
  > 
  > # should not cause error
  > ee = qgcomp.zi.noboot(f=y ~ z + x1 + x2 | z, expnms = c('x1', 'x2'), data=dat, q=2, dist="negbin")
  > qgcomp.zi.noboot(f=y ~ z + x1 + x2 | x1 + x2 + z, expnms = c('x1', 'x2'), data=dat, q=2, dist="negbin")
  Prob(Y ~ count):
  Scaled effect size (positive direction, sum of positive coefficients = 0.829)
  x2 
   1 
  
  Scaled effect size (negative direction, sum of negative coefficients = -0.0552)
  x1 
   1 
  
  Prob(Y ~ zero/count):
  Scaled effect size (positive direction, sum of positive coefficients = 10.7)
  x2 
   1 
  
  Scaled effect size (negative direction, sum of negative coefficients = -0.0138)
  x1 
   1 
  
  Mixture log(OR/RR) (Delta method CI):
  
  Prob(Y ~ count):
              Estimate Std. Error Lower CI  Upper CI Z value Pr(>|z|)
  (Intercept) -0.78158    0.37460 -1.51578 -0.047375 -2.0864  0.03694
  psi1         0.77419    0.51025 -0.22588  1.774258  1.5173  0.12920
  Prob(Y ~ zero):
              Estimate Std. Error Lower CI Upper CI Z value Pr(>|z|)
  (Intercept)  -11.575    125.694  -257.93   234.78 -0.0921   0.9266
  psi1          10.706    125.695  -235.65   257.06  0.0852   0.9321
  Warning message:
  In sqrt(diag(vc)[np]) : NaNs produced
  > 
  > qgcomp.zi.boot(f=y ~ z + x1 + x2 | z, expnms = c('x1', 'x2'), data=dat, q=4, B = 10, MCsize = 100, 
  +                dist="negbin")
  qgcomp.zi.boot function is still experimental. Please use with caution and be sure results are reasonable.
  
  Mixture log(OR/RR) (bootstrap CI):
  
  Prob(Y ~ count):
               Estimate Std. Error Lower CI Upper CI Z value Pr(>|z|)
  (Intercept) -0.709328   0.390252 -1.47421 0.055551 -1.8176  0.06912
  psi1         0.049487   0.245082 -0.43086 0.529838  0.2019  0.83998
  Prob(Y ~ zero):
              Estimate Std. Error Lower CI Upper CI Z value Pr(>|z|)
  (Intercept)  -2.5553     3.4677  -9.3519   4.2413 -0.7369   0.4612
  > qgcomp.zi.boot(f=y ~ z + x1 + x2 | z, expnms = c('x1', 'x2'), data=dat, q=4, B = 10, MCsize = 1000, 
  +                dist="negbin", msmcontrol = zimsm.fit.control(predmethod="catprobs"))
  qgcomp.zi.boot function is still experimental. Please use with caution and be sure results are reasonable.
  
  Error in solve.default(as.matrix(fit$hessian)) : 
    system is computationally singular: reciprocal condition number = 6.14272e-31
  Calls: qgcomp.zi.boot ... FUN -> zimsm.fit -> zeroinfl -> solve -> solve.default
  Execution halted