/ Stars of Shareware: Programmierung / SOURCE.mdf / programm / msdos / c / yamp / app2.tex

< prev

next >

MacOS 8.1

Win98

DOS

browse contents

view JSON data

view as text

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 5c 63 68 61 70 74 65 72 | 7b 74 65 73 74 72 65 67 |\chapter|{testreg|
|00000010| 2e 63 70 70 7d 0d 0a 54 | 68 69 73 20 69 73 20 74 |.cpp}..T|his is t|
|00000020| 68 65 20 74 65 73 74 20 | 70 72 6f 67 72 61 6d 20 |he test |program |
|00000030| 66 6f 72 20 74 68 65 20 | 72 65 67 72 65 73 73 69 |for the |regressi|
|00000040| 6f 6e 20 64 65 63 6f 6d | 70 6f 73 69 74 69 6f 6e |on decom|position|
|00000050| 73 20 0d 0a 69 6e 20 74 | 68 65 20 56 4d 61 74 72 |s ..in t|he VMatr|
|00000060| 69 78 20 73 79 73 74 65 | 6d 2e 20 54 68 65 20 73 |ix syste|m. The s|
|00000070| 69 6e 67 75 6c 61 72 20 | 76 61 6c 75 65 64 20 64 |ingular |valued d|
|00000080| 65 63 6f 6d 70 6f 73 69 | 74 69 6f 6e 20 0d 0a 69 |ecomposi|tion ..i|
|00000090| 6e 76 65 72 74 73 20 74 | 68 65 20 31 31 78 31 31 |nverts t|he 11x11|
|000000a0| 20 48 69 6c 62 65 72 74 | 20 6d 61 74 72 69 78 20 | Hilbert| matrix |
|000000b0| 77 69 74 68 20 61 20 6d | 61 78 69 6d 75 6d 20 65 |with a m|aximum e|
|000000c0| 72 72 6f 72 20 6f 66 20 | 30 2e 30 30 31 2e 0d 0a |rror of |0.001...|
|000000d0| 0d 0a 5c 62 65 67 69 6e | 7b 76 65 72 62 61 74 69 |..\begin|{verbati|
|000000e0| 6d 7d 0d 0a 0d 0a 0d 0a | 23 69 6e 63 6c 75 64 65 |m}......|#include|
|000000f0| 20 22 76 69 72 74 2e 68 | 22 0d 0a 0d 0a 2f 2f 72 | "virt.h|"....//r|
|00000100| 65 71 75 69 72 65 64 20 | 67 6c 6f 62 61 6c 20 64 |equired |global d|
|00000110| 65 63 6c 61 72 61 74 69 | 6f 6e 20 66 6f 72 20 74 |eclarati|on for t|
|00000120| 68 65 0d 0a 2f 2f 20 20 | 6d 61 74 72 69 78 20 73 |he..//  |matrix s|
|00000130| 74 61 63 6b 20 6f 62 6a | 65 63 74 0d 0a 0d 0a 75 |tack obj|ect....u|
|00000140| 6e 73 69 67 6e 65 64 20 | 69 6e 74 20 5f 73 74 6b |nsigned |int _stk|
|00000150| 6c 65 6e 20 3d 20 53 54 | 41 43 4b 4c 45 4e 47 54 |len = ST|ACKLENGT|
|00000160| 48 3b 0d 0a 4d 53 74 61 | 63 6b 20 2a 44 69 73 70 |H;..MSta|ck *Disp|
|00000170| 61 74 63 68 20 3d 20 6e | 65 77 20 4d 53 74 61 63 |atch = n|ew MStac|
|00000180| 6b 3b 0d 0a 0d 0a 56 4d | 61 74 72 69 78 26 20 45 |k;....VM|atrix& E|
|00000190| 78 61 63 74 48 69 6c 62 | 49 6e 76 28 20 69 6e 74 |xactHilb|Inv( int|
|000001a0| 20 6e 20 29 0d 0a 7b 0d | 0a 20 20 20 2f 2f 20 63 | n )..{.|.   // c|
|000001b0| 6f 6e 73 74 72 75 63 74 | 20 65 78 61 63 74 20 68 |onstruct| exact h|
|000001c0| 69 6c 62 65 72 74 20 6d | 61 74 72 69 78 20 69 6e |ilbert m|atrix in|
|000001d0| 76 65 72 73 65 0d 0a 20 | 20 20 44 69 73 70 61 74 |verse.. |  Dispat|
|000001e0| 63 68 2d 3e 49 6e 63 6c | 65 76 65 6c 28 29 3b 0d |ch->Incl|evel();.|
|000001f0| 0a 20 20 20 56 4d 61 74 | 72 69 78 20 48 69 28 22 |.   VMat|rix Hi("|
|00000200| 48 69 6c 62 65 72 74 20 | 49 6e 76 65 72 73 65 22 |Hilbert |Inverse"|
|00000210| 2c 6e 2c 6e 29 3b 0d 0a | 0d 0a 20 20 20 64 6f 75 |,n,n);..|..   dou|
|00000220| 62 6c 65 20 64 69 61 67 | 2c 20 64 6e 20 3d 20 6e |ble diag|, dn = n|
|00000230| 3b 0d 0a 20 20 20 66 6f | 72 28 20 69 6e 74 20 69 |;..   fo|r( int i|
|00000240| 3d 31 3b 20 69 3c 3d 6e | 3b 20 69 2b 2b 29 7b 0d |=1; i<=n|; i++){.|
|00000250| 0a 09 20 20 64 6f 75 62 | 6c 65 20 69 6d 31 20 3d |..  doub|le im1 =|
|00000260| 20 28 64 6f 75 62 6c 65 | 29 20 28 69 2d 31 2e 30 | (double|) (i-1.0|
|00000270| 29 3b 0d 0a 09 20 20 64 | 69 61 67 20 3d 20 28 20 |);...  d|iag = ( |
|00000280| 69 20 3d 3d 20 31 20 29 | 20 3f 20 64 6e 20 3a 20 |i == 1 )| ? dn : |
|00000290| 28 28 64 6e 2d 69 6d 31 | 29 2a 64 69 61 67 2a 28 |((dn-im1|)*diag*(|
|000002a0| 64 6e 2b 69 6d 31 29 29 | 2f 28 69 6d 31 2a 69 6d |dn+im1))|/(im1*im|
|000002b0| 31 29 3b 0d 0a 0d 0a 09 | 20 20 6c 6f 6e 67 20 64 |1);.....|  long d|
|000002c0| 6f 75 62 6c 65 20 72 65 | 73 74 20 3d 20 64 69 61 |ouble re|st = dia|
|000002d0| 67 2a 64 69 61 67 3b 0d | 0a 09 20 20 48 69 2e 4d |g*diag;.|..  Hi.M|
|000002e0| 28 69 2c 69 29 20 3d 20 | 72 65 73 74 2f 28 32 2e |(i,i) = |rest/(2.|
|000002f0| 30 2a 69 6d 31 2b 31 2e | 30 29 3b 0d 0a 09 20 20 |0*im1+1.|0);...  |
|00000300| 66 6f 72 28 20 69 6e 74 | 20 6a 3d 69 2b 31 3b 20 |for( int| j=i+1; |
|00000310| 6a 3c 3d 6e 3b 20 6a 2b | 2b 20 29 7b 0d 0a 09 09 |j<=n; j+|+ ){....|
|00000320| 20 64 6f 75 62 6c 65 20 | 6a 6d 31 20 3d 20 28 64 | double |jm1 = (d|
|00000330| 6f 75 62 6c 65 29 20 28 | 6a 2d 31 29 3b 0d 0a 09 |ouble) (|j-1);...|
|00000340| 09 20 72 65 73 74 20 3d | 20 2d 28 28 64 6e 2d 6a |. rest =| -((dn-j|
|00000350| 6d 31 29 2a 72 65 73 74 | 2a 28 64 6e 2b 6a 6d 31 |m1)*rest|*(dn+jm1|
|00000360| 29 29 2f 28 6a 6d 31 2a | 6a 6d 31 29 3b 0d 0a 09 |))/(jm1*|jm1);...|
|00000370| 09 20 48 69 2e 4d 28 69 | 2c 6a 29 20 3d 20 72 65 |. Hi.M(i|,j) = re|
|00000380| 73 74 2f 28 69 6d 31 2b | 6a 6d 31 2b 31 2e 30 29 |st/(im1+|jm1+1.0)|
|00000390| 3b 0d 0a 09 09 20 48 69 | 2e 4d 28 6a 2c 69 29 20 |;.... Hi|.M(j,i) |
|000003a0| 3d 20 48 69 2e 6d 28 69 | 2c 6a 29 3b 0d 0a 09 20 |= Hi.m(i|,j);... |
|000003b0| 20 7d 0d 0a 20 20 20 7d | 0d 0a 20 20 20 44 69 73 | }..   }|..   Dis|
|000003c0| 70 61 74 63 68 2d 3e 50 | 75 73 68 28 20 48 69 20 |patch->P|ush( Hi |
|000003d0| 29 3b 0d 0a 20 20 20 72 | 65 74 75 72 6e 20 44 69 |);..   r|eturn Di|
|000003e0| 73 70 61 74 63 68 2d 3e | 44 65 63 52 65 74 75 72 |spatch->|DecRetur|
|000003f0| 6e 28 29 3b 0d 0a 7d 0d | 0a 0d 0a 0d 0a 76 6f 69 |n();..}.|.....voi|
|00000400| 64 20 54 65 73 74 48 69 | 6c 62 28 20 69 6e 74 20 |d TestHi|lb( int |
|00000410| 68 69 6c 62 5f 69 6e 64 | 20 29 0d 0a 2f 2f 20 74 |hilb_ind| )..// t|
|00000420| 65 73 74 20 73 76 64 20 | 6f 6e 20 68 69 6c 62 65 |est svd |on hilbe|
|00000430| 72 74 20 6d 61 74 72 69 | 78 0d 0a 7b 0d 0a 20 20 |rt matri|x..{..  |
|00000440| 44 69 73 70 61 74 63 68 | 2d 3e 49 6e 63 6c 65 76 |Dispatch|->Inclev|
|00000450| 65 6c 28 29 3b 0d 0a 20 | 20 56 4d 61 74 72 69 78 |el();.. | VMatrix|
|00000460| 20 68 69 6c 62 20 3d 20 | 4b 72 6f 6e 28 20 46 69 | hilb = |Kron( Fi|
|00000470| 6c 6c 28 31 2c 68 69 6c | 62 5f 69 6e 64 2c 31 29 |ll(1,hil|b_ind,1)|
|00000480| 2c 20 49 6e 64 65 78 28 | 31 2c 20 68 69 6c 62 5f |, Index(|1, hilb_|
|00000490| 69 6e 64 29 29 3b 0d 0a | 20 20 2f 2f 20 68 28 69 |ind));..|  // h(i|
|000004a0| 2c 6a 29 20 3d 20 31 2e | 30 2f 28 69 2b 6a 2d 31 |,j) = 1.|0/(i+j-1|
|000004b0| 29 0d 0a 20 20 68 69 6c | 62 20 3d 20 31 2e 30 20 |)..  hil|b = 1.0 |
|000004c0| 2f 20 28 20 28 68 69 6c | 62 2b 54 72 61 6e 28 68 |/ ( (hil|b+Tran(h|
|000004d0| 69 6c 62 29 29 2d 31 2e | 30 29 3b 0d 0a 20 20 68 |ilb))-1.|0);..  h|
|000004e0| 69 6c 62 2e 4e 61 6d 65 | 69 74 28 22 48 69 6c 62 |ilb.Name|it("Hilb|
|000004f0| 65 72 74 20 4d 61 74 72 | 69 78 22 29 3b 0d 0a 20 |ert Matr|ix");.. |
|00000500| 20 68 69 6c 62 2e 44 69 | 73 70 6c 61 79 4d 61 74 | hilb.Di|splayMat|
|00000510| 28 29 3b 0d 0a 20 20 2f | 2f 20 75 73 65 20 66 6f |();..  /|/ use fo|
|00000520| 72 6d 75 6c 61 2e 20 4c | 6f 6f 73 65 73 20 61 63 |rmula. L|ooses ac|
|00000530| 63 75 72 61 63 79 20 6f | 66 20 30 2e 30 30 31 20 |curacy o|f 0.001 |
|00000540| 61 74 20 6e 3e 3d 31 31 | 0d 0a 20 20 28 45 78 61 |at n>=11|..  (Exa|
|00000550| 63 74 48 69 6c 62 49 6e | 76 28 20 68 69 6c 62 5f |ctHilbIn|v( hilb_|
|00000560| 69 6e 64 20 29 2a 68 69 | 6c 62 20 29 2e 44 69 73 |ind )*hi|lb ).Dis|
|00000570| 70 6c 61 79 4d 61 74 28 | 29 3b 0d 0a 0d 0a 20 20 |playMat(|);....  |
|00000580| 2f 2f 20 75 73 65 20 73 | 77 65 65 70 0d 0a 20 20 |// use s|weep..  |
|00000590| 56 4d 61 74 72 69 78 20 | 47 61 75 73 73 20 3d 20 |VMatrix |Gauss = |
|000005a0| 49 6e 76 28 68 69 6c 62 | 29 2a 68 69 6c 62 3b 0d |Inv(hilb|)*hilb;.|
|000005b0| 0a 20 20 47 61 75 73 73 | 2e 4e 61 6d 65 69 74 28 |.  Gauss|.Nameit(|
|000005c0| 22 47 61 75 73 73 69 61 | 6e 20 65 6c 69 6d 69 6e |"Gaussia|n elimin|
|000005d0| 61 74 69 6f 6e 22 29 3b | 0d 0a 20 20 47 61 75 73 |ation");|..  Gaus|
|000005e0| 73 2e 44 69 73 70 6c 61 | 79 4d 61 74 28 29 3b 0d |s.Displa|yMat();.|
|000005f0| 0a 0d 0a 20 20 2f 2f 20 | 75 73 65 20 73 69 6e 67 |...  // |use sing|
|00000600| 75 6c 61 72 20 76 61 6c | 75 65 64 20 64 65 63 6f |ular val|ued deco|
|00000610| 6d 70 6f 73 69 74 69 6f | 6e 0d 0a 20 20 56 4d 61 |mpositio|n..  VMa|
|00000620| 74 72 69 78 20 53 2c 20 | 56 2c 20 44 3b 0d 0a 20 |trix S, |V, D;.. |
|00000630| 20 53 76 64 28 20 68 69 | 6c 62 2c 20 53 2c 20 56 | Svd( hi|lb, S, V|
|00000640| 2c 20 44 29 3b 0d 0a 20 | 20 56 4d 61 74 72 69 78 |, D);.. | VMatrix|
|00000650| 20 74 20 3d 20 46 69 6c | 6c 28 56 2e 72 2c 20 56 | t = Fil|l(V.r, V|
|00000660| 2e 72 2c 20 30 2e 30 29 | 3b 0d 0a 20 20 66 6f 72 |.r, 0.0)|;..  for|
|00000670| 20 28 69 6e 74 20 69 20 | 3d 20 31 3b 20 69 20 3c | (int i |= 1; i <|
|00000680| 3d 20 56 2e 72 3b 20 69 | 2b 2b 29 20 7b 0d 0a 09 |= V.r; i|++) {...|
|00000690| 64 6f 75 62 6c 65 20 74 | 74 20 3d 20 56 2e 6d 28 |double t|t = V.m(|
|000006a0| 69 2c 20 31 29 3b 0d 0a | 09 74 2e 4d 28 69 2c 20 |i, 1);..|.t.M(i, |
|000006b0| 69 29 20 3d 20 31 2e 30 | 20 2f 20 74 74 3b 0d 0a |i) = 1.0| / tt;..|
|000006c0| 20 20 7d 0d 0a 0d 0a 20 | 20 56 4d 61 74 72 69 78 |  }.... | VMatrix|
|000006d0| 20 47 67 61 75 73 73 20 | 3d 20 28 44 2a 74 2a 54 | Ggauss |= (D*t*T|
|000006e0| 72 61 6e 28 53 29 29 2a | 68 69 6c 62 3b 0d 0a 20 |ran(S))*|hilb;.. |
|000006f0| 20 47 67 61 75 73 73 2e | 4e 61 6d 65 69 74 28 22 | Ggauss.|Nameit("|
|00000700| 47 65 6e 65 72 61 6c 69 | 7a 65 64 20 69 6e 76 65 |Generali|zed inve|
|00000710| 72 73 69 6f 6e 22 29 3b | 0d 0a 20 20 47 67 61 75 |rsion");|..  Ggau|
|00000720| 73 73 2e 44 69 73 70 6c | 61 79 4d 61 74 28 29 3b |ss.Displ|ayMat();|
|00000730| 0d 0a 0d 0a 20 20 44 69 | 73 70 61 74 63 68 2d 3e |....  Di|spatch->|
|00000740| 44 65 63 6c 65 76 65 6c | 28 29 3b 0d 0a 7d 0d 0a |Declevel|();..}..|
|00000750| 0d 0a 2f 2f 2f 2f 2f 2f | 2f 2f 2f 2f 2f 2f 2f 2f |..//////|////////|
|00000760| 2f 2f 2f 2f 2f 2f 2f 20 | 6e 6f 77 20 74 65 73 74 |/////// |now test|
|00000770| 20 72 65 67 72 65 73 73 | 69 6f 6e 73 20 6f 6e 20 | regress|ions on |
|00000780| 73 61 6d 70 6c 65 20 64 | 61 74 61 0d 0a 0d 0a 56 |sample d|ata....V|
|00000790| 4d 61 74 72 69 78 20 26 | 67 65 74 78 28 69 6e 74 |Matrix &|getx(int|
|000007a0| 20 4e 29 0d 0a 2f 2f 20 | 63 72 65 61 74 65 20 61 | N)..// |create a|
|000007b0| 6e 20 78 20 6d 61 74 72 | 69 78 0d 0a 7b 0d 0a 09 |n x matr|ix..{...|
|000007c0| 44 69 73 70 61 74 63 68 | 2d 3e 49 6e 63 6c 65 76 |Dispatch|->Inclev|
|000007d0| 65 6c 28 29 3b 0d 0a 0d | 0a 09 56 4d 61 74 72 69 |el();...|..VMatri|
|000007e0| 78 20 78 20 3d 20 49 6e | 64 65 78 28 20 31 2c 20 |x x = In|dex( 1, |
|000007f0| 4e 20 29 20 2d 20 4e 2f | 32 3b 0d 0a 09 56 4d 61 |N ) - N/|2;...VMa|
|00000800| 74 72 69 78 20 63 31 20 | 3d 20 46 69 6c 6c 28 4e |trix c1 |= Fill(N|
|00000810| 2c 20 31 2c 20 31 2e 30 | 29 3b 0d 0a 09 56 4d 61 |, 1, 1.0|);...VMa|
|00000820| 74 72 69 78 20 78 32 20 | 3d 20 78 20 25 20 78 3b |trix x2 |= x % x;|
|00000830| 0d 0a 09 78 20 3d 20 43 | 68 28 63 31 2c 20 43 68 |...x = C|h(c1, Ch|
|00000840| 28 78 2c 20 78 32 29 29 | 3b 0d 0a 0d 0a 09 2f 2f |(x, x2))|;.....//|
|00000850| 20 70 75 73 68 20 78 20 | 6f 6e 74 6f 20 74 68 65 | push x |onto the|
|00000860| 20 73 74 61 63 6b 0d 0a | 09 44 69 73 70 61 74 63 | stack..|.Dispatc|
|00000870| 68 2d 3e 50 75 73 68 28 | 78 29 3b 0d 0a 09 2f 2f |h->Push(|x);...//|
|00000880| 20 64 65 63 72 65 6d 65 | 6e 74 20 74 68 65 20 73 | decreme|nt the s|
|00000890| 75 62 72 6f 75 74 69 6e | 65 20 6e 65 73 74 69 6e |ubroutin|e nestin|
|000008a0| 67 20 6c 65 76 65 6c 0d | 0a 09 2f 2f 20 61 6e 64 |g level.|..// and|
|000008b0| 20 72 65 74 75 72 6e 20 | 74 68 65 20 73 74 61 63 | return |the stac|
|000008c0| 6b 20 74 6f 70 0d 0a 09 | 72 65 74 75 72 6e 20 44 |k top...|return D|
|000008d0| 69 73 70 61 74 63 68 2d | 3e 44 65 63 52 65 74 75 |ispatch-|>DecRetu|
|000008e0| 72 6e 28 29 3b 0d 0a 7d | 0d 0a 0d 0a 56 4d 61 74 |rn();..}|....VMat|
|000008f0| 72 69 78 20 26 67 65 74 | 79 28 56 4d 61 74 72 69 |rix &get|y(VMatri|
|00000900| 78 20 26 78 2c 20 56 4d | 61 74 72 69 78 20 26 62 |x &x, VM|atrix &b|
|00000910| 65 74 61 29 0d 0a 2f 2f | 20 63 72 65 61 74 65 20 |eta)..//| create |
|00000920| 61 20 79 20 6d 61 74 72 | 69 78 0d 0a 7b 0d 0a 09 |a y matr|ix..{...|
|00000930| 44 69 73 70 61 74 63 68 | 2d 3e 49 6e 63 6c 65 76 |Dispatch|->Inclev|
|00000940| 65 6c 28 29 3b 0d 0a 0d | 0a 09 56 4d 61 74 72 69 |el();...|..VMatri|
|00000950| 78 20 79 20 3d 20 78 2a | 62 65 74 61 3b 0d 0a 09 |x y = x*|beta;...|
|00000960| 73 72 61 6e 64 28 31 32 | 33 29 3b 0d 0a 09 66 6f |srand(12|3);...fo|
|00000970| 72 20 28 69 6e 74 20 69 | 20 3d 20 31 3b 20 69 20 |r (int i| = 1; i |
|00000980| 3c 3d 20 79 2e 72 3b 20 | 69 2b 2b 29 20 7b 0d 0a |<= y.r; |i++) {..|
|00000990| 09 09 2f 2f 20 75 73 65 | 20 73 75 6d 20 6f 66 20 |..// use| sum of |
|000009a0| 33 20 75 6e 69 66 6f 72 | 6d 73 20 66 6f 72 20 61 |3 unifor|ms for a|
|000009b0| 6e 20 61 70 70 72 6f 78 | 69 6d 61 74 65 0d 0a 09 |n approx|imate...|
|000009c0| 09 2f 2f 20 6e 6f 72 6d | 61 6c 20 72 61 6e 64 6f |.// norm|al rando|
|000009d0| 6d 20 76 61 72 69 61 62 | 6c 65 0d 0a 09 09 69 6e |m variab|le....in|
|000009e0| 74 20 75 20 3d 20 72 61 | 6e 64 6f 6d 28 31 30 30 |t u = ra|ndom(100|
|000009f0| 29 20 2b 20 72 61 6e 64 | 6f 6d 28 31 30 30 29 20 |) + rand|om(100) |
|00000a00| 2b 20 72 61 6e 64 6f 6d | 28 31 30 30 29 20 2b 20 |+ random|(100) + |
|00000a10| 33 3b 0d 0a 09 09 79 2e | 4d 28 69 2c 20 31 29 20 |3;....y.|M(i, 1) |
|00000a20| 3d 20 79 2e 6d 28 69 2c | 20 31 29 20 2b 20 28 28 |= y.m(i,| 1) + ((|
|00000a30| 64 6f 75 62 6c 65 29 20 | 28 75 20 2d 20 31 35 30 |double) |(u - 150|
|00000a40| 29 29 20 2f 20 33 30 30 | 2e 30 3b 0d 0a 09 7d 0d |)) / 300|.0;...}.|
|00000a50| 0a 09 44 69 73 70 61 74 | 63 68 2d 3e 50 75 73 68 |..Dispat|ch->Push|
|00000a60| 28 79 29 3b 0d 0a 09 72 | 65 74 75 72 6e 20 44 69 |(y);...r|eturn Di|
|00000a70| 73 70 61 74 63 68 2d 3e | 44 65 63 52 65 74 75 72 |spatch->|DecRetur|
|00000a80| 6e 28 29 3b 0d 0a 7d 0d | 0a 0d 0a 56 4d 61 74 72 |n();..}.|...VMatr|
|00000a90| 69 78 20 26 72 65 67 72 | 65 73 73 69 6f 6e 28 56 |ix &regr|ession(V|
|00000aa0| 4d 61 74 72 69 78 26 20 | 78 2c 20 56 4d 61 74 72 |Matrix& |x, VMatr|
|00000ab0| 69 78 26 20 79 29 0d 0a | 2f 2f 20 64 6f 20 61 20 |ix& y)..|// do a |
|00000ac0| 6d 75 6c 74 69 70 6c 65 | 20 6c 69 6e 65 61 72 20 |multiple| linear |
|00000ad0| 72 65 67 72 65 73 73 69 | 6f 6e 0d 0a 7b 0d 0a 09 |regressi|on..{...|
|00000ae0| 44 69 73 70 61 74 63 68 | 2d 3e 49 6e 63 6c 65 76 |Dispatch|->Inclev|
|00000af0| 65 6c 28 29 3b 0d 0a 0d | 0a 09 69 6e 74 20 4e 20 |el();...|..int N |
|00000b00| 3d 20 78 2e 72 2c 20 70 | 20 3d 20 78 2e 63 3b 0d |= x.r, p| = x.c;.|
|00000b10| 0a 09 2f 2f 20 73 6f 6c | 76 65 20 66 6f 72 20 72 |..// sol|ve for r|
|00000b20| 65 67 72 65 73 73 69 6f | 6e 20 70 61 72 61 6d 65 |egressio|n parame|
|00000b30| 74 65 72 73 20 75 73 69 | 6e 67 20 73 77 65 65 70 |ters usi|ng sweep|
|00000b40| 0d 0a 09 56 4d 61 74 72 | 69 78 20 79 78 20 3d 20 |...VMatr|ix yx = |
|00000b50| 43 68 28 79 2c 20 78 29 | 3b 0d 0a 09 56 4d 61 74 |Ch(y, x)|;...VMat|
|00000b60| 72 69 78 20 72 65 67 20 | 3d 20 53 77 65 65 70 28 |rix reg |= Sweep(|
|00000b70| 32 2c 20 70 20 2b 20 31 | 2c 20 54 72 61 6e 28 79 |2, p + 1|, Tran(y|
|00000b80| 78 29 20 2a 79 78 29 3b | 0d 0a 09 2f 2f 20 63 61 |x) *yx);|...// ca|
|00000b90| 6c 63 75 6c 61 74 65 20 | 6d 65 61 6e 20 73 71 75 |lculate |mean squ|
|00000ba0| 61 72 65 20 65 72 72 6f | 72 0d 0a 09 72 65 67 2e |are erro|r...reg.|
|00000bb0| 4d 28 31 2c 20 31 29 20 | 3d 20 72 65 67 2e 6d 28 |M(1, 1) |= reg.m(|
|00000bc0| 31 2c 20 31 29 20 2f 20 | 28 28 64 6f 75 62 6c 65 |1, 1) / |((double|
|00000bd0| 29 20 28 4e 20 2d 20 70 | 29 29 3b 0d 0a 09 72 65 |) (N - p|));...re|
|00000be0| 67 2e 44 69 73 70 6c 61 | 79 4d 61 74 28 29 3b 0d |g.Displa|yMat();.|
|00000bf0| 0a 0d 0a 0d 0a 09 2f 2f | 20 73 6f 6c 76 65 20 72 |......//| solve r|
|00000c00| 65 67 72 65 73 73 69 6f | 6e 20 75 73 69 6e 67 20 |egressio|n using |
|00000c10| 6e 6f 72 6d 61 6c 20 65 | 71 75 61 74 69 6f 6e 73 |normal e|quations|
|00000c20| 0d 0a 09 56 4d 61 74 72 | 69 78 20 62 65 74 61 68 |...VMatr|ix betah|
|00000c30| 61 74 20 3d 20 49 6e 76 | 28 54 72 61 6e 28 78 29 |at = Inv|(Tran(x)|
|00000c40| 20 2a 78 29 20 2a 54 72 | 61 6e 28 78 29 20 2a 79 | *x) *Tr|an(x) *y|
|00000c50| 3b 0d 0a 0d 0a 09 44 69 | 73 70 61 74 63 68 2d 3e |;.....Di|spatch->|
|00000c60| 50 75 73 68 28 62 65 74 | 61 68 61 74 29 3b 0d 0a |Push(bet|ahat);..|
|00000c70| 09 72 65 74 75 72 6e 20 | 44 69 73 70 61 74 63 68 |.return |Dispatch|
|00000c80| 2d 3e 44 65 63 52 65 74 | 75 72 6e 28 29 3b 0d 0a |->DecRet|urn();..|
|00000c90| 7d 0d 0a 56 4d 61 74 72 | 69 78 20 26 51 52 72 65 |}..VMatr|ix &QRre|
|00000ca0| 67 72 65 73 73 69 6f 6e | 28 56 4d 61 74 72 69 78 |gression|(VMatrix|
|00000cb0| 20 26 78 2c 20 56 4d 61 | 74 72 69 78 20 26 79 29 | &x, VMa|trix &y)|
|00000cc0| 0d 0a 7b 0d 0a 09 2f 2f | 20 75 73 65 20 51 52 20 |..{...//| use QR |
|00000cd0| 64 65 63 6f 6d 70 6f 73 | 69 74 69 6f 6e 20 74 6f |decompos|ition to|
|00000ce0| 20 73 6f 6c 76 65 20 72 | 65 67 72 65 73 73 69 6f | solve r|egressio|
|00000cf0| 6e 0d 0a 0d 0a 09 44 69 | 73 70 61 74 63 68 2d 3e |n.....Di|spatch->|
|00000d00| 49 6e 63 6c 65 76 65 6c | 28 29 3b 0d 0a 09 56 4d |Inclevel|();...VM|
|00000d10| 61 74 72 69 78 20 51 2c | 20 52 3b 0d 0a 0d 0a 09 |atrix Q,| R;.....|
|00000d20| 51 52 28 78 2c 20 51 2c | 20 52 29 3b 0d 0a 09 56 |QR(x, Q,| R);...V|
|00000d30| 4d 61 74 72 69 78 20 62 | 65 74 61 68 61 74 20 3d |Matrix b|etahat =|
|00000d40| 20 49 6e 76 28 54 72 61 | 6e 28 52 29 20 2a 52 29 | Inv(Tra|n(R) *R)|
|00000d50| 20 2a 54 72 61 6e 28 52 | 29 20 2a 54 72 61 6e 28 | *Tran(R|) *Tran(|
|00000d60| 51 29 20 2a 79 3b 0d 0a | 0d 0a 09 44 69 73 70 61 |Q) *y;..|...Dispa|
|00000d70| 74 63 68 2d 3e 50 75 73 | 68 28 62 65 74 61 68 61 |tch->Pus|h(betaha|
|00000d80| 74 29 3b 0d 0a 09 72 65 | 74 75 72 6e 20 44 69 73 |t);...re|turn Dis|
|00000d90| 70 61 74 63 68 2d 3e 44 | 65 63 52 65 74 75 72 6e |patch->D|ecReturn|
|00000da0| 28 29 3b 0d 0a 7d 0d 0a | 0d 0a 56 4d 61 74 72 69 |();..}..|..VMatri|
|00000db0| 78 20 26 47 69 6e 76 52 | 65 67 72 65 73 73 69 6f |x &GinvR|egressio|
|00000dc0| 6e 28 56 4d 61 74 72 69 | 78 20 26 78 2c 20 56 4d |n(VMatri|x &x, VM|
|00000dd0| 61 74 72 69 78 20 26 79 | 29 0d 0a 7b 0d 0a 09 2f |atrix &y|)..{.../|
|00000de0| 2f 20 75 73 65 20 47 69 | 6e 76 20 74 6f 20 73 6f |/ use Gi|nv to so|
|00000df0| 6c 76 65 20 6e 6f 72 6d | 61 6c 20 65 71 75 61 74 |lve norm|al equat|
|00000e00| 69 6f 6e 73 0d 0a 09 44 | 69 73 70 61 74 63 68 2d |ions...D|ispatch-|
|00000e10| 3e 49 6e 63 6c 65 76 65 | 6c 28 29 3b 0d 0a 0d 0a |>Incleve|l();....|
|00000e20| 09 56 4d 61 74 72 69 78 | 20 62 65 74 61 68 61 74 |.VMatrix| betahat|
|00000e30| 20 3d 20 47 69 6e 76 28 | 54 72 61 6e 28 78 29 20 | = Ginv(|Tran(x) |
|00000e40| 2a 78 29 20 2a 54 72 61 | 6e 28 78 29 20 2a 79 3b |*x) *Tra|n(x) *y;|
|00000e50| 0d 0a 0d 0a 09 44 69 73 | 70 61 74 63 68 2d 3e 50 |.....Dis|patch->P|
|00000e60| 75 73 68 28 62 65 74 61 | 68 61 74 29 3b 0d 0a 09 |ush(beta|hat);...|
|00000e70| 72 65 74 75 72 6e 20 44 | 69 73 70 61 74 63 68 2d |return D|ispatch-|
|00000e80| 3e 44 65 63 52 65 74 75 | 72 6e 28 29 3b 0d 0a 7d |>DecRetu|rn();..}|
|00000e90| 0d 0a 0d 0a 0d 0a 56 4d | 61 74 72 69 78 20 26 53 |......VM|atrix &S|
|00000ea0| 56 44 72 65 67 72 65 73 | 73 69 6f 6e 28 56 4d 61 |VDregres|sion(VMa|
|00000eb0| 74 72 69 78 20 26 78 2c | 20 56 4d 61 74 72 69 78 |trix &x,| VMatrix|
|00000ec0| 20 26 79 29 0d 0a 7b 0d | 0a 09 2f 2f 20 75 73 65 | &y)..{.|..// use|
|00000ed0| 20 53 56 44 20 74 6f 20 | 73 6f 6c 76 65 20 72 65 | SVD to |solve re|
|00000ee0| 67 72 65 73 73 69 6f 6e | 20 78 20 3d 20 53 20 44 |gression| x = S D|
|00000ef0| 69 61 67 28 56 29 20 54 | 72 61 6e 28 44 29 0d 0a |iag(V) T|ran(D)..|
|00000f00| 09 44 69 73 70 61 74 63 | 68 2d 3e 49 6e 63 6c 65 |.Dispatc|h->Incle|
|00000f10| 76 65 6c 28 29 3b 0d 0a | 09 56 4d 61 74 72 69 78 |vel();..|.VMatrix|
|00000f20| 20 53 2c 20 56 2c 20 44 | 3b 0d 0a 0d 0a 09 53 76 | S, V, D|;.....Sv|
|00000f30| 64 28 78 2c 20 53 2c 20 | 56 2c 20 44 29 3b 0d 0a |d(x, S, |V, D);..|
|00000f40| 09 56 4d 61 74 72 69 78 | 20 74 20 3d 20 46 69 6c |.VMatrix| t = Fil|
|00000f50| 6c 28 56 2e 72 2c 20 56 | 2e 72 2c 20 30 2e 30 29 |l(V.r, V|.r, 0.0)|
|00000f60| 3b 0d 0a 09 66 6f 72 20 | 28 69 6e 74 20 69 20 3d |;...for |(int i =|
|00000f70| 20 31 3b 20 69 20 3c 3d | 20 56 2e 72 3b 20 69 2b | 1; i <=| V.r; i+|
|00000f80| 2b 29 20 7b 0d 0a 09 09 | 64 6f 75 62 6c 65 20 74 |+) {....|double t|
|00000f90| 74 20 3d 20 56 2e 6d 28 | 69 2c 20 31 29 3b 0d 0a |t = V.m(|i, 1);..|
|00000fa0| 09 09 74 2e 4d 28 69 2c | 20 69 29 20 3d 20 28 66 |..t.M(i,| i) = (f|
|00000fb0| 61 62 73 28 74 74 29 20 | 3c 3d 20 30 2e 30 29 20 |abs(tt) |<= 0.0) |
|00000fc0| 3f 20 30 2e 30 20 3a 20 | 31 2e 30 20 2f 20 74 74 |? 0.0 : |1.0 / tt|
|00000fd0| 3b 0d 0a 09 7d 0d 0a 09 | 56 4d 61 74 72 69 78 20 |;...}...|VMatrix |
|00000fe0| 62 65 74 61 68 61 74 20 | 3d 20 44 2a 74 2a 54 72 |betahat |= D*t*Tr|
|00000ff0| 61 6e 28 53 29 2a 79 3b | 0d 0a 0d 0a 09 44 69 73 |an(S)*y;|.....Dis|
|00001000| 70 61 74 63 68 2d 3e 50 | 75 73 68 28 62 65 74 61 |patch->P|ush(beta|
|00001010| 68 61 74 29 3b 0d 0a 09 | 72 65 74 75 72 6e 20 44 |hat);...|return D|
|00001020| 69 73 70 61 74 63 68 2d | 3e 44 65 63 52 65 74 75 |ispatch-|>DecRetu|
|00001030| 72 6e 28 29 3b 0d 0a 7d | 0d 0a 56 4d 61 74 72 69 |rn();..}|..VMatri|
|00001040| 78 20 26 47 65 74 53 65 | 72 69 61 6c 43 6f 76 61 |x &GetSe|rialCova|
|00001050| 72 28 56 4d 61 74 72 69 | 78 20 26 52 29 0d 0a 7b |r(VMatri|x &R)..{|
|00001060| 0d 0a 20 20 20 2f 2f 20 | 50 61 72 61 6d 65 74 65 |..   // |Paramete|
|00001070| 72 73 20 74 6f 20 43 4f | 52 52 20 69 6e 20 54 69 |rs to CO|RR in Ti|
|00001080| 6d 65 73 6c 61 62 0d 0a | 20 20 20 2f 2f 20 63 6f |meslab..|   // co|
|00001090| 72 72 65 6c 6f 67 72 61 | 6d 20 3d 20 43 4f 52 52 |rrelogra|m = CORR|
|000010a0| 28 78 3d 52 2c 6e 3d 52 | 2e 72 2c 4d 3d 52 2e 72 |(x=R,n=R|.r,M=R.r|
|000010b0| 2d 31 2c 51 3d 32 2a 52 | 2e 72 2c 0d 0a 20 20 20 |-1,Q=2*R|.r,..   |
|000010c0| 2f 2f 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |//      |        |
|000010d0| 20 20 20 20 20 20 69 6f | 70 74 3d 31 2c 72 30 3d |      io|pt=1,r0=|
|000010e0| 72 30 2c 20 70 65 72 20 | 3d 20 73 70 65 63 74 64 |r0, per |= spectd|
|000010f0| 65 6e 73 29 0d 0a 20 20 | 20 2f 2f 20 63 6f 76 61 |ens)..  | // cova|
|00001100| 72 20 3d 20 63 6f 72 72 | 65 6c 6f 67 72 61 6d 2a |r = corr|elogram*|
|00001110| 72 30 0d 0a 20 20 20 44 | 69 73 70 61 74 63 68 2d |r0..   D|ispatch-|
|00001120| 3e 49 6e 63 6c 65 76 65 | 6c 28 29 3b 0d 0a 20 20 |>Incleve|l();..  |
|00001130| 20 56 4d 61 74 72 69 78 | 20 63 65 6e 74 65 72 65 | VMatrix| centere|
|00001140| 64 2c 20 7a 2c 20 63 6f | 76 61 72 2c 20 73 70 65 |d, z, co|var, spe|
|00001150| 63 74 64 65 6e 73 3b 0d | 0a 20 20 20 64 6f 75 62 |ctdens;.|.   doub|
|00001160| 6c 65 20 6e 20 3d 20 28 | 64 6f 75 62 6c 65 29 20 |le n = (|double) |
|00001170| 52 2e 72 3b 0d 0a 20 20 | 20 2f 2f 20 63 65 6e 74 |R.r;..  | // cent|
|00001180| 65 72 20 61 20 63 6f 6c | 75 6d 6e 20 76 65 63 74 |er a col|umn vect|
|00001190| 6f 72 0d 0a 20 20 20 63 | 65 6e 74 65 72 65 64 20 |or..   c|entered |
|000011a0| 3d 20 52 20 2d 20 53 75 | 6d 28 52 29 2e 6d 28 31 |= R - Su|m(R).m(1|
|000011b0| 2c 20 31 29 20 2f 20 6e | 3b 0d 0a 20 20 20 2f 2f |, 1) / n|;..   //|
|000011c0| 20 7a 65 72 6f 20 70 61 | 64 20 74 6f 20 6c 65 6e | zero pa|d to len|
|000011d0| 67 74 68 20 32 6e 3a 20 | 32 6e 20 70 65 72 69 6f |gth 2n: |2n perio|
|000011e0| 64 69 63 20 66 6f 72 20 | 66 75 6c 6c 0d 0a 20 20 |dic for |full..  |
|000011f0| 20 2f 2f 20 73 61 6d 70 | 6c 65 20 73 70 65 63 74 | // samp|le spect|
|00001200| 72 61 6c 20 64 65 6e 73 | 69 74 79 0d 0a 20 20 20 |ral dens|ity..   |
|00001210| 63 65 6e 74 65 72 65 64 | 20 3d 20 43 76 28 63 65 |centered| = Cv(ce|
|00001220| 6e 74 65 72 65 64 2c 20 | 46 69 6c 6c 28 52 2e 72 |ntered, |Fill(R.r|
|00001230| 2c 20 52 2e 63 2c 20 30 | 29 29 3b 0d 0a 20 20 20 |, R.c, 0|));..   |
|00001240| 2f 2f 20 74 61 6b 65 20 | 66 66 74 0d 0a 20 20 20 |// take |fft..   |
|00001250| 7a 20 3d 20 46 66 74 28 | 63 65 6e 74 65 72 65 64 |z = Fft(|centered|
|00001260| 29 3b 0d 0a 20 20 20 2f | 2f 20 74 61 6b 65 20 63 |);..   /|/ take c|
|00001270| 6f 6e 76 6f 6c 75 74 69 | 6f 6e 20 3a 20 67 69 76 |onvoluti|on : giv|
|00001280| 65 73 20 73 61 6d 70 6c | 65 20 73 70 65 63 74 72 |es sampl|e spectr|
|00001290| 61 6c 20 64 65 6e 73 69 | 74 79 0d 0a 20 20 20 73 |al densi|ty..   s|
|000012a0| 70 65 63 74 64 65 6e 73 | 20 3d 20 53 75 6d 28 7a |pectdens| = Sum(z|
|000012b0| 20 25 20 7a 20 2f 20 6e | 2c 20 43 4f 4c 55 4d 4e | % z / n|, COLUMN|
|000012c0| 53 29 3b 0d 0a 20 20 20 | 2f 2f 20 69 6e 76 65 72 |S);..   |// inver|
|000012d0| 73 65 20 66 66 74 20 66 | 6f 72 20 73 65 72 69 61 |se fft f|or seria|
|000012e0| 6c 20 63 6f 72 72 65 6c | 61 74 69 6f 6e 20 28 61 |l correl|ation (a|
|000012f0| 75 74 6f 63 6f 72 72 65 | 6c 61 74 69 6f 6e 29 0d |utocorre|lation).|
|00001300| 0a 20 20 20 63 6f 76 61 | 72 20 3d 20 46 66 74 28 |.   cova|r = Fft(|
|00001310| 73 70 65 63 74 64 65 6e | 73 2c 20 46 46 41 4c 53 |spectden|s, FFALS|
|00001320| 45 29 3b 0d 0a 20 20 20 | 2f 2f 20 74 68 72 6f 77 |E);..   |// throw|
|00001330| 20 61 77 61 79 20 6c 61 | 73 74 20 68 61 6c 66 2e | away la|st half.|
|00001340| 0d 0a 20 20 20 63 6f 76 | 61 72 20 3d 20 53 75 62 |..   cov|ar = Sub|
|00001350| 6d 61 74 28 63 6f 76 61 | 72 2c 20 52 2e 72 2c 20 |mat(cova|r, R.r, |
|00001360| 32 29 3b 0d 0a 20 20 20 | 44 69 73 70 61 74 63 68 |2);..   |Dispatch|
|00001370| 2d 3e 50 75 73 68 28 63 | 6f 76 61 72 29 3b 0d 0a |->Push(c|ovar);..|
|00001380| 20 20 20 72 65 74 75 72 | 6e 20 44 69 73 70 61 74 |   retur|n Dispat|
|00001390| 63 68 2d 3e 44 65 63 52 | 65 74 75 72 6e 28 29 3b |ch->DecR|eturn();|
|000013a0| 0d 0a 7d 0d 0a 0d 0a 6d | 61 69 6e 28 29 0d 0a 7b |..}....m|ain()..{|
|000013b0| 0d 0a 09 44 69 73 70 61 | 74 63 68 2d 3e 49 6e 63 |...Dispa|tch->Inc|
|000013c0| 6c 65 76 65 6c 28 29 3b | 0d 0a 09 56 4d 61 74 72 |level();|...VMatr|
|000013d0| 69 78 20 78 2c 20 62 65 | 74 61 28 22 62 65 74 61 |ix x, be|ta("beta|
|000013e0| 22 2c 20 33 2c 20 31 29 | 2c 20 79 2c 20 62 65 74 |", 3, 1)|, y, bet|
|000013f0| 61 68 61 74 2c 20 72 65 | 73 69 64 73 2c 20 73 65 |ahat, re|sids, se|
|00001400| 72 69 61 6c 3b 0d 0a 0d | 0a 09 54 65 73 74 48 69 |rial;...|..TestHi|
|00001410| 6c 62 28 20 31 31 20 29 | 3b 0d 0a 0d 0a 09 53 65 |lb( 11 )|;.....Se|
|00001420| 74 77 69 64 28 31 35 29 | 3b 0d 0a 09 53 65 74 64 |twid(15)|;...Setd|
|00001430| 65 63 28 31 30 29 3b 0d | 0a 0d 0a 0d 0a 09 62 65 |ec(10);.|......be|
|00001440| 74 61 2e 4d 28 31 2c 20 | 31 29 20 3d 20 31 3b 0d |ta.M(1, |1) = 1;.|
|00001450| 0a 09 62 65 74 61 2e 4d | 28 32 2c 20 31 29 20 3d |..beta.M|(2, 1) =|
|00001460| 20 30 2e 35 3b 0d 0a 09 | 62 65 74 61 2e 4d 28 33 | 0.5;...|beta.M(3|
|00001470| 2c 20 31 29 20 3d 20 30 | 2e 32 35 3b 0d 0a 0d 0a |, 1) = 0|.25;....|
|00001480| 09 78 20 3d 20 67 65 74 | 78 28 35 35 29 3b 0d 0a |.x = get|x(55);..|
|00001490| 09 79 20 3d 20 67 65 74 | 79 28 78 2c 20 62 65 74 |.y = get|y(x, bet|
|000014a0| 61 29 3b 0d 0a 0d 0a 09 | 62 65 74 61 68 61 74 20 |a);.....|betahat |
|000014b0| 3d 20 72 65 67 72 65 73 | 73 69 6f 6e 28 78 2c 20 |= regres|sion(x, |
|000014c0| 79 29 3b 0d 0a 09 62 65 | 74 61 68 61 74 2e 4e 61 |y);...be|tahat.Na|
|000014d0| 6d 65 69 74 28 22 54 65 | 78 74 20 62 6f 6f 6b 20 |meit("Te|xt book |
|000014e0| 62 65 74 61 68 61 74 22 | 29 3b 0d 0a 09 28 54 72 |betahat"|);...(Tr|
|000014f0| 61 6e 28 62 65 74 61 29 | 29 2e 44 69 73 70 6c 61 |an(beta)|).Displa|
|00001500| 79 4d 61 74 28 29 3b 0d | 0a 09 28 54 72 61 6e 28 |yMat();.|..(Tran(|
|00001510| 62 65 74 61 68 61 74 29 | 29 2e 44 69 73 70 6c 61 |betahat)|).Displa|
|00001520| 79 4d 61 74 28 29 3b 0d | 0a 0d 0a 09 62 65 74 61 |yMat();.|....beta|
|00001530| 68 61 74 20 3d 20 51 52 | 72 65 67 72 65 73 73 69 |hat = QR|regressi|
|00001540| 6f 6e 28 78 2c 20 79 29 | 3b 0d 0a 09 62 65 74 61 |on(x, y)|;...beta|
|00001550| 68 61 74 2e 4e 61 6d 65 | 69 74 28 22 51 52 20 62 |hat.Name|it("QR b|
|00001560| 65 74 61 68 61 74 22 29 | 3b 0d 0a 09 28 54 72 61 |etahat")|;...(Tra|
|00001570| 6e 28 62 65 74 61 29 29 | 2e 44 69 73 70 6c 61 79 |n(beta))|.Display|
|00001580| 4d 61 74 28 29 3b 0d 0a | 09 28 54 72 61 6e 28 62 |Mat();..|.(Tran(b|
|00001590| 65 74 61 68 61 74 29 29 | 2e 44 69 73 70 6c 61 79 |etahat))|.Display|
|000015a0| 4d 61 74 28 29 3b 0d 0a | 0d 0a 09 62 65 74 61 68 |Mat();..|...betah|
|000015b0| 61 74 20 3d 20 47 69 6e | 76 52 65 67 72 65 73 73 |at = Gin|vRegress|
|000015c0| 69 6f 6e 28 78 2c 20 79 | 29 3b 0d 0a 09 62 65 74 |ion(x, y|);...bet|
|000015d0| 61 68 61 74 2e 4e 61 6d | 65 69 74 28 22 47 69 6e |ahat.Nam|eit("Gin|
|000015e0| 76 20 62 65 74 61 68 61 | 74 22 29 3b 0d 0a 09 28 |v betaha|t");...(|
|000015f0| 54 72 61 6e 28 62 65 74 | 61 29 29 2e 44 69 73 70 |Tran(bet|a)).Disp|
|00001600| 6c 61 79 4d 61 74 28 29 | 3b 0d 0a 09 28 54 72 61 |layMat()|;...(Tra|
|00001610| 6e 28 62 65 74 61 68 61 | 74 29 29 2e 44 69 73 70 |n(betaha|t)).Disp|
|00001620| 6c 61 79 4d 61 74 28 29 | 3b 0d 0a 0d 0a 09 62 65 |layMat()|;.....be|
|00001630| 74 61 68 61 74 20 3d 20 | 53 56 44 72 65 67 72 65 |tahat = |SVDregre|
|00001640| 73 73 69 6f 6e 28 78 2c | 20 79 29 3b 0d 0a 09 62 |ssion(x,| y);...b|
|00001650| 65 74 61 68 61 74 2e 4e | 61 6d 65 69 74 28 22 53 |etahat.N|ameit("S|
|00001660| 56 44 20 72 65 67 72 65 | 73 73 69 6f 6e 22 29 3b |VD regre|ssion");|
|00001670| 0d 0a 09 28 54 72 61 6e | 28 62 65 74 61 29 29 2e |...(Tran|(beta)).|
|00001680| 44 69 73 70 6c 61 79 4d | 61 74 28 29 3b 0d 0a 09 |DisplayM|at();...|
|00001690| 28 54 72 61 6e 28 62 65 | 74 61 68 61 74 29 29 2e |(Tran(be|tahat)).|
|000016a0| 44 69 73 70 6c 61 79 4d | 61 74 28 29 3b 0d 0a 0d |DisplayM|at();...|
|000016b0| 0a 09 72 65 73 69 64 73 | 20 3d 20 79 20 2d 20 78 |..resids| = y - x|
|000016c0| 2a 62 65 74 61 68 61 74 | 3b 0d 0a 09 73 65 72 69 |*betahat|;...seri|
|000016d0| 61 6c 20 3d 20 47 65 74 | 53 65 72 69 61 6c 43 6f |al = Get|SerialCo|
|000016e0| 76 61 72 28 72 65 73 69 | 64 73 29 3b 0d 0a 09 28 |var(resi|ds);...(|
|000016f0| 54 72 61 6e 28 53 75 62 | 6d 61 74 28 73 65 72 69 |Tran(Sub|mat(seri|
|00001700| 61 6c 2c 20 35 2c 20 31 | 29 29 29 2e 44 69 73 70 |al, 5, 1|))).Disp|
|00001710| 6c 61 79 4d 61 74 28 29 | 3b 0d 0a 0d 0a 09 53 65 |layMat()|;.....Se|
|00001720| 74 77 69 64 28 36 29 3b | 0d 0a 09 53 65 74 64 65 |twid(6);|...Setde|
|00001730| 63 28 33 29 3b 0d 0a 0d | 0a 09 76 63 6c 6f 73 65 |c(3);...|..vclose|
|00001740| 28 29 3b 0d 0a 7d 0d 0a | 0d 0a 5c 65 6e 64 7b 76 |();..}..|..\end{v|
|00001750| 65 72 62 61 74 69 6d 7d |                         |erbatim}|        |
+--------+-------------------------+-------------------------+--------+--------+