PoissonBlocRnd {bikm1}

PoissonBlocRnd function for contingency data simulation

Description

Produce a simulated data matrix generated under the Poisson Latent Block Model.

Usage

PoissonBlocRnd(J,K,theta)

Arguments

J

a positive integer specifying the number of expected rows.

K

a positive integer specifying the number of expected columns.

theta

a list specifying the model parameters:

rho_h: a vector specifying the row mixing proportions.

tau_l: a vector specifying the column mixing proportions.

gamma_hl: a matrix specifying the distribution parameter.

Value

a list including the arguments:

x: simulated contingency data matrix.

xrow: numeric vector specifying row partition.

xcol: numeric vector specifying column partition.

Examples

require(bikm1)
J=200
K=120
h=3
l=2
theta=list()
theta$rho_h=1/h *matrix(1,h,1)
theta$tau_l=1/l *matrix(1,l,1)
theta$gamma_hl=matrix(c(1, 6,4, 1, 7, 1),ncol=2)
data=PoissonBlocRnd(J,K,theta)