Skip to content

Interpretation / Usage of davies.test

2 messages · Katharina Mersmann, vito muggeo

#
dear Katie,
Since you are looking for exactly 1 breakpoint (namely you know the 
number of breakpoints), I suggest to use bootstrap restarting (default 
in segmented) with the rough value from davies.test() as a starting 
value, namely

o<-davies.test(reg1.2,~lagBYmean)
start.psi<-as.numeric(o$statistic)
a<-segmented(reg1.2,seg.Z=~lagBYmean ,psi=start.psi, 
control=seg.control(n.boot=100))

If you need some other feedbacks please contact me off list

best,
vito


Il 17/02/2014 15.27, Katharina Mersmann ha scritto: