# Strauss

##### The Strauss Point Process Model

Creates an instance of the Strauss point process model which can then be fitted to point pattern data.

##### Usage

`Strauss(r)`

##### Arguments

- r
- The interaction radius of the Strauss process

##### Details

The (stationary) Strauss process with interaction radius $r$ and parameters $\beta$ and $\gamma$ is the pairwise interaction point process in which each point contributes a factor $\beta$ to the probability density of the point pattern, and each pair of points closer than $r$ units apart contributes a factor $\gamma$ to the density.

Thus the probability density is $$f(x_1,\ldots,x_n) = \alpha \beta^{n(x)} \gamma^{s(x)}$$ where $x_1,\ldots,x_n$ represent the points of the pattern, $n(x)$ is the number of points in the pattern, $s(x)$ is the number of distinct unordered pairs of points that are closer than $r$ units apart, and $\alpha$ is the normalising constant.

The interaction parameter $\gamma$ must be less than
or equal to $1$
so that this model describes an ``ordered'' or ``inhibitive'' pattern.
The nonstationary Strauss process is similar except that
the contribution of each individual point $x_i$
is a function $\beta(x_i)$
of location, rather than a constant beta.
The function `ppm()`

, which fits point process models to
point pattern data, requires an argument
of class `"interact"`

describing the interpoint interaction
structure of the model to be fitted.
The appropriate description of the Strauss process pairwise interaction is
yielded by the function `Strauss()`

. See the examples below.
Note the only argument is the interaction radius `r`

.
When `r`

is fixed, the model becomes an exponential family.
The canonical parameters $\log(\beta)$
and $\log(\gamma)$
are estimated by `ppm()`

, not fixed in
`Strauss()`

.

##### Value

- An object of class
`"interact"`

describing the interpoint interaction structure of the Strauss process with interaction radius $r$.

##### See Also

##### Examples

```
Strauss(r=0.1)
# prints a sensible description of itself
data(cells)
ppm(cells, ~1, Strauss(r=0.07))
# fit the stationary Strauss process to `cells'
ppm(cells, ~polynom(x,y,3), Strauss(r=0.07))
# fit a nonstationary Strauss process with log-cubic polynomial trend
```

*Documentation reproduced from package spatstat, version 1.15-2, License: GPL (>= 2)*