mhtboot (version 1.3.3)

plotchange: plotchange

Description

Plot the change function that is maximized to find the change point.

Usage

plotchange(zvec, rbuff = 25, h = 30, ...)

Arguments

zvec

vector of transformed order statistic of p values

rbuff

right buffer

h

window size

...

any graphical parameters passed to the plot function

Value

Nothing

Details

Currently there are two types of change functions supported. The difference between first difference series and the difference in gradients at each point. Both of these functions should have a theoretical maximum at the change point. We plot these two series side by side along with indicating the change point.

Examples

  
## Not run:       
# n = 50;m = 250;m0 = 20;
# sigeff = 1;
# Sigma <- 0.25*diag(m)
# X <- datgen(n,m,m0,sigeff,Sigma = Sigma)
# porder <- pboot.1sample(X=X,B=100,ncpus = 1)
# plotchange(porder[1,])
# ## End(Not run)