## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----setup--------------------------------------------------------------------
library(arkhaia)

## -----------------------------------------------------------------------------
x1 <- c(2, 0, 10, 11, 5)
x2 <- c(1, 1, 17, 23, 3)
x3 <- c(2, 81, 11)
x4 <- c(5, 0, 1)
T1 <- matrix(c(x1, x2), ncol = 2)
T2 <- matrix(c(x3, x4), ncol = 2)
colnames(T1) <- c("S1T1", "S2T1")
rownames(T1) <- LETTERS[1:5]

colnames(T2) <- c("S1T2", "S2T2")
rownames(T2) <- LETTERS[3:5]

T1
T2

## -----------------------------------------------------------------------------
T1 <- t(T1)
T2 <- t(T2)

pois_rcens(T1, omit_zero = TRUE) # omit_zero = TRUE is default
pois_rcens(T1, omit_zero = FALSE)

## -----------------------------------------------------------------------------
set.seed(9)
trunc_pois(T1) # omitting zeros
trunc_pois(T2) # omitting zeros

## -----------------------------------------------------------------------------
VB_T1 <- VB_trunc_pois(T1, n_iter = 10^2)
VB_T2 <- VB_trunc_pois(T2, n_iter = 10^2)

## -----------------------------------------------------------------------------
VB_T1 <- VB_T1[VB_T1 > 0]
VB_T2 <- VB_T2[VB_T2 > 0]

VB_T1 <- sample(VB_T1, 10^3, replace = TRUE)
VB_T2 <- sample(VB_T2, 10^3, replace = TRUE)

## ----fig.height = 4, fig.width = 5, fig.align = 'center'----------------------
D <- VB_T1 - VB_T2
hist(D)