Ia s- J1 M7 P Rl WT \( d hz j k$ nk p q qj rn Ip Ie IgDATE_VECT P J It IvSCALE_VECT P I3 J JSCRATCH_VECT P I9 I8 SUM_X I> I= SUM_Y ID IC SUM_X_SQ IJ II SUM_X_Y IP IO NUM IU IT MEAN I[ IZ INTERCEPT I I` SLOPE I% I I XS P H I) I*YS P H H STRIDE H H POINTS H$ H# X_POS H) H( Y_POS H/ H. LEN_S H4 H3 CNT H: H9 JCNT H@ H> X_SCALE HF HD Y_SCALE HL HJ X_AXIS HR HP Y_AXIS HW HV RANGE H] H\ X_ORIGIN Hc Hb Y_ORIGIN Hi Hh X_LENGTH Ho Hn Y_LENGTH Hu Ht POINT Hz Hy TOP I H LINE_LEN I I Y_TOP I I
Y1 I I Y2 H M J J J J J% J+ Ia J X J J Y J J XL J J$ YL 0 J J* 0 XLABEL 0 J% J0 0 YLABEL J8 JI JZ Jk J| Kb Ki L, LY M HF J HL J
Hc J Hi J I K
HL Hi
K K K/ K4 K? KN KY 0 K P HF HL
0 K2 V Hc
KC Hi HF HL
H$ Ky HL L
L L Hi L
K} L !9 J+ P L H$ 9 L LABEL L1 L5 LK LO 0 L4 P L? H$ HF LC !y 0 LN T(B)[0,25] J+ L^ Lb M 0 La T[15,0] Lm M J% M HF Lw Hc L{ P HL M
!^ ! 0 M T(E) P
M% M+ M1 Ia M* MIN VAL M% M0 MAXVAL M+ M6 ARRAY N M{ N5 N: M< ( M% M+ H H HW 9 M> LABEL
Hc N H H N' M1 9 N) LABEL H/ NE H P NN Nq O O0 P Ho N] HW Ng P] Pd Q O| ( P P H H H] 9 O~ LABEL
Qj H H/ Hi R R$ R1 R R Hc RE P RH 9 RJ LABEL W> R[ Ra Rg Ia R MINVAL R[ Re MAXVAL Ra Rk ARRAY Rr S S. SH Sa T T VQ V^ Vt W W: HW R|
R[ S S R[ S HW S'
Ra 0 HR SW
Hu HW
Hc HR
Sz Rg 9 S| LABEL H/ HW T& T-
0 ! T> T^ Tr U$ VF Ho TM H4 TT
TQ H/ ! H HF U U H: U HW Ho
U( U; UH o Uz V
V? UY V UM Uj HR UV
p U] Un
P
t
0 Hz V Hi VL P V# V2 Hz V- ! H H
W7 9 W9 LABEL \VAL WY Wl X X1 XJ Xc [ [4 [G [w \ H] Wc
WH Wq X WH W{ H] X
WN 0 HR X@
Hu H] Hi HR
H/ H] Xz Y H4 Y
Y [
0 ! Y Y1 YD Yt [ Ho Y! H4 Y(
Y% H/ ! HS HL Yj
Y` H@ YW H] Ho
Z* ZO Z Z: HR Z'
p Z. Z>
P
t
0 Hz Zk Hc Zo P Zr Z Hz Z| ! H
[c H@ [Z H] Ho 9 \ LABEL d \ \ \" Ia \ MINVAL \ \! MAXVAL \ \' ARRAY \- \N \d 3 F ] j `i a a$ HR \G
\= \K \ \ 0f3 \ Ia \h \u ] HR \r
! \ ] Ip ] X 0 ] ] ]/ HR ]
!5 ]$ ], Ip ]) X P ]3 ]@ _. HR ]=
!q@ ]J ]Q Ip ]N X | ]T ]i ]^ ]f Ia ]c X 0 ]m ]z ^ HR ]w
"` ^ ^
P 0 ^ ^! ^6 HR ^
"6x ^+ ^3 Ip ^0 P P ^: ^G ^[ HR ^D
"Vx ^Q ^X Ip ^U P X ^` ^l _* HR ^i
"vl ^v ^~ Ip ^{ P x J Ia P 0 ' Ip J$ 0 0 HW = Ip H/ H `` H | 0 H: `y
Hc a `} HW J a8 aK ay c7 c< Ho aA Ia H an ad HF aZ H: Ho HW
b
bW b b b b% b= bC 0 b" P H b4 HL b8
0 0 b? V(W(P4I(C)))
bP bF bT Hi X b\ c
H bm bt c HL bq
! Ho c c x ! H c c! c4 HL c
!Q@ Ho c. c1 | p ! H4 cM Ia
cW cq Ia H4
cg c{ Ia H4
Ip H4
hj
d d$ d1 d> dK f f f2 fr g gO gT h he I3 d! I9 d. I> d; ID dH H/ f I d[ d~ e eK I3 dj I3 dt I H/ I9 e IJ e| I IP f) I3 IJ I[ fU fB f ID fK IP I9 I> fh IP I3 IU g g IJ I9 g I[ I3 I g7 g- IU HF I[ gA I 9 gC LABEL I gw gm IU gc I[ HF Hc g{ I 9 g} LABEL h h h= hB 0 h P h- I h# I[ I IU 0 h@ V hV I hL I[ I IU h hu Ia hy ARRAY H/ is hu i iP i i) hu H/ HL i3
H@ iF i< H] hu H/ k ix j j Ia i| X ix j Y j j
AXIS j jk H/ jg I j! jD j* j: ix H/ I H/ jM j] j H/ I% H/ j j{ HL j
0 nc k k
k k k Ia k MIN_VAL k k MAX_VAL k
k POINTS k k STRIDE k k# ORIGIN k+ kD l l5 lk m m> n n$ HR k:
k
k kI kU kb
k k_ 0 k ko kl ks !
kw HR
l l l l
l HR k
X k l. k l2 P l: lT l> lQ
lG HR k
| k ld k lh P k l{ k l
m k k m m$ k m! k m4 k k mB mO n 0
n HR k
n( nL n2 k
k n< k k k n\ k n` 0 p no o. oA nt o) H
0 0 p4 V(W(P4I(C)))
p9 Hz H) pR H
H) pi pl H p| H
! ! q q Ia q Y q# q3 qB q HL !
q HL !
qG qV q I ! q I ! qo qe Ia qi ARRAY H/ rd qe r rA r r qe H/ ! HF r# H: r7
r- HW qe H/ ! qe H/ ! ! s rh Ia rm ZERO rr rh s " HL s
H@ H]
! Entrypoints are:
! 0 -- draw box
! 2 -- scale X linear (and scale array)
! 3 -- scale Y linear
! 4 -- scale X log (and scale array)
! 5 -- scale Y log
! 6 -- scale X date (and scale array)
! 7 -- compute linear regression
! 8 -- scale entire Y array
! 9 -- return X and Y arrays
! 10 -- scale axis (internal)
! 11 -- print Y scale axis line
! 12 -- print X scale axis lie
Y value
! 14 -- scale entry Y array
ale 0 axis
DECLARE STRING STR
$DECLARE REAL STRIDE, POINTS, X_POS, Y_POS, LEN_S, CNT, JCNT, X_SCALE,
Y_SCALE, X_AXIS, Y_AXIS, RANGE, X_ORIGIN, Y_ORIGIN,
X_LENGTH, Y_LENGTH, POINT, TOP, LINE_LEN,
Y_TOP, Y1, Y2, ZERO
DECLARE REAL VECTOR XS, YS
DECLARE REAL SUM_X, SUM_Y, SUM_X_SQ, SUM_X_Y, NUM, MEAN,
INTERCEPT, SLOPE
DECLARE REAL VECTOR DATE_VECT, SCALE_VECT, SCRATCH_VECT
! Draw boxes
ENTRY 0 (X,Y, XL, YL, X_LABEL : STRING, Y_LABEL : STRING)
BEGIN
X_AXIS = X
Y_AXIS = Y
X_LENGTH = XL
Y_LENGTH = YL
Y_TOP = Y_AXIS - Y_LENGTH
# PRINT "P", LXY (X_AXIS, Y_AXIS), "V", RX(X_LENGTH), RY( - Y_LENGTH),
LX(X_AXIS), LY(Y_AXIS)
Y_POS = Y_AXIS - (Y_LENGTH + 25 * LENGTH (Y_LABEL)) / 2
PLOT LABEL 13 (Y_POS)
PRINT "P", LXY (X_AXIS - 100, Y_POS),
"T(B)[0,25]", QUOTE(Y_LABEL)
PRINT "T[15,0]",
CENTER (X_AXIS + X_LENGTH / 2, Y_AXIS + 60, X_LABEL, 15),
"T(E)"
END
! Scale X"s linearily
ENTRY 2 (MIN_VAL, MAX_VAL, ARRAY : VECTOR)
BEGIN
PLOT LABEL 10 (MIN_VAL, MAX_VAL, POINTS, STRIDE, X_ORIGIN)
X_SCALE = X_LENGTH / (POINTS * STRIDE)
PLOT LABEL 14 (ARRAY)
INCR CNT FROM 0 TO POINTS
BEGIN
POINT = X_ORIGIN + CNT * STRIDE
STR = CVT (POINT)
X_POS = X_AXIS + CNT * X_LENGTH / POINTS
IF CNT EQ 0 OR CNT EQ POINTS
THEN LINE_LEN = 0
ELSE LINE_LEN = Y_LENGTH - 2
PLOT LABEL 12
END
END
! Scale Y"s linearily
ENTRY 3 (MIN_VAL, MAX_VAL)
BEGIN
PLOT LABEL 10 (MIN_VAL, MAX_VAL, POINTS, STRIDE, Y_ORIGIN)
Y_SCALE = Y_LENGTH / (POINTS * STRIDE)
INCR CNT FROM 0 TO POINTS
BEGIN
POINT = Y_ORIGIN + CNT * STRIDE
STR = CVT (POINT)
Y_POS = Y_AXIS - CNT * Y_LENGTH / POINTS
IF CNT EQ 0 OR CNT EQ POINTS
THEN LINE_LEN = 0
ELSE LINE_LEN = X_LENGTH - 2
PLOT LABEL 11
END
END
! Scale LOG (X)
ENTRY 4 (MIN_VAL, MAX_VAL, ARRAY : VECTOR)
BEGIN
INVAL)
IF MIN_VAL < 0
X_ORIGIN = X_ORIGIN - 1
XVAL) + 1
RANGE = TOP - X_ORIGIN
X_SCALE = X_LENGTH / RANGE
PLOT LABEL 14 (ARRAY)
INCR CNT FROM X_ORIGIN TO TOP - 1
INCR JCNT FROM 1 TO 9
BEGIN
POINT = JCNT * 10 ** CNT
STR = CVT (POINT)
X_POS = X_AXIS + (LOG (POINT) - X_ORIGIN) * X_SCALE
IF LOG (POINT) EQ X_ORIGIN
THEN
LINE_LEN = 0
ELSE
! IF (RANGE < 4 AND (JCNT EQ 2 OR JCNT EQ 5)) OR JCNT EQ 1
THEN LINE_LEN = Y_LENGTH - 2
ELSE BEGIN
LINE_LEN = 8
STR = " "
END
PLOT LABEL 12
END
LINE_LEN = 0
POINT = 10 ** TOP
STR = CVT (POINT)
X_POS = X_AXIS + (LOG (POINT) - X_ORIGIN) * X_SCALE
PLOT LABEL 12
END
! Scale (LOG Y) -- NOTE: Parameters are log of values
ENTRY 5 (MIN_VAL, MAX_VAL)
BEGIN
INVAL)
IF MIN_VAL < 0
Y_ORIGIN = Y_ORIGIN - 1
XVAL) + 1
RANGE = TOP - Y_ORIGIN
Y_SCALE = Y_LENGTH / RANGE
INCR CNT FROM Y_ORIGIN TO TOP - 1
INCR JCNT FROM 1 TO 9
BEGIN
POINT = JCNT * 10 ** CNT
STR = CVT (POINT)
Y_POS = Y_AXIS - (LOG (POINT) - Y_ORIGIN) * Y_SCALE
IF LOG (POINT) EQ Y_ORIGIN
THEN
LINE_LEN = 0
ELSE
! IF (RANGE < 4 AND (JCNT EQ 2 OR JCNT EQ 5)) OR JCNT EQ 1
THEN LINE_LEN = X_LENGTH - 2
ELSE BEGIN
LINE_LEN = 8
STR = " "
END
PLOT LABEL 11
END
POINT = 10 ** TOP
STR = CVT (POINT)
Y_POS = Y_AXIS - (LOG (POINT) - Y_ORIGIN) * Y_SCALE
LINE_LEN = 0
PLOT LABEL 11
END
! Scale dates
ENTRY 6 (MIN_VAL, MAX_VAL, ARRAY : VECTOR)
BEGIN
RANGE = (MAX_VAL - MIN_VAL) * 1.05
DATE_EXP (MIN_VAL, DATE_VECT)
IF RANGE < 11
THEN SCALE_VECT (3) = 1
ELSE IF RANGE < 21
ECT (3) = 2
ELSE IF RANGE < 70
ECT (3) = 7
ELSE BEGIN
DATE_VECT (3) = 1
IF RANGE < 190
THEN SCALE_VECT (2) = 1
ELSE IF RANGE < 367
THEN SCALE_VECT (2) = 2
ELSE IF RANGE < 734
THEN SCALE_VECT (2) = 3
ELSE IF RANGE <1462
THEN SCALE_VECT (2) = 6
ELSE BEGIN
DATE_VECT (2) = 1
SCALE_VECT (1) = 1
END
END
X_ORIGIN = DATE_COM (DATE_VECT)
DATE_EXP (X_ORIGIN, SCRATCH_VECT)
POINTS = 0
WHILE DATE_COM (SCRATCH_VECT) < MAX_VAL
BEGIN
INCR CNT OVER SCRATCH_VECT
! SCRATCH_VECT (CNT) = SCRATCH_VECT (CNT) + SCALE_VECT (CNT)
POINTS = POINTS + 1
END
X_SCALE = X_LENGTH / (DATE_COM (SCRATCH_VECT) - X_ORIGIN)
PLOT LABEL 14 (ARRAY)
INCR CNT FROM 0 TO POINTS
BEGIN
POINT = DATE_COM (DATE_VECT)
X_POS = (POINT - X_ORIGIN) * X_SCALE + X_AXIS
IF CNT NE 0 AND CNT NE POINTS
PRINT "P", LXY (X_POS, Y_AXIS - 1), "V(W(P4I(C)))",
RY( - Y_LENGTH + 3)
PRINT CENTER (X_POS, Y_AXIS + 15, DATE_STR (POINT,0,6), 9),
CENTER (X_POS, Y_AXIS + 35, DATE_STR (POINT,7,4), 9)
INCR JCNT OVER DATE_VECT
DATE_VECT (JCNT) = DATE_VECT (JCNT) + SCALE_VECT (JCNT)
END
END
! Compute linear regression (least squares)
ENTRY 7
BEGIN
SUM_X = 0
SUM_Y = 0
SUM_X_SQ = 0
SUM_X_Y = 0
INCR CNT OVER XS
BEGIN
XS(CNT)
YS(CNT)
SUM_X_SQ = SUM_X_SQ + XS(CNT) ** 2
SUM_X_Y = SUM_X_Y + XS(CNT) * YS(CNT)
END
NUM = SIZE (XS)
MEAN = SUM_X / NUM
" SLOPE = (SUM_X_Y - MEAN * SUM_Y) / (SUM_X_SQ - MEAN * SUM_X)
INTERCEPT = (SUM_Y - SLOPE * SUM_X) / NUM
Y1 = X_AXIS * SLOPE + INTERCEPT
PLOT LABEL 13 (Y1)
Y2 = (X_AXIS + X_LENGTH) * SLOPE + INTERCEPT
PLOT LABEL 13 (Y2)
PRINT "P", LXY ((Y1 - INTERCEPT) / SLOPE, Y1),
"V", LXY ((Y2 - INTERCEPT) / SLOPE, Y2)
END
! Scale Y array
VECTOR)
INCR CNT OVER ARRAY
BEGIN
ARRAY (CNT) = Y_AXIS - Y_SCALE * (ARRAY (CNT) - Y_ORIGIN)
YS (CNT) = ARRAY (CNT)
END
! Return X and Y arrays
ENTRY 9 (X : VECTOR, Y : VECTOR, AXIS)
BEGIN
INCR CNT OVER XS
BEGIN
X(CNT)=XS(CNT)
Y(CNT)=YS(CNT)
END
AXIS = Y_AXIS - 1
END
for axis
ENTRY 10 (MIN_VAL, MAX_VAL, POINTS, STRIDE, ORIGIN)
BEGIN
RANGE = MAX_VAL - MIN_VAL
IF RANGE EQ 0 THEN
STRIDE = 1 ELSE
STRIDE = 10 ** INT (LOG (RANGE))
IF INT (RANGE / STRIDE) < 3
THEN STRIDE = STRIDE / 2
IF INT (RANGE / STRIDE) > 7
THEN STRIDE = STRIDE * 2
ORIGIN = STRIDE * (INT (MIN_VAL / STRIDE))
IF MIN_VAL < 0 THEN ORIGIN = ORIGIN - STRIDE
IF MIN_VAL NE 0
THEN IF MIN_VAL EQ ORIGIN THEN ORIGIN = ORIGIN - STRIDE
POINTS = 1 + INT (RANGE / STRIDE)
IF ORIGIN + POINTS * STRIDE LE MAX_VAL
THEN POINTS = POINTS + 1
END
! Print a scale line and label
ENTRY 11
BEGIN
IF LINE_LEN NE 0
# PRINT "P", LXY (X_AXIS + 1, Y_POS), "V(W(P4I(C)))", RX(LINE_LEN)
LEN_S = LENGTH (STR)
IF LEN_S NE 0
PRINT CENTER (X_AXIS - 10, Y_POS - 5, STR, 18)
END
! Print a scale line and label
ENTRY 12
BEGIN
IF LINE_LEN NE 0
% PRINT "P", LXY (X_POS, Y_AXIS - 1), "V(W(P4I(C)))", RY( - LINE_LEN)
LEN_S = LENGTH (STR)
IF LEN_S NE 0
PRINT CENTER (X_POS, Y_AXIS + 15, STR, 9)
END
! Clip a Y value to fit into (onto) box
ENTRY 13 (Y)
IF Y GT Y_AXIS
THEN Y = Y_AXIS
ELSE IF Y LT Y_TOP
THEN Y = Y_TOP
! Scale entire X array
: VECTOR)
INCR CNT OVER ARRAY
BEGIN
ARRAY (CNT) = X_AXIS + X_SCALE * (ARRAY (CNT) - X_ORIGIN)
XS (CNT) = ARRAY (CNT)
END
! Scale 0 axis
ENTRY 15 (ZERO)
BEGIN
ZERO = Y_AXIS + Y_SCALE * Y_ORIGIN
END
END_PLOT