Initialize the ABC algorithm using the method of Sedki et al. (2013)
initSedki(y, neighbors, blocks, param = list(npart = 10000), priors = NULL)
| y | A vector of observed pixel data. |
|---|---|
| neighbors | A matrix of all neighbours in the lattice, one row per pixel. |
| blocks | A list of pixel indices, dividing the lattice into independent blocks. |
| param | A list of options for the ABC-SMC algorithm. |
| priors | A list of priors for the parameters of the model. |
A matrix containing SMC samples for the parameters of the Potts model.
Sedki, M.; Pudlo, P.; Marin, J.-M.; Robert, C. P. & Cornuet, J.-M. (2013) "Efficient learning in ABC algorithms" arXiv:1210.1388