prgmNEW1

Disp "DATA IN COLS","OF [A]"
dim([A]) ¿θ1:Lθ1(1)¿N:Lθ1(2) ¿L:L-1¿K
{N,1}¿dim([B]):Fill(1,[B])
augment([B],[A]) ¿ [B]
[B]T
[B] ¿ [D]
seq([D](I,1),I,2,L+1)/N¿Lθ1
{L,L}¿dim([C]):Ans¿dim([E])
{N,1}¿dim([C])
Input "RESPONSE COL=",R
For(I,1,N):[A](I,R) ¿ [C](I,1):End
If R≤K:Then
For(J,R+1,L):For(I,1,N)
[B](I,J+1) ¿ [B](I,J)
End:End
[B]T
[B] ¿ [D]
End
{L,L}¿dim([D]):[D]—¿ [D]
{N,L}¿dim([B])
[B] T [C]:[D]Ans¿ [E]
Matr►list([E],Lθ1)
Disp "COEF Lθ1=":Pause Lθ1
[B][E] ¿ [B]
Matr►list([C],LY)
Matr►list([B],LYP)
mean(LY) ¿Y
sum((LY-Y)2
)¿T
sum((LY-LYP) 2¿E
LYP¿L5
LY-LYP¿L6
DelVar LY:DelVar LYP
N-L¿M:T-E¿R:R/K¿Q:E/M¿D:√(D) ¿S:Q/D¿F
S√(seq([D](I,I),I,1,L)) ¿Lθ2
Disp "STDEV Lθ2=":Pause Lθ2
Lθ1/Lθ2¿Lθ3
Disp "T-RATIO Lθ3=":Pause Lθ3
Disp "COEF P Lθ4="
1-2seq(tcdf(0,abs(Lθ3(I)),M),I,1,L) ¿Lθ4
Pause Lθ4
1-(N-1)*D/T¿A
ClrHome
Output(1,1,"S="
Output(1,9,S
Output(2,1,"R2
="
Output(2,9,R/T
Output(3,1,"R2
ADJ="
Output(3,9,A
Output(4,1,"REG DF="
Output(4,9,K
Output(5,1,"ERR DF="
Output(5,9,M
Output(6,1,"TOT DF="
Output(6,9,N-1
Pause
ClrHome
Output(1,1,"SS REG="
Output(1,9,R
Output(2,1,"SS ERR="
Output(2,9,E
Output(3,1,"SS TOT="
Output(3,9,T
Output(4,1,"MS REG="
Output(4,9,Q
Output(5,1,"MS ERR="
Output(5,9,D
Output(6,1,"F= "
Output(6,9,F
Output(7,1,"P-VAL="
1-Fcdf(0,F,K,M) ¿P
round(P,5) ¿P
Output(7,9,P
Pause
ClrHome
Lbl C
Menu("PREDICT Y?","YES",A,"QUIT",B)
Lbl A
Input "X{}=",Lθ5
Input "CONF. LEVEL=",C
{L,1}¿dim([C])
1¿ [C](1,1)
For(I,1,K):Lθ5(I) ¿ [C](I+1,1):End
[E]T
[C]
Ans(1,1) ¿Y
ClrHome
round(Y,5) ¿Y
Output(1,1,"Y HAT="
Output(1,8,Y
[C]T
[D][C]
Ans(1,1) ¿V
round(S√(V),5) ¿T
Output(2,1,"S(FIT)="
Output(2,8,T
TInterval 0,√(M+1),M+1,C
upper¿T
Output(3,1,"C.I.="
Y+TS√(V) ¿D
Y-TS√(V) ¿E
Output(4,3,E
Output(5,3,D
Output(6,1,"P.I.="
Y-TS√(1+V) ¿E
Y+TS√(1+V) ¿D
Output(7,3,E
Output(8,3,D
Pause
Goto C
Lbl B
DelVar [C]:DelVar [D]:DelVar [E]:DelVar Lθ5
ClrHome