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

Description Usage Arguments Details Value See Also Examples

Generate random points according to the candidate probability distribution in the accept-reject method of generation of random variables applied to the distribution of the lambda parameter

1 2	rg(n, nMNO, nReg, fu, fv, flambda, relTol = 1e-06, nSim = 10000, nStrata = c(1, 100), verbose = FALSE)

`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
`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

The candidate distribution is a gamma distribution with parameters shape = nMNO + 1 and scale = λ^{*} / nMNO, where λ^{*} stands for the mode of the posterior distribution of the lambda parameter.

It is important to know that currently this function accepts only parameters for a single cell at a time. In case of interest for the candidate density function values for a set of cells, the user should program his/her own routine to apply this function to every cell.

rg generates n points according to the candidate distribution.

modeLambda, dlambda for related functions.

1
2
3

hist(rg(1e5, 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)), breaks = seq(1, 200, by = 2), main ='')