/ Stars of Shareware: Programmierung / SOURCE.mdf / programm / msdos / c / yamp / app3.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 67 72 61 |\chapter|{testgra|
|00000010| 66 2e 63 70 70 7d 0d 0a | 54 68 69 73 20 69 73 20 |f.cpp}..|This is |
|00000020| 74 68 65 20 74 65 73 74 | 20 73 75 69 74 65 20 66 |the test| suite f|
|00000030| 6f 72 20 74 68 65 20 67 | 72 61 70 68 69 63 73 20 |or the g|raphics |
|00000040| 72 6f 75 74 69 6e 65 73 | 2e 20 49 74 20 70 6c 6f |routines|. It plo|
|00000050| 74 73 20 74 68 65 0d 0a | 72 65 73 69 64 75 61 6c |ts the..|residual|
|00000060| 73 2c 20 61 6e 64 20 63 | 6f 72 72 65 6c 6f 67 72 |s, and c|orrelogr|
|00000070| 61 6d 20 61 6e 64 20 70 | 65 72 69 6f 64 6f 67 72 |am and p|eriodogr|
|00000080| 61 6d 20 6f 66 20 74 68 | 65 20 72 65 73 69 64 75 |am of th|e residu|
|00000090| 61 6c 73 2e 20 53 65 65 | 0d 0a 4e 65 77 74 6f 6e |als. See|..Newton|
|000000a0| 20 5c 63 69 74 65 7b 4e | 4a 3a 74 6c 7d 20 6f 72 | \cite{N|J:tl} or|
|000000b0| 20 50 72 65 73 73 20 65 | 74 2e 61 6c 2e 20 5c 63 | Press e|t.al. \c|
|000000c0| 69 74 65 7b 50 72 3a 6e | 72 63 7d 20 66 6f 72 20 |ite{Pr:n|rc} for |
|000000d0| 61 6e 0d 0a 65 78 70 6c | 61 6e 61 74 69 6f 6e 20 |an..expl|anation |
|000000e0| 6f 66 20 74 68 65 73 65 | 20 63 61 6c 63 75 6c 61 |of these| calcula|
|000000f0| 74 69 6f 6e 73 2e 20 54 | 68 65 79 20 69 64 65 6e |tions. T|hey iden|
|00000100| 74 69 66 79 20 73 65 72 | 69 61 6c 0d 0a 63 6f 72 |tify ser|ial..cor|
|00000110| 72 65 6c 61 74 69 6f 6e | 20 69 6e 20 74 68 65 20 |relation| in the |
|00000120| 72 65 73 69 64 75 61 6c | 73 20 6f 66 20 72 65 67 |residual|s of reg|
|00000130| 72 65 73 73 69 6f 6e 20 | 74 68 61 74 20 72 65 6d |ression |that rem|
|00000140| 6f 76 65 20 61 20 6c 69 | 6e 65 61 72 0d 0a 74 72 |ove a li|near..tr|
|00000150| 65 6e 64 20 66 72 6f 6d | 20 74 68 65 20 64 61 74 |end from| the dat|
|00000160| 61 2e 20 54 68 65 20 73 | 65 72 69 61 6c 20 63 6f |a. The s|erial co|
|00000170| 72 65 6c 61 74 69 6f 6e | 20 77 61 73 20 69 6e 74 |relation| was int|
|00000180| 65 6e 74 69 6f 6e 61 6c | 6c 79 0d 0a 70 6c 61 63 |entional|ly..plac|
|00000190| 65 64 20 6f 6e 20 74 68 | 65 20 72 65 73 69 64 75 |ed on th|e residu|
|000001a0| 61 6c 73 2e 20 0d 0a 0d | 0a 5c 62 65 67 69 6e 7b |als. ...|.\begin{|
|000001b0| 76 65 72 62 61 74 69 6d | 7d 0d 0a 0d 0a 0d 0a 23 |verbatim|}......#|
|000001c0| 69 6e 63 6c 75 64 65 20 | 22 76 69 72 74 2e 68 22 |include |"virt.h"|
|000001d0| 0d 0a 0d 0a 2f 2f 72 65 | 71 75 69 72 65 64 20 67 |....//re|quired g|
|000001e0| 6c 6f 62 61 6c 20 64 65 | 63 6c 61 72 61 74 69 6f |lobal de|claratio|
|000001f0| 6e 20 66 6f 72 20 74 68 | 65 0d 0a 2f 2f 20 20 6d |n for th|e..//  m|
|00000200| 61 74 72 69 78 20 73 74 | 61 63 6b 20 6f 62 6a 65 |atrix st|ack obje|
|00000210| 63 74 0d 0a 0d 0a 2f 2f | 20 75 6e 73 69 67 6e 65 |ct....//| unsigne|
|00000220| 64 20 69 6e 74 20 5f 73 | 74 6b 6c 65 6e 20 3d 20 |d int _s|tklen = |
|00000230| 53 54 41 43 4b 4c 45 4e | 47 54 48 3b 0d 0a 0d 0a |STACKLEN|GTH;....|
|00000240| 4d 53 74 61 63 6b 20 2a | 44 69 73 70 61 74 63 68 |MStack *|Dispatch|
|00000250| 20 3d 20 6e 65 77 20 4d | 53 74 61 63 6b 3b 0d 0a | = new M|Stack;..|
|00000260| 0d 0a 0d 0a 56 4d 61 74 | 72 69 78 20 26 67 65 74 |....VMat|rix &get|
|00000270| 78 28 20 69 6e 74 20 4e | 20 29 0d 0a 20 20 2f 2f |x( int N| )..  //|
|00000280| 20 63 72 65 61 74 65 20 | 61 6e 20 78 20 6d 61 74 | create |an x mat|
|00000290| 72 69 78 0d 0a 20 20 7b | 0d 0a 20 20 20 20 20 20 |rix..  {|..      |
|000002a0| 44 69 73 70 61 74 63 68 | 2d 3e 49 6e 63 6c 65 76 |Dispatch|->Inclev|
|000002b0| 65 6c 28 29 3b 0d 0a 20 | 20 20 20 20 20 56 4d 61 |el();.. |     VMa|
|000002c0| 74 72 69 78 20 78 2c 20 | 63 31 2c 20 78 32 3b 0d |trix x, |c1, x2;.|
|000002d0| 0a 0d 0a 20 20 20 20 20 | 20 63 31 20 3d 20 46 69 |...     | c1 = Fi|
|000002e0| 6c 6c 28 4e 2c 31 2c 31 | 2e 30 29 3b 0d 0a 20 20 |ll(N,1,1|.0);..  |
|000002f0| 20 20 20 20 78 20 3d 20 | 49 6e 64 65 78 28 20 31 |    x = |Index( 1|
|00000300| 2c 20 4e 20 29 20 2d 20 | 28 28 64 6f 75 62 6c 65 |, N ) - |((double|
|00000310| 29 4e 29 2a 30 2e 35 3b | 0d 0a 20 20 20 20 20 20 |)N)*0.5;|..      |
|00000320| 78 20 3d 20 43 68 28 20 | 63 31 2c 78 20 29 3b 0d |x = Ch( |c1,x );.|
|00000330| 0a 0d 0a 20 20 20 20 20 | 20 2f 2f 20 70 75 73 68 |...     | // push|
|00000340| 20 78 20 6f 6e 74 6f 20 | 74 68 65 20 73 74 61 63 | x onto |the stac|
|00000350| 6b 0d 0a 20 20 20 20 20 | 20 44 69 73 70 61 74 63 |k..     | Dispatc|
|00000360| 68 2d 3e 50 75 73 68 28 | 78 29 3b 0d 0a 20 20 20 |h->Push(|x);..   |
|00000370| 20 20 20 2f 2f 20 64 65 | 63 72 65 6d 65 6e 74 20 |   // de|crement |
|00000380| 74 68 65 20 73 75 62 72 | 6f 75 74 69 6e 65 20 6e |the subr|outine n|
|00000390| 65 73 74 69 6e 67 20 6c | 65 76 65 6c 0d 0a 20 20 |esting l|evel..  |
|000003a0| 20 20 20 20 2f 2f 20 61 | 6e 64 20 72 65 74 75 72 |    // a|nd retur|
|000003b0| 6e 20 74 68 65 20 73 74 | 61 63 6b 20 74 6f 70 0d |n the st|ack top.|
|000003c0| 0a 20 20 20 20 20 20 72 | 65 74 75 72 6e 20 44 69 |.      r|eturn Di|
|000003d0| 73 70 61 74 63 68 2d 3e | 44 65 63 52 65 74 75 72 |spatch->|DecRetur|
|000003e0| 6e 28 29 3b 0d 0a 20 20 | 7d 0d 0a 0d 0a 56 4d 61 |n();..  |}....VMa|
|000003f0| 74 72 69 78 20 26 67 65 | 74 79 28 20 56 4d 61 74 |trix &ge|ty( VMat|
|00000400| 72 69 78 20 26 78 2c 20 | 56 4d 61 74 72 69 78 20 |rix &x, |VMatrix |
|00000410| 26 62 65 74 61 29 0d 0a | 20 20 2f 2f 20 63 72 65 |&beta)..|  // cre|
|00000420| 61 74 65 20 61 20 79 20 | 6d 61 74 72 69 78 0d 0a |ate a y |matrix..|
|00000430| 20 20 7b 0d 0a 20 20 20 | 20 20 20 44 69 73 70 61 |  {..   |   Dispa|
|00000440| 74 63 68 2d 3e 49 6e 63 | 6c 65 76 65 6c 28 29 3b |tch->Inc|level();|
|00000450| 0d 0a 20 20 20 20 20 20 | 56 4d 61 74 72 69 78 20 |..      |VMatrix |
|00000460| 79 3b 0d 0a 0d 0a 20 20 | 20 20 20 20 79 20 3d 20 |y;....  |    y = |
|00000470| 78 2a 62 65 74 61 3b 0d | 0a 20 20 20 20 20 20 73 |x*beta;.|.      s|
|00000480| 72 61 6e 64 28 31 32 33 | 29 3b 0d 0a 20 20 20 20 |rand(123|);..    |
|00000490| 20 20 66 6f 72 28 69 6e | 74 20 69 3d 31 3b 20 69 |  for(in|t i=1; i|
|000004a0| 3c 3d 79 2e 72 3b 20 69 | 2b 2b 29 20 7b 0d 0a 20 |<=y.r; i|++) {.. |
|000004b0| 20 20 20 20 20 20 20 2f | 2f 20 75 73 65 20 73 75 |       /|/ use su|
|000004c0| 6d 20 6f 66 20 33 20 75 | 6e 69 66 6f 72 6d 73 20 |m of 3 u|niforms |
|000004d0| 66 6f 72 20 61 6e 20 61 | 70 70 72 6f 78 69 6d 61 |for an a|pproxima|
|000004e0| 74 65 0d 0a 20 20 20 20 | 20 20 20 20 2f 2f 20 6e |te..    |    // n|
|000004f0| 6f 72 6d 61 6c 20 72 61 | 6e 64 6f 6d 20 76 61 72 |ormal ra|ndom var|
|00000500| 69 61 62 6c 65 0d 0a 20 | 20 20 20 20 20 20 20 69 |iable.. |       i|
|00000510| 6e 74 20 75 20 3d 20 72 | 61 6e 64 6f 6d 28 31 30 |nt u = r|andom(10|
|00000520| 30 29 2b 72 61 6e 64 6f | 6d 28 31 30 30 29 2b 72 |0)+rando|m(100)+r|
|00000530| 61 6e 64 6f 6d 28 31 30 | 30 29 2b 33 3b 0d 0a 20 |andom(10|0)+3;.. |
|00000540| 20 20 20 20 20 20 20 79 | 2e 4d 28 69 2c 31 29 20 |       y|.M(i,1) |
|00000550| 3d 20 79 2e 6d 28 69 2c | 31 29 20 2b 20 35 2e 30 |= y.m(i,|1) + 5.0|
|00000560| 2a 73 69 6e 28 20 33 2e | 31 34 2a 28 28 64 6f 75 |*sin( 3.|14*((dou|
|00000570| 62 6c 65 29 20 28 69 25 | 38 29 29 2f 37 2e 30 20 |ble) (i%|8))/7.0 |
|00000580| 29 0d 0a 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |)..     |        |
|00000590| 20 20 20 20 2b 20 28 28 | 64 6f 75 62 6c 65 29 20 |    + ((|double) |
|000005a0| 28 75 2d 31 35 30 29 29 | 2f 33 30 30 2e 30 3b 0d |(u-150))|/300.0;.|
|000005b0| 0a 20 20 20 20 20 20 7d | 0d 0a 20 20 20 20 20 20 |.      }|..      |
|000005c0| 44 69 73 70 61 74 63 68 | 2d 3e 50 75 73 68 28 79 |Dispatch|->Push(y|
|000005d0| 29 3b 0d 0a 20 20 20 20 | 20 20 72 65 74 75 72 6e |);..    |  return|
|000005e0| 20 44 69 73 70 61 74 63 | 68 2d 3e 44 65 63 52 65 | Dispatc|h->DecRe|
|000005f0| 74 75 72 6e 28 29 3b 0d | 0a 20 20 7d 0d 0a 0d 0a |turn();.|.  }....|
|00000600| 56 4d 61 74 72 69 78 20 | 26 72 65 67 72 65 73 73 |VMatrix |&regress|
|00000610| 69 6f 6e 28 20 56 4d 61 | 74 72 69 78 26 20 78 2c |ion( VMa|trix& x,|
|00000620| 20 56 4d 61 74 72 69 78 | 26 20 79 20 29 0d 0a 2f | VMatrix|& y )../|
|00000630| 2f 20 64 6f 20 61 20 6d | 75 6c 74 69 70 6c 65 20 |/ do a m|ultiple |
|00000640| 6c 69 6e 65 61 72 20 72 | 65 67 72 65 73 73 69 6f |linear r|egressio|
|00000650| 6e 0d 0a 20 20 20 20 20 | 20 7b 0d 0a 20 20 20 20 |n..     | {..    |
|00000660| 20 20 44 69 73 70 61 74 | 63 68 2d 3e 49 6e 63 6c |  Dispat|ch->Incl|
|00000670| 65 76 65 6c 28 29 3b 0d | 0a 20 20 20 20 20 20 56 |evel();.|.      V|
|00000680| 4d 61 74 72 69 78 20 79 | 78 2c 20 72 65 67 2c 20 |Matrix y|x, reg, |
|00000690| 62 65 74 61 68 61 74 3b | 0d 0a 20 20 20 20 20 20 |betahat;|..      |
|000006a0| 69 6e 74 20 4e 3d 78 2e | 72 2c 20 70 3d 78 2e 63 |int N=x.|r, p=x.c|
|000006b0| 3b 0d 0a 0d 0a 20 20 20 | 20 20 20 2f 2f 20 73 6f |;....   |   // so|
|000006c0| 6c 76 65 20 66 6f 72 20 | 72 65 67 72 65 73 73 69 |lve for |regressi|
|000006d0| 6f 6e 20 70 61 72 61 6d | 65 74 65 72 73 20 75 73 |on param|eters us|
|000006e0| 69 6e 67 20 73 77 65 65 | 70 0d 0a 20 20 20 20 20 |ing swee|p..     |
|000006f0| 20 79 78 20 3d 20 43 68 | 28 79 2c 78 29 3b 0d 0a | yx = Ch|(y,x);..|
|00000700| 20 20 20 20 20 20 72 65 | 67 20 3d 20 53 77 65 65 |      re|g = Swee|
|00000710| 70 28 20 32 2c 70 2b 31 | 2c 20 54 72 61 6e 28 79 |p( 2,p+1|, Tran(y|
|00000720| 78 29 2a 79 78 29 3b 0d | 0a 20 20 20 20 20 20 2f |x)*yx);.|.      /|
|00000730| 2f 20 63 61 6c 63 75 6c | 61 74 65 20 6d 65 61 6e |/ calcul|ate mean|
|00000740| 20 73 71 75 61 72 65 20 | 65 72 72 6f 72 0d 0a 20 | square |error.. |
|00000750| 20 20 20 20 20 72 65 67 | 2e 4d 28 31 2c 31 29 20 |     reg|.M(1,1) |
|00000760| 3d 20 72 65 67 2e 6d 28 | 31 2c 31 29 2f 28 28 64 |= reg.m(|1,1)/((d|
|00000770| 6f 75 62 6c 65 20 29 28 | 20 4e 2d 70 29 29 3b 0d |ouble )(| N-p));.|
|00000780| 0a 20 20 20 20 20 20 72 | 65 67 2e 44 69 73 70 6c |.      r|eg.Displ|
|00000790| 61 79 4d 61 74 28 29 3b | 0d 0a 0d 0a 0d 0a 20 20 |ayMat();|......  |
|000007a0| 20 20 20 20 2f 2f 20 73 | 6f 6c 76 65 20 72 65 67 |    // s|olve reg|
|000007b0| 72 65 73 73 69 6f 6e 20 | 75 73 69 6e 67 20 6e 6f |ression |using no|
|000007c0| 72 6d 61 6c 20 65 71 75 | 61 74 69 6f 6e 73 0d 0a |rmal equ|ations..|
|000007d0| 20 20 20 20 20 20 62 65 | 74 61 68 61 74 20 3d 20 |      be|tahat = |
|000007e0| 49 6e 76 28 54 72 61 6e | 28 78 29 2a 78 29 2a 54 |Inv(Tran|(x)*x)*T|
|000007f0| 72 61 6e 28 78 29 2a 79 | 3b 0d 0a 0d 0a 20 20 20 |ran(x)*y|;....   |
|00000800| 20 20 20 44 69 73 70 61 | 74 63 68 2d 3e 50 75 73 |   Dispa|tch->Pus|
|00000810| 68 28 20 62 65 74 61 68 | 61 74 20 29 3b 0d 0a 20 |h( betah|at );.. |
|00000820| 20 20 20 20 20 72 65 74 | 75 72 6e 20 44 69 73 70 |     ret|urn Disp|
|00000830| 61 74 63 68 2d 3e 44 65 | 63 52 65 74 75 72 6e 28 |atch->De|cReturn(|
|00000840| 29 3b 0d 0a 7d 0d 0a 76 | 6f 69 64 20 70 6c 6f 74 |);..}..v|oid plot|
|00000850| 52 65 73 69 64 75 61 6c | 73 28 20 56 4d 61 74 72 |Residual|s( VMatr|
|00000860| 69 78 20 26 72 65 73 69 | 64 73 20 29 0d 0a 20 20 |ix &resi|ds )..  |
|00000870| 7b 0d 0a 20 20 20 20 44 | 69 73 70 61 74 63 68 2d |{..    D|ispatch-|
|00000880| 3e 49 6e 63 6c 65 76 65 | 6c 28 29 3b 0d 0a 20 20 |>Incleve|l();..  |
|00000890| 20 20 56 4d 61 74 72 69 | 78 20 67 72 66 20 3d 20 |  VMatri|x grf = |
|000008a0| 43 68 28 20 49 6e 64 65 | 78 28 20 31 2c 20 72 65 |Ch( Inde|x( 1, re|
|000008b0| 73 69 64 73 2e 72 20 29 | 2c 20 72 65 73 69 64 73 |sids.r )|, resids|
|000008c0| 20 29 3b 0d 0a 20 20 20 | 20 47 4d 61 74 72 69 78 | );..   | GMatrix|
|000008d0| 20 41 67 72 61 70 68 28 | 20 67 72 66 2c 20 2d 27 | Agraph(| grf, -'|
|000008e0| 25 27 20 29 3b 0d 0a 20 | 20 20 20 2a 41 67 72 61 |%' );.. |   *Agra|
|000008f0| 70 68 2e 50 61 74 68 54 | 6f 44 72 69 76 65 72 20 |ph.PathT|oDriver |
|00000900| 3d 20 22 43 3a 5c 5c 74 | 63 5c 5c 62 67 69 22 3b |= "C:\\t|c\\bgi";|
|00000910| 0d 0a 20 20 20 20 2a 41 | 67 72 61 70 68 2e 74 69 |..    *A|graph.ti|
|00000920| 74 6c 65 20 3d 20 22 52 | 65 73 69 64 75 61 6c 73 |tle = "R|esiduals|
|00000930| 20 66 6f 72 20 64 61 74 | 61 22 3b 0d 0a 20 20 20 | for dat|a";..   |
|00000940| 20 2a 41 67 72 61 70 68 | 2e 74 69 74 6c 65 32 3d | *Agraph|.title2=|
|00000950| 20 22 77 69 74 68 20 73 | 65 72 69 61 6c 20 63 6f | "with s|erial co|
|00000960| 72 72 65 6c 61 74 69 6f | 6e 73 20 77 69 74 68 20 |rrelatio|ns with |
|00000970| 66 72 65 71 75 65 6e 63 | 79 20 30 2e 32 35 22 3b |frequenc|y 0.25";|
|00000980| 0d 0a 20 20 20 20 2a 41 | 67 72 61 70 68 2e 79 6e |..    *A|graph.yn|
|00000990| 61 6d 65 20 3d 20 22 52 | 65 73 69 64 75 61 6c 73 |ame = "R|esiduals|
|000009a0| 22 3b 0d 0a 20 20 20 20 | 2a 41 67 72 61 70 68 2e |";..    |*Agraph.|
|000009b0| 78 6e 61 6d 65 20 3d 20 | 22 49 6e 64 65 78 22 3b |xname = |"Index";|
|000009c0| 0d 0a 20 20 20 20 41 67 | 72 61 70 68 2e 48 72 65 |..    Ag|raph.Hre|
|000009d0| 66 28 20 30 2e 30 20 29 | 3b 0d 0a 20 20 20 20 41 |f( 0.0 )|;..    A|
|000009e0| 67 72 61 70 68 2e 53 68 | 6f 77 28 29 3b 0d 0a 20 |graph.Sh|ow();.. |
|000009f0| 20 20 20 44 69 73 70 61 | 74 63 68 2d 3e 43 6c 65 |   Dispa|tch->Cle|
|00000a00| 61 6e 73 74 61 63 6b 28 | 29 3b 0d 0a 20 20 20 20 |anstack(|);..    |
|00000a10| 44 69 73 70 61 74 63 68 | 2d 3e 44 65 63 6c 65 76 |Dispatch|->Declev|
|00000a20| 65 6c 28 29 3b 0d 0a 20 | 20 7d 0d 0a 0d 0a 56 4d |el();.. | }....VM|
|00000a30| 61 74 72 69 78 20 26 47 | 65 74 53 65 72 69 61 6c |atrix &G|etSerial|
|00000a40| 43 6f 76 61 72 28 56 4d | 61 74 72 69 78 20 26 52 |Covar(VM|atrix &R|
|00000a50| 2c 20 56 4d 61 74 72 69 | 78 20 26 73 70 65 63 74 |, VMatri|x &spect|
|00000a60| 64 65 6e 73 29 0d 0a 7b | 0d 0a 20 20 20 2f 2f 20 |dens)..{|..   // |
|00000a70| 50 61 72 61 6d 65 74 65 | 72 73 20 74 6f 20 43 4f |Paramete|rs to CO|
|00000a80| 52 52 20 69 6e 20 54 69 | 6d 65 73 6c 61 62 0d 0a |RR in Ti|meslab..|
|00000a90| 20 20 20 2f 2f 20 63 6f | 72 72 65 6c 6f 67 72 61 |   // co|rrelogra|
|00000aa0| 6d 20 3d 20 43 4f 52 52 | 28 78 3d 52 2c 6e 3d 52 |m = CORR|(x=R,n=R|
|00000ab0| 2e 72 2c 4d 3d 52 2e 72 | 2d 31 2c 51 3d 32 2a 52 |.r,M=R.r|-1,Q=2*R|
|00000ac0| 2e 72 2c 0d 0a 20 20 20 | 2f 2f 20 20 20 20 20 20 |.r,..   |//      |
|00000ad0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 69 6f |        |      io|
|00000ae0| 70 74 3d 31 2c 72 30 3d | 72 30 2c 20 70 65 72 20 |pt=1,r0=|r0, per |
|00000af0| 3d 20 73 70 65 63 74 64 | 65 6e 73 29 0d 0a 20 20 |= spectd|ens)..  |
|00000b00| 20 2f 2f 20 63 6f 76 61 | 72 20 3d 20 63 6f 72 72 | // cova|r = corr|
|00000b10| 65 6c 6f 67 72 61 6d 2a | 72 30 0d 0a 20 20 20 44 |elogram*|r0..   D|
|00000b20| 69 73 70 61 74 63 68 2d | 3e 49 6e 63 6c 65 76 65 |ispatch-|>Incleve|
|00000b30| 6c 28 29 3b 0d 0a 20 20 | 20 56 4d 61 74 72 69 78 |l();..  | VMatrix|
|00000b40| 20 63 65 6e 74 65 72 65 | 64 2c 20 7a 2c 20 63 6f | centere|d, z, co|
|00000b50| 76 61 72 3b 0d 0a 20 20 | 20 64 6f 75 62 6c 65 20 |var;..  | double |
|00000b60| 6e 20 3d 20 28 64 6f 75 | 62 6c 65 29 20 52 2e 72 |n = (dou|ble) R.r|
|00000b70| 3b 0d 0a 20 20 20 2f 2f | 20 63 65 6e 74 65 72 20 |;..   //| center |
|00000b80| 61 20 63 6f 6c 75 6d 6e | 20 76 65 63 74 6f 72 0d |a column| vector.|
|00000b90| 0a 20 20 20 63 65 6e 74 | 65 72 65 64 20 3d 20 52 |.   cent|ered = R|
|00000ba0| 20 2d 20 53 75 6d 28 52 | 29 2e 6d 28 31 2c 20 31 | - Sum(R|).m(1, 1|
|00000bb0| 29 20 2f 20 6e 3b 0d 0a | 20 20 20 2f 2f 20 7a 65 |) / n;..|   // ze|
|00000bc0| 72 6f 20 70 61 64 20 74 | 6f 20 6c 65 6e 67 74 68 |ro pad t|o length|
|00000bd0| 20 32 6e 3a 20 32 6e 20 | 70 65 72 69 6f 64 69 63 | 2n: 2n |periodic|
|00000be0| 20 66 6f 72 20 66 75 6c | 6c 0d 0a 20 20 20 2f 2f | for ful|l..   //|
|00000bf0| 20 73 61 6d 70 6c 65 20 | 73 70 65 63 74 72 61 6c | sample |spectral|
|00000c00| 20 64 65 6e 73 69 74 79 | 0d 0a 20 20 20 63 65 6e | density|..   cen|
|00000c10| 74 65 72 65 64 20 3d 20 | 43 76 28 63 65 6e 74 65 |tered = |Cv(cente|
|00000c20| 72 65 64 2c 20 46 69 6c | 6c 28 52 2e 72 2c 20 52 |red, Fil|l(R.r, R|
|00000c30| 2e 63 2c 20 30 29 29 3b | 0d 0a 20 20 20 2f 2f 20 |.c, 0));|..   // |
|00000c40| 74 61 6b 65 20 66 66 74 | 0d 0a 20 20 20 7a 20 3d |take fft|..   z =|
|00000c50| 20 46 66 74 28 63 65 6e | 74 65 72 65 64 29 3b 0d | Fft(cen|tered);.|
|00000c60| 0a 20 20 20 2f 2f 20 74 | 61 6b 65 20 63 6f 6e 76 |.   // t|ake conv|
|00000c70| 6f 6c 75 74 69 6f 6e 20 | 3a 20 67 69 76 65 73 20 |olution |: gives |
|00000c80| 73 61 6d 70 6c 65 20 73 | 70 65 63 74 72 61 6c 20 |sample s|pectral |
|00000c90| 64 65 6e 73 69 74 79 0d | 0a 20 20 20 73 70 65 63 |density.|.   spec|
|00000ca0| 74 64 65 6e 73 20 3d 20 | 53 75 6d 28 7a 20 25 20 |tdens = |Sum(z % |
|00000cb0| 7a 20 2f 20 6e 2c 20 43 | 4f 4c 55 4d 4e 53 29 3b |z / n, C|OLUMNS);|
|00000cc0| 0d 0a 20 20 20 2f 2f 20 | 69 6e 76 65 72 73 65 20 |..   // |inverse |
|00000cd0| 66 66 74 20 66 6f 72 20 | 73 65 72 69 61 6c 20 63 |fft for |serial c|
|00000ce0| 6f 72 72 65 6c 61 74 69 | 6f 6e 20 28 61 75 74 6f |orrelati|on (auto|
|00000cf0| 63 6f 72 72 65 6c 61 74 | 69 6f 6e 29 0d 0a 20 20 |correlat|ion)..  |
|00000d00| 20 63 6f 76 61 72 20 3d | 20 46 66 74 28 73 70 65 | covar =| Fft(spe|
|00000d10| 63 74 64 65 6e 73 2c 20 | 46 46 41 4c 53 45 29 3b |ctdens, |FFALSE);|
|00000d20| 0d 0a 20 20 20 2f 2f 20 | 74 68 72 6f 77 20 61 77 |..   // |throw aw|
|00000d30| 61 79 20 6c 61 73 74 20 | 68 61 6c 66 2e 0d 0a 20 |ay last |half... |
|00000d40| 20 20 63 6f 76 61 72 20 | 3d 20 53 75 62 6d 61 74 |  covar |= Submat|
|00000d50| 28 63 6f 76 61 72 2c 20 | 52 2e 72 2c 20 32 29 3b |(covar, |R.r, 2);|
|00000d60| 0d 0a 20 20 20 44 69 73 | 70 61 74 63 68 2d 3e 50 |..   Dis|patch->P|
|00000d70| 75 73 68 28 63 6f 76 61 | 72 29 3b 0d 0a 20 20 20 |ush(cova|r);..   |
|00000d80| 72 65 74 75 72 6e 20 44 | 69 73 70 61 74 63 68 2d |return D|ispatch-|
|00000d90| 3e 44 65 63 52 65 74 75 | 72 6e 28 29 3b 0d 0a 7d |>DecRetu|rn();..}|
|00000da0| 0d 0a 0d 0a 76 6f 69 64 | 20 70 6c 6f 74 43 6f 72 |....void| plotCor|
|00000db0| 72 65 6c 6f 67 72 61 6d | 28 56 4d 61 74 72 69 78 |relogram|(VMatrix|
|00000dc0| 20 26 73 65 72 69 61 6c | 29 0d 0a 7b 0d 0a 20 20 | &serial|)..{..  |
|00000dd0| 20 44 69 73 70 61 74 63 | 68 2d 3e 49 6e 63 6c 65 | Dispatc|h->Incle|
|00000de0| 76 65 6c 28 29 3b 0d 0a | 20 20 20 0d 0a 20 20 20 |vel();..|   ..   |
|00000df0| 69 6e 74 20 6e 20 3d 20 | 73 65 72 69 61 6c 2e 72 |int n = |serial.r|
|00000e00| 3b 0d 0a 20 20 20 64 6f | 75 62 6c 65 20 73 69 67 |;..   do|uble sig|
|00000e10| 6d 61 20 3d 20 73 65 72 | 69 61 6c 2e 6d 28 31 2c |ma = ser|ial.m(1,|
|00000e20| 20 31 29 3b 0d 0a 20 20 | 20 56 4d 61 74 72 69 78 | 1);..  | VMatrix|
|00000e30| 20 43 6f 72 72 65 6c 6f | 67 72 61 6d 20 3d 20 73 | Correlo|gram = s|
|00000e40| 65 72 69 61 6c 20 2f 20 | 73 69 67 6d 61 3b 0d 0a |erial / |sigma;..|
|00000e50| 20 20 20 43 6f 72 72 65 | 6c 6f 67 72 61 6d 20 3d |   Corre|logram =|
|00000e60| 20 53 75 62 6d 61 74 28 | 43 6f 72 72 65 6c 6f 67 | Submat(|Correlog|
|00000e70| 72 61 6d 2c 20 6e 2c 20 | 31 2c 20 32 2c 20 31 29 |ram, n, |1, 2, 1)|
|00000e80| 3b 0d 0a 20 20 20 56 4d | 61 74 72 69 78 20 67 72 |;..   VM|atrix gr|
|00000e90| 61 66 20 3d 20 43 68 28 | 49 6e 64 65 78 28 31 2c |af = Ch(|Index(1,|
|00000ea0| 20 43 6f 72 72 65 6c 6f | 67 72 61 6d 2e 72 29 2c | Correlo|gram.r),|
|00000eb0| 20 43 6f 72 72 65 6c 6f | 67 72 61 6d 29 3b 0d 0a | Correlo|gram);..|
|00000ec0| 20 20 20 0d 0a 20 20 20 | 47 4d 61 74 72 69 78 20 |   ..   |GMatrix |
|00000ed0| 41 67 72 61 70 68 28 67 | 72 61 66 29 3b 0d 0a 20 |Agraph(g|raf);.. |
|00000ee0| 20 20 2a 41 67 72 61 70 | 68 2e 50 61 74 68 54 6f |  *Agrap|h.PathTo|
|00000ef0| 44 72 69 76 65 72 20 3d | 20 22 43 3a 5c 5c 74 63 |Driver =| "C:\\tc|
|00000f00| 5c 5c 62 67 69 22 3b 0d | 0a 20 20 20 2a 41 67 72 |\\bgi";.|.   *Agr|
|00000f10| 61 70 68 2e 74 69 74 6c | 65 20 3d 20 22 53 65 72 |aph.titl|e = "Ser|
|00000f20| 69 61 6c 20 43 6f 72 72 | 65 6c 61 74 69 6f 6e 73 |ial Corr|elations|
|00000f30| 22 3b 0d 0a 20 20 20 2a | 41 67 72 61 70 68 2e 74 |";..   *|Agraph.t|
|00000f40| 69 74 6c 65 32 20 3d 20 | 22 66 6f 72 20 73 61 6d |itle2 = |"for sam|
|00000f50| 70 6c 65 20 72 65 73 69 | 64 75 61 6c 73 22 3b 0d |ple resi|duals";.|
|00000f60| 0a 20 20 20 2a 41 67 72 | 61 70 68 2e 79 6e 61 6d |.   *Agr|aph.ynam|
|00000f70| 65 20 3d 20 22 53 65 72 | 69 61 6c 20 63 6f 72 72 |e = "Ser|ial corr|
|00000f80| 65 6c 61 74 69 6f 6e 73 | 22 3b 0d 0a 20 20 20 2a |elations|";..   *|
|00000f90| 41 67 72 61 70 68 2e 78 | 6e 61 6d 65 20 3d 20 22 |Agraph.x|name = "|
|00000fa0| 4c 61 67 73 22 3b 0d 0a | 20 20 20 41 67 72 61 70 |Lags";..|   Agrap|
|00000fb0| 68 2e 48 72 65 66 28 30 | 2e 30 29 3b 0d 0a 20 20 |h.Href(0|.0);..  |
|00000fc0| 20 41 67 72 61 70 68 2e | 53 68 6f 77 28 29 3b 0d | Agraph.|Show();.|
|00000fd0| 0a 20 20 20 0d 0a 20 20 | 20 44 69 73 70 61 74 63 |.   ..  | Dispatc|
|00000fe0| 68 2d 3e 43 6c 65 61 6e | 73 74 61 63 6b 28 29 3b |h->Clean|stack();|
|00000ff0| 0d 0a 20 20 20 44 69 73 | 70 61 74 63 68 2d 3e 44 |..   Dis|patch->D|
|00001000| 65 63 6c 65 76 65 6c 28 | 29 3b 0d 0a 7d 0d 0a 0d |eclevel(|);..}...|
|00001010| 0a 76 6f 69 64 20 70 6c | 6f 74 50 65 72 69 6f 64 |.void pl|otPeriod|
|00001020| 6f 67 72 61 6d 28 56 4d | 61 74 72 69 78 20 26 73 |ogram(VM|atrix &s|
|00001030| 70 65 63 74 64 65 6e 73 | 29 0d 0a 7b 0d 0a 20 20 |pectdens|)..{..  |
|00001040| 20 2f 2f 20 63 61 6c 63 | 75 6c 61 74 65 20 61 20 | // calc|ulate a |
|00001050| 73 74 61 6e 64 61 72 64 | 69 7a 65 64 20 70 65 72 |standard|ized per|
|00001060| 69 6f 64 6f 67 72 61 6d | 20 6f 6e 20 74 68 65 20 |iodogram| on the |
|00001070| 6c 6f 67 20 73 63 61 6c | 65 0d 0a 20 20 20 44 69 |log scal|e..   Di|
|00001080| 73 70 61 74 63 68 2d 3e | 49 6e 63 6c 65 76 65 6c |spatch->|Inclevel|
|00001090| 28 29 3b 0d 0a 20 20 20 | 69 6e 74 20 6e 20 3d 20 |();..   |int n = |
|000010a0| 73 70 65 63 74 64 65 6e | 73 2e 72 3b 0d 0a 20 20 |spectden|s.r;..  |
|000010b0| 20 2f 2f 20 74 68 69 73 | 20 77 6f 72 6b 73 20 62 | // this| works b|
|000010c0| 65 63 61 75 73 65 20 64 | 61 74 61 20 69 73 20 61 |ecause d|ata is a|
|000010d0| 6c 72 65 61 64 79 20 63 | 65 6e 74 65 72 65 64 2c |lready c|entered,|
|000010e0| 20 77 68 69 63 68 0d 0a | 20 20 20 2f 2f 20 66 6f | which..|   // fo|
|000010f0| 72 63 65 73 20 73 70 65 | 63 74 64 65 6e 73 2e 6d |rces spe|ctdens.m|
|00001100| 28 31 2c 31 29 20 3d 20 | 30 2e 30 3b 0d 0a 20 20 |(1,1) = |0.0;..  |
|00001110| 20 64 6f 75 62 6c 65 20 | 73 69 67 6d 61 20 3d 20 | double |sigma = |
|00001120| 53 75 6d 28 73 70 65 63 | 74 64 65 6e 73 29 2e 6d |Sum(spec|tdens).m|
|00001130| 28 31 2c 20 31 29 20 2f | 20 28 28 64 6f 75 62 6c |(1, 1) /| ((doubl|
|00001140| 65 29 20 6e 29 3b 0d 0a | 20 20 20 56 4d 61 74 72 |e) n);..|   VMatr|
|00001150| 69 78 20 50 65 72 69 6f | 64 6f 67 72 61 6d 20 3d |ix Perio|dogram =|
|00001160| 20 73 70 65 63 74 64 65 | 6e 73 20 2f 20 73 69 67 | spectde|ns / sig|
|00001170| 6d 61 3b 0d 0a 20 20 20 | 50 65 72 69 6f 64 6f 67 |ma;..   |Periodog|
|00001180| 72 61 6d 20 3d 20 4d 6c | 6f 67 28 53 75 62 6d 61 |ram = Ml|og(Subma|
|00001190| 74 28 50 65 72 69 6f 64 | 6f 67 72 61 6d 2c 20 6e |t(Period|ogram, n|
|000011a0| 20 2f 20 32 20 2b 20 31 | 2c 20 31 2c 20 32 2c 20 | / 2 + 1|, 1, 2, |
|000011b0| 31 29 29 3b 0d 0a 20 20 | 20 0d 0a 20 20 20 2f 2f |1));..  | ..   //|
|000011c0| 20 66 72 65 71 75 65 6e | 63 69 65 73 0d 0a 20 20 | frequen|cies..  |
|000011d0| 20 64 6f 75 62 6c 65 20 | 64 6e 20 3d 28 64 6f 75 | double |dn =(dou|
|000011e0| 62 6c 65 29 20 28 6e 20 | 2f 20 32 20 2b 20 31 29 |ble) (n |/ 2 + 1)|
|000011f0| 3b 0d 0a 20 20 20 56 4d | 61 74 72 69 78 20 67 72 |;..   VM|atrix gr|
|00001200| 61 66 20 3d 20 43 68 28 | 49 6e 64 65 78 28 31 2c |af = Ch(|Index(1,|
|00001210| 20 50 65 72 69 6f 64 6f | 67 72 61 6d 2e 72 29 20 | Periodo|gram.r) |
|00001220| 2f 20 64 6e 2c 20 50 65 | 72 69 6f 64 6f 67 72 61 |/ dn, Pe|riodogra|
|00001230| 6d 29 3b 0d 0a 20 20 20 | 0d 0a 20 20 20 47 4d 61 |m);..   |..   GMa|
|00001240| 74 72 69 78 20 41 67 72 | 61 70 68 28 67 72 61 66 |trix Agr|aph(graf|
|00001250| 29 3b 0d 0a 20 20 20 2a | 41 67 72 61 70 68 2e 50 |);..   *|Agraph.P|
|00001260| 61 74 68 54 6f 44 72 69 | 76 65 72 20 3d 20 22 43 |athToDri|ver = "C|
|00001270| 3a 5c 5c 74 63 5c 5c 62 | 67 69 22 3b 0d 0a 20 20 |:\\tc\\b|gi";..  |
|00001280| 20 2a 41 67 72 61 70 68 | 2e 74 69 74 6c 65 20 3d | *Agraph|.title =|
|00001290| 20 22 50 65 72 69 6f 64 | 6f 67 72 61 6d 22 3b 0d | "Period|ogram";.|
|000012a0| 0a 20 20 20 2a 41 67 72 | 61 70 68 2e 74 69 74 6c |.   *Agr|aph.titl|
|000012b0| 65 32 20 3d 20 22 66 6f | 72 20 73 61 6d 70 6c 65 |e2 = "fo|r sample|
|000012c0| 20 72 65 73 69 64 75 61 | 6c 73 22 3b 0d 0a 20 20 | residua|ls";..  |
|000012d0| 20 2a 41 67 72 61 70 68 | 2e 79 6e 61 6d 65 20 3d | *Agraph|.yname =|
|000012e0| 20 22 4c 6f 67 20 70 65 | 72 69 6f 64 6f 67 72 61 | "Log pe|riodogra|
|000012f0| 6d 22 3b 0d 0a 20 20 20 | 2a 41 67 72 61 70 68 2e |m";..   |*Agraph.|
|00001300| 78 6e 61 6d 65 20 3d 20 | 22 46 72 65 71 75 65 6e |xname = |"Frequen|
|00001310| 63 69 65 73 22 3b 0d 0a | 20 20 20 41 67 72 61 70 |cies";..|   Agrap|
|00001320| 68 2e 56 72 65 66 28 30 | 2e 32 35 29 3b 0d 0a 20 |h.Vref(0|.25);.. |
|00001330| 20 20 41 67 72 61 70 68 | 2e 53 68 6f 77 28 29 3b |  Agraph|.Show();|
|00001340| 0d 0a 20 20 20 0d 0a 20 | 20 20 44 69 73 70 61 74 |..   .. |  Dispat|
|00001350| 63 68 2d 3e 43 6c 65 61 | 6e 73 74 61 63 6b 28 29 |ch->Clea|nstack()|
|00001360| 3b 0d 0a 20 20 20 44 69 | 73 70 61 74 63 68 2d 3e |;..   Di|spatch->|
|00001370| 44 65 63 6c 65 76 65 6c | 28 29 3b 0d 0a 7d 0d 0a |Declevel|();..}..|
|00001380| 0d 0a 0d 0a 6d 61 69 6e | 28 29 0d 0a 7b 0d 0a 20 |....main|()..{.. |
|00001390| 20 20 20 20 20 44 69 73 | 70 61 74 63 68 2d 3e 49 |     Dis|patch->I|
|000013a0| 6e 63 6c 65 76 65 6c 28 | 29 3b 0d 0a 20 20 20 20 |nclevel(|);..    |
|000013b0| 20 20 56 4d 61 74 72 69 | 78 20 78 2c 20 62 65 74 |  VMatri|x x, bet|
|000013c0| 61 28 22 62 65 74 61 22 | 2c 32 2c 31 29 2c 20 79 |a("beta"|,2,1), y|
|000013d0| 2c 20 62 65 74 61 68 61 | 74 2c 20 72 65 73 69 64 |, betaha|t, resid|
|000013e0| 73 2c 20 73 65 72 69 61 | 6c 3b 0d 0a 20 20 20 20 |s, seria|l;..    |
|000013f0| 20 20 53 65 74 77 69 64 | 28 31 35 29 3b 0d 0a 20 |  Setwid|(15);.. |
|00001400| 20 20 20 20 20 53 65 74 | 64 65 63 28 31 30 29 3b |     Set|dec(10);|
|00001410| 0d 0a 0d 0a 20 20 20 20 | 20 20 62 65 74 61 2e 4d |....    |  beta.M|
|00001420| 28 31 2c 31 29 20 3d 20 | 31 3b 0d 0a 20 20 20 20 |(1,1) = |1;..    |
|00001430| 20 20 62 65 74 61 2e 4d | 28 32 2c 31 29 20 3d 20 |  beta.M|(2,1) = |
|00001440| 30 2e 35 3b 0d 0a 0d 0a | 20 20 20 20 20 20 78 20 |0.5;....|      x |
|00001450| 3d 20 67 65 74 78 28 20 | 31 32 38 20 29 3b 0d 0a |= getx( |128 );..|
|00001460| 20 20 20 20 20 20 79 20 | 3d 20 67 65 74 79 28 78 |      y |= gety(x|
|00001470| 2c 62 65 74 61 29 3b 0d | 0a 0d 0a 20 20 20 20 20 |,beta);.|...     |
|00001480| 20 62 65 74 61 68 61 74 | 20 3d 20 72 65 67 72 65 | betahat| = regre|
|00001490| 73 73 69 6f 6e 28 78 2c | 79 29 3b 0d 0a 20 20 20 |ssion(x,|y);..   |
|000014a0| 20 20 20 62 65 74 61 68 | 61 74 2e 4e 61 6d 65 69 |   betah|at.Namei|
|000014b0| 74 28 20 22 54 65 78 74 | 20 62 6f 6f 6b 20 62 65 |t( "Text| book be|
|000014c0| 74 61 68 61 74 22 29 3b | 0d 0a 20 20 20 20 20 20 |tahat");|..      |
|000014d0| 28 54 72 61 6e 28 62 65 | 74 61 29 29 2e 44 69 73 |(Tran(be|ta)).Dis|
|000014e0| 70 6c 61 79 4d 61 74 28 | 29 3b 0d 0a 20 20 20 20 |playMat(|);..    |
|000014f0| 20 20 28 54 72 61 6e 28 | 62 65 74 61 68 61 74 29 |  (Tran(|betahat)|
|00001500| 29 2e 44 69 73 70 6c 61 | 79 4d 61 74 28 29 3b 0d |).Displa|yMat();.|
|00001510| 0a 0d 0a 20 20 20 20 20 | 20 72 65 73 69 64 73 20 |...     | resids |
|00001520| 3d 20 79 20 2d 20 78 2a | 62 65 74 61 68 61 74 3b |= y - x*|betahat;|
|00001530| 0d 0a 20 20 20 20 20 20 | 70 6c 6f 74 52 65 73 69 |..      |plotResi|
|00001540| 64 75 61 6c 73 28 20 72 | 65 73 69 64 73 20 29 3b |duals( r|esids );|
|00001550| 0d 0a 0d 0a 20 20 20 20 | 20 20 56 4d 61 74 72 69 |....    |  VMatri|
|00001560| 78 20 73 70 65 63 74 64 | 65 6e 73 3b 0d 0a 20 20 |x spectd|ens;..  |
|00001570| 20 20 20 20 73 65 72 69 | 61 6c 20 3d 20 47 65 74 |    seri|al = Get|
|00001580| 53 65 72 69 61 6c 43 6f | 76 61 72 28 20 72 65 73 |SerialCo|var( res|
|00001590| 69 64 73 2c 20 73 70 65 | 63 74 64 65 6e 73 20 29 |ids, spe|ctdens )|
|000015a0| 3b 0d 0a 0d 0a 20 20 20 | 20 20 20 70 6c 6f 74 43 |;....   |   plotC|
|000015b0| 6f 72 72 65 6c 6f 67 72 | 61 6d 28 20 73 65 72 69 |orrelogr|am( seri|
|000015c0| 61 6c 20 29 3b 0d 0a 20 | 20 20 20 20 20 70 6c 6f |al );.. |     plo|
|000015d0| 74 50 65 72 69 6f 64 6f | 67 72 61 6d 28 20 73 70 |tPeriodo|gram( sp|
|000015e0| 65 63 74 64 65 6e 73 29 | 3b 0d 0a 0d 0a 20 20 20 |ectdens)|;....   |
|000015f0| 76 63 6c 6f 73 65 28 29 | 3b 0d 0a 7d 0d 0a 0d 0a |vclose()|;..}....|
|00001600| 5c 65 6e 64 7b 76 65 72 | 62 61 74 69 6d 7d       |\end{ver|batim}  |
+--------+-------------------------+-------------------------+--------+--------+