Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/BinSeg_one_func_minseglen.R

Implements the Binary Segmentation method for identifying changepoints in a given set of summary statistics for a specified cost function and penalty.

This function is called by `cpt.mean`

, `cpt.var`

and `cpt.meanvar`

when `method="BinSeg"`

. This is not intended for use by regular users of the package. It is exported for developers to call directly for speed increases or to fit alternative cost functions.

WARNING: No checks on arguments are performed!

1 |

`sumstat` |
A matrix containing the summary statistics of data within which you wish to find a changepoint. Currently assumes 3 columns and uses the number of rows as the length of the data +1 (initial value of 0). |

`pen` |
Default choice is 0, this should be evaluated elsewhere and a numerical value entered. This should be positive - this isn't checked but results are meaningless if it isn't. |

`cost_func` |
The friendly name of the cost function to be called in C. If using your own cost function, this must be the name of the C function to use. |

`shape` |
Only required for cost_func="Gamma",default is 1. Must be a positive value, this isn't checked. |

`minseglen` |
Positive integer giving the minimum segment length (no. of observations between changes), default is 2. No checks are performed on the input value so it could be larger than feasible to have changes in the data. |

`Q` |
The maximum number of changepoints to search for (positive integer). No checks are performed and so a number larger than allowed can be input. |

This function is used as a wrapper function to implement the Binary Segmentation algorithm in C. It simply creates the necessary worker vectors, ensures all inputs are the correct type, and passes everything to the C function.

This function is exported for developer use only. It does not perform any checks on inputs (other than type coersion) and is simply a wrapper function for the C code.

A list is returned with elements:

`cps` |
2xQ Matrix containing the changepoint positions on the first row and the test statistic on the second row in the order identified. |

`cpts` |
Ordered list of optimal number of changepoints ending with n. |

`op.cpts` |
The optimal number changepoint locations for the penalty supplied. |

`pen` |
Penalty used to find the optimal number of changepoints. |

Rebecca Killick

Binary Segmentation: Scott, A. J. and Knott, M. (1974) A Cluster Analysis Method for Grouping Means in the Analysis of Variance, *Biometrics* **30(3)**, 507–512

`cpt.mean`

,`cpt.meanvar`

,`plot-methods`

,`cpt`

1 | ```
#This function should only be used by developers, see its use in cpt.mean, cpt.var and cpt.meanvar.
``` |

```
Loading required package: zoo
Attaching package: 'zoo'
The following objects are masked from 'package:base':
as.Date, as.Date.numeric
Successfully loaded changepoint package version 2.2.2
NOTE: Predefined penalty values changed in version 2.2. Previous penalty values with a postfix 1 i.e. SIC1 are now without i.e. SIC and previous penalties without a postfix i.e. SIC are now with a postfix 0 i.e. SIC0. See NEWS and help files for further details.
```

changepoint documentation built on May 29, 2017, 3:19 p.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.