GetPMEst: Retrieve the Surface of Posterior Means

In sppmix: Modeling Spatial Poisson and Related Point Processes

Description

The function first calculates the posterior means of the parameters of the components of the mixture intensity, based on a DAMCMC or BDMCMC fit. Then the surface of posterior means is calculated using the posterior means of the parameters. For a BDMCMC fit, the number of components should be specified, and all realizations with that number of components are gathered to calculate the posterior intensity surface.

For examples see

http://faculty.missouri.edu/~micheasa/sppmix/sppmix_all_examples.html#GetPMEst

Usage

GetPMEst(fit, num_comp = 1, burnin = floor(fit$L/10))

Arguments

fit	Object of class damcmc_res or bdmcmc_res.
num_comp	Number of components requested. The posterior will be calculated only based on the posterior realizations that have this many mixture components. If missing the realizations corresponding to the MAP number of components are returned. This parameter is ignored if fit is of class damcmc_res.
burnin	Number of initial realizations to discard. By default, it is 1/10 of the total number of iterations.

Value

An object of class intensity_surface.

Author(s)

Jiaxun Chen, Sakis Micheas, Yuchen Wang

See Also

est_mix_damcmc,est_mix_bdmcmc

Examples

fit <- est_mix_damcmc(pp = spatstat::redwood, m = 3)
post_intsurf <- GetPMEst(fit, burnin = 1000)
plot(post_intsurf)
fit <- est_mix_bdmcmc(pp = spatstat::redwood, m = 5)
post_intsurf <- GetPMEst(fit, num_comp = 4, burnin = 1000)
plot(post_intsurf)
post_fixed = FixLS_da(fit,approx=FALSE, plot_result = TRUE)
plot(GetPMEst(post_fixed))

sppmix documentation built on Jan. 13, 2021, 10:04 p.m.