\start83P\ \comment=Program file 09/02/03, 17:43 \protected=FALSE \name=POWER \file=C:\MYTIDATA\BACKUPS\POWER\POWER.TXT If O\!=\9999 Then {\(-)\2,2,.25}\->\\L\DELTA 9999\->\O End Disp "SERIES LENGTH" Prompt N Disp "INTERVENTION" Prompt T Disp "TEST TYPE" Disp "1: LOWER TAIL" Disp "2: TWO TAIL" Disp "3: UPPER TAIL" Input A Disp "LEVEL (EG .05)" Input B {A,B,N,T}\->\TEST Menu("NOISE TERM","AR",I,"IMA",J) Lbl I {0}\->\IMA Disp "PARAMETER" Prompt P N*(1-P)\^2\\->\A (N-T)*(1-P)\^2\+1\->\C C-P\->\B A/(AC-B\^2\)\->\S \root\(S*(1-P\^2\))\->\S S\^-1\\->\S {1,P,S}\->\AR Goto K Lbl J {0}\->\AR Disp "PARAMETER" Prompt P If P=1 Then N/((N+1-T)(T-1))\->\S {1,P,S}\->\IMA Goto K End N/(1-P)\^2\\->\A N-T\->\M (1+sum(seq(P^K,K,1,M)))\->\B B/(1-P)\->\B sum(seq(P^(2K),K,1,M))+1\->\C A/(AC-B\^2\)\->\S \root\((1+P\^2\)/S)\->\S {1,P,S}\->\IMA Lbl K Func RectGC:CoordOn GridOn:AxesOn LabelOff ExprOff If dim(\L\AR)=3 Then \L\AR(3)\->\S Else \L\IMA(3)\->\S End If \L\TEST(1)=1 Then 1-\L\TEST(2)\->\A invNorm(A)\->\Z "normalcdf(\(-)\1\EE\99,\(-)\Z-S*X)"\->\Y\1\ Else If \L\TEST(1)=2 Then 1-\L\TEST(2)/2\->\A invNorm(A)\->\Z "normalcdf(\(-)\1\EE\99,\(-)\Z-S*X)+1-normalcdf(\(-)\1\EE\99,Z-S*X)"\->\Y\1\ Else 1-\L\TEST(2)\->\A invNorm(A)\->\Z "1-normalcdf(\(-)\1\EE\99,Z-S*X)"\->\Y\1\ End End FnOff :FnOn 1 \L\DELTA(1)\->\Xmin \L\DELTA(2)\->\Xmax 0\->\Ymin:1\->\Ymax \L\DELTA(3)\->\Xscl .1\->\Yscl \L\DELTA(1)\->\TblStart \L\DELTA(3)\->\\Delta\Tbl DispTable DispGraph \stop83P\