R: Intersection between model region and a given interval.

intersect.modal.region {circular} | R Documentation

## Intersection between model region and a given interval.

### Description

Find an estimates of the probability of the intersection between a modal region and a given interval.

### Usage

    intersect.modal.region(x, ...)
    ## Default S3 method:
    intersect.modal.region(x, ...)
    ## S3 method for class 'circular'
    intersect.modal.region(x, breaks, z = NULL, q = 0.95, bw,
      adjust = 1, type = c("K", "L"), kernel = c("vonmises", "wrappednormal"),
      na.rm = FALSE, step = 0.01, eps.lower = 10^(-4), eps.upper = 10^(-4), ...)

### Arguments

`x` | numeric or an object of class `circular`.
`breaks` | a matrix with two columns. Each row specifies a sub-interval.
`z` | numeric or object of class `circular`. The grid were the kernel density estimate will be evaluated. If `NULL` equally spaced points in the interval [0,2*pi) with step `step`.
`q` | numeric in the interval [0,1]. The quantile of the modal region.
`bw` | the smoothing bandwidth to be used. When the `kernel` is `vonmises` the bandwidth is equal to the concentration parameter.
`adjust` | the bandwidth used is actually `adjust*bw`. This makes it easy to specify values like “half the default bandwidth”.
`type` | Not Yet Used.
`kernel` | a character string giving the smoothing kernel to be used. This must be one of `"vonmises"` or `"wrappednormal"`, that are kernels of `type` `"K"`.
`na.rm` | logical; if `TRUE`, missing values are removed from `x`. If `FALSE` any missing values cause an error.
`step` | numeric. Used in the construction of the regular grid `z`.
`eps.lower,eps.upper` | the cut point in the density is searched in the interval [min(density)*(1+eps.lower),max(density)*(1-eps.upper)].
`...` | further arguments passed to the next methods.

### Details

Only the version for circular data is actually implemented.

### Value

For the circular method a list with the following three components

`tot` | the total area.
`areas` | information for each subinterval.
`breaks` | the extremes of each subinterval.

### Author(s)

Claudio Agostinelli

### See Also

`modal.region`

### Examples

      x <- rvonmises(100, circular(pi), 10)
      res <- intersect.modal.region(x, breaks=circular(matrix(c(pi,pi+pi/12,
        pi-pi/12, pi), ncol=2, byrow=TRUE)), bw=50)
      res$tot

      x <- rvonmises(100, circular(0), 10)
      res <- intersect.modal.region(x, breaks=circular(matrix(c(pi,pi+pi/12),
        ncol=2)), bw=50)
      res$tot

      res <- intersect.modal.region(x, breaks=circular(matrix(c(pi/12,
        2*pi-pi/12), ncol=2, byrow=TRUE)), bw=50)
      res$tot

[Package circular version 0.4-93 Index]