;;
;;  SIMWN.MAC: macro to test BARTTEST and QTEST
;;
;;  INPUT: nreps(number of white noise realizations), 
;;         n (length of realizations), 
;;         m (for QTEST) 
;;         seed (seed for WN)
;;
;;  OUTPUT: pvb, pvq: vectors of length nreps with
;;                    pvalues for the 2 tests for the
;;                    nreps samples
;;
pause
;start

promptoff
x=wn(seed,n)
pvb=line(nreps,0,0)
pvq=line(nreps,0,0)

rep=1
;stloop
if(rep.le.nreps)

  list(rep)
  x=wn(0,n)

  rho=corr(x,n,0,n,1,r0,per)
  cper=cumsp(per,n)
  pvb[rep]=barttest(cper,n,n,b)

  rho=corr(x,n,m,0,1,r0,per)
  pvq[rep]=qtest(rho,m,0,0,n,qt)

  rep=rep+1
  goto(stloop)
endif

pvb=sort(pvb,nreps)
label(pvb)='Sorted pvalues for BARTTEST (n=#n#)'
pvq=sort(pvq,nreps)
label(pvq)='Sorted pvalues for QTEST (n=#n#)'

u=line(nreps,{-1./nreps},{1./nreps})
label(u)=' '
plot2(pvb,u,nreps,nreps,1,1,0,nreps,0,1)
plot2(pvq,u,nreps,nreps,1,1,0,nreps,0,1)