TSS.RESTREND (version 0.3.1)

CHOW: Chow test on detected breakpoints

Description

Takes the breakpoints detected by VPR.BFAST, finds the most significant one then tests it in both the residuals and the VPR.

Usage

CHOW(anu.VI, acu.RF, VI.index, breakpoints, acu.TM = NULL, sig = 0.05)

Arguments

anu.VI

The annual (Growing season) max VI. Must be a object of class 'ts' without NA's. if anu.VI=FALSE, it will be calculated from the CTSR.VI using AnMaxVI.

acu.RF

The optimal accumulated rainfall for anu.VI. Must be a object of class 'ts' without NA's and be of equal length and temporal range to anu.VI. if anu.RF=FALSE, it will be calculated from ACP.table usingthe AnnualClim.Cal

VI.index

the index of the CTSR.VI ts that the anu.VI values occur at. Must be the same length as anu.VI. NOTE. R indexs from 1 rather than 0. if VI.index=FALSE, it will be calculated from the CTSR.VI using AnMaxVI.

breakpoints

vector containing the breakpoints detected by VPR.BFAST (bkps)

acu.TM

The optimal accumulated rainfall for anu.TM. Must be a object of class 'ts' without NA's and be of equal length and temporal range to anu.TM. if anu.TM=FALSE, it will be calculated from ACT.table usingthe AnnualClim.Cal

sig

Significance of all the functions. defualt sig=0.05

Value

n.Method The method that the ts should be tested with. TSSRESTREND internal communication.

bp.summary Summary of the most signifcant breakpoint in the residuals and VPR. see sctest

allbp.index the Annual index of every breakpoint. Used by plot.TSSRESTREND

bpRESID.chow Chow test in the VPR residuals. See sctest

bpVPR.chow Chow test in the VPR. See sctest

Examples

# NOT RUN {
#Test the complete time series for breakpoints
VPRBFdem <- VPR.BFAST(segVPRCTSR$cts.NDVI, segVPRCTSR$cts.precip)
bp <- as.numeric(VPRBFdem$bkps)
#test the significance of the breakpoints
reschow <- CHOW(segVPR$max.NDVI, segVPR$acum.RF, segVPR$index, bp)
print(reschow)
# }
# NOT RUN {
# }