rN0: Generation of random deviates of the posterior distribution...
In MobilePhoneESSnetBigData/pestim: Population Estimations Using Mobile Phone Data

Description Usage Arguments Details Value See Also Examples

Generate random points according to the posterior probability distribution of the number of individuals in the hierarchical model.

1
2
3

rN0(n, nMNO, nReg, fu, fv, flambda, scale = 1, relTol = 1e-06,
  nSim = 10000, nStrata = c(1, 100), verbose = FALSE,
  nThreads = RcppParallel::defaultNumThreads())

`n`	number of values to generate
`nMNO,`	nReg non-negative integer vectors with the number of individuals detected in each cell according to the network operator and the register
`fu,`	fv named lists with the prior marginal distributions of the two-dimensional points for the Monte Carlo integration
`flambda`	named list with the prior distribution of the lambda parameter
`scale`	numeric vector with the scale to count the number of individuals. Default value is 1
`relTol`	relative tolerance in the computation of the `kummer` function. Default value is `1e-6`
`nSim`	number of two-dimensional points to generate to compute the integral. Default value is `1e4`
`nStrata`	integer vector of length 2 with the number of strata in each dimension. Default value is `c(1, 1e2)`
`verbose`	logical (default `FALSE`) to report progress of the computation
`nThreads`	number (default the number of all cores, including logical cores) to use for computation

The posterior distribution is a Poisson distribution with parameter lambda * scale, where the values of lambda are generated with the function rlambda.

The prior distributions are specified as named lists where the first component of each list must be the name of distribution ('unif', 'triang', 'degen', 'gamma') and the rest components must be named according to the name of the parameters of the random generator of the corresponding distribution according to:

unif: xMin, xMax for the minimum, maximum of the sampled interval.
degen: x0 for the degenerate value of the random variable.
triang: xMin, xMax, xMode for minimum, maximum and mode (see qtriang).
gamma: scale and shape with the same meaning as in rgamma.

rN0 generates n points according to the posterior distribution. The function returns a data.table with these points (under the column N0) together with the additional variables:

The common length of nMNO and nReg identifies the number of territorial cells in which the number of individuals detected by the telecommunication network and official data. The column cellID identifies these territorial cells.
The different values of the generated values of lambda are returned under the column lambda.
The inputs nMNO and nReg are also included in the output data.table in columns under the same name.

rlambda, rg, rNt for related functions.

# It takes a couple of minutes
hist(rN0(500, nMNO = 20, nReg = 115, fu = list('unif', xMin = 0.3, xMax = 0.5),
        fv = list('unif', xMin = 100, xMax = 120),
        flambda = list('gamma', shape = 11, scale = 12))$N0,
        breaks = seq(1, 200, by = 1), main ='', xlab = 'number of individuals')