/ Chip 2002 June / Chip_2002-06_cd1.bin / zkuste / derive / download / Setup.exe / %MAINDIR% / Users / NonlinearRegression.dfw

< prev

next >

view JSON data

view as text

open on a Mac

open on a PC

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 44 45 52 49 56 45 20 66 | 6f 72 20 57 69 6e 64 6f |DERIVE f|or Windo|
|00000010| 77 73 20 76 65 72 73 69 | 6f 6e 20 35 2e 78 78 20 |ws versi|on 5.xx |
|00000020| 44 66 57 20 66 69 6c 65 | 20 73 61 76 65 64 20 6f |DfW file| saved o|
|00000030| 6e 20 30 32 20 4e 6f 76 | 20 32 30 30 31 0d 0a 1a |n 02 Nov| 2001...|
|00000040| 4c 00 00 00 41 44 4a 4f | 49 4e 54 28 61 29 3a 3d |L...ADJO|INT(a):=|
|00000050| 56 45 43 54 4f 52 28 56 | 45 43 54 4f 52 28 43 4f |VECTOR(V|ECTOR(CO|
|00000060| 46 41 43 54 4f 52 28 61 | 2c 20 6d 5f 2c 20 6e 5f |FACTOR(a|, m_, n_|
|00000070| 29 2c 20 6d 5f 2c 20 44 | 49 4d 28 61 29 29 2c 20 |), m_, D|IM(a)), |
|00000080| 6e 5f 2c 20 44 49 4d 28 | 61 29 29 0d 0a 41 44 4a |n_, DIM(|a))..ADJ|
|00000090| 4f 49 4e 5f 45 4c 45 4d | 45 4e 54 28 65 2c 20 76 |OIN_ELEM|ENT(e, v|
|000000a0| 29 3a 3d 41 44 4a 4f 49 | 4e 28 65 2c 20 76 29 0d |):=ADJOI|N(e, v).|
|000000b0| 0a 41 50 50 45 4e 44 5f | 43 4f 4c 55 4d 4e 53 28 |.APPEND_|COLUMNS(|
|000000c0| 61 2c 20 62 29 3a 3d 41 | 50 50 45 4e 44 28 61 60 |a, b):=A|PPEND(a`|
|000000d0| 2c 20 62 60 29 60 0d 0a | 41 50 50 52 4f 58 5f 45 |, b`)`..|APPROX_E|
|000000e0| 49 47 45 4e 56 45 43 54 | 4f 52 28 61 2c 20 b5 29 |IGENVECT|OR(a, .)|
|000000f0| 3a 3d 53 49 47 4e 28 52 | 4f 57 5f 52 45 44 55 43 |:=SIGN(R|OW_REDUC|
|00000100| 45 28 61 20 2d 20 b5 b7 | 49 44 45 4e 54 49 54 59 |E(a - ..|IDENTITY|
|00000110| 5f 4d 41 54 52 49 58 28 | 44 49 4d 28 61 29 29 2c |_MATRIX(|DIM(a)),|
|00000120| 20 56 45 43 54 4f 52 28 | 33 30 30 30 b7 43 4f 53 | VECTOR(|3000.COS|
|00000130| 28 b9 b7 6d 5f 29 2f 6d | 5f 2c 20 6d 5f 2c 20 44 |(..m_)/m|_, m_, D|
|00000140| 49 4d 28 61 29 29 29 60 | 99 28 44 49 4d 28 61 29 |IM(a)))`|.(DIM(a)|
|00000150| 20 2b 20 31 29 29 0d 0a | 42 49 4e 4f 4d 49 41 4c | + 1))..|BINOMIAL|
|00000160| 5f 44 45 4e 53 49 54 59 | 28 6b 2c 20 6e 2c 20 70 |_DENSITY|(k, n, p|
|00000170| 29 3a 3d 70 5e 6b b7 28 | 31 20 2d 20 70 29 5e 28 |):=p^k.(|1 - p)^(|
|00000180| 6e 20 2d 20 6b 29 b7 43 | 4f 4d 42 28 6e 2c 20 6b |n - k).C|OMB(n, k|
|00000190| 29 0d 0a 42 49 4e 4f 4d | 49 41 4c 5f 44 49 53 54 |)..BINOM|IAL_DIST|
|000001a0| 52 49 42 55 54 49 4f 4e | 28 6b 2c 20 6e 2c 20 70 |RIBUTION|(k, n, p|
|000001b0| 29 3a 3d a4 28 42 49 4e | 4f 4d 49 41 4c 5f 44 45 |):=.(BIN|OMIAL_DE|
|000001c0| 4e 53 49 54 59 28 6d 5f | 2c 20 6e 2c 20 70 29 2c |NSITY(m_|, n, p),|
|000001d0| 20 6d 5f 2c 20 30 2c 20 | 4d 49 4e 28 6b 2c 20 6e | m_, 0, |MIN(k, n|
|000001e0| 29 29 0d 0a 43 48 49 5f | 53 51 55 41 52 45 28 78 |))..CHI_|SQUARE(x|
|000001f0| 32 2c 20 76 29 3a 3d 49 | 4e 43 4f 4d 50 4c 45 54 |2, v):=I|NCOMPLET|
|00000200| 45 5f 47 41 4d 4d 41 28 | 76 2f 32 2c 20 78 32 2f |E_GAMMA(|v/2, x2/|
|00000210| 32 29 0d 0a 43 4f 46 41 | 43 54 4f 52 28 61 2c 20 |2)..COFA|CTOR(a, |
|00000220| 69 2c 20 6a 29 3a 3d 43 | 4f 53 28 b9 b7 28 69 20 |i, j):=C|OS(..(i |
|00000230| 2b 20 6a 29 29 b7 44 45 | 54 28 4d 49 4e 4f 52 28 |+ j)).DE|T(MINOR(|
|00000240| 61 2c 20 69 2c 20 6a 29 | 29 0d 0a 43 4f 56 41 52 |a, i, j)|)..COVAR|
|00000250| 49 41 4e 54 5f 4d 45 54 | 52 49 43 5f 54 45 4e 53 |IANT_MET|RIC_TENS|
|00000260| 4f 52 28 61 29 3a 3d 61 | 60 20 95 20 61 0d 0a 44 |OR(a):=a|` . a..D|
|00000270| 49 47 41 4d 4d 41 5f 50 | 53 49 28 7a 29 3a 3d 4c |IGAMMA_P|SI(z):=L|
|00000280| 4e 28 7a 29 20 2d 20 31 | 2f 28 32 b7 7a 29 20 2d |N(z) - 1|/(2.z) -|
|00000290| 20 32 b7 49 4e 54 28 74 | 5f 2f 28 28 74 5f 5e 32 | 2.INT(t|_/((t_^2|
|000002a0| 20 2b 20 7a 5e 32 29 b7 | 28 ea 5e 28 32 b7 b9 b7 | + z^2).|(.^(2...|
|000002b0| 74 5f 29 20 2d 20 31 29 | 29 2c 20 74 5f 2c 20 30 |t_) - 1)|), t_, 0|
|000002c0| 2c 20 96 29 0d 0a 45 55 | 4c 45 52 5f 42 45 54 41 |, .)..EU|LER_BETA|
|000002d0| 28 7a 2c 20 77 29 3a 3d | 9b 28 77 29 b7 9b 28 7a |(z, w):=|.(w)..(z|
|000002e0| 29 2f 9b 28 7a 20 2b 20 | 77 29 0d 0a 45 58 41 43 |)/.(z + |w)..EXAC|
|000002f0| 54 5f 45 49 47 45 4e 56 | 45 43 54 4f 52 28 61 2c |T_EIGENV|ECTOR(a,|
|00000300| 20 b5 29 3a 3d 53 4f 4c | 55 54 49 4f 4e 53 28 56 | .):=SOL|UTIONS(V|
|00000310| 45 43 54 4f 52 28 a4 28 | 41 50 50 45 4e 44 28 22 |ECTOR(.(|APPEND("|
|00000320| 78 22 2c 20 6e 5f 29 b7 | 28 61 99 6d 5f 99 6e 5f |x", n_).|(a.m_.n_|
|00000330| 20 2d 20 b5 b7 4b 52 4f | 4e 45 43 4b 45 52 28 6d | - ..KRO|NECKER(m|
|00000340| 5f 2c 20 6e 5f 29 29 2c | 20 6e 5f 2c 20 31 2c 20 |_, n_)),| n_, 1, |
|00000350| 44 49 4d 28 61 29 29 2c | 20 6d 5f 2c 20 44 49 4d |DIM(a)),| m_, DIM|
|00000360| 28 61 29 29 2c 20 56 45 | 43 54 4f 52 28 41 50 50 |(a)), VE|CTOR(APP|
|00000370| 45 4e 44 28 22 78 22 2c | 20 6e 5f 29 2c 20 6e 5f |END("x",| n_), n_|
|00000380| 2c 20 31 2c 20 44 49 4d | 28 61 29 29 29 0d 0a 46 |, 1, DIM|(a)))..F|
|00000390| 4f 52 43 45 30 28 61 2c | 20 69 2c 20 6a 2c 20 70 |ORCE0(a,| i, j, p|
|000003a0| 29 3a 3d 56 45 43 54 4f | 52 28 49 46 28 6d 5f 20 |):=VECTO|R(IF(m_ |
|000003b0| 3d 20 69 2c 20 61 99 69 | 20 2d 20 61 99 70 b7 61 |= i, a.i| - a.p.a|
|000003c0| 99 69 99 6a 2f 61 99 70 | 99 6a 2c 20 61 99 6d 5f |.i.j/a.p|.j, a.m_|
|000003d0| 29 2c 20 6d 5f 2c 20 44 | 49 4d 28 61 29 29 0d 0a |), m_, D|IM(a))..|
|000003e0| 46 5f 44 49 53 54 52 49 | 42 55 54 49 4f 4e 28 66 |F_DISTRI|BUTION(f|
|000003f0| 5f 2c 20 76 31 2c 20 76 | 32 29 3a 3d 49 4e 43 4f |_, v1, v|2):=INCO|
|00000400| 4d 50 4c 45 54 45 5f 42 | 45 54 41 28 76 32 2f 28 |MPLETE_B|ETA(v2/(|
|00000410| 76 32 20 2b 20 66 5f b7 | 76 31 29 2c 20 76 32 2f |v2 + f_.|v1), v2/|
|00000420| 32 2c 20 76 31 2f 32 29 | 0d 0a 47 41 55 53 53 5f |2, v1/2)|..GAUSS_|
|00000430| 4e 45 57 54 4f 4e 28 65 | 71 2c 20 70 61 72 6d 2c |NEWTON(e|q, parm,|
|00000440| 20 df 30 5f 2c 20 64 61 | 74 61 2c 20 6e 5f 20 3a | .0_, da|ta, n_ :|
|00000450| 3d 20 33 30 2c 20 73 6e | 5f 20 3a 3d 20 31 30 2c |= 30, sn|_ := 10,|
|00000460| 20 63 63 5f 20 3a 3d 20 | 31 30 5e 28 2d 38 29 29 | cc_ := |10^(-8))|
|00000470| 3a 3d 50 52 4f 47 28 65 | 71 71 20 3a 3d 20 52 48 |:=PROG(e|qq := RH|
|00000480| 53 28 65 71 29 2c 20 64 | 76 61 72 20 3a 3d 20 4c |S(eq), d|var := L|
|00000490| 48 53 28 65 71 29 2c 20 | 69 76 61 72 73 32 20 3a |HS(eq), |ivars2 :|
|000004a0| 3d 20 56 41 52 49 41 42 | 4c 45 53 28 65 71 71 29 |= VARIAB|LES(eqq)|
|000004b0| 2c 20 69 76 61 72 73 20 | 3a 3d 20 53 45 4c 45 43 |, ivars |:= SELEC|
|000004c0| 54 28 ac 20 4d 45 4d 42 | 45 52 3f 28 69 5f 2c 20 |T(. MEMB|ER?(i_, |
|000004d0| 70 61 72 6d 29 2c 20 69 | 5f 2c 20 69 76 61 72 73 |parm), i|_, ivars|
|000004e0| 32 29 2c 20 69 76 61 72 | 73 31 20 3a 3d 20 69 76 |2), ivar|s1 := iv|
|000004f0| 61 72 73 2c 20 70 61 72 | 6d 31 20 3a 3d 20 70 61 |ars, par|m1 := pa|
|00000500| 72 6d 2c 20 6f 62 73 5f | 20 3a 3d 20 44 49 4d 28 |rm, obs_| := DIM(|
|00000510| 64 61 74 61 29 20 2d 20 | 31 2c 20 76 61 72 73 20 |data) - |1, vars |
|00000520| 3a 3d 20 64 61 74 61 99 | 31 2c 20 64 61 74 61 31 |:= data.|1, data1|
|00000530| 20 3a 3d 20 64 61 74 61 | 99 5b 32 2c 20 2e 2e 2e | := data|.[2, ...|
|00000540| 2c 20 6f 62 73 5f 20 2b | 20 31 5d 2c 20 6b 5f 20 |, obs_ +| 1], k_ |
|00000550| 3a 3d 20 44 49 4d 28 70 | 61 72 6d 29 2c 20 70 6f |:= DIM(p|arm), po|
|00000560| 73 5f 20 3a 3d 20 70 6f | 73 69 74 69 6f 6e 28 64 |s_ := po|sition(d|
|00000570| 76 61 72 2c 20 76 61 72 | 73 29 2c 20 49 46 28 70 |var, var|s), IF(p|
|00000580| 6f 73 5f 20 3d 20 30 2c | 20 52 45 54 55 52 4e 20 |os_ = 0,| RETURN |
|00000590| 22 55 6e 64 65 66 69 6e | 65 64 20 64 65 70 65 6e |"Undefin|ed depen|
|000005a0| 64 65 6e 74 20 76 61 72 | 69 61 62 6c 65 21 22 2c |dent var|iable!",|
|000005b0| 20 79 5f 20 3a 3d 20 64 | 61 74 61 31 99 99 70 6f | y_ := d|ata1..po|
|000005c0| 73 5f 29 2c 20 70 6f 73 | 5f 20 3a 3d 20 56 45 43 |s_), pos|_ := VEC|
|000005d0| 54 4f 52 28 70 6f 73 69 | 74 69 6f 6e 28 69 76 61 |TOR(posi|tion(iva|
|000005e0| 72 73 99 69 5f 2c 20 76 | 61 72 73 29 2c 20 69 5f |rs.i_, v|ars), i_|
|000005f0| 2c 20 31 2c 20 44 49 4d | 28 69 76 61 72 73 29 29 |, 1, DIM|(ivars))|
|00000600| 2c 20 49 46 28 4d 45 4d | 42 45 52 3f 28 30 2c 20 |, IF(MEM|BER?(0, |
|00000610| 70 6f 73 5f 29 2c 20 52 | 45 54 55 52 4e 20 22 55 |pos_), R|ETURN "U|
|00000620| 6e 64 65 66 69 6e 65 64 | 20 69 6e 64 65 70 65 6e |ndefined| indepen|
|00000630| 64 65 6e 74 20 76 61 72 | 69 61 62 6c 65 28 73 29 |dent var|iable(s)|
|00000640| 21 22 2c 20 6d 5f 20 3a | 3d 20 64 61 74 61 31 99 |!", m_ :|= data1.|
|00000650| 99 70 6f 73 5f 29 2c 20 | 79 63 70 74 20 3a 3d 20 |.pos_), |ycpt := |
|00000660| 49 46 28 4d 45 4d 42 45 | 52 3f 28 31 2c 20 47 52 |IF(MEMBE|R?(1, GR|
|00000670| 41 44 28 65 71 71 2c 20 | 70 61 72 6d 29 29 2c 20 |AD(eqq, |parm)), |
|00000680| 31 2c 20 30 2c 20 30 29 | 2c 20 6b 5f 20 3a 2d 20 |1, 0, 0)|, k_ :- |
|00000690| 79 63 70 74 2c 20 df 5f | 20 3a 3d 20 df 30 5f 2c |ycpt, ._| := .0_,|
|000006a0| 20 65 71 31 20 3a 3d 20 | 56 45 43 54 4f 52 28 53 | eq1 := |VECTOR(S|
|000006b0| 55 42 53 54 28 65 71 71 | 2c 20 69 76 61 72 73 2c |UBST(eqq|, ivars,|
|000006c0| 20 6d 5f 99 69 29 2c 20 | 69 2c 20 6f 62 73 5f 29 | m_.i), |i, obs_)|
|000006d0| 2c 20 6a 5f 20 3a 3d 20 | 47 52 41 44 28 65 71 31 |, j_ := |GRAD(eq1|
|000006e0| 2c 20 70 61 72 6d 29 60 | 2c 20 70 5f 20 3a 3d 20 |, parm)`|, p_ := |
|000006f0| 53 55 42 53 54 28 65 71 | 31 2c 20 70 61 72 6d 2c |SUBST(eq|1, parm,|
|00000700| 20 df 5f 29 2c 20 72 5f | 20 3a 3d 20 79 5f 20 2d | ._), r_| := y_ -|
|00000710| 20 70 5f 2c 20 73 73 65 | 5f 20 3a 3d 20 72 5f 20 | p_, sse|_ := r_ |
|00000720| 95 20 72 5f 2c 20 74 5f | 20 3a 3d 20 49 46 28 79 |. r_, t_| := IF(y|
|00000730| 63 70 74 20 3d 20 31 2c | 20 79 5f 20 2d 20 56 45 |cpt = 1,| y_ - VE|
|00000740| 43 54 4f 52 28 41 56 45 | 52 41 47 45 28 79 5f 29 |CTOR(AVE|RAGE(y_)|
|00000750| 2c 20 69 2c 20 6f 62 73 | 5f 29 2c 20 79 5f 29 2c |, i, obs|_), y_),|
|00000760| 20 73 73 74 5f 20 3a 3d | 20 74 5f 20 95 20 74 5f | sst_ :=| t_ . t_|
|00000770| 2c 20 69 74 65 72 20 3a | 3d 20 5b 41 50 50 45 4e |, iter :|= [APPEN|
|00000780| 44 28 5b 30 5d 2c 20 df | 5f 2c 20 5b 73 73 65 5f |D([0], .|_, [sse_|
|00000790| 5d 29 5d 2c 20 65 70 73 | 5f 20 3a 3d 20 31 2c 20 |])], eps|_ := 1, |
|000007a0| 69 74 65 72 5f 20 3a 3d | 20 31 2c 20 4c 4f 4f 50 |iter_ :=| 1, LOOP|
|000007b0| 28 49 46 28 69 74 65 72 | 5f 20 3e 20 6e 5f 20 90 |(IF(iter|_ > n_ .|
|000007c0| 20 65 70 73 5f 20 3c 20 | 63 63 5f 2c 20 65 78 69 | eps_ < |cc_, exi|
|000007d0| 74 29 2c 20 78 5f 20 3a | 3d 20 53 55 42 53 54 28 |t), x_ :|= SUBST(|
|000007e0| 6a 5f 2c 20 70 61 72 6d | 2c 20 df 5f 29 2c 20 6c |j_, parm|, ._), l|
|000007f0| 73 73 65 5f 20 3a 3d 20 | 73 73 65 5f 2c 20 78 70 |sse_ := |sse_, xp|
|00000800| 78 69 20 3a 3d 20 31 2f | 28 78 5f 60 20 95 20 78 |xi := 1/|(x_` . x|
|00000810| 5f 29 2c 20 ab 5f 20 3a | 3d 20 78 70 78 69 20 95 |_), ._ :|= xpxi .|
|00000820| 20 78 5f 60 20 95 20 72 | 5f 2c 20 6f 6c 64 5f 20 | x_` . r|_, old_ |
|00000830| 3a 3d 20 df 5f 2c 20 df | 5f 20 3a 2b 20 ab 5f 2c |:= ._, .|_ :+ ._,|
|00000840| 20 70 5f 20 3a 3d 20 53 | 55 42 53 54 28 65 71 31 | p_ := S|UBST(eq1|
|00000850| 2c 20 70 61 72 6d 2c 20 | df 5f 29 2c 20 72 5f 20 |, parm, |._), r_ |
|00000860| 3a 3d 20 79 5f 20 2d 20 | 70 5f 2c 20 73 73 65 5f |:= y_ - |p_, sse_|
|00000870| 20 3a 3d 20 72 5f 20 95 | 20 72 5f 2c 20 69 74 65 | := r_ .| r_, ite|
|00000880| 72 20 3a 3d 20 41 50 50 | 45 4e 44 28 69 74 65 72 |r := APP|END(iter|
|00000890| 2c 20 5b 41 50 50 45 4e | 44 28 5b 69 74 65 72 5f |, [APPEN|D([iter_|
|000008a0| 5d 2c 20 df 5f 2c 20 5b | 73 73 65 5f 5d 29 5d 29 |], ._, [|sse_])])|
|000008b0| 2c 20 73 75 62 69 74 5f | 20 3a 3d 20 31 2c 20 4c |, subit_| := 1, L|
|000008c0| 4f 4f 50 28 49 46 28 73 | 73 65 5f 20 93 20 6c 73 |OOP(IF(s|se_ . ls|
|000008d0| 73 65 5f 20 90 20 73 75 | 62 69 74 5f 20 3e 20 73 |se_ . su|bit_ > s|
|000008e0| 6e 5f 2c 20 65 78 69 74 | 29 2c 20 ab 5f 20 3a 2f |n_, exit|), ._ :/|
|000008f0| 20 32 2c 20 df 5f 20 3a | 3d 20 6f 6c 64 5f 20 2b | 2, ._ :|= old_ +|
|00000900| 20 ab 5f 2c 20 70 5f 20 | 3a 3d 20 53 55 42 53 54 | ._, p_ |:= SUBST|
|00000910| 28 65 71 31 2c 20 70 61 | 72 6d 2c 20 df 5f 29 2c |(eq1, pa|rm, ._),|
|00000920| 20 72 5f 20 3a 3d 20 79 | 5f 20 2d 20 70 5f 2c 20 | r_ := y|_ - p_, |
|00000930| 73 73 65 5f 20 3a 3d 20 | 72 5f 20 95 20 72 5f 2c |sse_ := |r_ . r_,|
|00000940| 20 69 74 65 72 20 3a 3d | 20 41 50 50 45 4e 44 28 | iter :=| APPEND(|
|00000950| 69 74 65 72 2c 20 5b 41 | 50 50 45 4e 44 28 5b 69 |iter, [A|PPEND([i|
|00000960| 74 65 72 5f 20 2b 20 73 | 75 62 69 74 5f 2f 31 30 |ter_ + s|ubit_/10|
|00000970| 30 5d 2c 20 df 5f 2c 20 | 5b 73 73 65 5f 5d 29 5d |0], ._, |[sse_])]|
|00000980| 29 2c 20 73 75 62 69 74 | 5f 20 3a 2b 20 31 29 2c |), subit|_ :+ 1),|
|00000990| 20 49 46 28 73 75 62 69 | 74 5f 20 3e 20 73 6e 5f | IF(subi|t_ > sn_|
|000009a0| 2c 20 65 78 69 74 29 2c | 20 65 70 73 5f 20 3a 3d |, exit),| eps_ :=|
|000009b0| 20 41 42 53 28 28 6c 73 | 73 65 5f 20 2d 20 73 73 | ABS((ls|se_ - ss|
|000009c0| 65 5f 29 2f 28 73 73 65 | 5f 20 2b 20 31 30 5e 28 |e_)/(sse|_ + 10^(|
|000009d0| 2d 36 29 29 29 2c 20 69 | 74 65 72 5f 20 3a 2b 20 |-6))), i|ter_ :+ |
|000009e0| 31 29 2c 20 49 46 28 69 | 74 65 72 5f 20 3e 20 6e |1), IF(i|ter_ > n|
|000009f0| 5f 2c 20 6d 73 67 20 3a | 3d 20 41 50 50 45 4e 44 |_, msg :|= APPEND|
|00000a00| 28 22 43 6f 6e 76 65 72 | 67 65 6e 63 65 20 66 61 |("Conver|gence fa|
|00000a10| 69 6c 65 64 20 61 66 74 | 65 72 20 22 2c 20 6e 5f |iled aft|er ", n_|
|00000a20| 2c 20 22 20 69 74 65 72 | 61 74 69 6f 6e 73 21 22 |, " iter|ations!"|
|00000a30| 29 2c 20 49 46 28 73 75 | 62 69 74 5f 20 3e 20 73 |), IF(su|bit_ > s|
|00000a40| 6e 5f 2c 20 6d 73 67 20 | 3a 3d 20 22 53 53 45 20 |n_, msg |:= "SSE |
|00000a50| 64 69 64 20 6e 6f 74 20 | 69 6d 70 72 6f 76 65 20 |did not |improve |
|00000a60| 61 66 74 65 72 20 31 30 | 20 68 61 6c 76 69 6e 67 |after 10| halving|
|00000a70| 73 21 22 2c 20 6d 73 67 | 20 3a 3d 20 22 43 6f 6e |s!", msg| := "Con|
|00000a80| 76 65 72 67 65 6e 63 65 | 20 63 72 69 74 65 72 69 |vergence| criteri|
|00000a90| 61 20 6d 65 74 21 22 29 | 29 2c 20 64 66 65 5f 20 |a met!")|), dfe_ |
|00000aa0| 3a 3d 20 6f 62 73 5f 20 | 2d 20 6b 5f 20 2d 20 79 |:= obs_ |- k_ - y|
|00000ab0| 63 70 74 2c 20 6d 73 65 | 5f 20 3a 3d 20 73 73 65 |cpt, mse|_ := sse|
|00000ac0| 5f 2f 64 66 65 5f 2c 20 | 73 65 5f 20 3a 3d 20 8b |_/dfe_, |se_ := .|
|00000ad0| 6d 73 65 5f 2c 20 73 74 | 64 5f 20 3a 3d 20 56 45 |mse_, st|d_ := VE|
|00000ae0| 43 54 4f 52 28 8b 28 78 | 70 78 69 99 69 5f 99 69 |CTOR(.(x|pxi.i_.i|
|00000af0| 5f b7 6d 73 65 5f 29 2c | 20 69 5f 2c 20 6b 5f 20 |_.mse_),| i_, k_ |
|00000b00| 2b 20 79 63 70 74 29 2c | 20 74 5f 20 3a 3d 20 56 |+ ycpt),| t_ := V|
|00000b10| 45 43 54 4f 52 28 df 5f | 99 69 5f 2f 73 74 64 5f |ECTOR(._|.i_/std_|
|00000b20| 99 69 5f 2c 20 69 5f 2c | 20 6b 5f 20 2b 20 79 63 |.i_, i_,| k_ + yc|
|00000b30| 70 74 29 2c 20 70 72 6f | 62 5f 20 3a 3d 20 56 45 |pt), pro|b_ := VE|
|00000b40| 43 54 4f 52 28 31 20 2d | 20 53 54 55 44 45 4e 54 |CTOR(1 -| STUDENT|
|00000b50| 28 74 5f 99 69 5f 2c 20 | 64 66 65 5f 29 2c 20 69 |(t_.i_, |dfe_), i|
|00000b60| 5f 2c 20 6b 5f 20 2b 20 | 79 63 70 74 29 2c 20 73 |_, k_ + |ycpt), s|
|00000b70| 73 72 5f 20 3a 3d 20 73 | 73 74 5f 20 2d 20 73 73 |sr_ := s|st_ - ss|
|00000b80| 65 5f 2c 20 72 73 71 20 | 3a 3d 20 73 73 72 5f 2f |e_, rsq |:= ssr_/|
|00000b90| 73 73 74 5f 2c 20 61 72 | 73 71 20 3a 3d 20 31 20 |sst_, ar|sq := 1 |
|00000ba0| 2d 20 73 73 65 5f b7 28 | 6f 62 73 5f 20 2d 20 79 |- sse_.(|obs_ - y|
|00000bb0| 63 70 74 29 2f 28 28 6f | 62 73 5f 20 2d 20 6b 5f |cpt)/((o|bs_ - k_|
|00000bc0| 20 2d 20 79 63 70 74 29 | b7 73 73 74 5f 29 2c 20 | - ycpt)|.sst_), |
|00000bd0| 66 64 5f 20 3a 3d 20 46 | 5f 44 49 53 54 52 49 42 |fd_ := F|_DISTRIB|
|00000be0| 55 54 49 4f 4e 28 73 73 | 72 5f b7 64 66 65 5f 2f |UTION(ss|r_.dfe_/|
|00000bf0| 28 73 73 65 5f b7 6b 5f | 29 2c 20 6b 5f 2c 20 64 |(sse_.k_|), k_, d|
|00000c00| 66 65 5f 29 2c 20 61 31 | 5f 20 3a 3d 20 5b 22 53 |fe_), a1|_ := ["S|
|00000c10| 6f 75 72 63 65 22 2c 20 | 22 44 46 22 2c 20 22 53 |ource", |"DF", "S|
|00000c20| 53 22 2c 20 22 4d 53 22 | 2c 20 22 46 22 2c 20 22 |S", "MS"|, "F", "|
|00000c30| 50 72 6f 62 28 46 29 22 | 5d 2c 20 61 32 5f 20 3a |Prob(F)"|], a2_ :|
|00000c40| 3d 20 5b 22 52 65 67 22 | 2c 20 6b 5f 2c 20 73 73 |= ["Reg"|, k_, ss|
|00000c50| 72 5f 2c 20 73 73 72 5f | 2f 6b 5f 2c 20 73 73 72 |r_, ssr_|/k_, ssr|
|00000c60| 5f b7 64 66 65 5f 2f 28 | 73 73 65 5f b7 6b 5f 29 |_.dfe_/(|sse_.k_)|
|00000c70| 2c 20 49 46 28 66 64 5f | 20 3c 20 30 2c 20 30 2c |, IF(fd_| < 0, 0,|
|00000c80| 20 66 64 5f 29 5d 2c 20 | 61 33 5f 20 3a 3d 20 5b | fd_)], |a3_ := [|
|00000c90| 22 45 72 72 6f 72 22 2c | 20 64 66 65 5f 2c 20 73 |"Error",| dfe_, s|
|00000ca0| 73 65 5f 2c 20 73 73 65 | 5f 2f 64 66 65 5f 2c 20 |se_, sse|_/dfe_, |
|00000cb0| 22 20 22 2c 20 22 20 22 | 5d 2c 20 61 34 5f 20 3a |" ", " "|], a4_ :|
|00000cc0| 3d 20 5b 49 46 28 79 63 | 70 74 20 3d 20 31 2c 20 |= [IF(yc|pt = 1, |
|00000cd0| 22 43 6f 72 72 65 63 74 | 65 64 20 54 6f 74 61 6c |"Correct|ed Total|
|00000ce0| 22 2c 20 22 55 6e 63 6f | 72 72 65 63 74 65 64 20 |", "Unco|rrected |
|00000cf0| 54 6f 74 61 6c 22 29 2c | 20 64 66 65 5f 20 2b 20 |Total"),| dfe_ + |
|00000d00| 6b 5f 2c 20 73 73 74 5f | 2c 20 22 20 22 2c 20 22 |k_, sst_|, " ", "|
|00000d10| 20 22 2c 20 22 20 22 5d | 2c 20 61 6e 6f 76 61 20 | ", " "]|, anova |
|00000d20| 3a 3d 20 41 50 50 45 4e | 44 28 5b 61 31 5f 5d 2c |:= APPEN|D([a1_],|
|00000d30| 20 5b 61 32 5f 5d 2c 20 | 5b 61 33 5f 5d 2c 20 5b | [a2_], |[a3_], [|
|00000d40| 61 34 5f 5d 29 2c 20 66 | 65 71 20 3a 3d 20 64 76 |a4_]), f|eq := dv|
|00000d50| 61 72 20 3d 20 53 55 42 | 53 54 28 65 71 71 2c 20 |ar = SUB|ST(eqq, |
|00000d60| 70 61 72 6d 2c 20 df 5f | 29 2c 20 69 74 65 72 20 |parm, ._|), iter |
|00000d70| 3a 3d 20 41 50 50 45 4e | 44 28 5b 41 50 50 45 4e |:= APPEN|D([APPEN|
|00000d80| 44 28 5b 22 49 74 65 72 | 22 5d 2c 20 70 61 72 6d |D(["Iter|"], parm|
|00000d90| 2c 20 5b 22 53 53 45 22 | 5d 29 5d 2c 20 69 74 65 |, ["SSE"|])], ite|
|00000da0| 72 29 2c 20 74 69 74 6c | 65 20 3a 3d 20 5b 22 50 |r), titl|e := ["P|
|00000db0| 61 72 6d 22 3b 20 22 56 | 61 6c 75 65 22 3b 20 22 |arm"; "V|alue"; "|
|00000dc0| 53 54 44 22 3b 20 41 50 | 50 45 4e 44 28 22 74 28 |STD"; AP|PEND("t(|
|00000dd0| 22 2c 20 64 66 65 5f 2c | 20 22 29 22 29 3b 20 22 |", dfe_,| ")"); "|
|00000de0| 50 72 6f 62 28 74 29 22 | 5d 2c 20 6f 75 74 70 75 |Prob(t)"|], outpu|
|00000df0| 74 20 3a 3d 20 41 50 50 | 45 4e 44 5f 43 4f 4c 55 |t := APP|END_COLU|
|00000e00| 4d 4e 53 28 74 69 74 6c | 65 2c 20 41 50 50 45 4e |MNS(titl|e, APPEN|
|00000e10| 44 28 5b 70 61 72 6d 5d | 2c 20 5b df 5f 5d 2c 20 |D([parm]|, [._], |
|00000e20| 5b 73 74 64 5f 5d 2c 20 | 5b 74 5f 5d 2c 20 5b 70 |[std_], |[t_], [p|
|00000e30| 72 6f 62 5f 5d 29 29 60 | 2c 20 73 74 61 74 73 20 |rob_]))`|, stats |
|00000e40| 3a 3d 20 5b 22 53 45 22 | 2c 20 22 52 5e 32 22 2c |:= ["SE"|, "R^2",|
|00000e50| 20 22 41 64 6a 52 5e 32 | 22 3b 20 73 65 5f 2c 20 | "AdjR^2|"; se_, |
|00000e60| 72 73 71 2c 20 61 72 73 | 71 5d 2c 20 52 45 54 55 |rsq, ars|q], RETU|
|00000e70| 52 4e 20 5b 22 47 61 75 | 73 73 5f 4e 65 77 74 6f |RN ["Gau|ss_Newto|
|00000e80| 6e 20 4d 65 74 68 6f 64 | 22 3b 20 22 20 22 3b 20 |n Method|"; " "; |
|00000e90| 6d 73 67 3b 20 22 20 22 | 3b 20 6f 75 74 70 75 74 |msg; " "|; output|
|00000ea0| 3b 20 22 20 22 3b 20 61 | 6e 6f 76 61 3b 20 22 20 |; " "; a|nova; " |
|00000eb0| 22 3b 20 73 74 61 74 73 | 3b 20 22 20 22 3b 20 66 |"; stats|; " "; f|
|00000ec0| 65 71 5d 29 0d 0a 47 45 | 4f 4d 45 54 52 59 5f 4d |eq])..GE|OMETRY_M|
|00000ed0| 41 54 52 49 58 28 a9 2c | 20 67 29 3a 3d 5b a9 2c |ATRIX(.,| g):=[.,|
|00000ee0| 20 56 45 43 54 4f 52 28 | 8b 28 67 99 6d 5f 99 6d | VECTOR(|.(g.m_.m|
|00000ef0| 5f 29 2c 20 6d 5f 2c 20 | 44 49 4d 28 67 29 29 5d |_), m_, |DIM(g))]|
|00000f00| 0d 0a 48 59 50 45 52 47 | 45 4f 4d 45 54 52 49 43 |..HYPERG|EOMETRIC|
|00000f10| 5f 44 45 4e 53 49 54 59 | 28 6b 2c 20 6e 2c 20 6d |_DENSITY|(k, n, m|
|00000f20| 2c 20 6a 29 3a 3d 43 4f | 4d 42 28 6d 2c 20 6b 29 |, j):=CO|MB(m, k)|
|00000f30| b7 43 4f 4d 42 28 6a 20 | 2d 20 6d 2c 20 6e 20 2d |.COMB(j |- m, n -|
|00000f40| 20 6b 29 2f 43 4f 4d 42 | 28 6a 2c 20 6e 29 0d 0a | k)/COMB|(j, n)..|
|00000f50| 48 59 50 45 52 47 45 4f | 4d 45 54 52 49 43 5f 44 |HYPERGEO|METRIC_D|
|00000f60| 49 53 54 52 49 42 55 54 | 49 4f 4e 28 6b 2c 20 6e |ISTRIBUT|ION(k, n|
|00000f70| 2c 20 6d 2c 20 6a 29 3a | 3d a4 28 48 59 50 45 52 |, m, j):|=.(HYPER|
|00000f80| 47 45 4f 4d 45 54 52 49 | 43 5f 44 45 4e 53 49 54 |GEOMETRI|C_DENSIT|
|00000f90| 59 28 6c 5f 2c 20 6e 2c | 20 6d 2c 20 6a 29 2c 20 |Y(l_, n,| m, j), |
|00000fa0| 6c 5f 2c 20 4d 41 58 28 | 30 2c 20 6e 20 2d 20 6a |l_, MAX(|0, n - j|
|00000fb0| 20 2b 20 6d 29 2c 20 4d | 49 4e 28 6b 2c 20 6e 2c | + m), M|IN(k, n,|
|00000fc0| 20 6d 29 29 0d 0a 49 4e | 43 4f 4d 50 4c 45 54 45 | m))..IN|COMPLETE|
|00000fd0| 5f 42 45 54 41 28 78 2c | 20 7a 2c 20 77 29 3a 3d |_BETA(x,| z, w):=|
|00000fe0| 28 78 5e 7a 2f 7a 20 2b | 20 28 31 20 2d 20 28 31 |(x^z/z +| (1 - (1|
|00000ff0| 20 2d 20 78 29 5e 77 29 | 2f 77 20 2b 20 49 4e 54 | - x)^w)|/w + INT|
|00001000| 28 74 5f 5e 28 7a 20 2d | 20 31 29 b7 28 28 31 20 |(t_^(z -| 1).((1 |
|00001010| 2d 20 74 5f 29 5e 28 77 | 20 2d 20 31 29 20 2d 20 |- t_)^(w| - 1) - |
|00001020| 31 29 20 2d 20 28 31 20 | 2d 20 74 5f 29 5e 28 77 |1) - (1 |- t_)^(w|
|00001030| 20 2d 20 31 29 2c 20 74 | 5f 2c 20 30 2c 20 78 29 | - 1), t|_, 0, x)|
|00001040| 29 2f 45 55 4c 45 52 5f | 42 45 54 41 28 7a 2c 20 |)/EULER_|BETA(z, |
|00001050| 77 29 0d 0a 49 4e 43 4f | 4d 50 4c 45 54 45 5f 47 |w)..INCO|MPLETE_G|
|00001060| 41 4d 4d 41 28 7a 2c 20 | 77 29 3a 3d 28 77 5e 7a |AMMA(z, |w):=(w^z|
|00001070| 2f 7a 20 2b 20 49 4e 54 | 28 74 5f 5e 28 7a 20 2d |/z + INT|(t_^(z -|
|00001080| 20 31 29 b7 28 ea 5e 28 | 2d 74 5f 29 20 2d 20 31 | 1).(.^(|-t_) - 1|
|00001090| 29 2c 20 74 5f 2c 20 30 | 2c 20 77 29 29 2f 9b 28 |), t_, 0|, w))/.(|
|000010a0| 7a 29 0d 0a 49 4e 43 4f | 4d 50 4c 45 54 45 5f 47 |z)..INCO|MPLETE_G|
|000010b0| 41 4d 4d 41 5f 53 45 52 | 49 45 53 28 7a 2c 20 77 |AMMA_SER|IES(z, w|
|000010c0| 2c 20 6d 29 3a 3d ea 5e | 28 2d 77 29 b7 77 5e 7a |, m):=.^|(-w).w^z|
|000010d0| b7 a4 28 77 5e 28 6d 20 | 2d 20 6e 5f 29 2f 9b 28 |..(w^(m |- n_)/.(|
|000010e0| 7a 20 2b 20 6d 20 2d 20 | 6e 5f 20 2b 20 31 29 2c |z + m - |n_ + 1),|
|000010f0| 20 6e 5f 2c 20 30 2c 20 | 6d 29 0d 0a 4a 41 43 4f | n_, 0, |m)..JACO|
|00001100| 42 49 41 4e 28 75 2c 20 | a9 29 3a 3d 56 45 43 54 |BIAN(u, |.):=VECT|
|00001110| 4f 52 28 47 52 41 44 28 | 75 99 6d 5f 2c 20 a9 29 |OR(GRAD(|u.m_, .)|
|00001120| 2c 20 6d 5f 2c 20 44 49 | 4d 28 75 29 29 0d 0a 4b |, m_, DI|M(u))..K|
|00001130| 52 4f 4e 45 43 4b 45 52 | 28 69 2c 20 6a 29 3a 3d |RONECKER|(i, j):=|
|00001140| 49 46 28 69 20 3d 20 6a | 2c 20 31 2c 20 30 29 0d |IF(i = j|, 1, 0).|
|00001150| 0a 4d 41 52 51 55 41 52 | 44 54 28 65 71 2c 20 70 |.MARQUAR|DT(eq, p|
|00001160| 61 72 6d 2c 20 df 30 5f | 2c 20 64 61 74 61 2c 20 |arm, .0_|, data, |
|00001170| 6e 5f 20 3a 3d 20 35 30 | 2c 20 63 63 5f 20 3a 3d |n_ := 50|, cc_ :=|
|00001180| 20 31 30 5e 28 2d 38 29 | 29 3a 3d 50 52 4f 47 28 | 10^(-8)|):=PROG(|
|00001190| 65 71 71 20 3a 3d 20 52 | 48 53 28 65 71 29 2c 20 |eqq := R|HS(eq), |
|000011a0| 64 76 61 72 20 3a 3d 20 | 4c 48 53 28 65 71 29 2c |dvar := |LHS(eq),|
|000011b0| 20 69 76 61 72 73 32 20 | 3a 3d 20 56 41 52 49 41 | ivars2 |:= VARIA|
|000011c0| 42 4c 45 53 28 65 71 71 | 29 2c 20 69 76 61 72 73 |BLES(eqq|), ivars|
|000011d0| 20 3a 3d 20 53 45 4c 45 | 43 54 28 ac 20 4d 45 4d | := SELE|CT(. MEM|
|000011e0| 42 45 52 3f 28 69 5f 2c | 20 70 61 72 6d 29 2c 20 |BER?(i_,| parm), |
|000011f0| 69 5f 2c 20 69 76 61 72 | 73 32 29 2c 20 70 61 72 |i_, ivar|s2), par|
|00001200| 6d 31 20 3a 3d 20 70 61 | 72 6d 2c 20 69 76 61 72 |m1 := pa|rm, ivar|
|00001210| 73 31 20 3a 3d 20 69 76 | 61 72 73 2c 20 6f 62 73 |s1 := iv|ars, obs|
|00001220| 5f 20 3a 3d 20 44 49 4d | 28 64 61 74 61 29 20 2d |_ := DIM|(data) -|
|00001230| 20 31 2c 20 76 61 72 73 | 20 3a 3d 20 64 61 74 61 | 1, vars| := data|
|00001240| 99 31 2c 20 64 61 74 61 | 31 20 3a 3d 20 64 61 74 |.1, data|1 := dat|
|00001250| 61 99 5b 32 2c 20 2e 2e | 2e 2c 20 6f 62 73 5f 20 |a.[2, ..|., obs_ |
|00001260| 2b 20 31 5d 2c 20 6b 5f | 20 3a 3d 20 44 49 4d 28 |+ 1], k_| := DIM(|
|00001270| 70 61 72 6d 29 2c 20 70 | 6f 73 5f 20 3a 3d 20 70 |parm), p|os_ := p|
|00001280| 6f 73 69 74 69 6f 6e 28 | 64 76 61 72 2c 20 76 61 |osition(|dvar, va|
|00001290| 72 73 29 2c 20 49 46 28 | 70 6f 73 5f 20 3d 20 30 |rs), IF(|pos_ = 0|
|000012a0| 2c 20 52 45 54 55 52 4e | 20 22 55 6e 64 65 66 69 |, RETURN| "Undefi|
|000012b0| 6e 65 64 20 64 65 70 65 | 6e 64 65 6e 74 20 76 61 |ned depe|ndent va|
|000012c0| 72 69 61 62 6c 65 21 22 | 2c 20 79 5f 20 3a 3d 20 |riable!"|, y_ := |
|000012d0| 64 61 74 61 31 99 99 70 | 6f 73 5f 29 2c 20 70 6f |data1..p|os_), po|
|000012e0| 73 5f 20 3a 3d 20 56 45 | 43 54 4f 52 28 70 6f 73 |s_ := VE|CTOR(pos|
|000012f0| 69 74 69 6f 6e 28 69 76 | 61 72 73 99 69 5f 2c 20 |ition(iv|ars.i_, |
|00001300| 76 61 72 73 29 2c 20 69 | 5f 2c 20 31 2c 20 44 49 |vars), i|_, 1, DI|
|00001310| 4d 28 69 76 61 72 73 29 | 29 2c 20 49 46 28 4d 45 |M(ivars)|), IF(ME|
|00001320| 4d 42 45 52 3f 28 30 2c | 20 70 6f 73 5f 29 2c 20 |MBER?(0,| pos_), |
|00001330| 52 45 54 55 52 4e 20 22 | 55 6e 64 65 66 69 6e 65 |RETURN "|Undefine|
|00001340| 64 20 69 6e 64 65 70 65 | 6e 64 65 6e 74 20 76 61 |d indepe|ndent va|
|00001350| 72 69 61 62 6c 65 28 73 | 29 21 22 2c 20 6d 5f 20 |riable(s|)!", m_ |
|00001360| 3a 3d 20 64 61 74 61 31 | 99 99 70 6f 73 5f 29 2c |:= data1|..pos_),|
|00001370| 20 79 63 70 74 20 3a 3d | 20 49 46 28 4d 45 4d 42 | ycpt :=| IF(MEMB|
|00001380| 45 52 3f 28 31 2c 20 47 | 52 41 44 28 65 71 71 2c |ER?(1, G|RAD(eqq,|
|00001390| 20 70 61 72 6d 29 29 2c | 20 31 2c 20 30 2c 20 30 | parm)),| 1, 0, 0|
|000013a0| 29 2c 20 6b 5f 20 3a 2d | 20 79 63 70 74 2c 20 df |), k_ :-| ycpt, .|
|000013b0| 5f 20 3a 3d 20 df 30 5f | 2c 20 65 71 31 20 3a 3d |_ := .0_|, eq1 :=|
|000013c0| 20 56 45 43 54 4f 52 28 | 53 55 42 53 54 28 65 71 | VECTOR(|SUBST(eq|
|000013d0| 71 2c 20 69 76 61 72 73 | 2c 20 6d 5f 99 69 29 2c |q, ivars|, m_.i),|
|000013e0| 20 69 2c 20 6f 62 73 5f | 29 2c 20 6a 5f 20 3a 3d | i, obs_|), j_ :=|
|000013f0| 20 47 52 41 44 28 65 71 | 31 2c 20 70 61 72 6d 29 | GRAD(eq|1, parm)|
|00001400| 60 2c 20 70 5f 20 3a 3d | 20 53 55 42 53 54 28 65 |`, p_ :=| SUBST(e|
|00001410| 71 31 2c 20 70 61 72 6d | 2c 20 df 5f 29 2c 20 72 |q1, parm|, ._), r|
|00001420| 5f 20 3a 3d 20 79 5f 20 | 2d 20 70 5f 2c 20 73 73 |_ := y_ |- p_, ss|
|00001430| 65 5f 20 3a 3d 20 72 5f | 20 95 20 72 5f 2c 20 74 |e_ := r_| . r_, t|
|00001440| 5f 20 3a 3d 20 49 46 28 | 79 63 70 74 20 3d 20 31 |_ := IF(|ycpt = 1|
|00001450| 2c 20 79 5f 20 2d 20 56 | 45 43 54 4f 52 28 41 56 |, y_ - V|ECTOR(AV|
|00001460| 45 52 41 47 45 28 79 5f | 29 2c 20 69 2c 20 6f 62 |ERAGE(y_|), i, ob|
|00001470| 73 5f 29 2c 20 79 5f 29 | 2c 20 73 73 74 5f 20 3a |s_), y_)|, sst_ :|
|00001480| 3d 20 74 5f 20 95 20 74 | 5f 2c 20 bf 5f 20 3a 3d |= t_ . t|_, ._ :=|
|00001490| 20 30 2e 30 30 31 2c 20 | 69 74 65 72 20 3a 3d 20 | 0.001, |iter := |
|000014a0| 5b 41 50 50 45 4e 44 28 | 5b 30 5d 2c 20 5b 22 20 |[APPEND(|[0], [" |
|000014b0| 22 5d 2c 20 df 5f 2c 20 | 5b 73 73 65 5f 5d 29 5d |"], ._, |[sse_])]|
|000014c0| 2c 20 65 70 73 5f 20 3a | 3d 20 31 2c 20 69 74 65 |, eps_ :|= 1, ite|
|000014d0| 72 5f 20 3a 3d 20 31 2c | 20 4c 4f 4f 50 28 49 46 |r_ := 1,| LOOP(IF|
|000014e0| 28 69 74 65 72 5f 20 3e | 20 6e 5f 20 90 20 65 70 |(iter_ >| n_ . ep|
|000014f0| 73 5f 20 3c 20 63 63 5f | 2c 20 65 78 69 74 29 2c |s_ < cc_|, exit),|
|00001500| 20 78 5f 20 3a 3d 20 53 | 55 42 53 54 28 6a 5f 2c | x_ := S|UBST(j_,|
|00001510| 20 70 61 72 6d 2c 20 df | 5f 29 2c 20 6c 73 73 65 | parm, .|_), lsse|
|00001520| 5f 20 3a 3d 20 73 73 65 | 5f 2c 20 78 70 78 20 3a |_ := sse|_, xpx :|
|00001530| 3d 20 78 5f 60 20 95 20 | 78 5f 2c 20 64 69 61 67 |= x_` . |x_, diag|
|00001540| 5f 20 3a 3d 20 56 45 43 | 54 4f 52 28 56 45 43 54 |_ := VEC|TOR(VECT|
|00001550| 4f 52 28 49 46 28 71 5f | 20 3d 20 73 5f 2c 20 78 |OR(IF(q_| = s_, x|
|00001560| 70 78 99 71 5f 99 73 5f | 2c 20 30 29 2c 20 73 5f |px.q_.s_|, 0), s_|
|00001570| 2c 20 6b 5f 20 2b 20 79 | 63 70 74 29 2c 20 71 5f |, k_ + y|cpt), q_|
|00001580| 2c 20 6b 5f 20 2b 20 79 | 63 70 74 29 2c 20 78 70 |, k_ + y|cpt), xp|
|00001590| 78 69 20 3a 3d 20 31 2f | 28 78 70 78 20 2b 20 bf |xi := 1/|(xpx + .|
|000015a0| 5f b7 64 69 61 67 5f 29 | 2c 20 ab 5f 20 3a 3d 20 |_.diag_)|, ._ := |
|000015b0| 78 70 78 69 20 95 20 78 | 5f 60 20 95 20 72 5f 2c |xpxi . x|_` . r_,|
|000015c0| 20 df 5f 20 3a 2b 20 ab | 5f 2c 20 70 5f 20 3a 3d | ._ :+ .|_, p_ :=|
|000015d0| 20 53 55 42 53 54 28 65 | 71 31 2c 20 70 61 72 6d | SUBST(e|q1, parm|
|000015e0| 2c 20 df 5f 29 2c 20 72 | 5f 20 3a 3d 20 79 5f 20 |, ._), r|_ := y_ |
|000015f0| 2d 20 70 5f 2c 20 73 73 | 65 5f 20 3a 3d 20 72 5f |- p_, ss|e_ := r_|
|00001600| 20 95 20 72 5f 2c 20 69 | 74 65 72 20 3a 3d 20 41 | . r_, i|ter := A|
|00001610| 50 50 45 4e 44 28 69 74 | 65 72 2c 20 5b 41 50 50 |PPEND(it|er, [APP|
|00001620| 45 4e 44 28 5b 69 74 65 | 72 5f 5d 2c 20 5b bf 5f |END([ite|r_], [._|
|00001630| 5d 2c 20 df 5f 2c 20 5b | 73 73 65 5f 5d 29 5d 29 |], ._, [|sse_])])|
|00001640| 2c 20 49 46 28 73 73 65 | 5f 20 93 20 6c 73 73 65 |, IF(sse|_ . lsse|
|00001650| 5f 2c 20 50 52 4f 47 28 | 65 70 73 5f 20 3a 3d 20 |_, PROG(|eps_ := |
|00001660| 41 42 53 28 28 6c 73 73 | 65 5f 20 2d 20 73 73 65 |ABS((lss|e_ - sse|
|00001670| 5f 29 2f 28 73 73 65 5f | 20 2b 20 31 30 5e 28 2d |_)/(sse_| + 10^(-|
|00001680| 36 29 29 29 2c 20 bf 5f | 20 3a 2f 20 31 30 29 2c |6))), ._| :/ 10),|
|00001690| 20 50 52 4f 47 28 df 5f | 20 3a 2d 20 ab 5f 2c 20 | PROG(._| :- ._, |
|000016a0| bf 5f 20 3a 2a 20 31 30 | 29 2c 20 52 45 54 55 52 |._ :* 10|), RETUR|
|000016b0| 4e 20 22 43 61 6e 6e 6f | 74 20 72 65 73 6f 6c 76 |N "Canno|t resolv|
|000016c0| 65 20 53 53 45 20 74 65 | 73 74 21 22 29 2c 20 69 |e SSE te|st!"), i|
|000016d0| 74 65 72 5f 20 3a 2b 20 | 31 29 2c 20 49 46 28 69 |ter_ :+ |1), IF(i|
|000016e0| 74 65 72 5f 20 3e 20 6e | 5f 2c 20 6d 73 67 20 3a |ter_ > n|_, msg :|
|000016f0| 3d 20 41 50 50 45 4e 44 | 28 22 43 6f 6e 76 65 72 |= APPEND|("Conver|
|00001700| 67 65 6e 63 65 20 66 61 | 69 6c 65 64 20 61 66 74 |gence fa|iled aft|
|00001710| 65 72 20 22 2c 20 6e 5f | 2c 20 22 20 69 74 65 72 |er ", n_|, " iter|
|00001720| 61 74 69 6f 6e 73 21 22 | 29 2c 20 6d 73 67 20 3a |ations!"|), msg :|
|00001730| 3d 20 22 43 6f 6e 76 65 | 72 67 65 6e 63 65 20 63 |= "Conve|rgence c|
|00001740| 72 69 74 65 72 69 61 20 | 6d 65 74 21 22 29 2c 20 |riteria |met!"), |
|00001750| 64 66 65 5f 20 3a 3d 20 | 6f 62 73 5f 20 2d 20 6b |dfe_ := |obs_ - k|
|00001760| 5f 20 2d 20 79 63 70 74 | 2c 20 6d 73 65 5f 20 3a |_ - ycpt|, mse_ :|
|00001770| 3d 20 73 73 65 5f 2f 64 | 66 65 5f 2c 20 73 65 5f |= sse_/d|fe_, se_|
|00001780| 20 3a 3d 20 8b 6d 73 65 | 5f 2c 20 73 74 64 5f 20 | := .mse|_, std_ |
|00001790| 3a 3d 20 56 45 43 54 4f | 52 28 8b 28 78 70 78 69 |:= VECTO|R(.(xpxi|
|000017a0| 99 69 5f 99 69 5f b7 6d | 73 65 5f 29 2c 20 69 5f |.i_.i_.m|se_), i_|
|000017b0| 2c 20 6b 5f 20 2b 20 79 | 63 70 74 29 2c 20 74 5f |, k_ + y|cpt), t_|
|000017c0| 20 3a 3d 20 56 45 43 54 | 4f 52 28 df 5f 99 69 5f | := VECT|OR(._.i_|
|000017d0| 2f 73 74 64 5f 99 69 5f | 2c 20 69 5f 2c 20 6b 5f |/std_.i_|, i_, k_|
|000017e0| 20 2b 20 79 63 70 74 29 | 2c 20 70 72 6f 62 5f 20 | + ycpt)|, prob_ |
|000017f0| 3a 3d 20 56 45 43 54 4f | 52 28 31 20 2d 20 53 54 |:= VECTO|R(1 - ST|
|00001800| 55 44 45 4e 54 28 74 5f | 99 69 5f 2c 20 64 66 65 |UDENT(t_|.i_, dfe|
|00001810| 5f 29 2c 20 69 5f 2c 20 | 6b 5f 20 2b 20 79 63 70 |_), i_, |k_ + ycp|
|00001820| 74 29 2c 20 73 73 72 5f | 20 3a 3d 20 73 73 74 5f |t), ssr_| := sst_|
|00001830| 20 2d 20 73 73 65 5f 2c | 20 72 73 71 20 3a 3d 20 | - sse_,| rsq := |
|00001840| 73 73 72 5f 2f 73 73 74 | 5f 2c 20 61 72 73 71 20 |ssr_/sst|_, arsq |
|00001850| 3a 3d 20 31 20 2d 20 73 | 73 65 5f b7 28 6f 62 73 |:= 1 - s|se_.(obs|
|00001860| 5f 20 2d 20 79 63 70 74 | 29 2f 28 28 6f 62 73 5f |_ - ycpt|)/((obs_|
|00001870| 20 2d 20 6b 5f 20 2d 20 | 79 63 70 74 29 b7 73 73 | - k_ - |ycpt).ss|
|00001880| 74 5f 29 2c 20 66 64 5f | 20 3a 3d 20 46 5f 44 49 |t_), fd_| := F_DI|
|00001890| 53 54 52 49 42 55 54 49 | 4f 4e 28 73 73 72 5f b7 |STRIBUTI|ON(ssr_.|
|000018a0| 64 66 65 5f 2f 28 73 73 | 65 5f b7 6b 5f 29 2c 20 |dfe_/(ss|e_.k_), |
|000018b0| 6b 5f 2c 20 64 66 65 5f | 29 2c 20 61 31 5f 20 3a |k_, dfe_|), a1_ :|
|000018c0| 3d 20 5b 22 53 6f 75 72 | 63 65 22 2c 20 22 44 46 |= ["Sour|ce", "DF|
|000018d0| 22 2c 20 22 53 53 22 2c | 20 22 4d 53 22 2c 20 22 |", "SS",| "MS", "|
|000018e0| 46 22 2c 20 22 50 72 6f | 62 28 46 29 22 5d 2c 20 |F", "Pro|b(F)"], |
|000018f0| 61 32 5f 20 3a 3d 20 5b | 22 52 65 67 22 2c 20 6b |a2_ := [|"Reg", k|
|00001900| 5f 2c 20 73 73 72 5f 2c | 20 73 73 72 5f 2f 6b 5f |_, ssr_,| ssr_/k_|
|00001910| 2c 20 73 73 72 5f b7 64 | 66 65 5f 2f 28 73 73 65 |, ssr_.d|fe_/(sse|
|00001920| 5f b7 6b 5f 29 2c 20 49 | 46 28 66 64 5f 20 3c 20 |_.k_), I|F(fd_ < |
|00001930| 30 2c 20 30 2c 20 66 64 | 5f 29 5d 2c 20 61 33 5f |0, 0, fd|_)], a3_|
|00001940| 20 3a 3d 20 5b 22 45 72 | 72 6f 72 22 2c 20 64 66 | := ["Er|ror", df|
|00001950| 65 5f 2c 20 73 73 65 5f | 2c 20 73 73 65 5f 2f 64 |e_, sse_|, sse_/d|
|00001960| 66 65 5f 2c 20 22 20 22 | 2c 20 22 20 22 5d 2c 20 |fe_, " "|, " "], |
|00001970| 61 34 5f 20 3a 3d 20 5b | 49 46 28 79 63 70 74 20 |a4_ := [|IF(ycpt |
|00001980| 3d 20 31 2c 20 22 43 6f | 72 72 65 63 74 65 64 20 |= 1, "Co|rrected |
|00001990| 54 6f 74 61 6c 22 2c 20 | 22 55 6e 63 6f 72 72 65 |Total", |"Uncorre|
|000019a0| 63 74 65 64 20 54 6f 74 | 61 6c 22 29 2c 20 64 66 |cted Tot|al"), df|
|000019b0| 65 5f 20 2b 20 6b 5f 2c | 20 73 73 74 5f 2c 20 22 |e_ + k_,| sst_, "|
|000019c0| 20 22 2c 20 22 20 22 2c | 20 22 20 22 5d 2c 20 61 | ", " ",| " "], a|
|000019d0| 6e 6f 76 61 20 3a 3d 20 | 41 50 50 45 4e 44 28 5b |nova := |APPEND([|
|000019e0| 61 31 5f 5d 2c 20 5b 61 | 32 5f 5d 2c 20 5b 61 33 |a1_], [a|2_], [a3|
|000019f0| 5f 5d 2c 20 5b 61 34 5f | 5d 29 2c 20 66 65 71 20 |_], [a4_|]), feq |
|00001a00| 3a 3d 20 64 76 61 72 20 | 3d 20 53 55 42 53 54 28 |:= dvar |= SUBST(|
|00001a10| 65 71 71 2c 20 70 61 72 | 6d 2c 20 df 5f 29 2c 20 |eqq, par|m, ._), |
|00001a20| 69 74 65 72 20 3a 3d 20 | 41 50 50 45 4e 44 28 5b |iter := |APPEND([|
|00001a30| 41 50 50 45 4e 44 28 5b | 22 49 74 65 72 22 5d 2c |APPEND([|"Iter"],|
|00001a40| 20 5b 22 bf 22 5d 2c 20 | 70 61 72 6d 2c 20 5b 22 | ["."], |parm, ["|
|00001a50| 53 53 45 22 5d 29 5d 2c | 20 69 74 65 72 29 2c 20 |SSE"])],| iter), |
|00001a60| 74 69 74 6c 65 20 3a 3d | 20 5b 22 50 61 72 6d 22 |title :=| ["Parm"|
|00001a70| 3b 20 22 56 61 6c 75 65 | 22 3b 20 22 53 54 44 22 |; "Value|"; "STD"|
|00001a80| 3b 20 41 50 50 45 4e 44 | 28 22 74 28 22 2c 20 64 |; APPEND|("t(", d|
|00001a90| 66 65 5f 2c 20 22 29 22 | 29 3b 20 22 50 72 6f 62 |fe_, ")"|); "Prob|
|00001aa0| 28 74 29 22 5d 2c 20 6f | 75 74 70 75 74 20 3a 3d |(t)"], o|utput :=|
|00001ab0| 20 41 50 50 45 4e 44 5f | 43 4f 4c 55 4d 4e 53 28 | APPEND_|COLUMNS(|
|00001ac0| 74 69 74 6c 65 2c 20 41 | 50 50 45 4e 44 28 5b 70 |title, A|PPEND([p|
|00001ad0| 61 72 6d 5d 2c 20 5b df | 5f 5d 2c 20 5b 73 74 64 |arm], [.|_], [std|
|00001ae0| 5f 5d 2c 20 5b 74 5f 5d | 2c 20 5b 70 72 6f 62 5f |_], [t_]|, [prob_|
|00001af0| 5d 29 29 60 2c 20 73 74 | 61 74 73 20 3a 3d 20 5b |]))`, st|ats := [|
|00001b00| 22 53 45 22 2c 20 22 52 | 5e 32 22 2c 20 22 41 64 |"SE", "R|^2", "Ad|
|00001b10| 6a 52 5e 32 22 3b 20 73 | 65 5f 2c 20 72 73 71 2c |jR^2"; s|e_, rsq,|
|00001b20| 20 61 72 73 71 5d 2c 20 | 52 45 54 55 52 4e 20 5b | arsq], |RETURN [|
|00001b30| 22 4d 61 72 71 75 61 72 | 64 74 20 4d 65 74 68 6f |"Marquar|dt Metho|
|00001b40| 64 22 3b 20 22 20 22 3b | 20 6d 73 67 3b 20 22 20 |d"; " ";| msg; " |
|00001b50| 22 3b 20 6f 75 74 70 75 | 74 3b 20 22 20 22 3b 20 |"; outpu|t; " "; |
|00001b60| 61 6e 6f 76 61 3b 20 22 | 20 22 3b 20 73 74 61 74 |anova; "| "; stat|
|00001b70| 73 3b 20 22 20 22 3b 20 | 66 65 71 5d 29 0d 0a 4d |s; " "; |feq])..M|
|00001b80| 41 54 50 52 4f 44 28 61 | 2c 20 62 2c 20 69 2c 20 |ATPROD(a|, b, i, |
|00001b90| 6a 29 3a 3d a4 28 61 99 | 69 99 6e 5f b7 62 99 6e |j):=.(a.|i.n_.b.n|
|00001ba0| 5f 99 6a 2c 20 6e 5f 2c | 20 31 2c 20 44 49 4d 28 |_.j, n_,| 1, DIM(|
|00001bb0| 62 29 29 0d 0a 4d 49 4e | 4f 52 28 61 2c 20 69 2c |b))..MIN|OR(a, i,|
|00001bc0| 20 6a 29 3a 3d 44 45 4c | 45 54 45 28 44 45 4c 45 | j):=DEL|ETE(DELE|
|00001bd0| 54 45 28 61 2c 20 69 29 | 60 2c 20 6a 29 60 0d 0a |TE(a, i)|`, j)`..|
|00001be0| 4f 55 54 45 52 28 76 2c | 20 77 29 3a 3d 56 45 43 |OUTER(v,| w):=VEC|
|00001bf0| 54 4f 52 28 5b 76 99 6e | 5f 5d 2c 20 6e 5f 2c 20 |TOR([v.n|_], n_, |
|00001c00| 44 49 4d 28 76 29 29 20 | 95 20 5b 77 5d 0d 0a 50 |DIM(v)) |. [w]..P|
|00001c10| 41 52 54 49 54 49 4f 4e | 28 76 2c 20 6e 2c 20 64 |ARTITION|(v, n, d|
|00001c20| 29 3a 3d 49 46 28 4e 55 | 4d 42 45 52 28 64 29 2c |):=IF(NU|MBER(d),|
|00001c30| 20 56 45 43 54 4f 52 28 | 56 45 43 54 4f 52 28 76 | VECTOR(|VECTOR(v|
|00001c40| 99 6d 5f 2c 20 6d 5f 2c | 20 6e 5f 2c 20 6e 5f 20 |.m_, m_,| n_, n_ |
|00001c50| 2b 20 6e 20 2d 20 31 29 | 2c 20 6e 5f 2c 20 31 2c |+ n - 1)|, n_, 1,|
|00001c60| 20 31 20 2b 20 64 b7 46 | 4c 4f 4f 52 28 44 49 4d | 1 + d.F|LOOR(DIM|
|00001c70| 28 76 29 20 2d 20 6e 2c | 20 64 29 2c 20 64 29 2c |(v) - n,| d), d),|
|00001c80| 20 50 41 52 54 49 54 49 | 4f 4e 28 76 2c 20 6e 2c | PARTITI|ON(v, n,|
|00001c90| 20 6e 29 29 0d 0a 50 49 | 56 4f 54 28 61 2c 20 69 | n))..PI|VOT(a, i|
|00001ca0| 2c 20 6a 29 3a 3d 56 45 | 43 54 4f 52 28 49 46 28 |, j):=VE|CTOR(IF(|
|00001cb0| 6d 5f 20 93 20 69 2c 20 | 61 99 6d 5f 2c 20 61 99 |m_ . i, |a.m_, a.|
|00001cc0| 6d 5f 20 2d 20 61 99 69 | b7 61 99 6d 5f 99 6a 2f |m_ - a.i|.a.m_.j/|
|00001cd0| 61 99 69 99 6a 29 2c 20 | 6d 5f 2c 20 44 49 4d 28 |a.i.j), |m_, DIM(|
|00001ce0| 61 29 29 0d 0a 50 4f 43 | 48 48 41 4d 4d 45 52 28 |a))..POC|HHAMMER(|
|00001cf0| 61 2c 20 78 29 3a 3d 50 | 45 52 4d 28 78 20 2b 20 |a, x):=P|ERM(x + |
|00001d00| 61 20 2d 20 31 2c 20 78 | 29 0d 0a 50 4f 49 53 53 |a - 1, x|)..POISS|
|00001d10| 4f 4e 5f 44 45 4e 53 49 | 54 59 28 6b 2c 20 74 29 |ON_DENSI|TY(k, t)|
|00001d20| 3a 3d ea 5e 28 2d 74 29 | b7 74 5e 6b 2f 6b 21 0d |:=.^(-t)|.t^k/k!.|
|00001d30| 0a 50 4f 49 53 53 4f 4e | 5f 44 49 53 54 52 49 42 |.POISSON|_DISTRIB|
|00001d40| 55 54 49 4f 4e 28 6b 2c | 20 74 29 3a 3d a4 28 50 |UTION(k,| t):=.(P|
|00001d50| 4f 49 53 53 4f 4e 5f 44 | 45 4e 53 49 54 59 28 6d |OISSON_D|ENSITY(m|
|00001d60| 5f 2c 20 74 29 2c 20 6d | 5f 2c 20 30 2c 20 6b 29 |_, t), m|_, 0, k)|
|00001d70| 0d 0a 50 4f 4c 59 47 41 | 4d 4d 41 28 6e 2c 20 7a |..POLYGA|MMA(n, z|
|00001d80| 29 3a 3d 49 46 28 6e 20 | 3d 20 2d 31 2c 20 4c 4e |):=IF(n |= -1, LN|
|00001d90| 28 9b 28 7a 29 29 2c 20 | 49 46 28 6e 20 3d 20 30 |(.(z)), |IF(n = 0|
|00001da0| 2c 20 44 49 47 41 4d 4d | 41 28 7a 29 2c 20 28 2d |, DIGAMM|A(z), (-|
|00001db0| 31 29 5e 28 6e 20 2b 20 | 31 29 b7 6e 21 b7 ae 28 |1)^(n + |1).n!..(|
|00001dc0| 6e 20 2b 20 31 2c 20 7a | 29 29 29 0d 0a 50 72 65 |n + 1, z|)))..Pre|
|00001dd0| 64 69 63 74 65 64 5f 56 | 61 6c 75 65 73 28 78 30 |dicted_V|alues(x0|
|00001de0| 2c 20 63 69 20 3a 3d 20 | 30 2e 39 35 29 3a 3d 50 |, ci := |0.95):=P|
|00001df0| 52 4f 47 28 70 72 64 20 | 3a 3d 20 53 55 42 53 54 |ROG(prd |:= SUBST|
|00001e00| 28 52 48 53 28 66 65 71 | 29 2c 20 69 76 61 72 73 |(RHS(feq|), ivars|
|00001e10| 31 2c 20 78 30 29 2c 20 | 67 5f 20 3a 3d 20 47 52 |1, x0), |g_ := GR|
|00001e20| 41 44 28 65 71 71 2c 20 | 70 61 72 6d 31 29 2c 20 |AD(eqq, |parm1), |
|00001e30| 67 76 20 3a 3d 20 53 55 | 42 53 54 28 67 5f 2c 20 |gv := SU|BST(g_, |
|00001e40| 41 50 50 45 4e 44 28 70 | 61 72 6d 31 2c 20 69 76 |APPEND(p|arm1, iv|
|00001e50| 61 72 73 31 29 2c 20 41 | 50 50 45 4e 44 28 df 5f |ars1), A|PPEND(._|
|00001e60| 2c 20 78 30 29 29 2c 20 | 69 5f 20 3a 3d 20 67 76 |, x0)), |i_ := gv|
|00001e70| 20 95 20 78 70 78 69 20 | 95 20 67 76 60 2c 20 73 | . xpxi |. gv`, s|
|00001e80| 65 70 20 3a 3d 20 8b 69 | 5f b7 73 65 5f 2c 20 73 |ep := .i|_.se_, s|
|00001e90| 65 66 20 3a 3d 20 73 65 | 5f b7 8b 28 31 20 2b 20 |ef := se|_..(1 + |
|00001ea0| 69 5f 29 2c 20 50 52 4f | 47 28 50 72 65 63 69 73 |i_), PRO|G(Precis|
|00001eb0| 69 6f 6e 44 69 67 69 74 | 73 20 3a 3d 20 36 2c 20 |ionDigit|s := 6, |
|00001ec0| 69 6e 76 74 20 3a 3d 20 | 52 48 53 28 4e 53 4f 4c |invt := |RHS(NSOL|
|00001ed0| 56 45 28 53 54 55 44 45 | 4e 54 28 74 2c 20 64 66 |VE(STUDE|NT(t, df|
|00001ee0| 65 5f 29 20 2d 20 63 69 | 2c 20 74 2c 20 30 2c 20 |e_) - ci|, t, 0, |
|00001ef0| 96 29 29 2c 20 50 72 65 | 63 69 73 69 6f 6e 44 69 |.)), Pre|cisionDi|
|00001f00| 67 69 74 73 20 3a 3d 20 | 31 30 2c 20 4e 6f 74 61 |gits := |10, Nota|
|00001f10| 74 69 6f 6e 44 69 67 69 | 74 73 20 3a 3d 20 38 29 |tionDigi|ts := 8)|
|00001f20| 2c 20 63 79 68 31 20 3a | 3d 20 70 72 64 20 2d 20 |, cyh1 :|= prd - |
|00001f30| 69 6e 76 74 b7 73 65 70 | 2c 20 63 79 68 32 20 3a |invt.sep|, cyh2 :|
|00001f40| 3d 20 70 72 64 20 2b 20 | 69 6e 76 74 b7 73 65 70 |= prd + |invt.sep|
|00001f50| 2c 20 63 79 30 31 20 3a | 3d 20 70 72 64 20 2d 20 |, cy01 :|= prd - |
|00001f60| 69 6e 76 74 b7 73 65 66 | 2c 20 63 79 30 32 20 3a |invt.sef|, cy02 :|
|00001f70| 3d 20 70 72 64 20 2b 20 | 69 6e 76 74 b7 73 65 66 |= prd + |invt.sef|
|00001f80| 2c 20 74 69 74 6c 65 20 | 3a 3d 20 41 50 50 45 4e |, title |:= APPEN|
|00001f90| 44 28 53 54 52 49 4e 47 | 28 31 30 30 b7 63 69 29 |D(STRING|(100.ci)|
|00001fa0| 2c 20 22 25 20 43 6f 6e | 66 69 64 65 6e 63 65 20 |, "% Con|fidence |
|00001fb0| 49 6e 74 65 72 76 61 6c | 22 29 2c 20 5b 74 69 74 |Interval|"), [tit|
|00001fc0| 6c 65 3b 20 5b 22 56 61 | 6c 75 65 22 2c 20 70 72 |le; ["Va|lue", pr|
|00001fd0| 64 2c 20 22 20 22 3b 20 | 22 53 65 5f 79 68 61 74 |d, " "; |"Se_yhat|
|00001fe0| 2f 53 65 5f 59 30 22 2c | 20 73 65 70 2c 20 73 65 |/Se_Y0",| sep, se|
|00001ff0| 66 3b 20 22 43 49 5f 79 | 68 61 74 22 2c 20 63 79 |f; "CI_y|hat", cy|
|00002000| 68 31 2c 20 63 79 68 32 | 3b 20 22 43 49 5f 59 30 |h1, cyh2|; "CI_Y0|
|00002010| 22 2c 20 63 79 30 31 2c | 20 63 79 30 32 5d 5d 29 |", cy01,| cy02]])|
|00002020| 0d 0a 53 43 41 4c 45 5f | 45 4c 45 4d 45 4e 54 28 |..SCALE_|ELEMENT(|
|00002030| 76 2c 20 69 2c 20 73 29 | 3a 3d 56 45 43 54 4f 52 |v, i, s)|:=VECTOR|
|00002040| 28 49 46 28 6d 5f 20 3d | 20 69 2c 20 73 b7 76 99 |(IF(m_ =| i, s.v.|
|00002050| 69 2c 20 76 99 6d 5f 29 | 2c 20 6d 5f 2c 20 44 49 |i, v.m_)|, m_, DI|
|00002060| 4d 28 76 29 29 0d 0a 53 | 54 55 44 45 4e 54 28 74 |M(v))..S|TUDENT(t|
|00002070| 2c 20 76 29 3a 3d 31 20 | 2d 20 49 4e 43 4f 4d 50 |, v):=1 |- INCOMP|
|00002080| 4c 45 54 45 5f 42 45 54 | 41 28 76 2f 28 76 20 2b |LETE_BET|A(v/(v +|
|00002090| 20 74 5e 32 29 2c 20 76 | 2f 32 2c 20 31 2f 32 29 | t^2), v|/2, 1/2)|
|000020a0| 0d 0a 53 55 42 54 52 41 | 43 54 5f 45 4c 45 4d 45 |..SUBTRA|CT_ELEME|
|000020b0| 4e 54 53 28 76 2c 20 69 | 2c 20 6a 2c 20 73 20 3a |NTS(v, i|, j, s :|
|000020c0| 3d 20 31 29 3a 3d 56 45 | 43 54 4f 52 28 49 46 28 |= 1):=VE|CTOR(IF(|
|000020d0| 6d 5f 20 3d 20 69 2c 20 | 76 99 69 20 2d 20 73 b7 |m_ = i, |v.i - s.|
|000020e0| 76 99 6a 2c 20 76 99 6d | 5f 29 2c 20 6d 5f 2c 20 |v.j, v.m|_), m_, |
|000020f0| 44 49 4d 28 76 29 29 0d | 0a 53 57 41 50 5f 45 4c |DIM(v)).|.SWAP_EL|
|00002100| 45 4d 45 4e 54 53 28 76 | 2c 20 69 2c 20 6a 29 3a |EMENTS(v|, i, j):|
|00002110| 3d 56 45 43 54 4f 52 28 | 49 46 28 6d 5f 20 3d 20 |=VECTOR(|IF(m_ = |
|00002120| 69 2c 20 76 99 6a 2c 20 | 49 46 28 6d 5f 20 3d 20 |i, v.j, |IF(m_ = |
|00002130| 6a 2c 20 76 99 69 2c 20 | 76 99 6d 5f 29 29 2c 20 |j, v.i, |v.m_)), |
|00002140| 6d 5f 2c 20 44 49 4d 28 | 76 29 29 0d 0a 70 6f 73 |m_, DIM(|v))..pos|
|00002150| 69 74 69 6f 6e 28 65 2c | 20 75 29 3a 3d 50 52 4f |ition(e,| u):=PRO|
|00002160| 47 28 64 20 3a 3d 20 44 | 49 4d 28 75 29 2c 20 69 |G(d := D|IM(u), i|
|00002170| 5f 20 3a 3d 20 31 2c 20 | 4c 4f 4f 50 28 49 46 28 |_ := 1, |LOOP(IF(|
|00002180| 65 20 3d 20 75 99 69 5f | 2c 20 52 45 54 55 52 4e |e = u.i_|, RETURN|
|00002190| 20 69 5f 29 2c 20 49 46 | 28 69 5f 20 3e 20 64 2c | i_), IF|(i_ > d,|
|000021a0| 20 52 45 54 55 52 4e 20 | 30 29 2c 20 69 5f 20 3a | RETURN |0), i_ :|
|000021b0| 2b 20 31 29 29 0d 0a 52 | 65 73 69 64 75 61 6c 73 |+ 1))..R|esiduals|
|000021c0| 3a 3d 56 45 43 54 4f 52 | 28 5b 69 2c 20 72 5f 99 |:=VECTOR|([i, r_.|
|000021d0| 69 5d 2c 20 69 2c 20 31 | 2c 20 44 49 4d 28 72 5f |i], i, 1|, DIM(r_|
|000021e0| 29 29 0d 0a 61 6e 6f 76 | 61 3a 3d 5b 22 53 6f 75 |))..anov|a:=["Sou|
|000021f0| 72 63 65 22 2c 20 22 44 | 46 22 2c 20 22 53 53 22 |rce", "D|F", "SS"|
|00002200| 2c 20 22 4d 53 22 2c 20 | 22 46 22 2c 20 22 50 72 |, "MS", |"F", "Pr|
|00002210| 6f 62 28 46 29 22 3b 20 | 22 52 65 67 22 2c 20 33 |ob(F)"; |"Reg", 3|
|00002220| 2c 20 31 2e 36 34 32 32 | 37 38 39 b7 31 30 5e 35 |, 1.6422|789.10^5|
|00002230| 2c 20 35 2e 34 37 34 32 | 36 33 33 b7 31 30 5e 34 |, 5.4742|633.10^4|
|00002240| 2c 20 34 39 33 38 2e 31 | 36 39 34 2c 20 38 2e 39 |, 4938.1|694, 8.9|
|00002250| 30 33 30 31 38 34 b7 31 | 30 5e 28 2d 31 31 29 3b |030184.1|0^(-11);|
|00002260| 20 22 45 72 72 6f 72 22 | 2c 20 31 36 2c 20 31 37 | "Error"|, 16, 17|
|00002270| 37 2e 33 36 39 38 30 2c | 20 31 31 2e 30 38 35 36 |7.36980,| 11.0856|
|00002280| 31 32 2c 20 22 20 22 2c | 20 22 20 22 3b 20 22 55 |12, " ",| " "; "U|
|00002290| 6e 63 6f 72 72 65 63 74 | 65 64 20 54 6f 74 61 6c |ncorrect|ed Total|
|000022a0| 22 2c 20 31 39 2c 20 31 | 2e 36 34 34 30 35 32 36 |", 19, 1|.6440526|
|000022b0| b7 31 30 5e 35 2c 20 22 | 20 22 2c 20 22 20 22 2c |.10^5, "| ", " ",|
|000022c0| 20 22 20 22 5d 0d 0a 61 | 72 73 71 3a 3d 30 2e 39 | " "]..a|rsq:=0.9|
|000022d0| 39 38 37 31 38 38 35 0d | 0a 63 63 5f 3a 3d 0d 0a |9871885.|.cc_:=..|
|000022e0| 63 69 3a 3d 0d 0a 64 61 | 74 61 3a 3d 0d 0a 64 66 |ci:=..da|ta:=..df|
|000022f0| 65 5f 3a 3d 31 36 0d 0a | 64 76 61 72 3a 3d 70 0d |e_:=16..|dvar:=p.|
|00002300| 0a 65 71 3a 3d 0d 0a 65 | 71 71 3a 3d 30 2e 35 b7 |.eq:=..e|qq:=0.5.|
|00002310| 63 b7 28 45 52 46 28 30 | 2e 37 30 37 31 30 36 37 |c.(ERF(0|.7071067|
|00002320| 38 b7 28 62 b7 79 20 2b | 20 61 20 2d 20 31 37 39 |8.(b.y +| a - 179|
|00002330| 30 b7 62 29 29 20 2b 20 | 31 29 0d 0a 66 5f 3a 3d |0.b)) + |1)..f_:=|
|00002340| 0d 0a 66 65 71 3a 3d 70 | 20 3d 20 32 30 33 2e 35 |..feq:=p| = 203.5|
|00002350| 35 31 39 36 b7 28 45 52 | 46 28 31 2e 33 32 34 30 |5196.(ER|F(1.3240|
|00002360| 30 37 38 b7 31 30 5e 28 | 2d 31 34 29 b7 28 36 2e |078.10^(|-14).(6.|
|00002370| 37 34 34 32 37 30 36 b7 | 31 30 5e 31 31 b7 79 20 |7442706.|10^11.y |
|00002380| 2d 20 31 2e 33 33 30 32 | 30 39 39 b7 31 30 5e 31 |- 1.3302|099.10^1|
|00002390| 35 29 29 20 2b 20 31 29 | 0d 0a 68 43 72 6f 73 73 |5)) + 1)|..hCross|
|000023a0| 3a 3d 41 50 50 52 4f 58 | 28 31 2e 39 39 34 30 34 |:=APPROX|(1.99404|
|000023b0| 37 36 29 0d 0a 69 74 65 | 72 3a 3d 5b 22 49 74 65 |76)..ite|r:=["Ite|
|000023c0| 72 22 2c 20 22 bf 22 2c | 20 61 2c 20 62 2c 20 63 |r", ".",| a, b, c|
|000023d0| 2c 20 22 53 53 45 22 3b | 20 30 2c 20 22 20 22 2c |, "SSE";| 0, " ",|
|000023e0| 20 2d 32 2e 34 2c 20 30 | 2e 30 31 32 2c 20 34 30 | -2.4, 0|.012, 40|
|000023f0| 30 2c 20 37 31 37 34 2e | 35 39 30 38 3b 20 31 2c |0, 7174.|5908; 1,|
|00002400| 20 30 2e 30 30 31 2c 20 | 2d 32 2e 32 37 36 39 39 | 0.001, |-2.27699|
|00002410| 37 32 2c 20 30 2e 30 31 | 32 34 35 38 39 36 37 2c |72, 0.01|2458967,|
|00002420| 20 34 31 32 2e 38 34 36 | 33 30 2c 20 32 32 32 2e | 412.846|30, 222.|
|00002430| 33 32 35 39 37 3b 20 32 | 2c 20 30 2e 30 30 30 31 |32597; 2|, 0.0001|
|00002440| 2c 20 2d 32 2e 33 30 32 | 34 39 38 32 2c 20 30 2e |, -2.302|4982, 0.|
|00002450| 30 31 32 36 34 33 36 30 | 32 2c 20 34 30 35 2e 39 |01264360|2, 405.9|
|00002460| 39 30 33 38 2c 20 31 37 | 37 2e 33 37 33 38 37 3b |9038, 17|7.37387;|
|00002470| 20 33 2c 20 31 30 5e 28 | 2d 35 29 2c 20 2d 32 2e | 3, 10^(|-5), -2.|
|00002480| 33 30 32 37 39 32 35 2c | 20 30 2e 30 31 32 36 32 |3027925,| 0.01262|
|00002490| 38 35 34 39 2c 20 34 30 | 37 2e 30 36 37 30 30 2c |8549, 40|7.06700,|
|000024a0| 20 31 37 37 2e 33 36 39 | 38 31 3b 20 34 2c 20 31 | 177.369|81; 4, 1|
|000024b0| 30 5e 28 2d 36 29 2c 20 | 2d 32 2e 33 30 32 38 31 |0^(-6), |-2.30281|
|000024c0| 37 38 2c 20 30 2e 30 31 | 32 36 32 38 31 37 33 2c |78, 0.01|2628173,|
|000024d0| 20 34 30 37 2e 31 30 33 | 39 32 2c 20 31 37 37 2e | 407.103|92, 177.|
|000024e0| 33 36 39 38 30 3b 20 35 | 2c 20 31 30 5e 28 2d 37 |36980; 5|, 10^(-7|
|000024f0| 29 2c 20 2d 32 2e 33 30 | 32 38 31 38 33 2c 20 30 |), -2.30|28183, 0|
|00002500| 2e 30 31 32 36 32 38 31 | 34 33 2c 20 34 30 37 2e |.0126281|43, 407.|
|00002510| 31 30 36 30 33 2c 20 31 | 37 37 2e 33 36 39 38 30 |10603, 1|77.36980|
|00002520| 3b 20 36 2c 20 31 30 5e | 28 2d 36 29 2c 20 2d 32 |; 6, 10^|(-6), -2|
|00002530| 2e 33 30 32 38 31 37 38 | 2c 20 30 2e 30 31 32 36 |.3028178|, 0.0126|
|00002540| 32 38 31 37 33 2c 20 34 | 30 37 2e 31 30 33 39 32 |28173, 4|07.10392|
|00002550| 2c 20 31 37 37 2e 33 36 | 39 38 30 5d 0d 0a 69 76 |, 177.36|980]..iv|
|00002560| 61 72 73 3a 3d 5b 79 5d | 0d 0a 69 76 61 72 73 31 |ars:=[y]|..ivars1|
|00002570| 3a 3d 5b 79 5d 0d 0a 6d | 73 67 3a 3d 22 43 6f 6e |:=[y]..m|sg:="Con|
|00002580| 76 65 72 67 65 6e 63 65 | 20 63 72 69 74 65 72 69 |vergence| criteri|
|00002590| 61 20 6d 65 74 21 22 0d | 0a 6e 5f 3a 3d 0d 0a 6f |a met!".|.n_:=..o|
|000025a0| 75 74 70 75 74 3a 3d 5b | 22 50 61 72 6d 22 2c 20 |utput:=[|"Parm", |
|000025b0| 22 56 61 6c 75 65 22 2c | 20 22 53 54 44 22 2c 20 |"Value",| "STD", |
|000025c0| 22 74 28 31 36 29 22 2c | 20 22 50 72 6f 62 28 74 |"t(16)",| "Prob(t|
|000025d0| 29 22 3b 20 61 2c 20 2d | 32 2e 33 30 32 38 31 37 |)"; a, -|2.302817|
|000025e0| 38 2c 20 30 2e 30 33 32 | 38 34 35 33 35 36 2c 20 |8, 0.032|845356, |
|000025f0| 2d 37 30 2e 31 31 30 39 | 31 30 2c 20 30 3b 20 62 |-70.1109|10, 0; b|
|00002600| 2c 20 30 2e 30 31 32 36 | 32 38 31 37 33 2c 20 30 |, 0.0126|28173, 0|
|00002610| 2e 30 30 30 39 34 31 35 | 32 39 36 37 2c 20 31 33 |.0009415|2967, 13|
|00002620| 2e 34 31 32 34 30 31 2c | 20 30 3b 20 63 2c 20 34 |.412401,| 0; c, 4|
|00002630| 30 37 2e 31 30 33 39 32 | 2c 20 36 30 2e 36 34 35 |07.10392|, 60.645|
|00002640| 30 32 35 2c 20 36 2e 37 | 31 32 38 39 38 39 2c 20 |025, 6.7|128989, |
|00002650| 34 2e 39 38 34 31 33 39 | 37 b7 31 30 5e 28 2d 36 |4.984139|7.10^(-6|
|00002660| 29 5d 0d 0a 70 61 72 6d | 3a 3d 0d 0a 70 61 72 6d |)]..parm|:=..parm|
|00002670| 31 3a 3d 5b 61 2c 20 62 | 2c 20 63 5d 0d 0a 70 6f |1:=[a, b|, c]..po|
|00002680| 70 3a 3d 5b 70 2c 20 79 | 3b 20 33 2e 39 32 39 2c |p:=[p, y|; 3.929,|
|00002690| 20 31 37 39 30 3b 20 35 | 2e 33 30 38 2c 20 31 38 | 1790; 5|.308, 18|
|000026a0| 30 30 3b 20 37 2e 32 33 | 39 2c 20 31 38 31 30 3b |00; 7.23|9, 1810;|
|000026b0| 20 39 2e 36 33 38 2c 20 | 31 38 32 30 3b 20 31 32 | 9.638, |1820; 12|
|000026c0| 2e 38 36 36 2c 20 31 38 | 33 30 3b 20 31 37 2e 30 |.866, 18|30; 17.0|
|000026d0| 36 39 2c 20 31 38 34 30 | 3b 20 32 33 2e 31 39 31 |69, 1840|; 23.191|
|000026e0| 2c 20 31 38 35 30 3b 20 | 33 31 2e 34 34 33 2c 20 |, 1850; |31.443, |
|000026f0| 31 38 36 30 3b 20 33 39 | 2e 38 31 38 2c 20 31 38 |1860; 39|.818, 18|
|00002700| 37 30 3b 20 35 30 2e 31 | 35 35 2c 20 31 38 38 30 |70; 50.1|55, 1880|
|00002710| 3b 20 36 32 2e 39 34 37 | 2c 20 31 38 39 30 3b 20 |; 62.947|, 1890; |
|00002720| 37 35 2e 39 39 34 2c 20 | 31 39 30 30 3b 20 39 31 |75.994, |1900; 91|
|00002730| 2e 39 37 32 2c 20 31 39 | 31 30 3b 20 31 30 35 2e |.972, 19|10; 105.|
|00002740| 37 31 2c 20 31 39 32 30 | 3b 20 31 32 32 2e 37 37 |71, 1920|; 122.77|
|00002750| 35 2c 20 31 39 33 30 3b | 20 31 33 31 2e 36 36 39 |5, 1930;| 131.669|
|00002760| 2c 20 31 39 34 30 3b 20 | 31 35 31 2e 33 32 35 2c |, 1940; |151.325,|
|00002770| 20 31 39 35 30 3b 20 31 | 37 39 2e 33 32 33 2c 20 | 1950; 1|79.323, |
|00002780| 31 39 36 30 3b 20 32 30 | 33 2e 32 31 31 2c 20 31 |1960; 20|3.211, 1|
|00002790| 39 37 30 5d 0d 0a 72 5f | 3a 3d 5b 2d 30 2e 34 30 |970]..r_|:=[-0.40|
|000027a0| 34 34 33 36 38 35 2c 20 | 2d 30 2e 36 39 39 38 37 |443685, |-0.69987|
|000027b0| 36 31 39 2c 20 2d 30 2e | 39 37 32 32 30 37 37 31 |619, -0.|97220771|
|000027c0| 2c 20 2d 31 2e 34 32 36 | 36 38 30 33 2c 20 2d 31 |, -1.426|6803, -1|
|000027d0| 2e 38 33 35 37 34 34 39 | 2c 20 2d 32 2e 31 39 35 |.8357449|, -2.195|
|000027e0| 33 33 38 34 2c 20 2d 31 | 2e 37 30 36 35 34 38 38 |3384, -1|.7065488|
|000027f0| 2c 20 2d 30 2e 32 39 39 | 37 30 39 32 37 2c 20 2d |, -0.299|70927, -|
|00002800| 30 2e 31 31 31 31 34 32 | 38 30 2c 20 30 2e 35 39 |0.111142|80, 0.59|
|00002810| 30 30 32 33 32 34 2c 20 | 32 2e 32 31 39 33 37 32 |002324, |2.219372|
|00002820| 35 2c 20 32 2e 35 33 39 | 32 39 31 30 2c 20 34 2e |5, 2.539|2910, 4.|
|00002830| 32 33 35 37 39 32 35 2c | 20 32 2e 32 30 31 32 33 |2357925,| 2.20123|
|00002840| 32 30 2c 20 32 2e 31 32 | 32 31 35 35 38 2c 20 2d |20, 2.12|21558, -|
|00002850| 37 2e 33 32 34 32 37 31 | 31 2c 20 2d 36 2e 39 37 |7.324271|1, -6.97|
|00002860| 38 35 33 37 37 2c 20 31 | 2e 30 30 39 33 30 31 39 |85377, 1|.0093019|
|00002870| 2c 20 34 2e 34 38 39 34 | 39 35 31 5d 0d 0a 72 73 |, 4.4894|951]..rs|
|00002880| 71 3a 3d 30 2e 39 39 38 | 39 32 31 31 34 0d 0a 73 |q:=0.998|92114..s|
|00002890| 65 5f 3a 3d 33 2e 33 32 | 39 35 30 36 33 0d 0a 73 |e_:=3.32|95063..s|
|000028a0| 6e 5f 3a 3d 0d 0a 76 31 | 3a 3d 0d 0a 76 32 3a 3d |n_:=..v1|:=..v2:=|
|000028b0| 0d 0a 76 43 72 6f 73 73 | 3a 3d 41 50 50 52 4f 58 |..vCross|:=APPROX|
|000028c0| 28 34 2e 30 31 39 39 39 | 39 39 29 0d 0a 76 61 72 |(4.01999|99)..var|
|000028d0| 73 3a 3d 5b 70 2c 20 79 | 5d 0d 0a 78 30 3a 3d 0d |s:=[p, y|]..x0:=.|
|000028e0| 0a 78 32 3a 3d 0d 0a 78 | 70 78 69 3a 3d 5b 39 2e |.x2:=..x|pxi:=[9.|
|000028f0| 37 33 31 36 38 39 39 b7 | 31 30 5e 28 2d 35 29 2c |7316899.|10^(-5),|
|00002900| 20 2d 20 39 2e 31 33 37 | 36 39 31 35 b7 31 30 5e | - 9.137|6915.10^|
|00002910| 28 2d 38 29 2c 20 2d 30 | 2e 30 33 35 39 30 36 35 |(-8), -0|.0359065|
|00002920| 37 34 3b 20 2d 20 39 2e | 31 33 37 36 39 31 35 b7 |74; - 9.|1376915.|
|00002930| 31 30 5e 28 2d 38 29 2c | 20 37 2e 39 39 36 36 35 |10^(-8),| 7.99665|
|00002940| 34 34 b7 31 30 5e 28 2d | 38 29 2c 20 2d 30 2e 30 |44.10^(-|8), -0.0|
|00002950| 30 35 30 30 32 34 32 38 | 39 3b 20 2d 30 2e 30 33 |05002428|9; -0.03|
|00002960| 35 39 30 36 35 37 34 2c | 20 2d 30 2e 30 30 35 30 |5906574,| -0.0050|
|00002970| 30 32 34 32 38 39 2c 20 | 33 33 31 2e 37 36 35 30 |024289, |331.7650|
|00002980| 37 5d 0d 0a a9 3a 3d 0d | 0a b5 3a 3d 0d 0a df 30 |7]...:=.|..:=...0|
|00002990| 5f 3a 3d 0d 0a df 5f 3a | 3d 5b 30 2e 30 39 37 31 |_:=..._:|=[0.0971|
|000029a0| 38 32 31 35 32 20 2d 20 | 32 2e 34 2c 20 30 2e 30 |82152 - |2.4, 0.0|
|000029b0| 31 32 36 32 38 31 37 33 | 2c 20 34 30 37 2e 31 30 |12628173|, 407.10|
|000029c0| 33 39 32 5d 0d 0a 1f 00 | ff ff 00 00 08 00 43 54 |392]....|......CT|
|000029d0| 65 78 74 4f 62 6a 08 00 | 00 00 0c 00 00 00 cd 03 |extObj..|........|
|000029e0| 00 00 6e 00 00 00 00 ff | 30 02 7b 5c 72 74 66 31 |..n.....|0.{\rtf1|
|000029f0| 5c 61 6e 73 69 5c 64 65 | 66 66 30 5c 64 65 66 74 |\ansi\de|ff0\deft|
|00002a00| 61 62 37 32 30 7b 5c 66 | 6f 6e 74 74 62 6c 7b 5c |ab720{\f|onttbl{\|
|00002a10| 66 30 5c 66 73 77 69 73 | 73 20 4d 53 20 53 61 6e |f0\fswis|s MS San|
|00002a20| 73 20 53 65 72 69 66 3b | 7d 7b 5c 66 31 5c 66 64 |s Serif;|}{\f1\fd|
|00002a30| 65 63 6f 72 5c 66 63 68 | 61 72 73 65 74 32 20 53 |ecor\fch|arset2 S|
|00002a40| 79 6d 62 6f 6c 3b 7d 7b | 5c 66 32 5c 66 73 77 69 |ymbol;}{|\f2\fswi|
|00002a50| 73 73 5c 66 70 72 71 32 | 20 53 79 73 74 65 6d 3b |ss\fprq2| System;|
|00002a60| 7d 7b 5c 66 33 5c 66 73 | 77 69 73 73 5c 66 70 72 |}{\f3\fs|wiss\fpr|
|00002a70| 71 32 20 41 72 69 61 6c | 3b 7d 7b 5c 66 34 5c 66 |q2 Arial|;}{\f4\f|
|00002a80| 6d 6f 64 65 72 6e 5c 66 | 63 68 61 72 73 65 74 32 |modern\f|charset2|
|00002a90| 20 44 66 57 35 20 50 72 | 69 6e 74 65 72 3b 7d 7d | DfW5 Pr|inter;}}|
|00002aa0| 0d 0a 7b 5c 63 6f 6c 6f | 72 74 62 6c 5c 72 65 64 |..{\colo|rtbl\red|
|00002ab0| 30 5c 67 72 65 65 6e 30 | 5c 62 6c 75 65 30 3b 5c |0\green0|\blue0;\|
|00002ac0| 72 65 64 32 35 35 5c 67 | 72 65 65 6e 30 5c 62 6c |red255\g|reen0\bl|
|00002ad0| 75 65 30 3b 5c 72 65 64 | 30 5c 67 72 65 65 6e 30 |ue0;\red|0\green0|
|00002ae0| 5c 62 6c 75 65 32 35 35 | 3b 7d 0d 0a 5c 64 65 66 |\blue255|;}..\def|
|00002af0| 6c 61 6e 67 31 30 33 33 | 5c 70 61 72 64 5c 71 63 |lang1033|\pard\qc|
|00002b00| 5c 70 6c 61 69 6e 5c 66 | 33 5c 66 73 32 34 5c 63 |\plain\f|3\fs24\c|
|00002b10| 66 31 20 4e 6f 6e 6c 69 | 6e 65 61 72 20 52 65 67 |f1 Nonli|near Reg|
|00002b20| 72 65 73 73 69 6f 6e 0d | 0a 5c 70 61 72 20 5c 70 |ression.|.\par \p|
|00002b30| 6c 61 69 6e 5c 66 33 5c | 66 73 32 30 5c 63 66 31 |lain\f3\|fs20\cf1|
|00002b40| 20 47 61 75 73 73 2d 4e | 65 77 74 6f 6e 20 28 57 | Gauss-N|ewton (W|
|00002b50| 69 74 68 20 53 74 65 70 | 2d 48 61 6c 76 69 6e 67 |ith Step|-Halving|
|00002b60| 29 20 26 20 4d 61 72 71 | 75 61 72 64 74 20 4d 65 |) & Marq|uardt Me|
|00002b70| 74 68 6f 64 73 0d 0a 5c | 70 61 72 20 5c 70 61 72 |thods..\|par \par|
|00002b80| 64 5c 71 72 5c 70 6c 61 | 69 6e 5c 66 33 5c 66 73 |d\qr\pla|in\f3\fs|
|00002b90| 32 30 5c 63 66 31 20 4d | 61 63 44 6f 6e 61 6c 64 |20\cf1 M|acDonald|
|00002ba0| 20 52 2e 20 50 68 69 6c | 6c 69 70 73 20 20 20 20 | R. Phil|lips    |
|00002bb0| 20 0d 0a 5c 70 61 72 20 | 50 68 69 6c 6c 69 70 73 | ..\par |Phillips|
|00002bc0| 4d 40 67 61 6f 2e 67 6f | 76 0d 0a 5c 70 61 72 20 |M@gao.go|v..\par |
|00002bd0| 64 6f 6e 70 68 69 6c 6c | 69 70 73 40 73 74 61 72 |donphill|ips@star|
|00002be0| 70 6f 77 65 72 2e 6e 65 | 74 0d 0a 5c 70 61 72 20 |power.ne|t..\par |
|00002bf0| 4e 6f 76 65 6d 62 65 72 | 20 32 30 30 31 5c 70 6c |November| 2001\pl|
|00002c00| 61 69 6e 5c 66 33 5c 66 | 73 32 30 5c 63 66 30 20 |ain\f3\f|s20\cf0 |
|00002c10| 0d 0a 5c 70 61 72 20 7d | 0d 0a ff ff 00 00 08 00 |..\par }|........|
|00002c20| 43 45 78 70 6e 4f 62 6a | 38 00 00 00 7a 00 00 00 |CExpnObj|8...z...|
|00002c30| c0 03 00 00 86 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00002c40| f0 bf 01 00 00 00 01 00 | 00 00 63 5b 49 6e 70 75 |........|..c[Inpu|
|00002c50| 74 4d 6f 64 65 3a 3d 57 | 6f 72 64 2c 50 72 65 63 |tMode:=W|ord,Prec|
|00002c60| 69 73 69 6f 6e 3a 3d 41 | 70 70 72 6f 78 69 6d 61 |ision:=A|pproxima|
|00002c70| 74 65 2c 50 72 65 63 69 | 73 69 6f 6e 44 69 67 69 |te,Preci|sionDigi|
|00002c80| 74 73 3a 3d 31 30 2c 4e | 6f 74 61 74 69 6f 6e 3a |ts:=10,N|otation:|
|00002c90| 3d 53 63 69 65 6e 74 69 | 66 69 63 2c 4e 6f 74 61 |=Scienti|fic,Nota|
|00002ca0| 74 69 6f 6e 44 69 67 69 | 74 73 3a 3d 38 5d 03 80 |tionDigi|ts:=8]..|
|00002cb0| 38 00 00 00 92 00 00 00 | c8 03 00 00 c2 00 00 00 |8.......|........|
|00002cc0| 00 00 00 00 00 00 00 00 | f0 bf 02 00 00 00 01 00 |........|........|
|00002cd0| 00 00 73 5b 6f 75 74 70 | 75 74 3a 3d 2c 69 74 65 |..s[outp|ut:=,ite|
|00002ce0| 72 3a 3d 2c 6d 73 67 3a | 3d 22 22 2c 66 65 71 3a |r:=,msg:|="",feq:|
|00002cf0| 3d 2c 72 73 71 3a 3d 2c | 61 72 73 71 3a 3d 2c 61 |=,rsq:=,|arsq:=,a|
|00002d00| 6e 6f 76 61 3a 3d 2c 62 | 65 74 61 5f 3a 3d 2c 78 |nova:=,b|eta_:=,x|
|00002d10| 70 78 69 3a 3d 2c 76 61 | 72 73 3a 3d 2c 70 61 72 |pxi:=,va|rs:=,par|
|00002d20| 6d 31 3a 3d 2c 69 76 61 | 72 73 31 3a 3d 2c 65 71 |m1:=,iva|rs1:=,eq|
|00002d30| 71 3a 3d 2c 73 65 5f 3a | 3d 2c 64 66 65 5f 3a 3d |q:=,se_:|=,dfe_:=|
|00002d40| 2c 72 5f 3a 3d 5d 01 80 | 08 00 00 00 ce 00 00 00 |,r_:=]..|........|
|00002d50| cd 03 00 00 be 02 00 00 | 00 ff 29 0f 7b 5c 72 74 |........|..).{\rt|
|00002d60| 66 31 5c 61 6e 73 69 5c | 64 65 66 66 30 5c 64 65 |f1\ansi\|deff0\de|
|00002d70| 66 74 61 62 37 32 30 7b | 5c 66 6f 6e 74 74 62 6c |ftab720{|\fonttbl|
|00002d80| 7b 5c 66 30 5c 66 73 77 | 69 73 73 20 4d 53 20 53 |{\f0\fsw|iss MS S|
|00002d90| 61 6e 73 20 53 65 72 69 | 66 3b 7d 7b 5c 66 31 5c |ans Seri|f;}{\f1\|
|00002da0| 66 64 65 63 6f 72 5c 66 | 63 68 61 72 73 65 74 32 |fdecor\f|charset2|
|00002db0| 20 53 79 6d 62 6f 6c 3b | 7d 7b 5c 66 32 5c 66 73 | Symbol;|}{\f2\fs|
|00002dc0| 77 69 73 73 5c 66 70 72 | 71 32 20 53 79 73 74 65 |wiss\fpr|q2 Syste|
|00002dd0| 6d 3b 7d 7b 5c 66 33 5c | 66 73 77 69 73 73 5c 66 |m;}{\f3\|fswiss\f|
|00002de0| 70 72 71 32 20 41 72 69 | 61 6c 3b 7d 7b 5c 66 34 |prq2 Ari|al;}{\f4|
|00002df0| 5c 66 6d 6f 64 65 72 6e | 5c 66 63 68 61 72 73 65 |\fmodern|\fcharse|
|00002e00| 74 32 20 44 66 57 35 20 | 50 72 69 6e 74 65 72 3b |t2 DfW5 |Printer;|
|00002e10| 7d 7d 0d 0a 7b 5c 63 6f | 6c 6f 72 74 62 6c 5c 72 |}}..{\co|lortbl\r|
|00002e20| 65 64 30 5c 67 72 65 65 | 6e 30 5c 62 6c 75 65 30 |ed0\gree|n0\blue0|
|00002e30| 3b 5c 72 65 64 32 35 35 | 5c 67 72 65 65 6e 30 5c |;\red255|\green0\|
|00002e40| 62 6c 75 65 30 3b 5c 72 | 65 64 30 5c 67 72 65 65 |blue0;\r|ed0\gree|
|00002e50| 6e 30 5c 62 6c 75 65 32 | 35 35 3b 7d 0d 0a 5c 64 |n0\blue2|55;}..\d|
|00002e60| 65 66 6c 61 6e 67 31 30 | 33 33 5c 70 61 72 64 5c |eflang10|33\pard\|
|00002e70| 70 6c 61 69 6e 5c 66 33 | 5c 66 73 32 30 5c 63 66 |plain\f3|\fs20\cf|
|00002e80| 32 20 47 61 75 73 73 5f | 4e 65 77 74 6f 6e 5c 70 |2 Gauss_|Newton\p|
|00002e90| 6c 61 69 6e 5c 66 33 5c | 66 73 32 30 5c 63 66 30 |lain\f3\|fs20\cf0|
|00002ea0| 20 20 61 6e 64 20 5c 70 | 6c 61 69 6e 5c 66 33 5c |  and \p|lain\f3\|
|00002eb0| 66 73 32 30 5c 63 66 32 | 20 4d 61 72 71 75 61 72 |fs20\cf2| Marquar|
|00002ec0| 64 74 5c 70 6c 61 69 6e | 5c 66 33 5c 66 73 32 30 |dt\plain|\f3\fs20|
|00002ed0| 5c 63 66 30 20 2c 20 74 | 77 6f 20 6d 65 74 68 6f |\cf0 , t|wo metho|
|00002ee0| 64 73 20 6f 66 20 73 6f | 6c 76 69 6e 67 20 6e 6f |ds of so|lving no|
|00002ef0| 6e 6c 69 6e 65 61 72 20 | 72 65 67 72 65 73 73 69 |nlinear |regressi|
|00002f00| 6f 6e 20 70 72 6f 62 6c | 65 6d 73 2c 20 61 72 65 |on probl|ems, are|
|00002f10| 20 6f 66 66 65 72 65 64 | 20 27 61 73 20 69 73 2e | offered| 'as is.|
|00002f20| 27 20 20 49 20 6d 61 6b | 65 20 6e 6f 20 63 6c 61 |'  I mak|e no cla|
|00002f30| 69 6d 20 74 68 61 74 20 | 74 68 65 79 20 61 72 65 |im that |they are|
|00002f40| 20 62 75 67 20 66 72 65 | 65 2c 20 61 6c 74 68 6f | bug fre|e, altho|
|00002f50| 75 67 68 20 49 20 62 65 | 6c 69 65 76 65 20 74 68 |ugh I be|lieve th|
|00002f60| 65 79 20 61 72 65 2e 20 | 20 49 66 2c 20 62 79 20 |ey are. | If, by |
|00002f70| 63 68 61 6e 63 65 2c 20 | 79 6f 75 20 64 6f 20 66 |chance, |you do f|
|00002f80| 69 6e 64 20 61 6e 79 20 | 62 75 67 73 2c 20 70 6c |ind any |bugs, pl|
|00002f90| 65 61 73 65 20 6c 65 74 | 20 6d 65 20 6b 6e 6f 77 |ease let| me know|
|00002fa0| 20 73 6f 20 74 68 61 74 | 20 49 20 63 61 6e 20 66 | so that| I can f|
|00002fb0| 69 78 20 74 68 65 6d 2e | 0d 0a 5c 70 61 72 20 0d |ix them.|..\par .|
|00002fc0| 0a 5c 70 61 72 20 54 68 | 65 20 70 72 6f 67 72 61 |.\par Th|e progra|
|00002fd0| 6d 73 20 66 69 74 20 64 | 61 74 61 20 74 6f 20 61 |ms fit d|ata to a|
|00002fe0| 6e 20 61 72 62 69 74 72 | 61 72 79 20 66 75 6e 63 |n arbitr|ary func|
|00002ff0| 74 69 6f 6e 2c 20 69 2e | 65 2e 2c 20 74 6f 20 66 |tion, i.|e., to f|
|00003000| 75 6e 63 74 69 6f 6e 73 | 20 74 68 61 74 20 61 72 |unctions| that ar|
|00003010| 65 20 6e 6f 74 20 6e 65 | 63 65 73 73 61 72 69 6c |e not ne|cessaril|
|00003020| 79 20 6c 69 6e 65 61 72 | 20 69 6e 20 74 68 65 69 |y linear| in thei|
|00003030| 72 20 70 61 72 61 6d 65 | 74 65 72 73 2e 20 20 46 |r parame|ters.  F|
|00003040| 69 74 74 69 6e 67 20 6e | 6f 6e 6c 69 6e 65 61 72 |itting n|onlinear|
|00003050| 20 66 75 6e 63 74 69 6f | 6e 73 20 74 6f 20 64 61 | functio|ns to da|
|00003060| 74 61 20 69 73 20 6d 6f | 72 65 20 6f 66 20 61 6e |ta is mo|re of an|
|00003070| 20 61 72 74 20 74 68 61 | 6e 20 61 20 73 63 69 65 | art tha|n a scie|
|00003080| 6e 63 65 2e 20 20 43 6f | 6e 76 65 72 67 65 6e 63 |nce.  Co|nvergenc|
|00003090| 65 20 74 6f 77 61 72 64 | 73 20 61 20 73 6f 6c 75 |e toward|s a solu|
|000030a0| 74 69 6f 6e 20 6d 61 79 | 20 64 65 70 65 6e 64 20 |tion may| depend |
|000030b0| 6f 6e 20 74 68 65 20 6d | 65 74 68 6f 64 20 75 73 |on the m|ethod us|
|000030c0| 65 64 20 61 6e 64 20 63 | 68 6f 69 63 65 20 6f 66 |ed and c|hoice of|
|000030d0| 20 73 74 61 72 74 69 6e | 67 20 76 61 6c 75 65 73 | startin|g values|
|000030e0| 2e 20 20 41 6e 64 20 74 | 68 65 72 65 20 6d 61 79 |.  And t|here may|
|000030f0| 20 62 65 20 6d 6f 72 65 | 20 74 68 61 6e 20 6f 6e | be more| than on|
|00003100| 65 20 73 6f 6c 75 74 69 | 6f 6e 20 6f 72 20 6c 6f |e soluti|on or lo|
|00003110| 63 61 6c 20 6d 69 6e 69 | 6d 75 6d 20 77 69 74 68 |cal mini|mum with|
|00003120| 69 6e 20 61 20 72 65 61 | 73 6f 6e 61 62 6c 65 20 |in a rea|sonable |
|00003130| 72 61 6e 67 65 20 6f 66 | 20 70 61 72 61 6d 65 74 |range of| paramet|
|00003140| 65 72 20 76 61 6c 75 65 | 73 2e 20 20 49 6e 20 61 |er value|s.  In a|
|00003150| 64 64 69 74 69 6f 6e 2c | 20 74 68 65 20 70 72 65 |ddition,| the pre|
|00003160| 63 69 73 69 6f 6e 20 73 | 65 74 20 66 6f 72 20 44 |cision s|et for D|
|00003170| 65 72 69 76 65 20 6d 61 | 79 20 69 6e 66 6c 75 65 |erive ma|y influe|
|00003180| 6e 63 65 20 74 68 65 20 | 70 61 74 68 20 74 6f 77 |nce the |path tow|
|00003190| 61 72 64 73 20 63 6f 6e | 76 65 72 67 65 6e 63 65 |ards con|vergence|
|000031a0| 2e 20 20 54 68 65 72 65 | 20 69 73 20 6e 6f 20 67 |.  There| is no g|
|000031b0| 75 61 72 61 6e 74 65 65 | 20 74 68 61 74 20 74 68 |uarantee| that th|
|000031c0| 65 73 65 20 72 6f 75 74 | 69 6e 65 73 20 77 69 6c |ese rout|ines wil|
|000031d0| 6c 20 62 65 20 61 62 6c | 65 20 74 6f 20 73 6f 6c |l be abl|e to sol|
|000031e0| 76 65 20 61 6e 79 20 72 | 65 67 72 65 73 73 69 6f |ve any r|egressio|
|000031f0| 6e 20 70 72 6f 62 6c 65 | 6d 20 79 6f 75 20 6d 61 |n proble|m you ma|
|00003200| 79 20 74 68 72 6f 77 20 | 61 74 20 74 68 65 6d 2e |y throw |at them.|
|00003210| 0d 0a 5c 70 61 72 20 0d | 0a 5c 70 61 72 20 42 6f |..\par .|.\par Bo|
|00003220| 74 68 20 6d 65 74 68 6f | 64 73 20 72 65 71 75 69 |th metho|ds requi|
|00003230| 72 65 20 61 20 64 61 74 | 61 20 6d 61 74 72 69 78 |re a dat|a matrix|
|00003240| 20 77 68 65 72 65 20 74 | 68 65 20 66 69 72 73 74 | where t|he first|
|00003250| 20 72 6f 77 20 6f 66 20 | 74 68 65 20 6d 61 74 72 | row of |the matr|
|00003260| 69 78 20 63 6f 6e 74 61 | 69 6e 73 20 74 68 65 20 |ix conta|ins the |
|00003270| 6e 61 6d 65 73 20 6f 66 | 20 74 68 65 20 76 61 72 |names of| the var|
|00003280| 69 61 62 6c 65 73 2e 20 | 20 54 68 65 20 64 61 74 |iables. | The dat|
|00003290| 61 20 6d 61 74 72 69 78 | 20 6d 61 79 20 63 6f 6e |a matrix| may con|
|000032a0| 74 61 69 6e 20 74 68 65 | 20 64 65 70 65 6e 64 65 |tain the| depende|
|000032b0| 6e 74 20 61 6e 64 20 69 | 6e 64 65 70 65 6e 64 61 |nt and i|ndependa|
|000032c0| 6e 74 20 76 61 72 69 61 | 62 6c 65 73 20 69 6e 20 |nt varia|bles in |
|000032d0| 61 6e 79 20 6f 72 64 65 | 72 2e 20 20 4e 6f 74 20 |any orde|r.  Not |
|000032e0| 61 6c 6c 20 6f 66 20 74 | 68 65 20 69 6e 64 65 70 |all of t|he indep|
|000032f0| 65 6e 64 65 6e 74 20 76 | 61 72 69 61 62 6c 65 73 |endent v|ariables|
|00003300| 20 6e 65 65 64 20 62 65 | 20 75 73 65 64 20 77 68 | need be| used wh|
|00003310| 65 6e 20 66 69 74 74 69 | 6e 67 20 64 61 74 61 20 |en fitti|ng data |
|00003320| 74 6f 20 61 6e 20 65 71 | 75 61 74 69 6f 6e 2e 20 |to an eq|uation. |
|00003330| 20 54 68 65 20 72 6f 75 | 74 69 6e 65 73 20 62 75 | The rou|tines bu|
|00003340| 69 6c 64 20 61 20 6e 65 | 77 20 64 61 74 61 20 6d |ild a ne|w data m|
|00003350| 61 74 72 69 78 20 66 72 | 6f 6d 20 74 68 65 20 76 |atrix fr|om the v|
|00003360| 61 72 69 61 62 6c 65 20 | 6e 61 6d 65 73 20 69 6e |ariable |names in|
|00003370| 20 74 68 65 20 65 71 75 | 61 74 69 6f 6e 2e 20 20 | the equ|ation.  |
|00003380| 54 68 69 73 20 65 6e 61 | 62 6c 65 73 20 79 6f 75 |This ena|bles you|
|00003390| 20 74 6f 20 66 69 74 20 | 6d 61 6e 79 20 64 69 66 | to fit |many dif|
|000033a0| 66 65 72 65 6e 74 20 65 | 71 75 61 74 69 6f 6e 73 |ferent e|quations|
|000033b0| 20 66 72 6f 6d 20 74 68 | 65 20 73 61 6d 65 20 64 | from th|e same d|
|000033c0| 61 74 61 73 65 74 20 77 | 69 74 68 6f 75 74 20 68 |ataset w|ithout h|
|000033d0| 61 76 69 6e 67 20 74 6f | 20 65 6e 74 65 72 20 61 |aving to| enter a|
|000033e0| 20 6e 65 77 20 64 61 74 | 61 73 65 74 20 66 6f 72 | new dat|aset for|
|000033f0| 20 65 61 63 68 20 70 72 | 6f 62 6c 65 6d 2e 0d 0a | each pr|oblem...|
|00003400| 5c 70 61 72 20 0d 0a 5c | 70 61 72 20 54 68 65 20 |\par ..\|par The |
|00003410| 72 65 71 75 69 72 65 64 | 20 61 72 67 75 65 6d 65 |required| argueme|
|00003420| 6e 74 73 20 74 6f 20 74 | 68 65 20 66 75 6e 63 74 |nts to t|he funct|
|00003430| 69 6f 6e 73 20 61 72 65 | 20 28 31 29 20 27 65 71 |ions are| (1) 'eq|
|00003440| 27 2c 20 74 68 65 20 65 | 71 75 61 74 69 6f 6e 20 |', the e|quation |
|00003450| 74 6f 20 62 65 20 66 69 | 74 20 69 6e 20 74 68 65 |to be fi|t in the|
|00003460| 20 66 6f 72 6d 20 27 64 | 65 70 20 76 61 72 20 3d | form 'd|ep var =|
|00003470| 20 66 28 69 6e 64 65 70 | 20 76 61 72 73 29 27 2c | f(indep| vars)',|
|00003480| 20 28 32 29 20 27 70 61 | 72 6d 27 2c 20 61 20 76 | (2) 'pa|rm', a v|
|00003490| 65 63 74 6f 72 20 6f 66 | 20 74 68 65 20 70 61 72 |ector of| the par|
|000034a0| 61 6d 65 74 65 72 73 20 | 6f 66 20 74 68 65 20 66 |ameters |of the f|
|000034b0| 75 6e 63 74 69 6f 6e 2c | 20 20 28 33 29 20 27 5c |unction,|  (3) '\|
|000034c0| 70 6c 61 69 6e 5c 66 34 | 5c 66 73 32 30 5c 63 66 |plain\f4|\fs20\cf|
|000034d0| 30 20 5c 27 64 66 5c 70 | 6c 61 69 6e 5c 66 33 5c |0 \'df\p|lain\f3\|
|000034e0| 66 73 32 30 5c 63 66 30 | 20 30 5f 27 2c 20 61 20 |fs20\cf0| 0_', a |
|000034f0| 76 65 63 74 6f 72 20 6f | 66 20 74 68 65 20 73 74 |vector o|f the st|
|00003500| 61 72 74 69 6e 67 20 76 | 61 6c 75 65 73 20 28 69 |arting v|alues (i|
|00003510| 6e 69 74 69 61 6c 20 67 | 75 65 73 73 65 73 29 20 |nitial g|uesses) |
|00003520| 66 6f 72 20 74 68 65 20 | 70 61 72 61 6d 65 74 65 |for the |paramete|
|00003530| 72 73 20 69 6e 20 74 68 | 65 20 6f 72 64 65 72 20 |rs in th|e order |
|00003540| 74 68 65 79 20 61 70 70 | 65 61 72 20 69 6e 20 74 |they app|ear in t|
|00003550| 68 65 20 27 70 61 72 6d | 27 20 76 65 63 74 6f 72 |he 'parm|' vector|
|00003560| 2c 20 61 6e 64 20 28 34 | 29 20 27 64 61 74 61 27 |, and (4|) 'data'|
|00003570| 2c 20 74 68 65 20 6e 61 | 6d 65 20 6f 66 20 74 68 |, the na|me of th|
|00003580| 65 20 64 61 74 61 20 6d | 61 74 72 69 78 2e 20 20 |e data m|atrix.  |
|00003590| 4f 70 74 69 6f 6e 61 6c | 20 61 72 67 75 65 6d 65 |Optional| argueme|
|000035a0| 6e 74 73 20 66 6f 72 20 | 5c 70 6c 61 69 6e 5c 66 |nts for |\plain\f|
|000035b0| 33 5c 66 73 32 30 5c 63 | 66 32 20 47 61 75 73 73 |3\fs20\c|f2 Gauss|
|000035c0| 2d 4e 65 77 74 6f 6e 5c | 70 6c 61 69 6e 5c 66 33 |-Newton\|plain\f3|
|000035d0| 5c 66 73 32 30 5c 63 66 | 30 20 20 61 72 65 20 28 |\fs20\cf|0  are (|
|000035e0| 61 29 20 27 6e 5f 27 2c | 20 74 68 65 20 6d 61 78 |a) 'n_',| the max|
|000035f0| 69 6d 75 6d 20 6e 75 6d | 62 65 72 20 6f 66 20 69 |imum num|ber of i|
|00003600| 74 65 72 61 74 69 6f 6e | 73 20 28 64 65 66 61 75 |teration|s (defau|
|00003610| 6c 74 20 69 73 20 33 30 | 29 2c 20 28 62 29 20 27 |lt is 30|), (b) '|
|00003620| 73 6e 5f 27 2c 20 74 68 | 65 20 6d 61 78 69 6d 75 |sn_', th|e maximu|
|00003630| 6d 20 6e 75 6d 62 65 72 | 20 6f 66 20 73 74 65 70 |m number| of step|
|00003640| 20 68 61 6c 76 69 6e 67 | 73 20 70 65 72 20 69 74 | halving|s per it|
|00003650| 65 72 61 74 69 6f 6e 20 | 28 64 65 66 61 75 6c 74 |eration |(default|
|00003660| 20 69 73 20 31 30 29 2c | 20 61 6e 64 20 28 63 29 | is 10),| and (c)|
|00003670| 20 27 63 63 5f 27 2c 20 | 74 68 65 20 63 6f 6e 76 | 'cc_', |the conv|
|00003680| 65 72 67 65 6e 63 65 20 | 63 72 69 74 65 72 69 61 |ergence |criteria|
|00003690| 20 28 74 68 65 20 64 65 | 66 61 75 6c 74 20 69 73 | (the de|fault is|
|000036a0| 20 31 30 5e 28 2d 38 29 | 29 2e 20 20 4f 70 74 69 | 10^(-8)|).  Opti|
|000036b0| 6f 6e 61 6c 20 61 72 67 | 75 65 6d 65 6e 74 73 20 |onal arg|uements |
|000036c0| 66 6f 72 20 5c 70 6c 61 | 69 6e 5c 66 33 5c 66 73 |for \pla|in\f3\fs|
|000036d0| 32 30 5c 63 66 32 20 4d | 61 72 71 75 61 72 64 74 |20\cf2 M|arquardt|
|000036e0| 5c 70 6c 61 69 6e 5c 66 | 33 5c 66 73 32 30 5c 63 |\plain\f|3\fs20\c|
|000036f0| 66 30 20 20 61 72 65 20 | 27 6e 5f 27 20 28 64 65 |f0  are |'n_' (de|
|00003700| 66 61 75 6c 74 20 69 73 | 20 35 30 29 20 61 6e 64 |fault is| 50) and|
|00003710| 20 27 63 63 5f 27 20 28 | 64 65 66 61 75 6c 74 20 | 'cc_' (|default |
|00003720| 69 73 20 31 30 5e 28 2d | 38 29 29 2e 0d 0a 5c 70 |is 10^(-|8))...\p|
|00003730| 61 72 20 0d 0a 5c 70 61 | 72 20 43 6f 6e 76 65 72 |ar ..\pa|r Conver|
|00003740| 67 65 6e 63 65 20 69 73 | 20 62 61 73 65 64 20 6f |gence is| based o|
|00003750| 6e 20 74 68 65 20 72 65 | 6c 61 74 69 76 65 20 63 |n the re|lative c|
|00003760| 68 61 6e 67 65 20 69 6e | 20 74 68 65 20 73 75 6d |hange in| the sum|
|00003770| 20 6f 66 20 73 71 75 61 | 72 65 64 20 65 72 72 6f | of squa|red erro|
|00003780| 72 73 20 28 73 73 65 29 | 2e 0d 0a 5c 70 61 72 20 |rs (sse)|...\par |
|00003790| 0d 0a 5c 70 61 72 20 54 | 68 65 20 72 6f 75 74 69 |..\par T|he routi|
|000037a0| 6e 65 73 20 61 75 74 6f | 6d 61 74 69 63 61 6c 6c |nes auto|maticall|
|000037b0| 79 20 64 65 74 65 63 74 | 20 77 68 65 74 68 65 72 |y detect| whether|
|000037c0| 20 74 68 65 20 65 71 75 | 61 74 69 6f 6e 20 74 6f | the equ|ation to|
|000037d0| 20 62 65 20 66 69 74 74 | 65 64 20 68 61 73 20 61 | be fitt|ed has a|
|000037e0| 20 79 2d 69 6e 74 65 72 | 63 65 70 74 20 61 6e 64 | y-inter|cept and|
|000037f0| 20 61 64 6a 75 73 74 73 | 20 74 68 65 20 63 61 6c | adjusts| the cal|
|00003800| 63 75 6c 61 74 69 6f 6e | 20 6f 66 20 74 68 65 20 |culation| of the |
|00003810| 74 6f 74 61 6c 20 61 6e | 64 20 72 65 67 72 65 73 |total an|d regres|
|00003820| 73 69 6f 6e 20 73 75 6d | 20 6f 66 20 73 71 75 61 |sion sum| of squa|
|00003830| 72 65 73 20 61 63 63 6f | 72 64 69 6e 67 6c 79 2e |res acco|rdingly.|
|00003840| 0d 0a 5c 70 61 72 20 0d | 0a 5c 70 61 72 20 54 68 |..\par .|.\par Th|
|00003850| 65 20 6f 75 74 70 75 74 | 20 69 73 20 61 20 6d 61 |e output| is a ma|
|00003860| 74 72 69 78 20 64 69 73 | 70 6c 61 79 69 6e 67 20 |trix dis|playing |
|00003870| 61 20 6d 65 73 73 61 67 | 65 20 72 65 67 61 72 64 |a messag|e regard|
|00003880| 69 6e 67 20 74 68 65 20 | 63 6f 6e 76 65 72 67 65 |ing the |converge|
|00003890| 6e 63 65 20 6f 72 20 6c | 61 63 6b 20 74 68 65 72 |nce or l|ack ther|
|000038a0| 65 6f 66 2c 20 61 20 73 | 75 62 6d 61 74 72 69 78 |eof, a s|ubmatrix|
|000038b0| 20 6f 66 20 74 68 65 20 | 70 61 72 61 6d 65 74 65 | of the |paramete|
|000038c0| 72 20 76 61 6c 75 65 73 | 2c 20 74 68 65 69 72 20 |r values|, their |
|000038d0| 73 74 61 6e 64 61 72 64 | 20 65 72 72 6f 72 73 2c |standard| errors,|
|000038e0| 20 74 20 76 61 6c 75 65 | 73 2c 20 61 6e 64 20 70 | t value|s, and p|
|000038f0| 72 6f 62 61 62 69 6c 69 | 74 69 65 73 2c 20 61 20 |robabili|ties, a |
|00003900| 73 75 62 6d 61 74 72 69 | 78 20 73 68 6f 77 69 6e |submatri|x showin|
|00003910| 67 20 74 68 65 20 41 4e | 4f 56 41 20 74 61 62 6c |g the AN|OVA tabl|
|00003920| 65 2c 20 61 20 76 65 63 | 74 6f 72 20 6f 66 20 74 |e, a vec|tor of t|
|00003930| 68 65 20 73 74 61 6e 64 | 61 72 64 20 65 72 72 6f |he stand|ard erro|
|00003940| 72 20 6f 66 20 72 65 67 | 72 65 73 73 69 6f 6e 2c |r of reg|ression,|
|00003950| 20 52 5e 32 20 61 6e 64 | 20 61 64 6a 75 73 74 65 | R^2 and| adjuste|
|00003960| 64 20 52 5e 32 20 76 61 | 6c 75 65 73 2c 20 61 6e |d R^2 va|lues, an|
|00003970| 64 20 66 69 6e 61 6c 6c | 79 20 74 68 65 20 66 69 |d finall|y the fi|
|00003980| 74 74 65 64 20 65 71 75 | 61 74 69 6f 6e 2e 20 20 |tted equ|ation.  |
|00003990| 41 75 74 68 6f 72 69 6e | 67 20 27 6d 73 67 27 2c |Authorin|g 'msg',|
|000039a0| 20 27 6f 75 74 70 75 74 | 27 2c 20 27 61 6e 6f 76 | 'output|', 'anov|
|000039b0| 61 27 2c 20 27 72 73 71 | 27 2c 20 27 61 72 73 71 |a', 'rsq|', 'arsq|
|000039c0| 27 2c 20 6f 72 20 27 66 | 65 71 27 20 77 69 6c 6c |', or 'f|eq' will|
|000039d0| 20 72 65 64 69 73 70 6c | 61 79 20 74 68 65 20 69 | redispl|ay the i|
|000039e0| 6e 64 69 76 69 64 75 61 | 6c 20 63 6f 6d 70 6f 6e |ndividua|l compon|
|000039f0| 65 6e 74 73 20 6f 66 20 | 74 68 65 20 6f 75 74 70 |ents of |the outp|
|00003a00| 75 74 2e 20 20 41 75 74 | 68 6f 72 69 6e 67 20 27 |ut.  Aut|horing '|
|00003a10| 69 74 65 72 27 20 77 69 | 6c 6c 20 64 69 73 70 6c |iter' wi|ll displ|
|00003a20| 61 79 20 61 20 6d 61 74 | 72 69 78 20 6f 66 20 74 |ay a mat|rix of t|
|00003a30| 68 65 20 69 74 65 72 61 | 74 69 6f 6e 73 20 74 68 |he itera|tions th|
|00003a40| 65 20 70 72 6f 67 72 61 | 6d 20 77 65 6e 74 20 74 |e progra|m went t|
|00003a50| 68 72 6f 75 67 68 20 74 | 6f 20 61 72 72 69 76 65 |hrough t|o arrive|
|00003a60| 20 61 74 20 74 68 65 20 | 72 65 73 75 6c 74 73 2e | at the |results.|
|00003a70| 0d 0a 5c 70 61 72 20 0d | 0a 5c 70 61 72 20 41 73 |..\par .|.\par As|
|00003a80| 20 74 6f 20 77 68 69 63 | 68 20 6d 65 74 68 6f 64 | to whic|h method|
|00003a90| 20 69 73 20 62 65 73 74 | 20 74 6f 20 75 73 65 20 | is best| to use |
|00003aa0| 66 6f 72 20 61 20 70 72 | 6f 62 6c 65 6d 2c 20 69 |for a pr|oblem, i|
|00003ab0| 74 20 69 73 20 73 61 69 | 64 20 74 68 61 74 20 74 |t is sai|d that t|
|00003ac0| 68 65 20 4d 61 72 71 75 | 61 72 64 74 20 6d 65 74 |he Marqu|ardt met|
|00003ad0| 68 6f 64 20 69 73 20 6d | 6f 72 65 20 66 6f 72 67 |hod is m|ore forg|
|00003ae0| 69 76 69 6e 67 20 6f 66 | 20 61 20 70 6f 6f 72 20 |iving of| a poor |
|00003af0| 63 68 6f 69 63 65 20 6f | 66 20 73 74 61 72 74 69 |choice o|f starti|
|00003b00| 6e 67 20 76 61 6c 75 65 | 73 20 61 6c 74 68 6f 75 |ng value|s althou|
|00003b10| 67 68 20 49 20 68 61 76 | 65 20 66 6f 75 6e 64 20 |gh I hav|e found |
|00003b20| 74 68 69 73 20 69 73 20 | 6e 6f 74 20 61 6c 77 61 |this is |not alwa|
|00003b30| 79 73 20 74 68 65 20 63 | 61 73 65 2e 20 20 54 68 |ys the c|ase.  Th|
|00003b40| 69 73 20 64 6f 65 73 20 | 6e 6f 74 20 6d 65 61 6e |is does |not mean|
|00003b50| 20 74 68 61 74 20 61 6e | 79 20 73 65 74 20 6f 66 | that an|y set of|
|00003b60| 20 73 74 61 72 74 69 6e | 67 20 76 61 6c 75 65 73 | startin|g values|
|00003b70| 20 6d 61 79 20 62 65 20 | 75 73 65 64 2e 0d 0a 5c | may be |used...\|
|00003b80| 70 61 72 20 0d 0a 5c 70 | 61 72 20 28 54 68 65 73 |par ..\p|ar (Thes|
|00003b90| 65 20 72 6f 75 74 69 6e | 65 73 20 6d 61 79 2c 20 |e routin|es may, |
|00003ba0| 6f 66 20 63 6f 75 72 73 | 65 2c 20 62 65 20 75 73 |of cours|e, be us|
|00003bb0| 65 64 20 74 6f 20 66 69 | 74 20 61 6e 79 20 74 79 |ed to fi|t any ty|
|00003bc0| 70 65 20 6f 66 20 6c 69 | 6e 65 61 72 20 72 65 67 |pe of li|near reg|
|00003bd0| 72 65 73 73 69 6f 6e 2e | 20 20 46 6f 72 20 74 68 |ression.|  For th|
|00003be0| 69 73 20 74 61 73 6b 2c | 20 5c 70 6c 61 69 6e 5c |is task,| \plain\|
|00003bf0| 66 33 5c 66 73 32 30 5c | 63 66 32 20 47 61 75 73 |f3\fs20\|cf2 Gaus|
|00003c00| 73 5f 4e 65 77 74 6f 6e | 5c 70 6c 61 69 6e 5c 66 |s_Newton|\plain\f|
|00003c10| 33 5c 66 73 32 30 5c 63 | 66 30 20 20 73 65 65 6d |3\fs20\c|f0  seem|
|00003c20| 73 20 74 6f 20 62 65 20 | 74 68 65 20 62 65 74 74 |s to be |the bett|
|00003c30| 65 72 20 6f 70 74 69 6f | 6e 2e 20 20 53 69 6d 70 |er optio|n.  Simp|
|00003c40| 6c 79 20 73 65 74 20 74 | 68 65 20 69 6e 69 74 69 |ly set t|he initi|
|00003c50| 61 6c 20 76 61 6c 75 65 | 73 20 6f 66 20 74 68 65 |al value|s of the|
|00003c60| 20 70 61 72 61 6d 65 74 | 65 72 20 65 73 74 69 6d | paramet|er estim|
|00003c70| 61 74 65 73 20 74 6f 20 | 31 2e 29 0d 0a 5c 70 61 |ates to |1.)..\pa|
|00003c80| 72 20 7d 0d 0a 03 80 38 | 00 00 00 ca 02 00 00 c8 |r }....8|........|
|00003c90| 00 00 00 42 03 00 00 00 | 00 00 00 00 00 00 00 f0 |...B....|........|
|00003ca0| bf 03 00 00 00 01 00 00 | 00 5d 70 6f 73 69 74 69 |........|.]positi|
|00003cb0| 6f 6e 28 65 2c 75 29 3a | 3d 50 52 4f 47 28 64 3a |on(e,u):|=PROG(d:|
|00003cc0| 3d 44 49 4d 28 75 29 2c | 69 5f 3a 3d 31 2c 4c 4f |=DIM(u),|i_:=1,LO|
|00003cd0| 4f 50 28 49 46 28 65 3d | 75 20 53 55 42 20 69 5f |OP(IF(e=|u SUB i_|
|00003ce0| 2c 52 45 54 55 52 4e 28 | 69 5f 29 29 2c 49 46 28 |,RETURN(|i_)),IF(|
|00003cf0| 69 5f 3e 64 2c 52 45 54 | 55 52 4e 28 30 29 29 2c |i_>d,RET|URN(0)),|
|00003d00| 69 5f 3a 2b 31 29 29 03 | 80 38 00 00 00 4e 03 00 |i_:+1)).|.8...N..|
|00003d10| 00 40 03 00 00 3e 07 00 | 00 00 04 55 73 65 72 00 |.@...>..|...User.|
|00003d20| 00 00 00 00 00 f0 bf 04 | 00 00 00 01 00 00 00 ff |........|........|
|00003d30| aa 09 47 41 55 53 53 5f | 4e 45 57 54 4f 4e 28 65 |..GAUSS_|NEWTON(e|
|00003d40| 71 2c 70 61 72 6d 2c 62 | 65 74 61 30 5f 2c 64 61 |q,parm,b|eta0_,da|
|00003d50| 74 61 2c 6e 5f 3a 3d 33 | 30 2c 73 6e 5f 3a 3d 31 |ta,n_:=3|0,sn_:=1|
|00003d60| 30 2c 63 63 5f 3a 3d 31 | 30 5e 28 2d 38 29 29 3a |0,cc_:=1|0^(-8)):|
|00003d70| 3d 50 52 4f 47 28 65 71 | 71 3a 3d 52 48 53 28 65 |=PROG(eq|q:=RHS(e|
|00003d80| 71 29 2c 64 76 61 72 3a | 3d 4c 48 53 28 65 71 29 |q),dvar:|=LHS(eq)|
|00003d90| 2c 69 76 61 72 73 32 3a | 3d 56 41 52 49 41 42 4c |,ivars2:|=VARIABL|
|00003da0| 45 53 28 65 71 71 29 2c | 69 76 61 72 73 3a 3d 53 |ES(eqq),|ivars:=S|
|00003db0| 45 4c 45 43 54 28 4e 4f | 54 28 4d 45 4d 42 45 52 |ELECT(NO|T(MEMBER|
|00003dc0| 3f 28 69 5f 2c 70 61 72 | 6d 29 29 2c 69 5f 2c 69 |?(i_,par|m)),i_,i|
|00003dd0| 76 61 72 73 32 29 2c 69 | 76 61 72 73 31 3a 3d 69 |vars2),i|vars1:=i|
|00003de0| 76 61 72 73 2c 70 61 72 | 6d 31 3a 3d 70 61 72 6d |vars,par|m1:=parm|
|00003df0| 2c 6f 62 73 5f 3a 3d 44 | 49 4d 28 64 61 74 61 29 |,obs_:=D|IM(data)|
|00003e00| 2d 31 2c 76 61 72 73 3a | 3d 64 61 74 61 20 53 55 |-1,vars:|=data SU|
|00003e10| 42 20 31 2c 64 61 74 61 | 31 3a 3d 64 61 74 61 20 |B 1,data|1:=data |
|00003e20| 53 55 42 20 5b 32 2c 2e | 2e 2e 2c 6f 62 73 5f 2b |SUB [2,.|..,obs_+|
|00003e30| 31 5d 2c 6b 5f 3a 3d 44 | 49 4d 28 70 61 72 6d 29 |1],k_:=D|IM(parm)|
|00003e40| 2c 70 6f 73 5f 3a 3d 70 | 6f 73 69 74 69 6f 6e 28 |,pos_:=p|osition(|
|00003e50| 64 76 61 72 2c 76 61 72 | 73 29 2c 49 46 28 70 6f |dvar,var|s),IF(po|
|00003e60| 73 5f 3d 30 2c 52 45 54 | 55 52 4e 28 22 55 6e 64 |s_=0,RET|URN("Und|
|00003e70| 65 66 69 6e 65 64 20 64 | 65 70 65 6e 64 65 6e 74 |efined d|ependent|
|00003e80| 20 76 61 72 69 61 62 6c | 65 21 22 29 2c 79 5f 3a | variabl|e!"),y_:|
|00003e90| 3d 64 61 74 61 31 20 53 | 55 42 20 20 53 55 42 20 |=data1 S|UB  SUB |
|00003ea0| 70 6f 73 5f 29 2c 70 6f | 73 5f 3a 3d 56 45 43 54 |pos_),po|s_:=VECT|
|00003eb0| 4f 52 28 70 6f 73 69 74 | 69 6f 6e 28 69 76 61 72 |OR(posit|ion(ivar|
|00003ec0| 73 20 53 55 42 20 69 5f | 2c 76 61 72 73 29 2c 69 |s SUB i_|,vars),i|
|00003ed0| 5f 2c 31 2c 44 49 4d 28 | 69 76 61 72 73 29 29 2c |_,1,DIM(|ivars)),|
|00003ee0| 49 46 28 4d 45 4d 42 45 | 52 3f 28 30 2c 70 6f 73 |IF(MEMBE|R?(0,pos|
|00003ef0| 5f 29 2c 52 45 54 55 52 | 4e 28 22 55 6e 64 65 66 |_),RETUR|N("Undef|
|00003f00| 69 6e 65 64 20 69 6e 64 | 65 70 65 6e 64 65 6e 74 |ined ind|ependent|
|00003f10| 20 76 61 72 69 61 62 6c | 65 28 73 29 21 22 29 2c | variabl|e(s)!"),|
|00003f20| 6d 5f 3a 3d 64 61 74 61 | 31 20 53 55 42 20 20 53 |m_:=data|1 SUB  S|
|00003f30| 55 42 20 70 6f 73 5f 29 | 2c 79 63 70 74 3a 3d 49 |UB pos_)|,ycpt:=I|
|00003f40| 46 28 4d 45 4d 42 45 52 | 3f 28 31 2c 47 52 41 44 |F(MEMBER|?(1,GRAD|
|00003f50| 28 65 71 71 2c 70 61 72 | 6d 29 29 2c 31 2c 30 2c |(eqq,par|m)),1,0,|
|00003f60| 30 29 2c 6b 5f 3a 2d 79 | 63 70 74 2c 62 65 74 61 |0),k_:-y|cpt,beta|
|00003f70| 5f 3a 3d 62 65 74 61 30 | 5f 2c 65 71 31 3a 3d 56 |_:=beta0|_,eq1:=V|
|00003f80| 45 43 54 4f 52 28 53 55 | 42 53 54 28 65 71 71 2c |ECTOR(SU|BST(eqq,|
|00003f90| 69 76 61 72 73 2c 6d 5f | 20 53 55 42 20 69 29 2c |ivars,m_| SUB i),|
|00003fa0| 69 2c 6f 62 73 5f 29 2c | 6a 5f 3a 3d 47 52 41 44 |i,obs_),|j_:=GRAD|
|00003fb0| 28 65 71 31 2c 70 61 72 | 6d 29 60 2c 70 5f 3a 3d |(eq1,par|m)`,p_:=|
|00003fc0| 53 55 42 53 54 28 65 71 | 31 2c 70 61 72 6d 2c 62 |SUBST(eq|1,parm,b|
|00003fd0| 65 74 61 5f 29 2c 72 5f | 3a 3d 79 5f 2d 70 5f 2c |eta_),r_|:=y_-p_,|
|00003fe0| 73 73 65 5f 3a 3d 72 5f | 20 2e 20 72 5f 2c 74 5f |sse_:=r_| . r_,t_|
|00003ff0| 3a 3d 49 46 28 79 63 70 | 74 3d 31 2c 79 5f 2d 56 |:=IF(ycp|t=1,y_-V|
|00004000| 45 43 54 4f 52 28 41 56 | 45 52 41 47 45 28 79 5f |ECTOR(AV|ERAGE(y_|
|00004010| 29 2c 69 2c 6f 62 73 5f | 29 2c 79 5f 29 2c 73 73 |),i,obs_|),y_),ss|
|00004020| 74 5f 3a 3d 74 5f 20 2e | 20 74 5f 2c 69 74 65 72 |t_:=t_ .| t_,iter|
|00004030| 3a 3d 5b 41 50 50 45 4e | 44 28 5b 30 5d 2c 62 65 |:=[APPEN|D([0],be|
|00004040| 74 61 5f 2c 5b 73 73 65 | 5f 5d 29 5d 2c 65 70 73 |ta_,[sse|_])],eps|
|00004050| 5f 3a 3d 31 2c 69 74 65 | 72 5f 3a 3d 31 2c 4c 4f |_:=1,ite|r_:=1,LO|
|00004060| 4f 50 28 49 46 28 69 74 | 65 72 5f 3e 6e 5f 20 4f |OP(IF(it|er_>n_ O|
|00004070| 52 20 65 70 73 5f 3c 63 | 63 5f 2c 65 78 69 74 29 |R eps_<c|c_,exit)|
|00004080| 2c 78 5f 3a 3d 53 55 42 | 53 54 28 6a 5f 2c 70 61 |,x_:=SUB|ST(j_,pa|
|00004090| 72 6d 2c 62 65 74 61 5f | 29 2c 6c 73 73 65 5f 3a |rm,beta_|),lsse_:|
|000040a0| 3d 73 73 65 5f 2c 78 70 | 78 69 3a 3d 28 78 5f 60 |=sse_,xp|xi:=(x_`|
|000040b0| 20 2e 20 78 5f 29 5e 28 | 2d 31 29 2c 64 65 6c 74 | . x_)^(|-1),delt|
|000040c0| 61 5f 3a 3d 78 70 78 69 | 20 2e 20 78 5f 60 20 2e |a_:=xpxi| . x_` .|
|000040d0| 20 72 5f 2c 6f 6c 64 5f | 3a 3d 62 65 74 61 5f 2c | r_,old_|:=beta_,|
|000040e0| 62 65 74 61 5f 3a 2b 64 | 65 6c 74 61 5f 2c 70 5f |beta_:+d|elta_,p_|
|000040f0| 3a 3d 53 55 42 53 54 28 | 65 71 31 2c 70 61 72 6d |:=SUBST(|eq1,parm|
|00004100| 2c 62 65 74 61 5f 29 2c | 72 5f 3a 3d 79 5f 2d 70 |,beta_),|r_:=y_-p|
|00004110| 5f 2c 73 73 65 5f 3a 3d | 72 5f 20 2e 20 72 5f 2c |_,sse_:=|r_ . r_,|
|00004120| 69 74 65 72 3a 3d 41 50 | 50 45 4e 44 28 69 74 65 |iter:=AP|PEND(ite|
|00004130| 72 2c 5b 41 50 50 45 4e | 44 28 5b 69 74 65 72 5f |r,[APPEN|D([iter_|
|00004140| 5d 2c 62 65 74 61 5f 2c | 5b 73 73 65 5f 5d 29 5d |],beta_,|[sse_])]|
|00004150| 29 2c 73 75 62 69 74 5f | 3a 3d 31 2c 4c 4f 4f 50 |),subit_|:=1,LOOP|
|00004160| 28 49 46 28 73 73 65 5f | 3c 3d 6c 73 73 65 5f 20 |(IF(sse_|<=lsse_ |
|00004170| 4f 52 20 73 75 62 69 74 | 5f 3e 73 6e 5f 2c 65 78 |OR subit|_>sn_,ex|
|00004180| 69 74 29 2c 64 65 6c 74 | 61 5f 3a 2f 32 2c 62 65 |it),delt|a_:/2,be|
|00004190| 74 61 5f 3a 3d 6f 6c 64 | 5f 2b 64 65 6c 74 61 5f |ta_:=old|_+delta_|
|000041a0| 2c 70 5f 3a 3d 53 55 42 | 53 54 28 65 71 31 2c 70 |,p_:=SUB|ST(eq1,p|
|000041b0| 61 72 6d 2c 62 65 74 61 | 5f 29 2c 72 5f 3a 3d 79 |arm,beta|_),r_:=y|
|000041c0| 5f 2d 70 5f 2c 73 73 65 | 5f 3a 3d 72 5f 20 2e 20 |_-p_,sse|_:=r_ . |
|000041d0| 72 5f 2c 69 74 65 72 3a | 3d 41 50 50 45 4e 44 28 |r_,iter:|=APPEND(|
|000041e0| 69 74 65 72 2c 5b 41 50 | 50 45 4e 44 28 5b 69 74 |iter,[AP|PEND([it|
|000041f0| 65 72 5f 2b 73 75 62 69 | 74 5f 2f 31 30 30 5d 2c |er_+subi|t_/100],|
|00004200| 62 65 74 61 5f 2c 5b 73 | 73 65 5f 5d 29 5d 29 2c |beta_,[s|se_])]),|
|00004210| 73 75 62 69 74 5f 3a 2b | 31 29 2c 49 46 28 73 75 |subit_:+|1),IF(su|
|00004220| 62 69 74 5f 3e 73 6e 5f | 2c 65 78 69 74 29 2c 65 |bit_>sn_|,exit),e|
|00004230| 70 73 5f 3a 3d 41 42 53 | 28 28 6c 73 73 65 5f 2d |ps_:=ABS|((lsse_-|
|00004240| 73 73 65 5f 29 2f 28 73 | 73 65 5f 2b 31 30 5e 28 |sse_)/(s|se_+10^(|
|00004250| 2d 36 29 29 29 2c 69 74 | 65 72 5f 3a 2b 31 29 2c |-6))),it|er_:+1),|
|00004260| 49 46 28 69 74 65 72 5f | 3e 6e 5f 2c 6d 73 67 3a |IF(iter_|>n_,msg:|
|00004270| 3d 41 50 50 45 4e 44 28 | 22 43 6f 6e 76 65 72 67 |=APPEND(|"Converg|
|00004280| 65 6e 63 65 20 66 61 69 | 6c 65 64 20 61 66 74 65 |ence fai|led afte|
|00004290| 72 20 22 2c 6e 5f 2c 22 | 20 69 74 65 72 61 74 69 |r ",n_,"| iterati|
|000042a0| 6f 6e 73 21 22 29 2c 49 | 46 28 73 75 62 69 74 5f |ons!"),I|F(subit_|
|000042b0| 3e 73 6e 5f 2c 6d 73 67 | 3a 3d 22 53 53 45 20 64 |>sn_,msg|:="SSE d|
|000042c0| 69 64 20 6e 6f 74 20 69 | 6d 70 72 6f 76 65 20 61 |id not i|mprove a|
|000042d0| 66 74 65 72 20 31 30 20 | 68 61 6c 76 69 6e 67 73 |fter 10 |halvings|
|000042e0| 21 22 2c 6d 73 67 3a 3d | 22 43 6f 6e 76 65 72 67 |!",msg:=|"Converg|
|000042f0| 65 6e 63 65 20 63 72 69 | 74 65 72 69 61 20 6d 65 |ence cri|teria me|
|00004300| 74 21 22 29 29 2c 64 66 | 65 5f 3a 3d 6f 62 73 5f |t!")),df|e_:=obs_|
|00004310| 2d 6b 5f 2d 79 63 70 74 | 2c 6d 73 65 5f 3a 3d 73 |-k_-ycpt|,mse_:=s|
|00004320| 73 65 5f 2f 64 66 65 5f | 2c 73 65 5f 3a 3d 53 51 |se_/dfe_|,se_:=SQ|
|00004330| 52 54 28 6d 73 65 5f 29 | 2c 73 74 64 5f 3a 3d 56 |RT(mse_)|,std_:=V|
|00004340| 45 43 54 4f 52 28 53 51 | 52 54 28 78 70 78 69 20 |ECTOR(SQ|RT(xpxi |
|00004350| 53 55 42 20 69 5f 20 53 | 55 42 20 69 5f 2a 6d 73 |SUB i_ S|UB i_*ms|
|00004360| 65 5f 29 2c 69 5f 2c 6b | 5f 2b 79 63 70 74 29 2c |e_),i_,k|_+ycpt),|
|00004370| 74 5f 3a 3d 56 45 43 54 | 4f 52 28 62 65 74 61 5f |t_:=VECT|OR(beta_|
|00004380| 20 53 55 42 20 69 5f 2f | 73 74 64 5f 20 53 55 42 | SUB i_/|std_ SUB|
|00004390| 20 69 5f 2c 69 5f 2c 6b | 5f 2b 79 63 70 74 29 2c | i_,i_,k|_+ycpt),|
|000043a0| 70 72 6f 62 5f 3a 3d 56 | 45 43 54 4f 52 28 31 2d |prob_:=V|ECTOR(1-|
|000043b0| 53 54 55 44 45 4e 54 28 | 74 5f 20 53 55 42 20 69 |STUDENT(|t_ SUB i|
|000043c0| 5f 2c 64 66 65 5f 29 2c | 69 5f 2c 6b 5f 2b 79 63 |_,dfe_),|i_,k_+yc|
|000043d0| 70 74 29 2c 73 73 72 5f | 3a 3d 73 73 74 5f 2d 73 |pt),ssr_|:=sst_-s|
|000043e0| 73 65 5f 2c 72 73 71 3a | 3d 73 73 72 5f 2f 73 73 |se_,rsq:|=ssr_/ss|
|000043f0| 74 5f 2c 61 72 73 71 3a | 3d 31 2d 73 73 65 5f 2f |t_,arsq:|=1-sse_/|
|00004400| 28 6f 62 73 5f 2d 6b 5f | 2d 79 63 70 74 29 2f 28 |(obs_-k_|-ycpt)/(|
|00004410| 73 73 74 5f 2f 28 6f 62 | 73 5f 2d 79 63 70 74 29 |sst_/(ob|s_-ycpt)|
|00004420| 29 2c 66 64 5f 3a 3d 46 | 5f 44 49 53 54 52 49 42 |),fd_:=F|_DISTRIB|
|00004430| 55 54 49 4f 4e 28 73 73 | 72 5f 2a 64 66 65 5f 2f |UTION(ss|r_*dfe_/|
|00004440| 28 73 73 65 5f 2a 6b 5f | 29 2c 6b 5f 2c 64 66 65 |(sse_*k_|),k_,dfe|
|00004450| 5f 29 2c 61 31 5f 3a 3d | 5b 22 53 6f 75 72 63 65 |_),a1_:=|["Source|
|00004460| 22 2c 22 44 46 22 2c 22 | 53 53 22 2c 22 4d 53 22 |","DF","|SS","MS"|
|00004470| 2c 22 46 22 2c 22 50 72 | 6f 62 28 46 29 22 5d 2c |,"F","Pr|ob(F)"],|
|00004480| 61 32 5f 3a 3d 5b 22 52 | 65 67 22 2c 6b 5f 2c 73 |a2_:=["R|eg",k_,s|
|00004490| 73 72 5f 2c 73 73 72 5f | 2f 6b 5f 2c 73 73 72 5f |sr_,ssr_|/k_,ssr_|
|000044a0| 2a 64 66 65 5f 2f 28 73 | 73 65 5f 2a 6b 5f 29 2c |*dfe_/(s|se_*k_),|
|000044b0| 49 46 28 66 64 5f 3c 30 | 2c 30 2c 66 64 5f 29 5d |IF(fd_<0|,0,fd_)]|
|000044c0| 2c 61 33 5f 3a 3d 5b 22 | 45 72 72 6f 72 22 2c 64 |,a3_:=["|Error",d|
|000044d0| 66 65 5f 2c 73 73 65 5f | 2c 73 73 65 5f 2f 64 66 |fe_,sse_|,sse_/df|
|000044e0| 65 5f 2c 22 20 22 2c 22 | 20 22 5d 2c 61 34 5f 3a |e_," ","| "],a4_:|
|000044f0| 3d 5b 49 46 28 79 63 70 | 74 3d 31 2c 22 43 6f 72 |=[IF(ycp|t=1,"Cor|
|00004500| 72 65 63 74 65 64 20 54 | 6f 74 61 6c 22 2c 22 55 |rected T|otal","U|
|00004510| 6e 63 6f 72 72 65 63 74 | 65 64 20 54 6f 74 61 6c |ncorrect|ed Total|
|00004520| 22 29 2c 64 66 65 5f 2b | 6b 5f 2c 73 73 74 5f 2c |"),dfe_+|k_,sst_,|
|00004530| 22 20 22 2c 22 20 22 2c | 22 20 22 5d 2c 61 6e 6f |" "," ",|" "],ano|
|00004540| 76 61 3a 3d 41 50 50 45 | 4e 44 28 5b 61 31 5f 5d |va:=APPE|ND([a1_]|
|00004550| 2c 5b 61 32 5f 5d 2c 5b | 61 33 5f 5d 2c 5b 61 34 |,[a2_],[|a3_],[a4|
|00004560| 5f 5d 29 2c 66 65 71 3a | 3d 64 76 61 72 3d 53 55 |_]),feq:|=dvar=SU|
|00004570| 42 53 54 28 65 71 71 2c | 70 61 72 6d 2c 62 65 74 |BST(eqq,|parm,bet|
|00004580| 61 5f 29 2c 69 74 65 72 | 3a 3d 41 50 50 45 4e 44 |a_),iter|:=APPEND|
|00004590| 28 5b 41 50 50 45 4e 44 | 28 5b 22 49 74 65 72 22 |([APPEND|(["Iter"|
|000045a0| 5d 2c 70 61 72 6d 2c 5b | 22 53 53 45 22 5d 29 5d |],parm,[|"SSE"])]|
|000045b0| 2c 69 74 65 72 29 2c 74 | 69 74 6c 65 3a 3d 5b 5b |,iter),t|itle:=[[|
|000045c0| 22 50 61 72 6d 22 5d 2c | 5b 22 56 61 6c 75 65 22 |"Parm"],|["Value"|
|000045d0| 5d 2c 5b 22 53 54 44 22 | 5d 2c 5b 41 50 50 45 4e |],["STD"|],[APPEN|
|000045e0| 44 28 22 74 28 22 2c 64 | 66 65 5f 2c 22 29 22 29 |D("t(",d|fe_,")")|
|000045f0| 5d 2c 5b 22 50 72 6f 62 | 28 74 29 22 5d 5d 2c 6f |],["Prob|(t)"]],o|
|00004600| 75 74 70 75 74 3a 3d 41 | 50 50 45 4e 44 5f 43 4f |utput:=A|PPEND_CO|
|00004610| 4c 55 4d 4e 53 28 74 69 | 74 6c 65 2c 41 50 50 45 |LUMNS(ti|tle,APPE|
|00004620| 4e 44 28 5b 70 61 72 6d | 5d 2c 5b 62 65 74 61 5f |ND([parm|],[beta_|
|00004630| 5d 2c 5b 73 74 64 5f 5d | 2c 5b 74 5f 5d 2c 5b 70 |],[std_]|,[t_],[p|
|00004640| 72 6f 62 5f 5d 29 29 60 | 2c 73 74 61 74 73 3a 3d |rob_]))`|,stats:=|
|00004650| 5b 5b 22 53 45 22 2c 22 | 52 5e 32 22 2c 22 41 64 |[["SE","|R^2","Ad|
|00004660| 6a 52 5e 32 22 5d 2c 5b | 73 65 5f 2c 72 73 71 2c |jR^2"],[|se_,rsq,|
|00004670| 61 72 73 71 5d 5d 2c 52 | 45 54 55 52 4e 28 5b 5b |arsq]],R|ETURN([[|
|00004680| 22 47 61 75 73 73 5f 4e | 65 77 74 6f 6e 20 4d 65 |"Gauss_N|ewton Me|
|00004690| 74 68 6f 64 22 5d 2c 5b | 22 20 22 5d 2c 5b 6d 73 |thod"],[|" "],[ms|
|000046a0| 67 5d 2c 5b 22 20 22 5d | 2c 5b 6f 75 74 70 75 74 |g],[" "]|,[output|
|000046b0| 5d 2c 5b 22 20 22 5d 2c | 5b 61 6e 6f 76 61 5d 2c |],[" "],|[anova],|
|000046c0| 5b 22 20 22 5d 2c 5b 73 | 74 61 74 73 5d 2c 5b 22 |[" "],[s|tats],["|
|000046d0| 20 22 5d 2c 5b 66 65 71 | 5d 5d 29 29 03 80 38 00 | "],[feq|]]))..8.|
|000046e0| 00 00 4a 07 00 00 40 03 | 00 00 0a 0b 00 00 00 04 |..J...@.|........|
|000046f0| 55 73 65 72 00 00 00 00 | 00 00 f0 bf 05 00 00 00 |User....|........|
|00004700| 01 00 00 00 ff 62 09 4d | 41 52 51 55 41 52 44 54 |.....b.M|ARQUARDT|
|00004710| 28 65 71 2c 70 61 72 6d | 2c 62 65 74 61 30 5f 2c |(eq,parm|,beta0_,|
|00004720| 64 61 74 61 2c 6e 5f 3a | 3d 35 30 2c 63 63 5f 3a |data,n_:|=50,cc_:|
|00004730| 3d 31 30 5e 28 2d 38 29 | 29 3a 3d 50 52 4f 47 28 |=10^(-8)|):=PROG(|
|00004740| 65 71 71 3a 3d 52 48 53 | 28 65 71 29 2c 64 76 61 |eqq:=RHS|(eq),dva|
|00004750| 72 3a 3d 4c 48 53 28 65 | 71 29 2c 69 76 61 72 73 |r:=LHS(e|q),ivars|
|00004760| 32 3a 3d 56 41 52 49 41 | 42 4c 45 53 28 65 71 71 |2:=VARIA|BLES(eqq|
|00004770| 29 2c 69 76 61 72 73 3a | 3d 53 45 4c 45 43 54 28 |),ivars:|=SELECT(|
|00004780| 4e 4f 54 28 4d 45 4d 42 | 45 52 3f 28 69 5f 2c 70 |NOT(MEMB|ER?(i_,p|
|00004790| 61 72 6d 29 29 2c 69 5f | 2c 69 76 61 72 73 32 29 |arm)),i_|,ivars2)|
|000047a0| 2c 70 61 72 6d 31 3a 3d | 70 61 72 6d 2c 69 76 61 |,parm1:=|parm,iva|
|000047b0| 72 73 31 3a 3d 69 76 61 | 72 73 2c 6f 62 73 5f 3a |rs1:=iva|rs,obs_:|
|000047c0| 3d 44 49 4d 28 64 61 74 | 61 29 2d 31 2c 76 61 72 |=DIM(dat|a)-1,var|
|000047d0| 73 3a 3d 64 61 74 61 20 | 53 55 42 20 31 2c 64 61 |s:=data |SUB 1,da|
|000047e0| 74 61 31 3a 3d 64 61 74 | 61 20 53 55 42 20 5b 32 |ta1:=dat|a SUB [2|
|000047f0| 2c 2e 2e 2e 2c 6f 62 73 | 5f 2b 31 5d 2c 6b 5f 3a |,...,obs|_+1],k_:|
|00004800| 3d 44 49 4d 28 70 61 72 | 6d 29 2c 70 6f 73 5f 3a |=DIM(par|m),pos_:|
|00004810| 3d 70 6f 73 69 74 69 6f | 6e 28 64 76 61 72 2c 76 |=positio|n(dvar,v|
|00004820| 61 72 73 29 2c 49 46 28 | 70 6f 73 5f 3d 30 2c 52 |ars),IF(|pos_=0,R|
|00004830| 45 54 55 52 4e 28 22 55 | 6e 64 65 66 69 6e 65 64 |ETURN("U|ndefined|
|00004840| 20 64 65 70 65 6e 64 65 | 6e 74 20 76 61 72 69 61 | depende|nt varia|
|00004850| 62 6c 65 21 22 29 2c 79 | 5f 3a 3d 64 61 74 61 31 |ble!"),y|_:=data1|
|00004860| 20 53 55 42 20 20 53 55 | 42 20 70 6f 73 5f 29 2c | SUB  SU|B pos_),|
|00004870| 70 6f 73 5f 3a 3d 56 45 | 43 54 4f 52 28 70 6f 73 |pos_:=VE|CTOR(pos|
|00004880| 69 74 69 6f 6e 28 69 76 | 61 72 73 20 53 55 42 20 |ition(iv|ars SUB |
|00004890| 69 5f 2c 76 61 72 73 29 | 2c 69 5f 2c 31 2c 44 49 |i_,vars)|,i_,1,DI|
|000048a0| 4d 28 69 76 61 72 73 29 | 29 2c 49 46 28 4d 45 4d |M(ivars)|),IF(MEM|
|000048b0| 42 45 52 3f 28 30 2c 70 | 6f 73 5f 29 2c 52 45 54 |BER?(0,p|os_),RET|
|000048c0| 55 52 4e 28 22 55 6e 64 | 65 66 69 6e 65 64 20 69 |URN("Und|efined i|
|000048d0| 6e 64 65 70 65 6e 64 65 | 6e 74 20 76 61 72 69 61 |ndepende|nt varia|
|000048e0| 62 6c 65 28 73 29 21 22 | 29 2c 6d 5f 3a 3d 64 61 |ble(s)!"|),m_:=da|
|000048f0| 74 61 31 20 53 55 42 20 | 20 53 55 42 20 70 6f 73 |ta1 SUB | SUB pos|
|00004900| 5f 29 2c 79 63 70 74 3a | 3d 49 46 28 4d 45 4d 42 |_),ycpt:|=IF(MEMB|
|00004910| 45 52 3f 28 31 2c 47 52 | 41 44 28 65 71 71 2c 70 |ER?(1,GR|AD(eqq,p|
|00004920| 61 72 6d 29 29 2c 31 2c | 30 2c 30 29 2c 6b 5f 3a |arm)),1,|0,0),k_:|
|00004930| 2d 79 63 70 74 2c 62 65 | 74 61 5f 3a 3d 62 65 74 |-ycpt,be|ta_:=bet|
|00004940| 61 30 5f 2c 65 71 31 3a | 3d 56 45 43 54 4f 52 28 |a0_,eq1:|=VECTOR(|
|00004950| 53 55 42 53 54 28 65 71 | 71 2c 69 76 61 72 73 2c |SUBST(eq|q,ivars,|
|00004960| 6d 5f 20 53 55 42 20 69 | 29 2c 69 2c 6f 62 73 5f |m_ SUB i|),i,obs_|
|00004970| 29 2c 6a 5f 3a 3d 47 52 | 41 44 28 65 71 31 2c 70 |),j_:=GR|AD(eq1,p|
|00004980| 61 72 6d 29 60 2c 70 5f | 3a 3d 53 55 42 53 54 28 |arm)`,p_|:=SUBST(|
|00004990| 65 71 31 2c 70 61 72 6d | 2c 62 65 74 61 5f 29 2c |eq1,parm|,beta_),|
|000049a0| 72 5f 3a 3d 79 5f 2d 70 | 5f 2c 73 73 65 5f 3a 3d |r_:=y_-p|_,sse_:=|
|000049b0| 72 5f 20 2e 20 72 5f 2c | 74 5f 3a 3d 49 46 28 79 |r_ . r_,|t_:=IF(y|
|000049c0| 63 70 74 3d 31 2c 79 5f | 2d 56 45 43 54 4f 52 28 |cpt=1,y_|-VECTOR(|
|000049d0| 41 56 45 52 41 47 45 28 | 79 5f 29 2c 69 2c 6f 62 |AVERAGE(|y_),i,ob|
|000049e0| 73 5f 29 2c 79 5f 29 2c | 73 73 74 5f 3a 3d 74 5f |s_),y_),|sst_:=t_|
|000049f0| 20 2e 20 74 5f 2c 6c 61 | 6d 62 64 61 5f 3a 3d 30 | . t_,la|mbda_:=0|
|00004a00| 2e 30 30 31 2c 69 74 65 | 72 3a 3d 5b 41 50 50 45 |.001,ite|r:=[APPE|
|00004a10| 4e 44 28 5b 30 5d 2c 5b | 22 20 22 5d 2c 62 65 74 |ND([0],[|" "],bet|
|00004a20| 61 5f 2c 5b 73 73 65 5f | 5d 29 5d 2c 65 70 73 5f |a_,[sse_|])],eps_|
|00004a30| 3a 3d 31 2c 69 74 65 72 | 5f 3a 3d 31 2c 4c 4f 4f |:=1,iter|_:=1,LOO|
|00004a40| 50 28 49 46 28 69 74 65 | 72 5f 3e 6e 5f 20 4f 52 |P(IF(ite|r_>n_ OR|
|00004a50| 20 65 70 73 5f 3c 63 63 | 5f 2c 65 78 69 74 29 2c | eps_<cc|_,exit),|
|00004a60| 78 5f 3a 3d 53 55 42 53 | 54 28 6a 5f 2c 70 61 72 |x_:=SUBS|T(j_,par|
|00004a70| 6d 2c 62 65 74 61 5f 29 | 2c 6c 73 73 65 5f 3a 3d |m,beta_)|,lsse_:=|
|00004a80| 73 73 65 5f 2c 78 70 78 | 3a 3d 78 5f 60 20 2e 20 |sse_,xpx|:=x_` . |
|00004a90| 78 5f 2c 64 69 61 67 5f | 3a 3d 56 45 43 54 4f 52 |x_,diag_|:=VECTOR|
|00004aa0| 28 56 45 43 54 4f 52 28 | 49 46 28 71 5f 3d 73 5f |(VECTOR(|IF(q_=s_|
|00004ab0| 2c 78 70 78 20 53 55 42 | 20 71 5f 20 53 55 42 20 |,xpx SUB| q_ SUB |
|00004ac0| 73 5f 2c 30 29 2c 73 5f | 2c 6b 5f 2b 79 63 70 74 |s_,0),s_|,k_+ycpt|
|00004ad0| 29 2c 71 5f 2c 6b 5f 2b | 79 63 70 74 29 2c 78 70 |),q_,k_+|ycpt),xp|
|00004ae0| 78 69 3a 3d 28 78 70 78 | 2b 6c 61 6d 62 64 61 5f |xi:=(xpx|+lambda_|
|00004af0| 2a 64 69 61 67 5f 29 5e | 28 2d 31 29 2c 64 65 6c |*diag_)^|(-1),del|
|00004b00| 74 61 5f 3a 3d 78 70 78 | 69 20 2e 20 78 5f 60 20 |ta_:=xpx|i . x_` |
|00004b10| 2e 20 72 5f 2c 62 65 74 | 61 5f 3a 2b 64 65 6c 74 |. r_,bet|a_:+delt|
|00004b20| 61 5f 2c 70 5f 3a 3d 53 | 55 42 53 54 28 65 71 31 |a_,p_:=S|UBST(eq1|
|00004b30| 2c 70 61 72 6d 2c 62 65 | 74 61 5f 29 2c 72 5f 3a |,parm,be|ta_),r_:|
|00004b40| 3d 79 5f 2d 70 5f 2c 73 | 73 65 5f 3a 3d 72 5f 20 |=y_-p_,s|se_:=r_ |
|00004b50| 2e 20 72 5f 2c 69 74 65 | 72 3a 3d 41 50 50 45 4e |. r_,ite|r:=APPEN|
|00004b60| 44 28 69 74 65 72 2c 5b | 41 50 50 45 4e 44 28 5b |D(iter,[|APPEND([|
|00004b70| 69 74 65 72 5f 5d 2c 5b | 6c 61 6d 62 64 61 5f 5d |iter_],[|lambda_]|
|00004b80| 2c 62 65 74 61 5f 2c 5b | 73 73 65 5f 5d 29 5d 29 |,beta_,[|sse_])])|
|00004b90| 2c 49 46 28 73 73 65 5f | 3c 3d 6c 73 73 65 5f 2c |,IF(sse_|<=lsse_,|
|00004ba0| 50 52 4f 47 28 65 70 73 | 5f 3a 3d 41 42 53 28 28 |PROG(eps|_:=ABS((|
|00004bb0| 6c 73 73 65 5f 2d 73 73 | 65 5f 29 2f 28 73 73 65 |lsse_-ss|e_)/(sse|
|00004bc0| 5f 2b 31 30 5e 28 2d 36 | 29 29 29 2c 6c 61 6d 62 |_+10^(-6|))),lamb|
|00004bd0| 64 61 5f 3a 2f 31 30 29 | 2c 50 52 4f 47 28 62 65 |da_:/10)|,PROG(be|
|00004be0| 74 61 5f 3a 2d 64 65 6c | 74 61 5f 2c 6c 61 6d 62 |ta_:-del|ta_,lamb|
|00004bf0| 64 61 5f 3a 2a 31 30 29 | 2c 52 45 54 55 52 4e 28 |da_:*10)|,RETURN(|
|00004c00| 22 43 61 6e 6e 6f 74 20 | 72 65 73 6f 6c 76 65 20 |"Cannot |resolve |
|00004c10| 53 53 45 20 74 65 73 74 | 21 22 29 29 2c 69 74 65 |SSE test|!")),ite|
|00004c20| 72 5f 3a 2b 31 29 2c 49 | 46 28 69 74 65 72 5f 3e |r_:+1),I|F(iter_>|
|00004c30| 6e 5f 2c 6d 73 67 3a 3d | 41 50 50 45 4e 44 28 22 |n_,msg:=|APPEND("|
|00004c40| 43 6f 6e 76 65 72 67 65 | 6e 63 65 20 66 61 69 6c |Converge|nce fail|
|00004c50| 65 64 20 61 66 74 65 72 | 20 22 2c 6e 5f 2c 22 20 |ed after| ",n_," |
|00004c60| 69 74 65 72 61 74 69 6f | 6e 73 21 22 29 2c 6d 73 |iteratio|ns!"),ms|
|00004c70| 67 3a 3d 22 43 6f 6e 76 | 65 72 67 65 6e 63 65 20 |g:="Conv|ergence |
|00004c80| 63 72 69 74 65 72 69 61 | 20 6d 65 74 21 22 29 2c |criteria| met!"),|
|00004c90| 64 66 65 5f 3a 3d 6f 62 | 73 5f 2d 6b 5f 2d 79 63 |dfe_:=ob|s_-k_-yc|
|00004ca0| 70 74 2c 6d 73 65 5f 3a | 3d 73 73 65 5f 2f 64 66 |pt,mse_:|=sse_/df|
|00004cb0| 65 5f 2c 73 65 5f 3a 3d | 53 51 52 54 28 6d 73 65 |e_,se_:=|SQRT(mse|
|00004cc0| 5f 29 2c 73 74 64 5f 3a | 3d 56 45 43 54 4f 52 28 |_),std_:|=VECTOR(|
|00004cd0| 53 51 52 54 28 78 70 78 | 69 20 53 55 42 20 69 5f |SQRT(xpx|i SUB i_|
|00004ce0| 20 53 55 42 20 69 5f 2a | 6d 73 65 5f 29 2c 69 5f | SUB i_*|mse_),i_|
|00004cf0| 2c 6b 5f 2b 79 63 70 74 | 29 2c 74 5f 3a 3d 56 45 |,k_+ycpt|),t_:=VE|
|00004d00| 43 54 4f 52 28 62 65 74 | 61 5f 20 53 55 42 20 69 |CTOR(bet|a_ SUB i|
|00004d10| 5f 2f 73 74 64 5f 20 53 | 55 42 20 69 5f 2c 69 5f |_/std_ S|UB i_,i_|
|00004d20| 2c 6b 5f 2b 79 63 70 74 | 29 2c 70 72 6f 62 5f 3a |,k_+ycpt|),prob_:|
|00004d30| 3d 56 45 43 54 4f 52 28 | 31 2d 53 54 55 44 45 4e |=VECTOR(|1-STUDEN|
|00004d40| 54 28 74 5f 20 53 55 42 | 20 69 5f 2c 64 66 65 5f |T(t_ SUB| i_,dfe_|
|00004d50| 29 2c 69 5f 2c 6b 5f 2b | 79 63 70 74 29 2c 73 73 |),i_,k_+|ycpt),ss|
|00004d60| 72 5f 3a 3d 73 73 74 5f | 2d 73 73 65 5f 2c 72 73 |r_:=sst_|-sse_,rs|
|00004d70| 71 3a 3d 73 73 72 5f 2f | 73 73 74 5f 2c 61 72 73 |q:=ssr_/|sst_,ars|
|00004d80| 71 3a 3d 31 2d 73 73 65 | 5f 2f 28 6f 62 73 5f 2d |q:=1-sse|_/(obs_-|
|00004d90| 6b 5f 2d 79 63 70 74 29 | 2f 28 73 73 74 5f 2f 28 |k_-ycpt)|/(sst_/(|
|00004da0| 6f 62 73 5f 2d 79 63 70 | 74 29 29 2c 66 64 5f 3a |obs_-ycp|t)),fd_:|
|00004db0| 3d 46 5f 44 49 53 54 52 | 49 42 55 54 49 4f 4e 28 |=F_DISTR|IBUTION(|
|00004dc0| 73 73 72 5f 2a 64 66 65 | 5f 2f 28 73 73 65 5f 2a |ssr_*dfe|_/(sse_*|
|00004dd0| 6b 5f 29 2c 6b 5f 2c 64 | 66 65 5f 29 2c 61 31 5f |k_),k_,d|fe_),a1_|
|00004de0| 3a 3d 5b 22 53 6f 75 72 | 63 65 22 2c 22 44 46 22 |:=["Sour|ce","DF"|
|00004df0| 2c 22 53 53 22 2c 22 4d | 53 22 2c 22 46 22 2c 22 |,"SS","M|S","F","|
|00004e00| 50 72 6f 62 28 46 29 22 | 5d 2c 61 32 5f 3a 3d 5b |Prob(F)"|],a2_:=[|
|00004e10| 22 52 65 67 22 2c 6b 5f | 2c 73 73 72 5f 2c 73 73 |"Reg",k_|,ssr_,ss|
|00004e20| 72 5f 2f 6b 5f 2c 73 73 | 72 5f 2a 64 66 65 5f 2f |r_/k_,ss|r_*dfe_/|
|00004e30| 28 73 73 65 5f 2a 6b 5f | 29 2c 49 46 28 66 64 5f |(sse_*k_|),IF(fd_|
|00004e40| 3c 30 2c 30 2c 66 64 5f | 29 5d 2c 61 33 5f 3a 3d |<0,0,fd_|)],a3_:=|
|00004e50| 5b 22 45 72 72 6f 72 22 | 2c 64 66 65 5f 2c 73 73 |["Error"|,dfe_,ss|
|00004e60| 65 5f 2c 73 73 65 5f 2f | 64 66 65 5f 2c 22 20 22 |e_,sse_/|dfe_," "|
|00004e70| 2c 22 20 22 5d 2c 61 34 | 5f 3a 3d 5b 49 46 28 79 |," "],a4|_:=[IF(y|
|00004e80| 63 70 74 3d 31 2c 22 43 | 6f 72 72 65 63 74 65 64 |cpt=1,"C|orrected|
|00004e90| 20 54 6f 74 61 6c 22 2c | 22 55 6e 63 6f 72 72 65 | Total",|"Uncorre|
|00004ea0| 63 74 65 64 20 54 6f 74 | 61 6c 22 29 2c 64 66 65 |cted Tot|al"),dfe|
|00004eb0| 5f 2b 6b 5f 2c 73 73 74 | 5f 2c 22 20 22 2c 22 20 |_+k_,sst|_," "," |
|00004ec0| 22 2c 22 20 22 5d 2c 61 | 6e 6f 76 61 3a 3d 41 50 |"," "],a|nova:=AP|
|00004ed0| 50 45 4e 44 28 5b 61 31 | 5f 5d 2c 5b 61 32 5f 5d |PEND([a1|_],[a2_]|
|00004ee0| 2c 5b 61 33 5f 5d 2c 5b | 61 34 5f 5d 29 2c 66 65 |,[a3_],[|a4_]),fe|
|00004ef0| 71 3a 3d 64 76 61 72 3d | 53 55 42 53 54 28 65 71 |q:=dvar=|SUBST(eq|
|00004f00| 71 2c 70 61 72 6d 2c 62 | 65 74 61 5f 29 2c 69 74 |q,parm,b|eta_),it|
|00004f10| 65 72 3a 3d 41 50 50 45 | 4e 44 28 5b 41 50 50 45 |er:=APPE|ND([APPE|
|00004f20| 4e 44 28 5b 22 49 74 65 | 72 22 5d 2c 5b 22 bf 22 |ND(["Ite|r"],["."|
|00004f30| 5d 2c 70 61 72 6d 2c 5b | 22 53 53 45 22 5d 29 5d |],parm,[|"SSE"])]|
|00004f40| 2c 69 74 65 72 29 2c 74 | 69 74 6c 65 3a 3d 5b 5b |,iter),t|itle:=[[|
|00004f50| 22 50 61 72 6d 22 5d 2c | 5b 22 56 61 6c 75 65 22 |"Parm"],|["Value"|
|00004f60| 5d 2c 5b 22 53 54 44 22 | 5d 2c 5b 41 50 50 45 4e |],["STD"|],[APPEN|
|00004f70| 44 28 22 74 28 22 2c 64 | 66 65 5f 2c 22 29 22 29 |D("t(",d|fe_,")")|
|00004f80| 5d 2c 5b 22 50 72 6f 62 | 28 74 29 22 5d 5d 2c 6f |],["Prob|(t)"]],o|
|00004f90| 75 74 70 75 74 3a 3d 41 | 50 50 45 4e 44 5f 43 4f |utput:=A|PPEND_CO|
|00004fa0| 4c 55 4d 4e 53 28 74 69 | 74 6c 65 2c 41 50 50 45 |LUMNS(ti|tle,APPE|
|00004fb0| 4e 44 28 5b 70 61 72 6d | 5d 2c 5b 62 65 74 61 5f |ND([parm|],[beta_|
|00004fc0| 5d 2c 5b 73 74 64 5f 5d | 2c 5b 74 5f 5d 2c 5b 70 |],[std_]|,[t_],[p|
|00004fd0| 72 6f 62 5f 5d 29 29 60 | 2c 73 74 61 74 73 3a 3d |rob_]))`|,stats:=|
|00004fe0| 5b 5b 22 53 45 22 2c 22 | 52 5e 32 22 2c 22 41 64 |[["SE","|R^2","Ad|
|00004ff0| 6a 52 5e 32 22 5d 2c 5b | 73 65 5f 2c 72 73 71 2c |jR^2"],[|se_,rsq,|
|00005000| 61 72 73 71 5d 5d 2c 52 | 45 54 55 52 4e 28 5b 5b |arsq]],R|ETURN([[|
|00005010| 22 4d 61 72 71 75 61 72 | 64 74 20 4d 65 74 68 6f |"Marquar|dt Metho|
|00005020| 64 22 5d 2c 5b 22 20 22 | 5d 2c 5b 6d 73 67 5d 2c |d"],[" "|],[msg],|
|00005030| 5b 22 20 22 5d 2c 5b 6f | 75 74 70 75 74 5d 2c 5b |[" "],[o|utput],[|
|00005040| 22 20 22 5d 2c 5b 61 6e | 6f 76 61 5d 2c 5b 22 20 |" "],[an|ova],[" |
|00005050| 22 5d 2c 5b 73 74 61 74 | 73 5d 2c 5b 22 20 22 5d |"],[stat|s],[" "]|
|00005060| 2c 5b 66 65 71 5d 5d 29 | 29 01 80 08 00 00 00 16 |,[feq]])|).......|
|00005070| 0b 00 00 cd 03 00 00 66 | 0b 00 00 00 ff e2 02 7b |.......f|.......{|
|00005080| 5c 72 74 66 31 5c 61 6e | 73 69 5c 64 65 66 66 30 |\rtf1\an|si\deff0|
|00005090| 5c 64 65 66 74 61 62 37 | 32 30 7b 5c 66 6f 6e 74 |\deftab7|20{\font|
|000050a0| 74 62 6c 7b 5c 66 30 5c | 66 73 77 69 73 73 20 4d |tbl{\f0\|fswiss M|
|000050b0| 53 20 53 61 6e 73 20 53 | 65 72 69 66 3b 7d 7b 5c |S Sans S|erif;}{\|
|000050c0| 66 31 5c 66 64 65 63 6f | 72 5c 66 63 68 61 72 73 |f1\fdeco|r\fchars|
|000050d0| 65 74 32 20 53 79 6d 62 | 6f 6c 3b 7d 7b 5c 66 32 |et2 Symb|ol;}{\f2|
|000050e0| 5c 66 73 77 69 73 73 5c | 66 70 72 71 32 20 53 79 |\fswiss\|fprq2 Sy|
|000050f0| 73 74 65 6d 3b 7d 7b 5c | 66 33 5c 66 73 77 69 73 |stem;}{\|f3\fswis|
|00005100| 73 5c 66 70 72 71 32 20 | 41 72 69 61 6c 3b 7d 7b |s\fprq2 |Arial;}{|
|00005110| 5c 66 34 5c 66 6d 6f 64 | 65 72 6e 5c 66 63 68 61 |\f4\fmod|ern\fcha|
|00005120| 72 73 65 74 32 20 44 66 | 57 35 20 50 72 69 6e 74 |rset2 Df|W5 Print|
|00005130| 65 72 3b 7d 7d 0d 0a 7b | 5c 63 6f 6c 6f 72 74 62 |er;}}..{|\colortb|
|00005140| 6c 5c 72 65 64 30 5c 67 | 72 65 65 6e 30 5c 62 6c |l\red0\g|reen0\bl|
|00005150| 75 65 30 3b 5c 72 65 64 | 32 35 35 5c 67 72 65 65 |ue0;\red|255\gree|
|00005160| 6e 30 5c 62 6c 75 65 30 | 3b 5c 72 65 64 30 5c 67 |n0\blue0|;\red0\g|
|00005170| 72 65 65 6e 30 5c 62 6c | 75 65 32 35 35 3b 7d 0d |reen0\bl|ue255;}.|
|00005180| 0a 5c 64 65 66 6c 61 6e | 67 31 30 33 33 5c 70 61 |.\deflan|g1033\pa|
|00005190| 72 64 5c 70 6c 61 69 6e | 5c 66 33 5c 66 73 32 30 |rd\plain|\f3\fs20|
|000051a0| 5c 63 66 30 20 0d 0a 5c | 70 61 72 20 4f 6e 63 65 |\cf0 ..\|par Once|
|000051b0| 20 61 20 72 65 67 72 65 | 73 73 69 6f 6e 20 68 61 | a regre|ssion ha|
|000051c0| 73 20 62 65 65 6e 20 72 | 75 6e 2c 20 65 73 74 69 |s been r|un, esti|
|000051d0| 6d 61 74 65 64 20 76 61 | 6c 75 65 73 20 61 6e 64 |mated va|lues and|
|000051e0| 20 74 68 65 69 72 20 63 | 6f 6e 66 69 64 65 6e 63 | their c|onfidenc|
|000051f0| 65 20 69 6e 74 65 72 76 | 61 6c 73 20 63 61 6e 20 |e interv|als can |
|00005200| 62 65 20 63 61 6c 63 75 | 6c 61 74 65 64 20 77 69 |be calcu|lated wi|
|00005210| 74 68 20 5c 70 6c 61 69 | 6e 5c 66 33 5c 66 73 32 |th \plai|n\f3\fs2|
|00005220| 30 5c 63 66 32 20 50 72 | 65 64 69 63 74 65 64 5f |0\cf2 Pr|edicted_|
|00005230| 56 61 6c 75 65 73 5c 70 | 6c 61 69 6e 5c 66 33 5c |Values\p|lain\f3\|
|00005240| 66 73 32 30 5c 63 66 30 | 20 2e 5c 70 6c 61 69 6e |fs20\cf0| .\plain|
|00005250| 5c 66 33 5c 66 73 32 30 | 5c 63 66 31 20 20 20 5c |\f3\fs20|\cf1   \|
|00005260| 70 6c 61 69 6e 5c 66 33 | 5c 66 73 32 30 5c 63 66 |plain\f3|\fs20\cf|
|00005270| 30 20 45 6e 74 65 72 20 | 61 20 76 65 63 74 6f 72 |0 Enter |a vector|
|00005280| 20 6f 66 20 58 30 20 76 | 61 6c 75 65 73 20 61 6e | of X0 v|alues an|
|00005290| 64 20 74 68 65 20 63 6f | 6e 66 69 64 65 6e 63 65 |d the co|nfidence|
|000052a0| 20 69 6e 74 65 72 76 61 | 6c 20 64 65 73 69 72 65 | interva|l desire|
|000052b0| 64 2e 20 20 45 2e 67 2e | 2c 20 69 66 20 79 6f 75 |d.  E.g.|, if you|
|000052c0| 20 77 61 6e 74 20 61 20 | 39 35 25 20 63 6f 6e 66 | want a |95% conf|
|000052d0| 69 64 65 6e 63 65 20 69 | 6e 74 65 72 76 61 6c 2c |idence i|nterval,|
|000052e0| 20 65 6e 74 65 72 20 2e | 39 35 2e 20 20 49 74 20 | enter .|95.  It |
|000052f0| 77 69 6c 6c 20 63 6f 6d | 70 75 74 65 20 61 6e 20 |will com|pute an |
|00005300| 69 6e 74 65 72 76 61 6c | 20 65 73 74 69 6d 61 74 |interval| estimat|
|00005310| 65 20 66 6f 72 20 74 68 | 65 20 6d 65 61 6e 20 76 |e for th|e mean v|
|00005320| 61 6c 75 65 20 6f 66 20 | 79 20 61 6e 64 20 6f 6e |alue of |y and on|
|00005330| 65 20 66 6f 72 20 74 68 | 65 20 69 6e 64 69 76 69 |e for th|e indivi|
|00005340| 64 75 61 6c 20 76 61 6c | 75 65 20 6f 66 20 79 2e |dual val|ue of y.|
|00005350| 0d 0a 5c 70 61 72 20 0d | 0a 5c 70 61 72 20 7d 0d |..\par .|.\par }.|
|00005360| 0a 03 80 38 00 00 00 72 | 0b 00 00 70 03 00 00 56 |...8...r|...p...V|
|00005370| 0c 00 00 00 00 00 00 00 | 00 00 00 f0 bf 06 00 00 |........|........|
|00005380| 00 01 00 00 00 ff 18 02 | 50 72 65 64 69 63 74 65 |........|Predicte|
|00005390| 64 5f 56 61 6c 75 65 73 | 28 78 30 2c 63 69 3a 3d |d_Values|(x0,ci:=|
|000053a0| 30 2e 39 35 29 3a 3d 50 | 52 4f 47 28 70 72 64 3a |0.95):=P|ROG(prd:|
|000053b0| 3d 53 55 42 53 54 28 52 | 48 53 28 66 65 71 29 2c |=SUBST(R|HS(feq),|
|000053c0| 69 76 61 72 73 31 2c 78 | 30 29 2c 67 5f 3a 3d 47 |ivars1,x|0),g_:=G|
|000053d0| 52 41 44 28 65 71 71 2c | 70 61 72 6d 31 29 2c 67 |RAD(eqq,|parm1),g|
|000053e0| 76 3a 3d 53 55 42 53 54 | 28 67 5f 2c 41 50 50 45 |v:=SUBST|(g_,APPE|
|000053f0| 4e 44 28 70 61 72 6d 31 | 2c 69 76 61 72 73 31 29 |ND(parm1|,ivars1)|
|00005400| 2c 41 50 50 45 4e 44 28 | 62 65 74 61 5f 2c 78 30 |,APPEND(|beta_,x0|
|00005410| 29 29 2c 69 5f 3a 3d 67 | 76 20 2e 20 78 70 78 69 |)),i_:=g|v . xpxi|
|00005420| 20 2e 20 67 76 60 2c 73 | 65 70 3a 3d 73 65 5f 2a | . gv`,s|ep:=se_*|
|00005430| 53 51 52 54 28 69 5f 29 | 2c 73 65 66 3a 3d 73 65 |SQRT(i_)|,sef:=se|
|00005440| 5f 2a 53 51 52 54 28 31 | 2b 69 5f 29 2c 50 52 4f |_*SQRT(1|+i_),PRO|
|00005450| 47 28 50 72 65 63 69 73 | 69 6f 6e 44 69 67 69 74 |G(Precis|ionDigit|
|00005460| 73 3a 3d 36 2c 69 6e 76 | 74 3a 3d 52 48 53 28 4e |s:=6,inv|t:=RHS(N|
|00005470| 53 4f 4c 56 45 28 53 54 | 55 44 45 4e 54 28 74 2c |SOLVE(ST|UDENT(t,|
|00005480| 64 66 65 5f 29 2d 63 69 | 2c 74 2c 30 2c 69 6e 66 |dfe_)-ci|,t,0,inf|
|00005490| 29 29 2c 50 72 65 63 69 | 73 69 6f 6e 44 69 67 69 |)),Preci|sionDigi|
|000054a0| 74 73 3a 3d 31 30 2c 4e | 6f 74 61 74 69 6f 6e 44 |ts:=10,N|otationD|
|000054b0| 69 67 69 74 73 3a 3d 38 | 29 2c 63 79 68 31 3a 3d |igits:=8|),cyh1:=|
|000054c0| 70 72 64 2d 73 65 70 2a | 69 6e 76 74 2c 63 79 68 |prd-sep*|invt,cyh|
|000054d0| 32 3a 3d 70 72 64 2b 73 | 65 70 2a 69 6e 76 74 2c |2:=prd+s|ep*invt,|
|000054e0| 63 79 30 31 3a 3d 70 72 | 64 2d 73 65 66 2a 69 6e |cy01:=pr|d-sef*in|
|000054f0| 76 74 2c 63 79 30 32 3a | 3d 70 72 64 2b 73 65 66 |vt,cy02:|=prd+sef|
|00005500| 2a 69 6e 76 74 2c 74 69 | 74 6c 65 3a 3d 41 50 50 |*invt,ti|tle:=APP|
|00005510| 45 4e 44 28 53 54 52 49 | 4e 47 28 63 69 2a 31 30 |END(STRI|NG(ci*10|
|00005520| 30 29 2c 22 25 20 43 6f | 6e 66 69 64 65 6e 63 65 |0),"% Co|nfidence|
|00005530| 20 49 6e 74 65 72 76 61 | 6c 22 29 2c 5b 5b 74 69 | Interva|l"),[[ti|
|00005540| 74 6c 65 5d 2c 5b 5b 5b | 22 56 61 6c 75 65 22 2c |tle],[[[|"Value",|
|00005550| 70 72 64 2c 22 20 22 5d | 2c 5b 22 53 65 5f 79 68 |prd," "]|,["Se_yh|
|00005560| 61 74 2f 53 65 5f 59 30 | 22 2c 73 65 70 2c 73 65 |at/Se_Y0|",sep,se|
|00005570| 66 5d 2c 5b 22 43 49 5f | 79 68 61 74 22 2c 63 79 |f],["CI_|yhat",cy|
|00005580| 68 31 2c 63 79 68 32 5d | 2c 5b 22 43 49 5f 59 30 |h1,cyh2]|,["CI_Y0|
|00005590| 22 2c 63 79 30 31 2c 63 | 79 30 32 5d 5d 5d 5d 29 |",cy01,c|y02]]]])|
|000055a0| 01 80 08 00 00 00 62 0c | 00 00 cd 03 00 00 b2 0c |......b.|........|
|000055b0| 00 00 00 ff 8f 02 7b 5c | 72 74 66 31 5c 61 6e 73 |......{\|rtf1\ans|
|000055c0| 69 5c 64 65 66 66 30 5c | 64 65 66 74 61 62 37 32 |i\deff0\|deftab72|
|000055d0| 30 7b 5c 66 6f 6e 74 74 | 62 6c 7b 5c 66 30 5c 66 |0{\fontt|bl{\f0\f|
|000055e0| 73 77 69 73 73 20 4d 53 | 20 53 61 6e 73 20 53 65 |swiss MS| Sans Se|
|000055f0| 72 69 66 3b 7d 7b 5c 66 | 31 5c 66 64 65 63 6f 72 |rif;}{\f|1\fdecor|
|00005600| 5c 66 63 68 61 72 73 65 | 74 32 20 53 79 6d 62 6f |\fcharse|t2 Symbo|
|00005610| 6c 3b 7d 7b 5c 66 32 5c | 66 73 77 69 73 73 5c 66 |l;}{\f2\|fswiss\f|
|00005620| 70 72 71 32 20 53 79 73 | 74 65 6d 3b 7d 7b 5c 66 |prq2 Sys|tem;}{\f|
|00005630| 33 5c 66 73 77 69 73 73 | 5c 66 70 72 71 32 20 41 |3\fswiss|\fprq2 A|
|00005640| 72 69 61 6c 3b 7d 7b 5c | 66 34 5c 66 6d 6f 64 65 |rial;}{\|f4\fmode|
|00005650| 72 6e 5c 66 63 68 61 72 | 73 65 74 32 20 44 66 57 |rn\fchar|set2 DfW|
|00005660| 35 20 50 72 69 6e 74 65 | 72 3b 7d 7d 0d 0a 7b 5c |5 Printe|r;}}..{\|
|00005670| 63 6f 6c 6f 72 74 62 6c | 5c 72 65 64 30 5c 67 72 |colortbl|\red0\gr|
|00005680| 65 65 6e 30 5c 62 6c 75 | 65 30 3b 5c 72 65 64 32 |een0\blu|e0;\red2|
|00005690| 35 35 5c 67 72 65 65 6e | 30 5c 62 6c 75 65 30 3b |55\green|0\blue0;|
|000056a0| 5c 72 65 64 30 5c 67 72 | 65 65 6e 30 5c 62 6c 75 |\red0\gr|een0\blu|
|000056b0| 65 32 35 35 3b 7d 0d 0a | 5c 64 65 66 6c 61 6e 67 |e255;}..|\deflang|
|000056c0| 31 30 33 33 5c 70 61 72 | 64 5c 70 6c 61 69 6e 5c |1033\par|d\plain\|
|000056d0| 66 33 5c 66 73 32 30 5c | 63 66 31 20 0d 0a 5c 70 |f3\fs20\|cf1 ..\p|
|000056e0| 61 72 20 4e 6f 74 65 3a | 20 57 68 69 6c 65 20 74 |ar Note:| While t|
|000056f0| 68 65 20 64 6f 63 75 6d | 65 6e 74 61 74 69 6f 6e |he docum|entation|
|00005700| 20 6f 6e 20 74 68 65 20 | 73 74 75 64 65 6e 74 27 | on the |student'|
|00005710| 73 20 74 20 64 69 73 74 | 72 69 62 75 74 69 6f 6e |s t dist|ribution|
|00005720| 20 73 61 79 73 20 69 74 | 20 67 69 76 65 73 20 74 | says it| gives t|
|00005730| 68 65 20 63 75 6d 75 6c | 61 74 69 76 65 20 70 72 |he cumul|ative pr|
|00005740| 6f 62 61 62 69 6c 69 74 | 79 2c 20 69 74 20 72 65 |obabilit|y, it re|
|00005750| 61 6c 6c 79 20 67 69 76 | 65 73 20 74 68 65 20 62 |ally giv|es the b|
|00005760| 65 74 77 65 65 6e 20 74 | 61 69 6c 73 20 70 72 6f |etween t|ails pro|
|00005770| 62 61 62 69 6c 69 74 79 | 20 66 6f 72 20 5c 70 6c |bability| for \pl|
|00005780| 61 69 6e 5c 66 34 5c 66 | 73 32 30 5c 63 66 31 20 |ain\f4\f|s20\cf1 |
|00005790| 5c 27 62 31 5c 70 6c 61 | 69 6e 5c 66 33 5c 66 73 |\'b1\pla|in\f3\fs|
|000057a0| 32 30 5c 63 66 31 20 74 | 2e 0d 0a 5c 70 61 72 20 |20\cf1 t|...\par |
|000057b0| 5c 70 6c 61 69 6e 5c 66 | 33 5c 66 73 32 30 5c 63 |\plain\f|3\fs20\c|
|000057c0| 66 30 20 0d 0a 5c 70 61 | 72 20 54 68 65 20 66 6f |f0 ..\pa|r The fo|
|000057d0| 6c 6c 6f 77 69 6e 67 20 | 66 75 6e 63 74 69 6f 6e |llowing |function|
|000057e0| 20 63 6f 6d 70 75 74 65 | 73 20 61 20 6d 61 74 72 | compute|s a matr|
|000057f0| 69 78 20 6f 66 20 72 65 | 73 69 64 75 61 6c 73 20 |ix of re|siduals |
|00005800| 74 68 61 74 20 6d 61 79 | 20 62 65 20 61 6e 61 6c |that may| be anal|
|00005810| 79 7a 65 64 20 61 6e 64 | 20 70 6c 6f 74 74 65 64 |yzed and| plotted|
|00005820| 2e 0d 0a 5c 70 61 72 20 | 5c 70 6c 61 69 6e 5c 66 |...\par |\plain\f|
|00005830| 33 5c 66 73 32 30 5c 63 | 66 31 20 0d 0a 5c 70 61 |3\fs20\c|f1 ..\pa|
|00005840| 72 20 7d 0d 0a 03 80 38 | 00 00 00 be 0c 00 00 98 |r }....8|........|
|00005850| 01 00 00 d6 0c 00 00 00 | 00 00 00 00 00 00 00 f0 |........|........|
|00005860| bf 07 00 00 00 01 00 00 | 00 2b 52 65 73 69 64 75 |........|.+Residu|
|00005870| 61 6c 73 3a 3d 56 45 43 | 54 4f 52 28 5b 69 2c 72 |als:=VEC|TOR([i,r|
|00005880| 5f 20 53 55 42 20 69 5d | 2c 69 2c 31 2c 44 49 4d |_ SUB i]|,i,1,DIM|
|00005890| 28 72 5f 29 29 01 80 08 | 00 00 00 e2 0c 00 00 cd |(r_))...|........|
|000058a0| 03 00 00 02 0d 00 00 00 | ff 1c 02 7b 5c 72 74 66 |........|...{\rtf|
|000058b0| 31 5c 61 6e 73 69 5c 64 | 65 66 66 30 5c 64 65 66 |1\ansi\d|eff0\def|
|000058c0| 74 61 62 37 32 30 7b 5c | 66 6f 6e 74 74 62 6c 7b |tab720{\|fonttbl{|
|000058d0| 5c 66 30 5c 66 73 77 69 | 73 73 20 4d 53 20 53 61 |\f0\fswi|ss MS Sa|
|000058e0| 6e 73 20 53 65 72 69 66 | 3b 7d 7b 5c 66 31 5c 66 |ns Serif|;}{\f1\f|
|000058f0| 64 65 63 6f 72 5c 66 63 | 68 61 72 73 65 74 32 20 |decor\fc|harset2 |
|00005900| 53 79 6d 62 6f 6c 3b 7d | 7b 5c 66 32 5c 66 73 77 |Symbol;}|{\f2\fsw|
|00005910| 69 73 73 5c 66 70 72 71 | 32 20 53 79 73 74 65 6d |iss\fprq|2 System|
|00005920| 3b 7d 7b 5c 66 33 5c 66 | 73 77 69 73 73 5c 66 70 |;}{\f3\f|swiss\fp|
|00005930| 72 71 32 20 41 72 69 61 | 6c 3b 7d 7b 5c 66 34 5c |rq2 Aria|l;}{\f4\|
|00005940| 66 6d 6f 64 65 72 6e 5c | 66 63 68 61 72 73 65 74 |fmodern\|fcharset|
|00005950| 32 20 44 66 57 35 20 50 | 72 69 6e 74 65 72 3b 7d |2 DfW5 P|rinter;}|
|00005960| 7d 0d 0a 7b 5c 63 6f 6c | 6f 72 74 62 6c 5c 72 65 |}..{\col|ortbl\re|
|00005970| 64 30 5c 67 72 65 65 6e | 30 5c 62 6c 75 65 30 3b |d0\green|0\blue0;|
|00005980| 5c 72 65 64 32 35 35 5c | 67 72 65 65 6e 30 5c 62 |\red255\|green0\b|
|00005990| 6c 75 65 30 3b 5c 72 65 | 64 30 5c 67 72 65 65 6e |lue0;\re|d0\green|
|000059a0| 30 5c 62 6c 75 65 32 35 | 35 3b 7d 0d 0a 5c 64 65 |0\blue25|5;}..\de|
|000059b0| 66 6c 61 6e 67 31 30 33 | 33 5c 70 61 72 64 5c 70 |flang103|3\pard\p|
|000059c0| 6c 61 69 6e 5c 66 33 5c | 66 73 32 30 5c 63 66 30 |lain\f3\|fs20\cf0|
|000059d0| 20 54 68 65 20 66 6f 6c | 6c 6f 77 69 6e 67 20 65 | The fol|lowing e|
|000059e0| 78 61 6d 70 6c 65 20 75 | 73 65 73 20 61 20 64 61 |xample u|ses a da|
|000059f0| 74 61 73 65 74 20 6f 66 | 20 74 68 65 20 55 2e 53 |taset of| the U.S|
|00005a00| 2e 20 70 6f 70 75 6c 61 | 74 69 6f 6e 20 66 72 6f |. popula|tion fro|
|00005a10| 6d 20 31 37 39 30 20 74 | 6f 20 31 39 39 30 2e 20 |m 1790 t|o 1990. |
|00005a20| 20 54 68 65 20 76 61 72 | 69 61 62 6c 65 20 70 20 | The var|iable p |
|00005a30| 69 73 20 74 68 65 20 70 | 6f 70 75 6c 61 74 69 6f |is the p|opulatio|
|00005a40| 6e 20 69 6e 20 6d 69 6c | 6c 69 6f 6e 73 2e 20 20 |n in mil|lions.  |
|00005a50| 54 68 65 20 76 61 72 69 | 61 62 6c 65 20 79 20 69 |The vari|able y i|
|00005a60| 73 20 74 68 65 20 79 65 | 61 72 2e 20 20 57 65 20 |s the ye|ar.  We |
|00005a70| 63 61 6e 20 75 73 65 20 | 74 68 69 73 20 64 61 74 |can use |this dat|
|00005a80| 61 20 74 6f 20 66 69 74 | 20 61 20 70 72 6f 62 69 |a to fit| a probi|
|00005a90| 74 20 6d 6f 64 65 6c 20 | 6f 66 20 74 68 65 20 66 |t model |of the f|
|00005aa0| 6f 72 6d 20 70 3d 63 2a | 4e 4f 52 4d 41 4c 28 61 |orm p=c*|NORMAL(a|
|00005ab0| 2b 62 2a 28 79 2d 31 37 | 39 30 29 29 2e 0d 0a 5c |+b*(y-17|90))...\|
|00005ac0| 70 61 72 20 7d 0d 0a 03 | 80 38 00 00 00 0e 0d 00 |par }...|.8......|
|00005ad0| 00 f8 00 00 00 e2 0e 00 | 00 00 04 55 73 65 72 00 |........|...User.|
|00005ae0| 00 00 00 00 00 f0 bf 08 | 00 00 00 01 00 00 00 ff |........|........|
|00005af0| 17 01 70 6f 70 3a 3d 5b | 5b 70 2c 79 5d 2c 5b 33 |..pop:=[|[p,y],[3|
|00005b00| 2e 39 32 39 2c 31 37 39 | 30 5d 2c 5b 35 2e 33 30 |.929,179|0],[5.30|
|00005b10| 38 2c 31 38 30 30 5d 2c | 5b 37 2e 32 33 39 2c 31 |8,1800],|[7.239,1|
|00005b20| 38 31 30 5d 2c 5b 39 2e | 36 33 38 2c 31 38 32 30 |810],[9.|638,1820|
|00005b30| 5d 2c 5b 31 32 2e 38 36 | 36 2c 31 38 33 30 5d 2c |],[12.86|6,1830],|
|00005b40| 5b 31 37 2e 30 36 39 2c | 31 38 34 30 5d 2c 5b 32 |[17.069,|1840],[2|
|00005b50| 33 2e 31 39 31 2c 31 38 | 35 30 5d 2c 5b 33 31 2e |3.191,18|50],[31.|
|00005b60| 34 34 33 2c 31 38 36 30 | 5d 2c 5b 33 39 2e 38 31 |443,1860|],[39.81|
|00005b70| 38 2c 31 38 37 30 5d 2c | 5b 35 30 2e 31 35 35 2c |8,1870],|[50.155,|
|00005b80| 31 38 38 30 5d 2c 5b 36 | 32 2e 39 34 37 2c 31 38 |1880],[6|2.947,18|
|00005b90| 39 30 5d 2c 5b 37 35 2e | 39 39 34 2c 31 39 30 30 |90],[75.|994,1900|
|00005ba0| 5d 2c 5b 39 31 2e 39 37 | 32 2c 31 39 31 30 5d 2c |],[91.97|2,1910],|
|00005bb0| 5b 31 30 35 2e 37 31 2c | 31 39 32 30 5d 2c 5b 31 |[105.71,|1920],[1|
|00005bc0| 32 32 2e 37 37 35 2c 31 | 39 33 30 5d 2c 5b 31 33 |22.775,1|930],[13|
|00005bd0| 31 2e 36 36 39 2c 31 39 | 34 30 5d 2c 5b 31 35 31 |1.669,19|40],[151|
|00005be0| 2e 33 32 35 2c 31 39 35 | 30 5d 2c 5b 31 37 39 2e |.325,195|0],[179.|
|00005bf0| 33 32 33 2c 31 39 36 30 | 5d 2c 5b 32 30 33 2e 32 |323,1960|],[203.2|
|00005c00| 31 31 2c 31 39 37 30 5d | 5d 03 80 38 00 00 00 ee |11,1970]|]..8....|
|00005c10| 0e 00 00 b8 02 00 00 fa | 0e 00 00 00 04 55 73 65 |........|.....Use|
|00005c20| 72 00 00 00 00 00 00 f0 | bf 09 00 00 00 01 00 00 |r.......|........|
|00005c30| 00 43 47 41 55 53 53 5f | 4e 45 57 54 4f 4e 28 70 |.CGAUSS_|NEWTON(p|
|00005c40| 3d 63 2a 4e 4f 52 4d 41 | 4c 28 61 2b 62 2a 28 79 |=c*NORMA|L(a+b*(y|
|00005c50| 2d 31 37 39 30 29 29 2c | 5b 61 2c 62 2c 63 5d 2c |-1790)),|[a,b,c],|
|00005c60| 5b 2d 32 2e 34 2c 30 2e | 30 31 32 2c 34 30 30 5d |[-2.4,0.|012,400]|
|00005c70| 2c 70 6f 70 29 03 80 a0 | 00 00 00 06 0f 00 00 58 |,pop)...|.......X|
|00005c80| 03 00 00 da 10 00 00 01 | 08 53 69 6d 70 28 23 39 |........|.Simp(#9|
|00005c90| 29 62 10 58 39 b4 c8 01 | 40 0a 00 00 00 01 00 00 |)b.X9...|@.......|
|00005ca0| 00 ff 5f 02 5b 5b 22 47 | 61 75 73 73 5f 4e 65 77 |.._.[["G|auss_New|
|00005cb0| 74 6f 6e 20 4d 65 74 68 | 6f 64 22 5d 2c 5b 22 20 |ton Meth|od"],[" |
|00005cc0| 22 5d 2c 5b 22 43 6f 6e | 76 65 72 67 65 6e 63 65 |"],["Con|vergence|
|00005cd0| 20 63 72 69 74 65 72 69 | 61 20 6d 65 74 21 22 5d | criteri|a met!"]|
|00005ce0| 2c 5b 22 20 22 5d 2c 5b | 5b 5b 22 50 61 72 6d 22 |,[" "],[|[["Parm"|
|00005cf0| 2c 22 56 61 6c 75 65 22 | 2c 22 53 54 44 22 2c 22 |,"Value"|,"STD","|
|00005d00| 74 28 31 36 29 22 2c 22 | 50 72 6f 62 28 74 29 22 |t(16)","|Prob(t)"|
|00005d10| 5d 2c 5b 61 2c 2d 32 2e | 33 30 32 38 31 38 33 2c |],[a,-2.|3028183,|
|00005d20| 30 2e 30 33 32 38 34 35 | 38 36 31 2c 2d 37 30 2e |0.032845|861,-70.|
|00005d30| 31 30 39 38 34 37 2c 30 | 5d 2c 5b 62 2c 30 2e 30 |109847,0|],[b,0.0|
|00005d40| 31 32 36 32 38 31 34 33 | 2c 30 2e 30 30 30 39 34 |12628143|,0.00094|
|00005d50| 31 38 32 34 33 36 2c 31 | 33 2e 34 30 38 31 37 32 |182436,1|3.408172|
|00005d60| 2c 30 5d 2c 5b 63 2c 34 | 30 37 2e 31 30 36 30 35 |,0],[c,4|07.10605|
|00005d70| 2c 36 30 2e 36 36 33 31 | 39 37 2c 36 2e 37 31 30 |,60.6631|97,6.710|
|00005d80| 39 32 33 31 2c 35 2e 30 | 30 31 38 32 35 36 2a 31 |9231,5.0|018256*1|
|00005d90| 30 5e 28 2d 36 29 5d 5d | 5d 2c 5b 22 20 22 5d 2c |0^(-6)]]|],[" "],|
|00005da0| 5b 5b 5b 22 53 6f 75 72 | 63 65 22 2c 22 44 46 22 |[[["Sour|ce","DF"|
|00005db0| 2c 22 53 53 22 2c 22 4d | 53 22 2c 22 46 22 2c 22 |,"SS","M|S","F","|
|00005dc0| 50 72 6f 62 28 46 29 22 | 5d 2c 5b 22 52 65 67 22 |Prob(F)"|],["Reg"|
|00005dd0| 2c 33 2c 31 2e 36 34 32 | 32 37 38 39 2a 31 30 5e |,3,1.642|2789*10^|
|00005de0| 35 2c 35 2e 34 37 34 32 | 36 33 33 2a 31 30 5e 34 |5,5.4742|633*10^4|
|00005df0| 2c 34 39 33 38 2e 31 36 | 39 34 2c 32 2e 38 32 35 |,4938.16|94,2.825|
|00005e00| 36 34 32 36 2a 31 30 5e | 28 2d 31 31 29 5d 2c 5b |6426*10^|(-11)],[|
|00005e10| 22 45 72 72 6f 72 22 2c | 31 36 2c 31 37 37 2e 33 |"Error",|16,177.3|
|00005e20| 36 39 38 30 2c 31 31 2e | 30 38 35 36 31 32 2c 22 |6980,11.|085612,"|
|00005e30| 20 22 2c 22 20 22 5d 2c | 5b 22 55 6e 63 6f 72 72 | "," "],|["Uncorr|
|00005e40| 65 63 74 65 64 20 54 6f | 74 61 6c 22 2c 31 39 2c |ected To|tal",19,|
|00005e50| 31 2e 36 34 34 30 35 32 | 36 2a 31 30 5e 35 2c 22 |1.644052|6*10^5,"|
|00005e60| 20 22 2c 22 20 22 2c 22 | 20 22 5d 5d 5d 2c 5b 22 | "," ","| "]]],["|
|00005e70| 20 22 5d 2c 5b 5b 5b 22 | 53 45 22 2c 22 52 5e 32 | "],[[["|SE","R^2|
|00005e80| 22 2c 22 41 64 6a 52 5e | 32 22 5d 2c 5b 33 2e 33 |","AdjR^|2"],[3.3|
|00005e90| 32 39 35 30 36 33 2c 30 | 2e 39 39 38 39 32 31 31 |295063,0|.9989211|
|00005ea0| 34 2c 30 2e 39 39 38 37 | 31 38 38 35 5d 5d 5d 2c |4,0.9987|1885]]],|
|00005eb0| 5b 22 20 22 5d 2c 5b 70 | 3d 32 30 33 2e 35 35 33 |[" "],[p|=203.553|
|00005ec0| 30 32 2a 28 45 52 46 28 | 32 2e 30 37 37 30 30 30 |02*(ERF(|2.077000|
|00005ed0| 32 2a 31 30 5e 28 2d 31 | 32 29 2a 28 34 2e 32 39 |2*10^(-1|2)*(4.29|
|00005ee0| 39 32 30 33 30 2a 31 30 | 5e 39 2a 79 2d 38 2e 34 |92030*10|^9*y-8.4|
|00005ef0| 37 39 35 35 39 31 2a 31 | 30 5e 31 32 29 29 2b 31 |795591*1|0^12))+1|
|00005f00| 29 5d 5d 01 80 08 00 00 | 00 e6 10 00 00 cd 03 00 |)]].....|........|
|00005f10| 00 f6 10 00 00 00 ff 48 | 01 7b 5c 72 74 66 31 5c |.......H|.{\rtf1\|
|00005f20| 61 6e 73 69 5c 64 65 66 | 66 30 5c 64 65 66 74 61 |ansi\def|f0\defta|
|00005f30| 62 37 32 30 7b 5c 66 6f | 6e 74 74 62 6c 7b 5c 66 |b720{\fo|nttbl{\f|
|00005f40| 30 5c 66 73 77 69 73 73 | 20 4d 53 20 53 61 6e 73 |0\fswiss| MS Sans|
|00005f50| 20 53 65 72 69 66 3b 7d | 7b 5c 66 31 5c 66 64 65 | Serif;}|{\f1\fde|
|00005f60| 63 6f 72 5c 66 63 68 61 | 72 73 65 74 32 20 53 79 |cor\fcha|rset2 Sy|
|00005f70| 6d 62 6f 6c 3b 7d 7b 5c | 66 32 5c 66 73 77 69 73 |mbol;}{\|f2\fswis|
|00005f80| 73 5c 66 70 72 71 32 20 | 53 79 73 74 65 6d 3b 7d |s\fprq2 |System;}|
|00005f90| 7b 5c 66 33 5c 66 73 77 | 69 73 73 5c 66 70 72 71 |{\f3\fsw|iss\fprq|
|00005fa0| 32 20 41 72 69 61 6c 3b | 7d 7b 5c 66 34 5c 66 6d |2 Arial;|}{\f4\fm|
|00005fb0| 6f 64 65 72 6e 5c 66 63 | 68 61 72 73 65 74 32 20 |odern\fc|harset2 |
|00005fc0| 44 66 57 35 20 50 72 69 | 6e 74 65 72 3b 7d 7d 0d |DfW5 Pri|nter;}}.|
|00005fd0| 0a 7b 5c 63 6f 6c 6f 72 | 74 62 6c 5c 72 65 64 30 |.{\color|tbl\red0|
|00005fe0| 5c 67 72 65 65 6e 30 5c | 62 6c 75 65 30 3b 5c 72 |\green0\|blue0;\r|
|00005ff0| 65 64 32 35 35 5c 67 72 | 65 65 6e 30 5c 62 6c 75 |ed255\gr|een0\blu|
|00006000| 65 30 3b 5c 72 65 64 30 | 5c 67 72 65 65 6e 30 5c |e0;\red0|\green0\|
|00006010| 62 6c 75 65 32 35 35 3b | 7d 0d 0a 5c 64 65 66 6c |blue255;|}..\defl|
|00006020| 61 6e 67 31 30 33 33 5c | 70 61 72 64 5c 70 6c 61 |ang1033\|pard\pla|
|00006030| 69 6e 5c 66 33 5c 66 73 | 32 30 5c 63 66 30 20 54 |in\f3\fs|20\cf0 T|
|00006040| 68 65 20 69 74 65 72 61 | 74 69 6f 6e 20 6d 61 74 |he itera|tion mat|
|00006050| 72 69 78 20 69 73 3a 0d | 0a 5c 70 61 72 20 7d 0d |rix is:.|.\par }.|
|00006060| 0a 03 80 38 00 00 00 02 | 11 00 00 58 00 00 00 0e |...8....|...X....|
|00006070| 11 00 00 00 04 55 73 65 | 72 00 00 00 00 00 00 f0 |.....Use|r.......|
|00006080| bf 0b 00 00 00 01 00 00 | 00 04 69 74 65 72 03 80 |........|..iter..|
|00006090| 20 01 00 00 1a 11 00 00 | d8 02 00 00 b6 11 00 00 | .......|........|
|000060a0| 01 09 53 69 6d 70 28 23 | 31 31 29 00 00 00 00 00 |..Simp(#|11).....|
|000060b0| 00 00 00 0c 00 00 00 01 | 00 00 00 ff 1e 01 5b 5b |........|......[[|
|000060c0| 22 49 74 65 72 22 2c 61 | 2c 62 2c 63 2c 22 53 53 |"Iter",a|,b,c,"SS|
|000060d0| 45 22 5d 2c 5b 30 2c 2d | 32 2e 34 2c 30 2e 30 31 |E"],[0,-|2.4,0.01|
|000060e0| 32 2c 34 30 30 2c 37 31 | 37 34 2e 35 39 30 38 5d |2,400,71|74.5908]|
|000060f0| 2c 5b 31 2c 2d 32 2e 32 | 37 32 35 37 34 33 2c 30 |,[1,-2.2|725743,0|
|00006100| 2e 30 31 32 36 30 38 36 | 31 34 2c 34 30 30 2e 32 |.0126086|14,400.2|
|00006110| 31 39 36 32 2c 32 30 38 | 2e 39 39 32 36 34 5d 2c |1962,208|.99264],|
|00006120| 5b 32 2c 2d 32 2e 33 30 | 32 34 38 38 32 2c 30 2e |[2,-2.30|24882,0.|
|00006130| 30 31 32 36 35 39 32 32 | 33 2c 34 30 34 2e 39 35 |01265922|3,404.95|
|00006140| 33 38 30 2c 31 37 37 2e | 33 38 36 37 31 5d 2c 5b |380,177.|38671],[|
|00006150| 33 2c 2d 32 2e 33 30 32 | 37 36 32 34 2c 30 2e 30 |3,-2.302|7624,0.0|
|00006160| 31 32 36 32 38 33 36 34 | 2c 34 30 37 2e 30 35 39 |12628364|,407.059|
|00006170| 30 39 2c 31 37 37 2e 33 | 36 39 39 31 5d 2c 5b 34 |09,177.3|6991],[4|
|00006180| 2c 2d 32 2e 33 30 32 38 | 31 38 31 2c 30 2e 30 31 |,-2.3028|181,0.01|
|00006190| 32 36 32 38 32 31 32 2c | 34 30 37 2e 31 30 31 36 |2628212,|407.1016|
|000061a0| 31 2c 31 37 37 2e 33 36 | 39 38 30 5d 2c 5b 35 2c |1,177.36|980],[5,|
|000061b0| 2d 32 2e 33 30 32 38 31 | 38 33 2c 30 2e 30 31 32 |-2.30281|83,0.012|
|000061c0| 36 32 38 31 34 33 2c 34 | 30 37 2e 31 30 36 30 35 |628143,4|07.10605|
|000061d0| 2c 31 37 37 2e 33 36 39 | 38 30 5d 5d 01 80 08 00 |,177.369|80]]....|
|000061e0| 00 00 c2 11 00 00 cd 03 | 00 00 d2 11 00 00 00 ff |........|........|
|000061f0| be 01 7b 5c 72 74 66 31 | 5c 61 6e 73 69 5c 64 65 |..{\rtf1|\ansi\de|
|00006200| 66 66 30 5c 64 65 66 74 | 61 62 37 32 30 7b 5c 66 |ff0\deft|ab720{\f|
|00006210| 6f 6e 74 74 62 6c 7b 5c | 66 30 5c 66 73 77 69 73 |onttbl{\|f0\fswis|
|00006220| 73 20 4d 53 20 53 61 6e | 73 20 53 65 72 69 66 3b |s MS San|s Serif;|
|00006230| 7d 7b 5c 66 31 5c 66 64 | 65 63 6f 72 5c 66 63 68 |}{\f1\fd|ecor\fch|
|00006240| 61 72 73 65 74 32 20 53 | 79 6d 62 6f 6c 3b 7d 7b |arset2 S|ymbol;}{|
|00006250| 5c 66 32 5c 66 73 77 69 | 73 73 5c 66 70 72 71 32 |\f2\fswi|ss\fprq2|
|00006260| 20 53 79 73 74 65 6d 3b | 7d 7b 5c 66 33 5c 66 73 | System;|}{\f3\fs|
|00006270| 77 69 73 73 5c 66 70 72 | 71 32 20 41 72 69 61 6c |wiss\fpr|q2 Arial|
|00006280| 3b 7d 7b 5c 66 34 5c 66 | 6d 6f 64 65 72 6e 5c 66 |;}{\f4\f|modern\f|
|00006290| 63 68 61 72 73 65 74 32 | 20 44 66 57 35 20 50 72 |charset2| DfW5 Pr|
|000062a0| 69 6e 74 65 72 3b 7d 7d | 0d 0a 7b 5c 63 6f 6c 6f |inter;}}|..{\colo|
|000062b0| 72 74 62 6c 5c 72 65 64 | 30 5c 67 72 65 65 6e 30 |rtbl\red|0\green0|
|000062c0| 5c 62 6c 75 65 30 3b 5c | 72 65 64 32 35 35 5c 67 |\blue0;\|red255\g|
|000062d0| 72 65 65 6e 30 5c 62 6c | 75 65 30 3b 5c 72 65 64 |reen0\bl|ue0;\red|
|000062e0| 30 5c 67 72 65 65 6e 30 | 5c 62 6c 75 65 32 35 35 |0\green0|\blue255|
|000062f0| 3b 7d 0d 0a 5c 64 65 66 | 6c 61 6e 67 31 30 33 33 |;}..\def|lang1033|
|00006300| 5c 70 61 72 64 5c 70 6c | 61 69 6e 5c 66 33 5c 66 |\pard\pl|ain\f3\f|
|00006310| 73 32 30 5c 63 66 30 20 | 54 68 65 20 65 78 70 65 |s20\cf0 |The expe|
|00006320| 63 74 65 64 20 70 6f 70 | 75 6c 61 74 69 6f 6e 20 |cted pop|ulation |
|00006330| 66 6f 72 20 74 68 65 20 | 79 65 61 72 20 32 30 30 |for the |year 200|
|00006340| 31 2c 20 61 6c 6f 6e 67 | 20 77 69 74 68 20 74 68 |1, along| with th|
|00006350| 65 20 73 74 61 6e 64 61 | 72 64 20 65 72 72 6f 72 |e standa|rd error|
|00006360| 73 20 66 6f 72 20 74 68 | 65 20 6d 65 61 6e 20 61 |s for th|e mean a|
|00006370| 6e 64 20 66 6f 72 65 63 | 61 73 74 20 76 61 6c 75 |nd forec|ast valu|
|00006380| 65 73 20 61 6e 64 20 74 | 68 65 69 72 20 63 6f 6e |es and t|heir con|
|00006390| 66 69 64 65 6e 63 65 20 | 69 6e 74 65 72 76 61 6c |fidence |interval|
|000063a0| 73 2c 20 69 73 3a 0d 0a | 5c 70 61 72 20 7d 0d 0a |s, is:..|\par }..|
|000063b0| 03 80 38 00 00 00 de 11 | 00 00 f8 00 00 00 ea 11 |..8.....|........|
|000063c0| 00 00 00 04 55 73 65 72 | 00 00 00 00 00 00 f0 bf |....User|........|
|000063d0| 0d 00 00 00 01 00 00 00 | 18 50 72 65 64 69 63 74 |........|.Predict|
|000063e0| 65 64 5f 56 61 6c 75 65 | 73 28 5b 32 30 30 31 5d |ed_Value|s([2001]|
|000063f0| 29 03 80 50 01 00 00 f6 | 11 00 00 a8 02 00 00 62 |)..P....|.......b|
+--------+-------------------------+-------------------------+--------+--------+