# EstimateIP: Parameter Estimation of the Inverse-Power Type Model In NScluster: Simulation and Estimation of the Neyman-Scott Type Spatial Cluster Models

## Description

Parameter estimation of the Inverse-power type model via numerical calculation of the Ripley's K-function.

## Usage

 ```1 2``` ``` EstimateIP(xy.points, pars, eps = 0.001, uplimit = 0.3, skip = 1, process.report = 0, plot = TRUE) ```

## Arguments

 `xy.points` a matrix containing the coordinates `(x,y)` of points in a unit square: W=[0,1]*[0,1]. `pars` a named vector of containing the initial guess of the model parameters (`mu`, `nu`, `p`, `c`), where `mu` is an intensity of parents, `nu` is an expected number of descendants for each parent, `p` is the decay order and `c` is the scaling parameter. `eps` the optimization procedure is iterated at most 1000 times until `process2\$stderr` becomes smaller than `eps`. `uplimit` upper limit value in place of ∞ in the integral in distribution function. `skip` the variable for the fast likelihood but rough approximation of the initial estimates. The `skip` calculate the Palm intensity function in the log-likelihood function for every `skip`-th r(i,j) in the ordered distances of the pairs i and j. `process.report` the level of reporting the process of minimizing. Allowed values are as follows: 0no report (default). 1output the process of minimizing the negative Palm log-likelihood function until the values converge to MPLEs. (`process1`) 2output the process of optimizing by the simplex with the normalized parameters. (`process2`) 3output both processes. `plot` logical. If `TRUE` (default), the process of optimizing by the simplex with the normalized parameters is plotted.

## Value

 `mple` MPLE (maximum Palm likelihood estimate). `process1` a list with following components. (Only returned if `process.report` = 1 or 3.) cflg1 (="update") or -1 (="testfn"), where "update" indicates that -log L value has attained the minimum so far, otherwise not. logl.palmthe minimized -log L in the process of minimizing the negative Palm log-likelihood function. mplescorresponding MPLEs. `process2` a list with following components. (Only returned if `process.report` = 2 or 3.) logl.simplexthe minimized -log L by the simplex method. stderrthe standard deviations. pa.normalthe normalized variables corresponding the initial estimates.

## Note

`EstimateIP` and `EstimateTypeA` have to use numerical integration and difference to compute the functions, which need very long c.p.u. time in the minimization procedure.

## References

U. Tanaka, Y. Ogata and K. Katsura, Simulation and estimation of the Neyman-Scott type spatial cluster models, Computer Science Monographs No.34, 2008, 1-44. The Institute of Statistical Mathematics.

## Examples

 ``` 1 2 3 4 5 6 7 8 9 10 11``` ```## simulation pars <- c(mu = 50.0, nu = 30.0, p = 1.5, c = 0.005) z <- SimulateIP(pars, seed = 353) ## estimation ## need very long c.p.u time in the minimization procedure ## Not run: init.pars <- c(mu = 55.0, nu = 35.0, p = 1.0, c = 0.01) EstimateIP(z\$offspring\$xy, init.pars, skip = 100) ## End(Not run) ```

NScluster documentation built on March 19, 2018, 9:03 a.m.