mimetic {ads} | R Documentation |

Simulates replicates of an observed univariate point pattern by stochastic optimization of its L-function properties.

mimetic(x,upto=NULL,by=NULL,prec=NULL,nsimax=3000,conv=50)

`x ` |
either a |

`upto ` |
(optional) maximum radius of the sample circles when |

`by ` |
(optional) interval length between successive sample circles radii when |

`prec ` |
precision of point coordinates generated during simulations when |

`nsimax ` |
maximum number of simulations allowed. By default the process stops after |

`conv ` |
maximum number of simulations without optimization gain (convergence criterion). |

Function `mimetic`

uses a stepwise depletion-replacement algorithm to generate a point pattern whose L-function is optimized with regards to an observed one, following the mimetic point process principle (Goreaud et al. 2004).
Four points are randomly deleted at each step of the process and replaced by new points that minimize the following cost function:||*Lobs(r) - Lsim (r)*||)^2. The simulation stops as soon as the cost function doesn't decrease
after `conv`

simulations or after a maximum of `nsimax`

simulations. The process apply to rectangular, circular or complex sampling windows (see `spp`

). There exist a `plot`

method that displays diagnostic
plots, i.e. the observed and simulated L-function, the simulated point pattern and the successive values of the cost function.

A list of class `"mimetic"`

with essentially the following components:

`call ` |
the function call. |

`fads ` |
an object of class |

`..r ` |
a vector of regularly spaced out distances corresponding to seq(by,upto,by). |

`..l ` |
a dataframe with 2 components: |

`.. ..obs` |
a vector of values of the L-function estimated for the initial observed pattern |

`.. ..sim` |
a vector of values of the L-function estimated for the simulated pattern |

`spp ` |
a object of class |

`theo ` |
a vector of theoretical values, i.e. Simpson |

`cost ` |
a vector of the successive values of the cost function. |

There are printing and plotting methods for `"mimetic"`

objects.

Goreaud F., Loussier, B., Ngo Bieng, M.-A. & Allain R. 2004. Simulating realistic spatial structure for forest stands: a mimetic point process. In Proceedings of Interdisciplinary Spatial Statistics Workshop, 2-3 December, 2004. Paris, France.

data(BPoirier) BP<-BPoirier ## Not run: performing point pattern analysis in a rectangle sampling window swr <- spp(BP$trees, win=BP$rect) plot(swr) ## Not run: performing the mimetic point process from "spp" object mimswr <- mimetic(swr, 20, 2) plot(mimswr) ## Not run: performing the mimetic point process from "fads" object mimkswr <- mimetic(kfun(swr, 20, 2)) plot(mimkswr)

[Package *ads* version 1.5-5 Index]