head	1.1;
access;
symbols;
locks
	wmc:1.1; strict;
comment	@# @;


1.1
date	2004.03.02.20.33.59;	author wmc;	state Exp;
branches;
next	;


desc
@@


1.1
log
@Initial revision
@
text
@pro lund_plot,y1,nop=nop,tis=tis,_extra=e

; Purpose: draw graphs of is-there-a-break-in-series a la Lund
; Category: statistics
; Method:
;		- call lund 
;		- plot the F_c function. Add on the line-of-significance (which may be off scale)
;		- plot the data (excluding, of course, missing data)
;	 	  add regression lines for the break		

fc=lund(y1,tis,a1=a1,mu1=mu1,a2=a2,mu2=mu2,c=c,ar=ar,mur=mur,t1=t1,t2=t2,tr=tr,y=y,_extra=e,sig=sig)

if (not keyword_set(nop)) then !p.multi=[0,1,2]

n=n_elements(y)

; fc is the f-type-statistic
plot,tr,fc,/xs
oplot,makerange(tr),[1,1]*lund_sig(n_elements(tr))
if (sig) then $
  print,'Max(F_c) of ',max(fc),' IS considered significant' $
else $
  print,'Max(F_c) of ',max(fc),' is NOT considered significant' 

; tr,y are the time and data filtered to remove missing data
plot,tr,y,/xs
; plots,psym=2,tr,y		; Add stars on the data points

print,a1,mu1
print,a2,mu2
print,'Break point index is after: ',c,' which is after time: ',tr(c)

oplot,tis,tis*ar+mur,th=3,li=2
oplot,t1,t1*a1+mu1,th=3
oplot,t2,t2*a2+mu2,th=3

end
@