GetDensityValues: Retrieve density values
In sppmix: Modeling Spatial Poisson and Related Point Processes

Description Usage Arguments Value Author(s) See Also Examples

This function operates on the point pattern and the realizations of a DAMCMC or BDMCMC fit (object damcmc_res or bdmcmc_res) and returns a plethora of information about the fit. When a bdmcmc_res is passed, only the realizations corresponding to the MAP number of components are used for calculations.

For examples see

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

1	GetDensityValues(fit)

fit

Object of class damcmc_res or bdmcmc_res.

A list containing the following components:

`Marginal`	the value of the Marginal (approximately)
`LogLikelihood`	the value of the LogLikelihood
`CompDensityAtXi`	the value of the component densities across all realizations and for each data point
`DensityAtXi`	the value of the mixture density across all realizations and for each data point
`EntropyMAP`	an approximation of the entropy of the distribution of the component indicators
`Density`	the joint density at each posterior realization, i.e., for each iteration of `ps`, `mus` and `sigmas`.

Sakis Micheas

normmix, est_mix_damcmc, est_mix_bdmcmc, rsppmix

# create the true mixture intensity surface
truesurf =normmix(ps=c(.2, .6,.2), mus=list(c(0.3, 0.3), c(0.7, 0.7),
 c(0.5, 0.5)),sigmas=list(.01*diag(2), .01*diag(2), .01*diag(2)),
 lambda=100,win=spatstat::square(1))
plot(truesurf)
# generate the point pattern, truncate=TRUE by default
genPP=rsppmix(truesurf,truncate=FALSE)
fit=est_mix_damcmc(pp = genPP, m = 3)
allvals=GetDensityValues(fit)
MAPest=GetMAPEst(fit,vals=allvals)
plot(MAPest,main="IPPP intensity surface of MAP estimates")