safer handling of posterior structural zi probability in 0/0 form
Defaults pseudo lambda LRT now uses use_cont
Lots of changes. Primarily analysis is now centered around a pvlrt object that is created when the function pvlrt() is used. Various summary, plotting, printing methods for pvlrt are created.
Changes in lrt_zi_poisson():
removes the argument skip_null_omega_estimation. Null omega estimation is not performed as it is not necessary.
Changes in lrt_zi_poisson():
used analytic expression for pseudo lrt statistic for lambda, and a reduced algebraic expression for llik of omega. Both speed up computations.
Replaced omega_est_vec by omega_vec. omega_est_vec is now deprecated: will throw a warning if used instead of omega_vec.
Added the argument test_drug_idx to specify which columns to test.
use_gamma_smooth_omega: determines whether or not to use a gamma prior (smoothing) while estimating the omegas. Defaults to FALSE. NOTE: this changes the default behavior of the function. In previous versions a gamma prior was always used, which in current version can be achieved by setting use_gamma_smooth_omega = FALSE
added 4 new data sets: statin (all pts), statin1491, statin46 and gbca
added function r_contin_table_zip() to generate random contingency tables for adverse effect (across rows) and drug (across columns) combinations given row and column marginal totals, embedded signals, and possibly with zero inflation
Instead of storing random null tables and LRT statistics through lapply, they are now computed inside a for loop, and p-value is updated on the fly.
Removed pbapply from import; added progress
The first development version of the package
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.