summary.bdmcmc_res: Summarize BDMCMC results

In sppmix: Modeling Spatial Poisson and Related Point Processes

Description

Prints a brief summary of the results of a BDMCMC fit.

For examples see

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

Usage

## S3 method for class 'bdmcmc_res'
summary(object, num_comp, burnin = floor(object$L/10),
  alpha = 0.05, dgt = 4, ...)

Arguments

object	Object of class bdmcmc_res.
num_comp	Number of components requested. Only the posterior realizations that have this many components will be returned. The function fails if the BDMCMC chain never visited this number of components. We can also pass a vector of integer values and present the posterior means summary. If this argument is missing, the MAP estimator is chosen by default.
burnin	Number of initial realizations to discard. By default, it is 1/10 of the total number of iterations.
alpha	Level alpha for the credible sets. Default is 0.05, for 95 sets of the mixture parameters.
dgt	Number of digits to use (formatting the output).
...	Additional arguments for the S3 method.

Author(s)

Sakis Micheas

See Also

rnormmix, to_int_surf, rsppmix, est_mix_bdmcmc, owin

Examples

# generate data
truemix<- rnormmix(m = 3, sig0 = .1, df = 5, xlim= c(0, 5), ylim = c(0, 5))
summary(truemix)
intsurf=to_int_surf(truemix, lambda = 100, win =spatstat::owin( c(0, 5),c(0, 5)))
pp1 <- rsppmix(intsurf = intsurf)# draw points
#Run BDMCMC and get posterior realizations
postfit=est_mix_bdmcmc(pp1,m=5)
#summary of the posterior results
summary(postfit)
summary(postfit, num_comp=2)
summary(postfit, num_comp=c(2,4))

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