Skip to contents

Simulate a spatial Poisson random variable with a specific mean and covariance structure.

Usage

sprpois(
  spcov_params,
  mean = 0,
  samples = 1,
  data,
  randcov_params,
  partition_factor,
  ...
)

Arguments

spcov_params

An spcov_params() object.

mean

A numeric vector representing the mean. mean must have length 1 (in which case it is recycled) or length equal to the number of rows in data. The default is 0.

samples

The number of independent samples to generate. The default is 1.

data

A data frame or sf object containing spatial information.

randcov_params

A randcov_params() object.

partition_factor

A formula indicating the partition factor.

...

Additional arguments passed to sprnorm().

Value

If samples is 1, a vector of random variables for each row of data

is returned. If samples is greater than one, a matrix of random variables is returned, where the rows correspond to each row of data and the columns correspond to independent samples.

Details

The values of spcov_params, mean, and randcov_params are assumed to be on the link scale. They are used to simulate a latent normal (Gaussian) response variable using sprnorm(). This latent variable is the conditional mean used with dispersion to simulate a Poisson random variable.

Examples

spcov_params_val <- spcov_params("exponential", de = 0.2, ie = 0.1, range = 1)
sprpois(spcov_params_val, data = caribou, xcoord = x, ycoord = y)
#>  [1] 1 0 1 0 1 0 0 0 1 1 1 0 1 2 1 1 1 0 0 0 1 0 2 0 2 0 0 2 0 3
sprpois(spcov_params_val, samples = 5, data = caribou, xcoord = x, ycoord = y)
#>       1 2 3 4 5
#>  [1,] 2 3 1 1 0
#>  [2,] 5 3 0 0 1
#>  [3,] 0 2 0 0 0
#>  [4,] 0 4 0 0 0
#>  [5,] 0 1 2 2 1
#>  [6,] 1 1 0 2 1
#>  [7,] 1 0 5 3 0
#>  [8,] 2 0 1 3 0
#>  [9,] 0 0 0 1 1
#> [10,] 0 0 0 0 0
#> [11,] 1 0 0 0 1
#> [12,] 1 1 2 2 0
#> [13,] 0 0 2 3 1
#> [14,] 2 0 1 2 0
#> [15,] 2 0 1 0 3
#> [16,] 1 1 1 4 0
#> [17,] 2 1 2 2 0
#> [18,] 0 0 2 2 1
#> [19,] 1 0 1 1 4
#> [20,] 0 4 0 4 1
#> [21,] 0 0 0 2 1
#> [22,] 3 3 1 1 1
#> [23,] 1 0 0 0 0
#> [24,] 1 1 2 1 1
#> [25,] 0 4 2 2 0
#> [26,] 0 0 1 4 1
#> [27,] 3 2 2 2 0
#> [28,] 1 1 1 0 1
#> [29,] 1 1 1 4 4
#> [30,] 0 0 2 0 0