home *** CD-ROM | disk | FTP | other *** search
open in: 
MacOS 8.1
|
Win98
|
DOS
view JSON data
|
view as text
This file was processed as: SHell self-extracting ARchive
(archive/shar).
You can browse this item here: 2398
| Confidence | Program | Detection | Match Type | Support
|
|---|
| 100%
| dexvert
| Newsgroup Content (archive/news)
| magic
| Supported |
| 100%
| dexvert
| SHell self-extracting ARchive (archive/shar)
| magic
| Supported |
| 100%
| dexvert
| Internet Message Format (text/imf)
| magic
| Supported |
| 1%
| dexvert
| Text File (text/txt)
| fallback
| Supported |
| 100%
| file
| news or mail, ASCII text
| default
| |
| 100%
| TrID
| E-Mail message (Var. 2)
| default
| |
| 100%
| checkBytes
| Printable ASCII
| default
| |
| 100%
| perlTextCheck
| Likely Text (Perl)
| default
| |
| 100%
| siegfried
| fmt/329 Shell Archive Format
| default
| |
| 100%
| detectItEasy
| Format: plain text[LF]
| default (weak)
| |
| 100%
| xdgMime
| message/rfc822
| default
|
|
hex view+--------+-------------------------+-------------------------+--------+--------+
|00000000| 46 72 6f 6d 3a 20 61 72 | 67 40 63 63 76 72 31 2e |From: ar|g@ccvr1.|
|00000010| 6e 63 73 75 2e 65 64 75 | 20 28 52 6f 6e 20 47 61 |ncsu.edu| (Ron Ga|
|00000020| 6c 6c 61 6e 74 29 0a 4e | 65 77 73 67 72 6f 75 70 |llant).N|ewsgroup|
|00000030| 73 3a 20 61 6c 74 2e 73 | 6f 75 72 63 65 73 0a 53 |s: alt.s|ources.S|
|00000040| 75 62 6a 65 63 74 3a 20 | 6e 6c 6d 64 6c 20 50 61 |ubject: |nlmdl Pa|
|00000050| 72 74 20 31 20 6f 66 20 | 36 0a 4d 65 73 73 61 67 |rt 1 of |6.Messag|
|00000060| 65 2d 49 44 3a 20 3c 31 | 39 39 30 44 65 63 32 32 |e-ID: <1|990Dec22|
|00000070| 2e 32 32 35 36 34 33 2e | 39 34 31 32 40 6e 63 73 |.225643.|9412@ncs|
|00000080| 75 76 78 2e 6e 63 73 75 | 2e 65 64 75 3e 0a 44 61 |uvx.ncsu|.edu>.Da|
|00000090| 74 65 3a 20 32 32 20 44 | 65 63 20 39 30 20 32 32 |te: 22 D|ec 90 22|
|000000a0| 3a 35 36 3a 34 33 20 47 | 4d 54 0a 0a 23 21 20 2f |:56:43 G|MT..#! /|
|000000b0| 62 69 6e 2f 73 68 0a 23 | 20 54 68 69 73 20 69 73 |bin/sh.#| This is|
|000000c0| 20 61 20 73 68 65 6c 6c | 20 61 72 63 68 69 76 65 | a shell| archive|
|000000d0| 2e 20 20 52 65 6d 6f 76 | 65 20 61 6e 79 74 68 69 |. Remov|e anythi|
|000000e0| 6e 67 20 62 65 66 6f 72 | 65 20 74 68 69 73 20 6c |ng befor|e this l|
|000000f0| 69 6e 65 2c 20 74 68 65 | 6e 20 75 6e 70 61 63 6b |ine, the|n unpack|
|00000100| 0a 23 20 69 74 20 62 79 | 20 73 61 76 69 6e 67 20 |.# it by| saving |
|00000110| 69 74 20 69 6e 74 6f 20 | 61 20 66 69 6c 65 20 61 |it into |a file a|
|00000120| 6e 64 20 74 79 70 69 6e | 67 20 22 73 68 20 66 69 |nd typin|g "sh fi|
|00000130| 6c 65 22 2e 20 20 54 6f | 20 6f 76 65 72 77 72 69 |le". To| overwri|
|00000140| 74 65 20 65 78 69 73 74 | 69 6e 67 0a 23 20 66 69 |te exist|ing.# fi|
|00000150| 6c 65 73 2c 20 74 79 70 | 65 20 22 73 68 20 66 69 |les, typ|e "sh fi|
|00000160| 6c 65 20 2d 63 22 2e 20 | 20 59 6f 75 20 63 61 6e |le -c". | You can|
|00000170| 20 61 6c 73 6f 20 66 65 | 65 64 20 74 68 69 73 20 | also fe|ed this |
|00000180| 61 73 20 73 74 61 6e 64 | 61 72 64 20 69 6e 70 75 |as stand|ard inpu|
|00000190| 74 20 76 69 61 0a 23 20 | 75 6e 73 68 61 72 2c 20 |t via.# |unshar, |
|000001a0| 6f 72 20 62 79 20 74 79 | 70 69 6e 67 20 22 73 68 |or by ty|ping "sh|
|000001b0| 20 3c 66 69 6c 65 22 2c | 20 65 2e 67 2e 2e 20 20 | <file",| e.g.. |
|000001c0| 49 66 20 74 68 69 73 20 | 61 72 63 68 69 76 65 20 |If this |archive |
|000001d0| 69 73 20 63 6f 6d 70 6c | 65 74 65 2c 20 79 6f 75 |is compl|ete, you|
|000001e0| 0a 23 20 77 69 6c 6c 20 | 73 65 65 20 74 68 65 20 |.# will |see the |
|000001f0| 66 6f 6c 6c 6f 77 69 6e | 67 20 6d 65 73 73 61 67 |followin|g messag|
|00000200| 65 20 61 74 20 74 68 65 | 20 65 6e 64 3a 0a 23 09 |e at the| end:.#.|
|00000210| 09 22 45 6e 64 20 6f 66 | 20 61 72 63 68 69 76 65 |."End of| archive|
|00000220| 20 31 20 28 6f 66 20 36 | 29 2e 22 0a 23 20 43 6f | 1 (of 6|).".# Co|
|00000230| 6e 74 65 6e 74 73 3a 20 | 20 4d 61 6b 65 66 69 6c |ntents: | Makefil|
|00000240| 65 20 52 45 41 44 4d 45 | 20 63 68 31 65 67 31 20 |e README| ch1eg1 |
|00000250| 63 68 31 65 67 31 2f 61 | 6d 73 74 61 74 2e 64 61 |ch1eg1/a|mstat.da|
|00000260| 74 20 63 68 31 65 67 31 | 2f 65 6e 64 69 6e 67 2e |t ch1eg1|/ending.|
|00000270| 62 61 6b 0a 23 20 20 20 | 63 68 31 65 67 31 2f 6d |bak.# |ch1eg1/m|
|00000280| 6f 64 65 6c 2e 63 63 20 | 63 68 31 65 67 31 2f 6d |odel.cc |ch1eg1/m|
|00000290| 6f 64 65 6c 2e 68 20 63 | 68 31 65 67 31 2f 73 74 |odel.h c|h1eg1/st|
|000002a0| 61 72 74 69 6e 67 2e 64 | 61 74 20 63 68 35 65 67 |arting.d|at ch5eg|
|000002b0| 31 0a 23 20 20 20 63 68 | 35 65 67 31 2f 65 6c 65 |1.# ch|5eg1/ele|
|000002c0| 63 74 62 2e 64 61 74 20 | 63 68 35 65 67 31 2f 6d |ctb.dat |ch5eg1/m|
|000002d0| 6f 64 65 6c 2e 68 20 63 | 68 35 65 67 31 2f 73 74 |odel.h c|h5eg1/st|
|000002e0| 61 72 74 69 6e 67 2e 64 | 61 74 20 63 68 36 65 67 |arting.d|at ch6eg|
|000002f0| 31 0a 23 20 20 20 63 68 | 36 65 67 31 2f 52 45 41 |1.# ch|6eg1/REA|
|00000300| 44 4d 45 20 63 68 36 65 | 67 31 2f 73 74 61 72 74 |DME ch6e|g1/start|
|00000310| 69 6e 67 2e 64 61 74 20 | 63 68 36 65 67 32 20 63 |ing.dat |ch6eg2 c|
|00000320| 68 36 65 67 32 2f 65 6e | 64 69 6e 67 2e 62 61 6b |h6eg2/en|ding.bak|
|00000330| 0a 23 20 20 20 63 68 36 | 65 67 32 2f 6d 6f 64 65 |.# ch6|eg2/mode|
|00000340| 6c 2e 63 63 20 63 68 36 | 65 67 32 2f 6d 6f 64 65 |l.cc ch6|eg2/mode|
|00000350| 6c 2e 68 20 63 68 36 65 | 67 32 2f 73 74 61 72 74 |l.h ch6e|g2/start|
|00000360| 33 6c 2e 64 61 74 0a 23 | 20 20 20 63 68 36 65 67 |3l.dat.#| ch6eg|
|00000370| 32 2f 73 74 61 72 74 69 | 6e 67 2e 64 61 74 20 64 |2/starti|ng.dat d|
|00000380| 63 6f 6e 64 2e 63 63 20 | 64 73 77 65 65 70 2e 63 |cond.cc |dsweep.c|
|00000390| 63 20 6e 6c 6d 64 6c 2e | 63 63 20 6e 6c 6d 64 6c |c nlmdl.|cc nlmdl|
|000003a0| 2e 68 20 73 74 61 72 62 | 6f 78 2e 63 63 0a 23 20 |.h starb|ox.cc.# |
|000003b0| 20 20 73 74 61 74 75 73 | 2e 68 20 74 6f 6f 6c 73 | status|.h tools|
|000003c0| 2e 68 20 75 73 75 61 6c | 2e 68 0a 23 20 57 72 61 |.h usual|.h.# Wra|
|000003d0| 70 70 65 64 20 62 79 20 | 61 72 67 40 73 70 61 72 |pped by |arg@spar|
|000003e0| 63 20 6f 6e 20 53 61 74 | 20 44 65 63 20 32 32 20 |c on Sat| Dec 22 |
|000003f0| 31 35 3a 32 35 3a 30 37 | 20 31 39 39 30 0a 50 41 |15:25:07| 1990.PA|
|00000400| 54 48 3d 2f 62 69 6e 3a | 2f 75 73 72 2f 62 69 6e |TH=/bin:|/usr/bin|
|00000410| 3a 2f 75 73 72 2f 75 63 | 62 20 3b 20 65 78 70 6f |:/usr/uc|b ; expo|
|00000420| 72 74 20 50 41 54 48 0a | 69 66 20 74 65 73 74 20 |rt PATH.|if test |
|00000430| 2d 66 20 27 4d 61 6b 65 | 66 69 6c 65 27 20 2d 61 |-f 'Make|file' -a|
|00000440| 20 22 24 7b 31 7d 22 20 | 21 3d 20 22 2d 63 22 20 | "${1}" |!= "-c" |
|00000450| 3b 20 74 68 65 6e 20 0a | 20 20 65 63 68 6f 20 73 |; then .| echo s|
|00000460| 68 61 72 3a 20 57 69 6c | 6c 20 6e 6f 74 20 63 6c |har: Wil|l not cl|
|00000470| 6f 62 62 65 72 20 65 78 | 69 73 74 69 6e 67 20 66 |obber ex|isting f|
|00000480| 69 6c 65 20 5c 22 27 4d | 61 6b 65 66 69 6c 65 27 |ile \"'M|akefile'|
|00000490| 5c 22 0a 65 6c 73 65 0a | 65 63 68 6f 20 73 68 61 |\".else.|echo sha|
|000004a0| 72 3a 20 45 78 74 72 61 | 63 74 69 6e 67 20 5c 22 |r: Extra|cting \"|
|000004b0| 27 4d 61 6b 65 66 69 6c | 65 27 5c 22 20 5c 28 31 |'Makefil|e'\" \(1|
|000004c0| 34 38 33 20 63 68 61 72 | 61 63 74 65 72 73 5c 29 |483 char|acters\)|
|000004d0| 0a 73 65 64 20 22 73 2f | 5e 58 2f 2f 22 20 3e 27 |.sed "s/|^X//" >'|
|000004e0| 4d 61 6b 65 66 69 6c 65 | 27 20 3c 3c 27 45 4e 44 |Makefile|' <<'END|
|000004f0| 5f 4f 46 5f 46 49 4c 45 | 27 0a 58 43 43 20 20 20 |_OF_FILE|'.XCC |
|00000500| 20 20 20 20 3d 20 67 2b | 2b 0a 58 43 46 4c 41 47 | = g+|+.XCFLAG|
|00000510| 53 20 20 20 3d 20 2d 4f | 0a 58 52 4d 44 49 52 20 |S = -O|.XRMDIR |
|00000520| 20 20 20 3d 20 2e 2e 0a | 58 4e 4c 44 49 52 20 20 | = ...|XNLDIR |
|00000530| 20 20 3d 20 2e 2e 0a 58 | 4d 44 4c 44 49 52 20 20 | = ...X|MDLDIR |
|00000540| 20 3d 20 2e 0a 58 4c 44 | 46 4c 41 47 53 20 20 3d | = ..XLD|FLAGS =|
|00000550| 20 2d 6c 6d 0a 58 0a 58 | 6e 6c 6d 64 6c 09 3a 20 | -lm.X.X|nlmdl.: |
|00000560| 6e 6c 6d 64 6c 2e 6f 20 | 73 74 61 74 75 73 2e 6f |nlmdl.o |status.o|
|00000570| 20 6d 6f 64 65 6c 2e 6f | 20 5c 0a 58 09 09 72 65 | model.o| \.X..re|
|00000580| 61 6c 6d 61 74 2e 6f 20 | 64 63 6f 6e 64 2e 6f 20 |almat.o |dcond.o |
|00000590| 64 73 77 65 65 70 2e 6f | 20 73 74 61 72 62 6f 78 |dsweep.o| starbox|
|000005a0| 2e 6f 0a 58 09 24 28 43 | 43 29 20 2d 6f 20 6e 6c |.o.X.$(C|C) -o nl|
|000005b0| 6d 64 6c 20 6e 6c 6d 64 | 6c 2e 6f 20 73 74 61 74 |mdl nlmd|l.o stat|
|000005c0| 75 73 2e 6f 20 6d 6f 64 | 65 6c 2e 6f 20 5c 0a 58 |us.o mod|el.o \.X|
|000005d0| 09 09 72 65 61 6c 6d 61 | 74 2e 6f 20 64 63 6f 6e |..realma|t.o dcon|
|000005e0| 64 2e 6f 20 64 73 77 65 | 65 70 2e 6f 20 73 74 61 |d.o dswe|ep.o sta|
|000005f0| 72 62 6f 78 2e 6f 20 24 | 28 4c 44 46 4c 41 47 53 |rbox.o $|(LDFLAGS|
|00000600| 29 20 0a 58 0a 58 6e 6c | 6d 64 6c 2e 6f 09 3a 20 |) .X.Xnl|mdl.o.: |
|00000610| 24 28 4e 4c 44 49 52 29 | 2f 6e 6c 6d 64 6c 2e 63 |$(NLDIR)|/nlmdl.c|
|00000620| 63 20 24 28 4e 4c 44 49 | 52 29 2f 6e 6c 6f 70 72 |c $(NLDI|R)/nlopr|
|00000630| 2e 63 63 20 5c 0a 58 09 | 09 24 28 4d 44 4c 44 49 |.cc \.X.|.$(MDLDI|
|00000640| 52 29 2f 6d 6f 64 65 6c | 2e 68 20 24 28 4e 4c 44 |R)/model|.h $(NLD|
|00000650| 49 52 29 2f 6e 6c 6d 64 | 6c 2e 68 20 24 28 4e 4c |IR)/nlmd|l.h $(NL|
|00000660| 44 49 52 29 2f 73 74 61 | 74 75 73 2e 68 20 20 5c |DIR)/sta|tus.h \|
|00000670| 0a 58 09 09 24 28 52 4d | 44 49 52 29 2f 75 73 75 |.X..$(RM|DIR)/usu|
|00000680| 61 6c 2e 68 20 24 28 52 | 4d 44 49 52 29 2f 72 65 |al.h $(R|MDIR)/re|
|00000690| 61 6c 6d 61 74 2e 68 20 | 24 28 52 4d 44 49 52 29 |almat.h |$(RMDIR)|
|000006a0| 2f 74 6f 6f 6c 73 2e 68 | 0a 58 09 24 28 43 43 29 |/tools.h|.X.$(CC)|
|000006b0| 20 24 28 43 46 4c 41 47 | 53 29 20 2d 63 20 2d 49 | $(CFLAG|S) -c -I|
|000006c0| 24 28 4d 44 4c 44 49 52 | 29 20 2d 49 24 28 4e 4c |$(MDLDIR|) -I$(NL|
|000006d0| 44 49 52 29 20 2d 49 24 | 28 52 4d 44 49 52 29 20 |DIR) -I$|(RMDIR) |
|000006e0| 24 28 4e 4c 44 49 52 29 | 2f 6e 6c 6d 64 6c 2e 63 |$(NLDIR)|/nlmdl.c|
|000006f0| 63 20 0a 58 0a 58 6d 6f | 64 65 6c 2e 6f 09 3a 20 |c .X.Xmo|del.o.: |
|00000700| 24 28 4d 44 4c 44 49 52 | 29 2f 6d 6f 64 65 6c 2e |$(MDLDIR|)/model.|
|00000710| 63 63 20 5c 0a 58 09 09 | 24 28 4d 44 4c 44 49 52 |cc \.X..|$(MDLDIR|
|00000720| 29 2f 6d 6f 64 65 6c 2e | 68 20 24 28 4e 4c 44 49 |)/model.|h $(NLDI|
|00000730| 52 29 2f 6e 6c 6d 64 6c | 2e 68 20 24 28 4e 4c 44 |R)/nlmdl|.h $(NLD|
|00000740| 49 52 29 2f 73 74 61 74 | 75 73 2e 68 20 20 5c 0a |IR)/stat|us.h \.|
|00000750| 58 09 09 24 28 52 4d 44 | 49 52 29 2f 75 73 75 61 |X..$(RMD|IR)/usua|
|00000760| 6c 2e 68 20 24 28 52 4d | 44 49 52 29 2f 72 65 61 |l.h $(RM|DIR)/rea|
|00000770| 6c 6d 61 74 2e 68 20 24 | 28 52 4d 44 49 52 29 2f |lmat.h $|(RMDIR)/|
|00000780| 74 6f 6f 6c 73 2e 68 0a | 58 09 24 28 43 43 29 20 |tools.h.|X.$(CC) |
|00000790| 24 28 43 46 4c 41 47 53 | 29 20 2d 63 20 2d 49 24 |$(CFLAGS|) -c -I$|
|000007a0| 28 4d 44 4c 44 49 52 29 | 20 2d 49 24 28 4e 4c 44 |(MDLDIR)| -I$(NLD|
|000007b0| 49 52 29 20 2d 49 24 28 | 52 4d 44 49 52 29 20 24 |IR) -I$(|RMDIR) $|
|000007c0| 28 4d 44 4c 44 49 52 29 | 2f 6d 6f 64 65 6c 2e 63 |(MDLDIR)|/model.c|
|000007d0| 63 20 0a 58 0a 58 73 74 | 61 74 75 73 2e 6f 09 3a |c .X.Xst|atus.o.:|
|000007e0| 20 24 28 4e 4c 44 49 52 | 29 2f 73 74 61 74 75 73 | $(NLDIR|)/status|
|000007f0| 2e 63 63 20 5c 0a 58 09 | 09 24 28 4e 4c 44 49 52 |.cc \.X.|.$(NLDIR|
|00000800| 29 2f 73 74 61 74 75 73 | 2e 68 20 5c 0a 58 09 09 |)/status|.h \.X..|
|00000810| 24 28 52 4d 44 49 52 29 | 2f 75 73 75 61 6c 2e 68 |$(RMDIR)|/usual.h|
|00000820| 20 24 28 52 4d 44 49 52 | 29 2f 72 65 61 6c 6d 61 | $(RMDIR|)/realma|
|00000830| 74 2e 68 20 24 28 52 4d | 44 49 52 29 2f 74 6f 6f |t.h $(RM|DIR)/too|
|00000840| 6c 73 2e 68 20 0a 58 09 | 24 28 43 43 29 20 24 28 |ls.h .X.|$(CC) $(|
|00000850| 43 46 4c 41 47 53 29 20 | 2d 63 20 2d 49 24 28 4d |CFLAGS) |-c -I$(M|
|00000860| 44 4c 44 49 52 29 20 2d | 49 24 28 4e 4c 44 49 52 |DLDIR) -|I$(NLDIR|
|00000870| 29 20 2d 49 24 28 52 4d | 44 49 52 29 20 24 28 4e |) -I$(RM|DIR) $(N|
|00000880| 4c 44 49 52 29 2f 73 74 | 61 74 75 73 2e 63 63 20 |LDIR)/st|atus.cc |
|00000890| 0a 58 0a 58 72 65 61 6c | 6d 61 74 2e 6f 09 3a 20 |.X.Xreal|mat.o.: |
|000008a0| 24 28 52 4d 44 49 52 29 | 2f 72 65 61 6c 6d 61 74 |$(RMDIR)|/realmat|
|000008b0| 2e 63 63 20 5c 0a 58 09 | 09 24 28 52 4d 44 49 52 |.cc \.X.|.$(RMDIR|
|000008c0| 29 2f 75 73 75 61 6c 2e | 68 20 24 28 52 4d 44 49 |)/usual.|h $(RMDI|
|000008d0| 52 29 2f 72 65 61 6c 6d | 61 74 2e 68 20 24 28 52 |R)/realm|at.h $(R|
|000008e0| 4d 44 49 52 29 2f 74 6f | 6f 6c 73 2e 68 0a 58 09 |MDIR)/to|ols.h.X.|
|000008f0| 24 28 43 43 29 20 24 28 | 43 46 4c 41 47 53 29 20 |$(CC) $(|CFLAGS) |
|00000900| 2d 63 20 2d 49 24 28 52 | 4d 44 49 52 29 20 24 28 |-c -I$(R|MDIR) $(|
|00000910| 52 4d 44 49 52 29 2f 72 | 65 61 6c 6d 61 74 2e 63 |RMDIR)/r|ealmat.c|
|00000920| 63 0a 58 0a 58 73 74 61 | 72 62 6f 78 2e 6f 09 3a |c.X.Xsta|rbox.o.:|
|00000930| 20 24 28 52 4d 44 49 52 | 29 2f 73 74 61 72 62 6f | $(RMDIR|)/starbo|
|00000940| 78 2e 63 63 20 5c 0a 58 | 09 09 24 28 52 4d 44 49 |x.cc \.X|..$(RMDI|
|00000950| 52 29 2f 75 73 75 61 6c | 2e 68 20 24 28 52 4d 44 |R)/usual|.h $(RMD|
|00000960| 49 52 29 2f 74 6f 6f 6c | 73 2e 68 20 0a 58 09 24 |IR)/tool|s.h .X.$|
|00000970| 28 43 43 29 20 24 28 43 | 46 4c 41 47 53 29 20 2d |(CC) $(C|FLAGS) -|
|00000980| 63 20 2d 49 24 28 52 4d | 44 49 52 29 20 24 28 52 |c -I$(RM|DIR) $(R|
|00000990| 4d 44 49 52 29 2f 73 74 | 61 72 62 6f 78 2e 63 63 |MDIR)/st|arbox.cc|
|000009a0| 0a 58 0a 58 64 63 6f 6e | 64 2e 6f 09 3a 20 24 28 |.X.Xdcon|d.o.: $(|
|000009b0| 52 4d 44 49 52 29 2f 64 | 63 6f 6e 64 2e 63 63 20 |RMDIR)/d|cond.cc |
|000009c0| 5c 0a 58 09 09 24 28 52 | 4d 44 49 52 29 2f 75 73 |\.X..$(R|MDIR)/us|
|000009d0| 75 61 6c 2e 68 20 24 28 | 52 4d 44 49 52 29 2f 74 |ual.h $(|RMDIR)/t|
|000009e0| 6f 6f 6c 73 2e 68 20 0a | 58 09 24 28 43 43 29 20 |ools.h .|X.$(CC) |
|000009f0| 24 28 43 46 4c 41 47 53 | 29 20 2d 63 20 2d 49 24 |$(CFLAGS|) -c -I$|
|00000a00| 28 52 4d 44 49 52 29 20 | 24 28 52 4d 44 49 52 29 |(RMDIR) |$(RMDIR)|
|00000a10| 2f 64 63 6f 6e 64 2e 63 | 63 0a 58 0a 58 64 73 77 |/dcond.c|c.X.Xdsw|
|00000a20| 65 65 70 2e 6f 09 3a 20 | 24 28 52 4d 44 49 52 29 |eep.o.: |$(RMDIR)|
|00000a30| 2f 64 73 77 65 65 70 2e | 63 63 20 5c 0a 58 09 09 |/dsweep.|cc \.X..|
|00000a40| 24 28 52 4d 44 49 52 29 | 2f 75 73 75 61 6c 2e 68 |$(RMDIR)|/usual.h|
|00000a50| 20 24 28 52 4d 44 49 52 | 29 2f 74 6f 6f 6c 73 2e | $(RMDIR|)/tools.|
|00000a60| 68 0a 58 09 24 28 43 43 | 29 20 24 28 43 46 4c 41 |h.X.$(CC|) $(CFLA|
|00000a70| 47 53 29 20 2d 63 20 2d | 49 24 28 52 4d 44 49 52 |GS) -c -|I$(RMDIR|
|00000a80| 29 20 24 28 52 4d 44 49 | 52 29 2f 64 73 77 65 65 |) $(RMDI|R)/dswee|
|00000a90| 70 2e 63 63 0a 58 0a 58 | 63 6c 65 61 6e 09 3a 0a |p.cc.X.X|clean.:.|
|00000aa0| 58 09 72 6d 20 2d 66 20 | 2a 2e 6f 20 63 6f 72 65 |X.rm -f |*.o core|
|00000ab0| 0a 58 0a 58 76 65 72 79 | 63 6c 65 61 6e 09 3a 0a |.X.Xvery|clean.:.|
|00000ac0| 58 09 72 6d 20 2d 66 20 | 2a 2e 6f 20 63 6f 72 65 |X.rm -f |*.o core|
|00000ad0| 20 6e 6c 6d 64 6c 20 64 | 65 74 61 69 6c 2e 64 61 | nlmdl d|etail.da|
|00000ae0| 74 20 65 6e 64 69 6e 67 | 2e 64 61 74 20 74 6d 70 |t ending|.dat tmp|
|00000af0| 2e 64 61 74 0a 45 4e 44 | 5f 4f 46 5f 46 49 4c 45 |.dat.END|_OF_FILE|
|00000b00| 0a 69 66 20 74 65 73 74 | 20 31 34 38 33 20 2d 6e |.if test| 1483 -n|
|00000b10| 65 20 60 77 63 20 2d 63 | 20 3c 27 4d 61 6b 65 66 |e `wc -c| <'Makef|
|00000b20| 69 6c 65 27 60 3b 20 74 | 68 65 6e 0a 20 20 20 20 |ile'`; t|hen. |
|00000b30| 65 63 68 6f 20 73 68 61 | 72 3a 20 5c 22 27 4d 61 |echo sha|r: \"'Ma|
|00000b40| 6b 65 66 69 6c 65 27 5c | 22 20 75 6e 70 61 63 6b |kefile'\|" unpack|
|00000b50| 65 64 20 77 69 74 68 20 | 77 72 6f 6e 67 20 73 69 |ed with |wrong si|
|00000b60| 7a 65 21 0a 66 69 0a 23 | 20 65 6e 64 20 6f 66 20 |ze!.fi.#| end of |
|00000b70| 27 4d 61 6b 65 66 69 6c | 65 27 0a 66 69 0a 69 66 |'Makefil|e'.fi.if|
|00000b80| 20 74 65 73 74 20 2d 66 | 20 27 52 45 41 44 4d 45 | test -f| 'README|
|00000b90| 27 20 2d 61 20 22 24 7b | 31 7d 22 20 21 3d 20 22 |' -a "${|1}" != "|
|00000ba0| 2d 63 22 20 3b 20 74 68 | 65 6e 20 0a 20 20 65 63 |-c" ; th|en . ec|
|00000bb0| 68 6f 20 73 68 61 72 3a | 20 57 69 6c 6c 20 6e 6f |ho shar:| Will no|
|00000bc0| 74 20 63 6c 6f 62 62 65 | 72 20 65 78 69 73 74 69 |t clobbe|r existi|
|00000bd0| 6e 67 20 66 69 6c 65 20 | 5c 22 27 52 45 41 44 4d |ng file |\"'READM|
|00000be0| 45 27 5c 22 0a 65 6c 73 | 65 0a 65 63 68 6f 20 73 |E'\".els|e.echo s|
|00000bf0| 68 61 72 3a 20 45 78 74 | 72 61 63 74 69 6e 67 20 |har: Ext|racting |
|00000c00| 5c 22 27 52 45 41 44 4d | 45 27 5c 22 20 5c 28 32 |\"'READM|E'\" \(2|
|00000c10| 38 37 37 20 63 68 61 72 | 61 63 74 65 72 73 5c 29 |877 char|acters\)|
|00000c20| 0a 73 65 64 20 22 73 2f | 5e 58 2f 2f 22 20 3e 27 |.sed "s/|^X//" >'|
|00000c30| 52 45 41 44 4d 45 27 20 | 3c 3c 27 45 4e 44 5f 4f |README' |<<'END_O|
|00000c40| 46 5f 46 49 4c 45 27 0a | 58 0a 58 6e 6c 6d 64 6c |F_FILE'.|X.Xnlmdl|
|00000c50| 20 69 73 20 61 20 43 2b | 2b 20 69 6d 70 6c 65 6d | is a C+|+ implem|
|00000c60| 65 6e 74 61 74 69 6f 6e | 20 6f 66 20 74 68 65 20 |entation| of the |
|00000c70| 73 74 61 74 69 73 74 69 | 63 61 6c 20 6d 65 74 68 |statisti|cal meth|
|00000c80| 6f 64 73 20 69 6e 20 41 | 2e 20 52 6f 6e 61 6c 64 |ods in A|. Ronald|
|00000c90| 20 47 61 6c 6c 61 6e 74 | 2c 20 0a 58 22 4e 6f 6e | Gallant|, .X"Non|
|00000ca0| 6c 69 6e 65 61 72 20 53 | 74 61 74 69 73 74 69 63 |linear S|tatistic|
|00000cb0| 61 6c 20 4d 6f 64 65 6c | 73 2c 22 20 4e 65 77 20 |al Model|s," New |
|00000cc0| 59 6f 72 6b 3a 20 4a 6f | 68 6e 20 57 69 6c 65 79 |York: Jo|hn Wiley|
|00000cd0| 20 61 6e 64 20 53 6f 6e | 73 2c 20 31 39 38 37 2c | and Son|s, 1987,|
|00000ce0| 20 0a 58 49 53 42 4e 20 | 30 2d 34 37 31 2d 38 30 | .XISBN |0-471-80|
|00000cf0| 32 36 30 2d 33 2c 20 75 | 73 69 6e 67 20 61 20 6d |260-3, u|sing a m|
|00000d00| 61 74 72 69 78 20 63 6c | 61 73 73 20 72 65 61 6c |atrix cl|ass real|
|00000d10| 6d 61 74 20 74 68 61 74 | 20 69 73 20 64 69 73 74 |mat that| is dist|
|00000d20| 72 69 62 75 74 65 64 20 | 77 69 74 68 20 69 74 2e |ributed |with it.|
|00000d30| 20 20 0a 58 54 68 65 20 | 6d 61 74 72 69 78 20 63 | .XThe |matrix c|
|00000d40| 6c 61 73 73 20 6d 61 79 | 20 62 65 20 6f 66 20 73 |lass may| be of s|
|00000d50| 6f 6d 65 20 69 6e 64 65 | 70 65 6e 64 65 6e 74 20 |ome inde|pendent |
|00000d60| 69 6e 74 65 72 65 73 74 | 2e 20 20 0a 58 0a 58 54 |interest|. .X.XT|
|00000d70| 68 65 20 70 72 6f 67 72 | 61 6d 20 63 6f 6d 70 75 |he progr|am compu|
|00000d80| 74 65 73 20 6c 65 61 73 | 74 20 73 71 75 61 72 65 |tes leas|t square|
|00000d90| 73 20 65 73 74 69 6d 61 | 74 65 73 20 66 6f 72 20 |s estima|tes for |
|00000da0| 61 20 75 6e 69 76 61 72 | 69 61 74 65 20 6e 6f 6e |a univar|iate non|
|00000db0| 6c 69 6e 65 61 72 20 0a | 58 72 65 67 72 65 73 73 |linear .|Xregress|
|00000dc0| 69 6f 6e 20 6d 6f 64 65 | 6c 20 61 6e 64 20 67 65 |ion mode|l and ge|
|00000dd0| 6e 65 72 61 6c 69 7a 65 | 64 20 6c 65 61 73 74 20 |neralize|d least |
|00000de0| 73 71 75 61 72 65 73 20 | 65 73 74 69 6d 61 74 65 |squares |estimate|
|00000df0| 73 20 66 6f 72 20 61 20 | 6d 75 6c 74 69 76 61 72 |s for a |multivar|
|00000e00| 69 61 74 65 20 0a 58 6e | 6f 6e 6c 69 6e 65 61 72 |iate .Xn|onlinear|
|00000e10| 20 72 65 67 72 65 73 73 | 69 6f 6e 20 6d 6f 64 65 | regress|ion mode|
|00000e20| 6c 2e 20 20 55 73 65 20 | 74 68 65 20 53 55 52 20 |l. Use |the SUR |
|00000e30| 6f 70 74 69 6f 6e 20 66 | 6f 72 20 62 6f 74 68 2e |option f|or both.|
|00000e40| 20 20 4f 74 68 65 72 20 | 6f 70 74 69 6f 6e 73 20 | Other |options |
|00000e50| 63 61 6e 20 62 65 20 0a | 58 73 65 74 20 74 6f 20 |can be .|Xset to |
|00000e60| 63 6f 72 72 65 63 74 20 | 76 61 72 69 61 6e 63 65 |correct |variance|
|00000e70| 20 65 73 74 69 6d 61 74 | 65 73 20 66 6f 72 20 68 | estimat|es for h|
|00000e80| 65 74 65 72 6f 73 6b 65 | 64 61 73 74 69 63 69 74 |eteroske|dasticit|
|00000e90| 79 20 61 6e 64 2f 6f 72 | 20 73 65 72 69 61 6c 20 |y and/or| serial |
|00000ea0| 0a 58 63 6f 72 72 65 6c | 61 74 69 6f 6e 2e 20 20 |.Xcorrel|ation. |
|00000eb0| 46 6f 72 20 73 69 6d 75 | 6c 74 61 6e 65 6f 75 73 |For simu|ltaneous|
|00000ec0| 20 6e 6f 6e 6c 69 6e 65 | 61 72 20 73 79 73 74 65 | nonline|ar syste|
|00000ed0| 6d 73 2c 20 70 6f 73 73 | 69 62 6c 79 20 69 6e 20 |ms, poss|ibly in |
|00000ee0| 69 6d 70 6c 69 63 69 74 | 20 66 6f 72 6d 2c 20 69 |implicit| form, i|
|00000ef0| 74 20 0a 58 63 6f 6d 70 | 75 74 65 73 20 74 68 72 |t .Xcomp|utes thr|
|00000f00| 65 65 2d 73 74 61 67 65 | 20 6c 65 61 73 74 2d 73 |ee-stage| least-s|
|00000f10| 71 75 61 72 65 73 20 65 | 73 74 69 6d 61 74 65 73 |quares e|stimates|
|00000f20| 2c 20 54 53 4c 53 20 6f | 70 74 69 6f 6e 2c 20 61 |, TSLS o|ption, a|
|00000f30| 6e 64 20 66 6f 72 20 6e | 6f 6e 6c 69 6e 65 61 72 |nd for n|onlinear|
|00000f40| 20 0a 58 64 79 6e 61 6d | 69 63 20 73 79 73 74 65 | .Xdynam|ic syste|
|00000f50| 6d 73 2c 20 70 6f 73 73 | 69 62 6c 79 20 69 6e 20 |ms, poss|ibly in |
|00000f60| 69 6d 70 6c 69 63 69 74 | 20 66 6f 72 6d 2c 20 67 |implicit| form, g|
|00000f70| 65 6e 65 72 61 6c 69 7a | 65 64 20 6d 65 74 68 6f |eneraliz|ed metho|
|00000f80| 64 2d 6f 66 2d 6d 6f 6d | 65 6e 74 73 20 0a 58 65 |d-of-mom|ents .Xe|
|00000f90| 73 74 69 6d 61 74 65 73 | 2c 20 47 4d 4d 20 6f 70 |stimates|, GMM op|
|00000fa0| 74 69 6f 6e 2e 20 20 0a | 58 0a 58 54 68 65 20 64 |tion. .|X.XThe d|
|00000fb0| 6f 63 75 6d 65 6e 74 61 | 74 69 6f 6e 20 63 6f 6e |ocumenta|tion con|
|00000fc0| 73 69 73 74 73 20 6f 66 | 20 74 68 69 73 20 66 69 |sists of| this fi|
|00000fd0| 6c 65 2c 20 74 68 65 20 | 74 77 6f 20 68 65 61 64 |le, the |two head|
|00000fe0| 65 72 20 66 69 6c 65 73 | 20 6e 6c 6d 64 6c 2e 68 |er files| nlmdl.h|
|00000ff0| 20 61 6e 64 20 0a 58 72 | 65 61 6c 6d 61 74 2e 68 | and .Xr|ealmat.h|
|00001000| 2c 20 61 6e 64 20 74 68 | 65 20 65 78 61 6d 70 6c |, and th|e exampl|
|00001010| 65 73 20 69 6e 20 74 68 | 65 20 64 69 72 65 63 74 |es in th|e direct|
|00001020| 6f 72 69 65 73 20 63 68 | 31 65 67 31 2c 20 63 68 |ories ch|1eg1, ch|
|00001030| 35 65 67 31 2c 20 63 68 | 36 65 67 31 2c 20 61 6e |5eg1, ch|6eg1, an|
|00001040| 64 20 0a 58 63 68 36 65 | 67 32 2e 20 20 0a 58 0a |d .Xch6e|g2. .X.|
|00001050| 58 41 73 20 64 69 73 74 | 72 69 62 75 74 65 64 2c |XAs dist|ributed,|
|00001060| 20 74 68 65 20 63 6f 64 | 65 20 69 73 20 73 65 74 | the cod|e is set|
|00001070| 20 75 70 20 66 6f 72 20 | 75 73 65 20 77 69 74 68 | up for |use with|
|00001080| 20 47 4e 55 27 73 20 67 | 2b 2b 20 63 6f 6d 70 69 | GNU's g|++ compi|
|00001090| 6c 65 72 20 6f 6e 20 61 | 20 55 6e 69 78 20 0a 58 |ler on a| Unix .X|
|000010a0| 62 6f 78 2e 20 20 54 6f | 20 63 68 65 63 6b 20 75 |box. To| check u|
|000010b0| 73 69 6e 67 20 6f 6e 65 | 20 6f 66 20 74 68 65 20 |sing one| of the |
|000010c0| 65 78 61 6d 70 6c 65 73 | 2c 20 65 2e 67 2e 20 63 |examples|, e.g. c|
|000010d0| 68 31 65 67 31 2c 20 63 | 64 20 74 6f 20 63 68 31 |h1eg1, c|d to ch1|
|000010e0| 65 67 31 20 61 6e 64 20 | 74 79 70 65 20 0a 58 6d |eg1 and |type .Xm|
|000010f0| 61 6b 65 20 2d 66 20 2e | 2e 2f 4d 61 6b 65 66 69 |ake -f .|./Makefi|
|00001100| 6c 65 2e 20 20 45 61 63 | 68 20 64 69 72 65 63 74 |le. Eac|h direct|
|00001110| 6f 72 79 20 68 61 73 20 | 61 20 66 69 6c 65 20 64 |ory has |a file d|
|00001120| 65 74 61 69 6c 2e 62 61 | 6b 20 63 72 65 61 74 65 |etail.ba|k create|
|00001130| 64 20 75 73 69 6e 67 20 | 0a 58 6e 6c 6d 64 6c 20 |d using |.Xnlmdl |
|00001140| 3e 20 64 65 74 61 69 6c | 2e 64 61 74 20 74 68 61 |> detail|.dat tha|
|00001150| 74 20 79 6f 75 20 63 61 | 6e 20 75 73 65 20 66 6f |t you ca|n use fo|
|00001160| 72 20 63 6f 6d 70 61 72 | 69 73 6f 6e 2e 20 20 0a |r compar|ison. .|
|00001170| 58 0a 58 41 6c 74 68 6f | 75 67 68 20 69 6e 74 65 |X.XAltho|ugh inte|
|00001180| 6e 64 65 64 20 66 6f 72 | 20 75 73 65 20 77 69 74 |nded for| use wit|
|00001190| 68 20 47 4e 55 27 73 20 | 67 2b 2b 20 63 6f 6d 70 |h GNU's |g++ comp|
|000011a0| 69 6c 65 72 2c 20 74 68 | 65 20 63 6f 64 65 20 63 |iler, th|e code c|
|000011b0| 61 6e 20 62 65 20 75 73 | 65 64 20 77 69 74 68 20 |an be us|ed with |
|000011c0| 0a 58 42 6f 72 6c 61 6e | 64 27 73 20 54 75 72 62 |.XBorlan|d's Turb|
|000011d0| 6f 20 43 2b 2b 20 63 6f | 6d 70 69 6c 65 72 2e 20 |o C++ co|mpiler. |
|000011e0| 20 54 6f 20 75 73 65 20 | 54 75 72 62 6f 20 43 2b | To use |Turbo C+|
|000011f0| 2b 2c 20 70 65 72 6d 75 | 74 65 20 74 68 65 20 23 |+, permu|te the #|
|00001200| 75 6e 64 65 66 20 61 6e | 64 20 23 64 65 66 69 6e |undef an|d #defin|
|00001210| 65 20 0a 58 66 6f 72 20 | 63 6f 6d 70 69 6c 65 72 |e .Xfor |compiler|
|00001220| 20 63 68 6f 69 63 65 20 | 69 6e 20 75 73 75 61 6c | choice |in usual|
|00001230| 2e 68 2c 20 63 6f 70 79 | 20 74 68 65 20 63 6f 6e |.h, copy| the con|
|00001240| 74 65 6e 74 73 20 6f 66 | 2c 20 65 2e 67 2e 2c 20 |tents of|, e.g., |
|00001250| 63 68 31 65 67 31 20 74 | 6f 20 74 68 65 20 0a 58 |ch1eg1 t|o the .X|
|00001260| 77 6f 72 6b 69 6e 67 20 | 64 69 72 65 63 74 6f 72 |working |director|
|00001270| 79 2c 20 61 6e 64 20 6c | 69 73 74 20 65 76 65 72 |y, and l|ist ever|
|00001280| 79 20 66 69 6c 65 20 77 | 69 74 68 20 61 20 2e 63 |y file w|ith a .c|
|00001290| 63 20 65 78 74 65 6e 73 | 69 6f 6e 20 69 6e 20 74 |c extens|ion in t|
|000012a0| 68 65 20 2e 50 52 4a 20 | 66 69 6c 65 20 0a 58 65 |he .PRJ |file .Xe|
|000012b0| 78 63 65 70 74 20 6e 6c | 6f 70 72 2e 63 63 2e 20 |xcept nl|opr.cc. |
|000012c0| 20 53 6d 61 6c 6c 20 6d | 6f 64 65 6c 73 2c 20 73 | Small m|odels, s|
|000012d0| 75 63 68 20 61 73 20 63 | 68 31 65 67 31 2c 20 65 |uch as c|h1eg1, e|
|000012e0| 78 65 63 75 74 65 20 63 | 6f 72 72 65 63 74 6c 79 |xecute c|orrectly|
|000012f0| 20 77 68 65 6e 20 63 6f | 6d 70 69 6c 65 64 20 0a | when co|mpiled .|
|00001300| 58 77 69 74 68 20 74 68 | 65 20 6c 61 72 67 65 20 |Xwith th|e large |
|00001310| 6d 65 6d 6f 72 79 20 6f | 70 74 69 6f 6e 20 62 75 |memory o|ption bu|
|00001320| 74 20 6c 61 72 67 65 72 | 20 6f 6e 65 73 2c 20 73 |t larger| ones, s|
|00001330| 75 63 68 20 61 73 20 63 | 68 35 65 67 32 2c 20 72 |uch as c|h5eg2, r|
|00001340| 75 6e 20 6f 75 74 20 6f | 66 20 6d 65 6d 6f 72 79 |un out o|f memory|
|00001350| 20 0a 58 64 75 72 69 6e | 67 20 65 78 65 63 75 74 | .Xdurin|g execut|
|00001360| 69 6f 6e 2e 20 20 49 20 | 68 61 76 65 6e 27 74 20 |ion. I |haven't |
|00001370| 74 72 69 65 64 20 54 75 | 72 62 6f 27 73 20 6f 76 |tried Tu|rbo's ov|
|00001380| 65 72 6c 61 79 20 61 6e | 64 20 6d 65 6d 6f 72 79 |erlay an|d memory|
|00001390| 20 65 78 70 61 6e 73 69 | 6f 6e 20 0a 58 66 65 61 | expansi|on .Xfea|
|000013a0| 74 75 72 65 73 20 74 6f | 20 73 65 65 20 69 66 20 |tures to| see if |
|000013b0| 6c 61 72 67 65 72 20 6d | 6f 64 65 6c 73 20 63 61 |larger m|odels ca|
|000013c0| 6e 20 62 65 20 61 63 63 | 6f 6d 6d 6f 64 61 74 65 |n be acc|ommodate|
|000013d0| 64 2e 20 20 54 68 65 20 | 63 6f 64 65 20 68 61 73 |d. The |code has|
|000013e0| 6e 27 74 20 62 65 65 6e | 20 0a 58 74 65 73 74 65 |n't been| .Xteste|
|000013f0| 64 20 77 69 74 68 20 6f | 74 68 65 72 20 43 2b 2b |d with o|ther C++|
|00001400| 20 63 6f 6d 70 69 6c 65 | 72 73 2e 20 20 0a 58 0a | compile|rs. .X.|
|00001410| 58 43 6f 70 79 72 69 67 | 68 74 20 28 43 29 20 31 |XCopyrig|ht (C) 1|
|00001420| 39 39 30 20 62 79 20 41 | 2e 20 52 6f 6e 61 6c 64 |990 by A|. Ronald|
|00001430| 20 47 61 6c 6c 61 6e 74 | 2e 0a 58 0a 58 50 65 72 | Gallant|..X.XPer|
|00001440| 6d 69 73 73 69 6f 6e 20 | 74 6f 20 75 73 65 2c 20 |mission |to use, |
|00001450| 63 6f 70 79 2c 20 6d 6f | 64 69 66 79 2c 20 61 6e |copy, mo|dify, an|
|00001460| 64 20 64 69 73 74 72 69 | 62 75 74 65 20 74 68 69 |d distri|bute thi|
|00001470| 73 20 73 6f 66 74 77 61 | 72 65 20 61 6e 64 20 69 |s softwa|re and i|
|00001480| 74 73 20 0a 58 64 6f 63 | 75 6d 65 6e 74 61 74 69 |ts .Xdoc|umentati|
|00001490| 6f 6e 20 66 6f 72 20 61 | 6e 79 20 70 75 72 70 6f |on for a|ny purpo|
|000014a0| 73 65 20 61 6e 64 20 77 | 69 74 68 6f 75 74 20 66 |se and w|ithout f|
|000014b0| 65 65 20 69 73 20 68 65 | 72 65 62 79 20 67 72 61 |ee is he|reby gra|
|000014c0| 6e 74 65 64 2c 20 70 72 | 6f 76 69 64 65 64 20 74 |nted, pr|ovided t|
|000014d0| 68 61 74 20 0a 58 74 68 | 65 20 61 62 6f 76 65 20 |hat .Xth|e above |
|000014e0| 63 6f 70 79 72 69 67 68 | 74 20 6e 6f 74 69 63 65 |copyrigh|t notice|
|000014f0| 20 61 70 70 65 61 72 20 | 69 6e 20 61 6c 6c 20 63 | appear |in all c|
|00001500| 6f 70 69 65 73 20 61 6e | 64 20 74 68 61 74 20 62 |opies an|d that b|
|00001510| 6f 74 68 20 74 68 61 74 | 20 63 6f 70 79 72 69 67 |oth that| copyrig|
|00001520| 68 74 20 0a 58 6e 6f 74 | 69 63 65 20 61 6e 64 20 |ht .Xnot|ice and |
|00001530| 74 68 69 73 20 70 65 72 | 6d 69 73 73 69 6f 6e 20 |this per|mission |
|00001540| 6e 6f 74 69 63 65 20 61 | 70 70 65 61 72 20 69 6e |notice a|ppear in|
|00001550| 20 73 75 70 70 6f 72 74 | 69 6e 67 20 64 6f 63 75 | support|ing docu|
|00001560| 6d 65 6e 74 61 74 69 6f | 6e 2e 20 20 0a 58 0a 58 |mentatio|n. .X.X|
|00001570| 54 68 69 73 20 73 6f 66 | 74 77 61 72 65 20 69 73 |This sof|tware is|
|00001580| 20 70 72 6f 76 69 64 65 | 64 20 22 61 73 20 69 73 | provide|d "as is|
|00001590| 22 20 77 69 74 68 6f 75 | 74 20 61 6e 79 20 65 78 |" withou|t any ex|
|000015a0| 70 72 65 73 73 65 64 20 | 6f 72 20 69 6d 70 6c 69 |pressed |or impli|
|000015b0| 65 64 20 77 61 72 72 61 | 6e 74 79 2e 0a 58 0a 58 |ed warra|nty..X.X|
|000015c0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000015d0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000015e0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000015f0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00001600| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 0a |--------|-------.|
|00001610| 58 50 72 6f 66 65 73 73 | 6f 72 20 41 2e 20 52 6f |XProfess|or A. Ro|
|00001620| 6e 61 6c 64 20 47 61 6c | 6c 61 6e 74 20 20 20 20 |nald Gal|lant |
|00001630| 20 20 20 20 50 68 6f 6e | 65 3a 20 20 39 31 39 2d | Phon|e: 919-|
|00001640| 37 33 37 2d 32 35 33 31 | 0a 58 44 65 70 61 72 74 |737-2531|.XDepart|
|00001650| 6d 65 6e 74 20 6f 66 20 | 53 74 61 74 69 73 74 69 |ment of |Statisti|
|00001660| 63 73 20 20 20 20 20 20 | 20 20 20 20 20 46 41 58 |cs | FAX|
|00001670| 3a 20 20 20 20 39 31 39 | 2d 37 33 37 2d 37 35 39 |: 919|-737-759|
|00001680| 31 0a 58 4e 6f 72 74 68 | 20 43 61 72 6f 6c 69 6e |1.XNorth| Carolin|
|00001690| 61 20 53 74 61 74 65 20 | 55 6e 69 76 65 72 73 69 |a State |Universi|
|000016a0| 74 79 20 20 20 20 65 2d | 6d 61 69 6c 3a 20 61 72 |ty e-|mail: ar|
|000016b0| 67 40 63 63 76 72 31 2e | 63 63 2e 6e 63 73 75 2e |g@ccvr1.|cc.ncsu.|
|000016c0| 65 64 75 0a 58 43 61 6d | 70 75 73 20 42 6f 78 20 |edu.XCam|pus Box |
|000016d0| 38 32 30 33 20 20 20 20 | 20 20 20 20 20 20 20 20 |8203 | |
|000016e0| 20 20 20 20 20 20 20 20 | 66 74 70 3a 20 20 20 20 | |ftp: |
|000016f0| 63 63 76 72 31 2e 63 63 | 2e 6e 63 73 75 2e 65 64 |ccvr1.cc|.ncsu.ed|
|00001700| 75 20 28 31 32 38 2e 31 | 30 39 2e 32 31 32 2e 32 |u (128.1|09.212.2|
|00001710| 30 29 2c 20 0a 58 52 61 | 6c 65 69 67 68 20 4e 43 |0), .XRa|leigh NC|
|00001720| 20 32 37 36 39 35 2d 38 | 32 30 33 20 55 53 41 20 | 27695-8|203 USA |
|00001730| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00001740| 20 75 73 65 72 3a 20 61 | 6e 6f 6e 79 6d 6f 75 73 | user: a|nonymous|
|00001750| 2c 20 64 69 72 65 63 74 | 6f 72 79 3a 20 70 75 62 |, direct|ory: pub|
|00001760| 2f 61 72 67 20 0a 58 2d | 2d 2d 2d 2d 2d 2d 2d 2d |/arg .X-|--------|
|00001770| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00001780| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|00001790| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000017a0| 2d 2d 2d 2d 2d 2d 2d 2d | 2d 2d 2d 2d 2d 2d 2d 2d |--------|--------|
|000017b0| 2d 2d 2d 2d 2d 2d 0a 45 | 4e 44 5f 4f 46 5f 46 49 |------.E|ND_OF_FI|
|000017c0| 4c 45 0a 69 66 20 74 65 | 73 74 20 32 38 37 37 20 |LE.if te|st 2877 |
|000017d0| 2d 6e 65 20 60 77 63 20 | 2d 63 20 3c 27 52 45 41 |-ne `wc |-c <'REA|
|000017e0| 44 4d 45 27 60 3b 20 74 | 68 65 6e 0a 20 20 20 20 |DME'`; t|hen. |
|000017f0| 65 63 68 6f 20 73 68 61 | 72 3a 20 5c 22 27 52 45 |echo sha|r: \"'RE|
|00001800| 41 44 4d 45 27 5c 22 20 | 75 6e 70 61 63 6b 65 64 |ADME'\" |unpacked|
|00001810| 20 77 69 74 68 20 77 72 | 6f 6e 67 20 73 69 7a 65 | with wr|ong size|
|00001820| 21 0a 66 69 0a 23 20 65 | 6e 64 20 6f 66 20 27 52 |!.fi.# e|nd of 'R|
|00001830| 45 41 44 4d 45 27 0a 66 | 69 0a 69 66 20 74 65 73 |EADME'.f|i.if tes|
|00001840| 74 20 21 20 2d 64 20 27 | 63 68 31 65 67 31 27 20 |t ! -d '|ch1eg1' |
|00001850| 3b 20 74 68 65 6e 0a 20 | 20 20 20 65 63 68 6f 20 |; then. | echo |
|00001860| 73 68 61 72 3a 20 43 72 | 65 61 74 69 6e 67 20 64 |shar: Cr|eating d|
|00001870| 69 72 65 63 74 6f 72 79 | 20 5c 22 27 63 68 31 65 |irectory| \"'ch1e|
|00001880| 67 31 27 5c 22 0a 20 20 | 20 20 6d 6b 64 69 72 20 |g1'\". | mkdir |
|00001890| 27 63 68 31 65 67 31 27 | 0a 66 69 0a 69 66 20 74 |'ch1eg1'|.fi.if t|
|000018a0| 65 73 74 20 2d 66 20 27 | 63 68 31 65 67 31 2f 61 |est -f '|ch1eg1/a|
|000018b0| 6d 73 74 61 74 2e 64 61 | 74 27 20 2d 61 20 22 24 |mstat.da|t' -a "$|
|000018c0| 7b 31 7d 22 20 21 3d 20 | 22 2d 63 22 20 3b 20 74 |{1}" != |"-c" ; t|
|000018d0| 68 65 6e 20 0a 20 20 65 | 63 68 6f 20 73 68 61 72 |hen . e|cho shar|
|000018e0| 3a 20 57 69 6c 6c 20 6e | 6f 74 20 63 6c 6f 62 62 |: Will n|ot clobb|
|000018f0| 65 72 20 65 78 69 73 74 | 69 6e 67 20 66 69 6c 65 |er exist|ing file|
|00001900| 20 5c 22 27 63 68 31 65 | 67 31 2f 61 6d 73 74 61 | \"'ch1e|g1/amsta|
|00001910| 74 2e 64 61 74 27 5c 22 | 0a 65 6c 73 65 0a 65 63 |t.dat'\"|.else.ec|
|00001920| 68 6f 20 73 68 61 72 3a | 20 45 78 74 72 61 63 74 |ho shar:| Extract|
|00001930| 69 6e 67 20 5c 22 27 63 | 68 31 65 67 31 2f 61 6d |ing \"'c|h1eg1/am|
|00001940| 73 74 61 74 2e 64 61 74 | 27 5c 22 20 5c 28 31 30 |stat.dat|'\" \(10|
|00001950| 32 30 20 63 68 61 72 61 | 63 74 65 72 73 5c 29 0a |20 chara|cters\).|
|00001960| 73 65 64 20 22 73 2f 5e | 58 2f 2f 22 20 3e 27 63 |sed "s/^|X//" >'c|
|00001970| 68 31 65 67 31 2f 61 6d | 73 74 61 74 2e 64 61 74 |h1eg1/am|stat.dat|
|00001980| 27 20 3c 3c 27 45 4e 44 | 5f 4f 46 5f 46 49 4c 45 |' <<'END|_OF_FILE|
|00001990| 27 0a 58 20 31 20 20 20 | 20 30 2e 39 38 36 31 30 |'.X 1 | 0.98610|
|000019a0| 20 20 20 20 31 20 20 20 | 20 20 31 20 20 20 20 20 | 1 | 1 |
|000019b0| 36 2e 32 38 0a 58 20 32 | 20 20 20 20 31 2e 30 33 |6.28.X 2| 1.03|
|000019c0| 38 34 38 20 20 20 20 30 | 20 20 20 20 20 31 20 20 |848 0| 1 |
|000019d0| 20 20 20 39 2e 38 36 0a | 58 20 33 20 20 20 20 30 | 9.86.|X 3 0|
|000019e0| 2e 39 35 34 38 32 20 20 | 20 20 31 20 20 20 20 20 |.95482 | 1 |
|000019f0| 31 20 20 20 20 20 39 2e | 31 31 0a 58 20 34 20 20 |1 9.|11.X 4 |
|00001a00| 20 20 31 2e 30 34 31 38 | 34 20 20 20 20 30 20 20 | 1.0418|4 0 |
|00001a10| 20 20 20 31 20 20 20 20 | 20 38 2e 34 33 0a 58 20 | 1 | 8.43.X |
|00001a20| 35 20 20 20 20 31 2e 30 | 32 33 32 34 20 20 20 20 |5 1.0|2324 |
|00001a30| 31 20 20 20 20 20 31 20 | 20 20 20 20 38 2e 31 31 |1 1 | 8.11|
|00001a40| 0a 58 20 36 20 20 20 20 | 30 2e 39 30 34 37 35 20 |.X 6 |0.90475 |
|00001a50| 20 20 20 30 20 20 20 20 | 20 31 20 20 20 20 20 31 | 0 | 1 1|
|00001a60| 2e 38 32 0a 58 20 37 20 | 20 20 20 30 2e 39 36 32 |.82.X 7 | 0.962|
|00001a70| 36 33 20 20 20 20 31 20 | 20 20 20 20 31 20 20 20 |63 1 | 1 |
|00001a80| 20 20 36 2e 35 38 0a 58 | 20 38 20 20 20 20 31 2e | 6.58.X| 8 1.|
|00001a90| 30 35 30 32 36 20 20 20 | 20 30 20 20 20 20 20 31 |05026 | 0 1|
|00001aa0| 20 20 20 20 20 35 2e 30 | 32 0a 58 20 39 20 20 20 | 5.0|2.X 9 |
|00001ab0| 20 30 2e 39 38 38 36 31 | 20 20 20 20 31 20 20 20 | 0.98861| 1 |
|00001ac0| 20 20 31 20 20 20 20 20 | 36 2e 35 32 0a 58 31 30 | 1 |6.52.X10|
|00001ad0| 20 20 20 20 31 2e 30 33 | 34 33 37 20 20 20 20 30 | 1.03|437 0|
|00001ae0| 20 20 20 20 20 31 20 20 | 20 20 20 33 2e 37 35 0a | 1 | 3.75.|
|00001af0| 58 31 31 20 20 20 20 30 | 2e 39 38 39 38 32 20 20 |X11 0|.98982 |
|00001b00| 20 20 31 20 20 20 20 20 | 31 20 20 20 20 20 39 2e | 1 |1 9.|
|00001b10| 38 36 0a 58 31 32 20 20 | 20 20 31 2e 30 31 32 31 |86.X12 | 1.0121|
|00001b20| 34 20 20 20 20 30 20 20 | 20 20 20 31 20 20 20 20 |4 0 | 1 |
|00001b30| 20 37 2e 33 31 0a 58 31 | 33 20 20 20 20 30 2e 36 | 7.31.X1|3 0.6|
|00001b40| 36 37 36 38 20 20 20 20 | 31 20 20 20 20 20 31 20 |6768 |1 1 |
|00001b50| 20 20 20 20 30 2e 34 37 | 0a 58 31 34 20 20 20 20 | 0.47|.X14 |
|00001b60| 30 2e 35 35 31 30 37 20 | 20 20 20 30 20 20 20 20 |0.55107 | 0 |
|00001b70| 20 31 20 20 20 20 20 30 | 2e 30 37 0a 58 31 35 20 | 1 0|.07.X15 |
|00001b80| 20 20 20 30 2e 39 36 38 | 32 32 20 20 20 20 31 20 | 0.968|22 1 |
|00001b90| 20 20 20 20 31 20 20 20 | 20 20 34 2e 30 37 0a 58 | 1 | 4.07.X|
|00001ba0| 31 36 20 20 20 20 30 2e | 39 38 38 32 33 20 20 20 |16 0.|98823 |
|00001bb0| 20 30 20 20 20 20 20 31 | 20 20 20 20 20 34 2e 36 | 0 1| 4.6|
|00001bc0| 31 0a 58 31 37 20 20 20 | 20 30 2e 35 39 37 35 39 |1.X17 | 0.59759|
|00001bd0| 20 20 20 20 31 20 20 20 | 20 20 31 20 20 20 20 20 | 1 | 1 |
|00001be0| 30 2e 31 37 0a 58 31 38 | 20 20 20 20 30 2e 39 39 |0.17.X18| 0.99|
|00001bf0| 34 31 38 20 20 20 20 30 | 20 20 20 20 20 31 20 20 |418 0| 1 |
|00001c00| 20 20 20 36 2e 39 39 0a | 58 31 39 20 20 20 20 31 | 6.99.|X19 1|
|00001c10| 2e 30 31 39 36 32 20 20 | 20 20 31 20 20 20 20 20 |.01962 | 1 |
|00001c20| 31 20 20 20 20 20 34 2e | 33 39 0a 58 32 30 20 20 |1 4.|39.X20 |
|00001c30| 20 20 30 2e 36 39 31 36 | 33 20 20 20 20 30 20 20 | 0.6916|3 0 |
|00001c40| 20 20 20 31 20 20 20 20 | 20 30 2e 33 39 0a 58 32 | 1 | 0.39.X2|
|00001c50| 31 20 20 20 20 31 2e 30 | 34 32 35 35 20 20 20 20 |1 1.0|4255 |
|00001c60| 31 20 20 20 20 20 31 20 | 20 20 20 20 34 2e 37 33 |1 1 | 4.73|
|00001c70| 0a 58 32 32 20 20 20 20 | 31 2e 30 34 33 34 33 20 |.X22 |1.04343 |
|00001c80| 20 20 20 30 20 20 20 20 | 20 31 20 20 20 20 20 39 | 0 | 1 9|
|00001c90| 2e 34 32 0a 58 32 33 20 | 20 20 20 30 2e 39 37 35 |.42.X23 | 0.975|
|00001ca0| 32 36 20 20 20 20 31 20 | 20 20 20 20 31 20 20 20 |26 1 | 1 |
|00001cb0| 20 20 38 2e 39 30 0a 58 | 32 34 20 20 20 20 31 2e | 8.90.X|24 1.|
|00001cc0| 30 34 39 36 39 20 20 20 | 20 30 20 20 20 20 20 31 |04969 | 0 1|
|00001cd0| 20 20 20 20 20 33 2e 30 | 32 0a 58 32 35 20 20 20 | 3.0|2.X25 |
|00001ce0| 20 30 2e 38 30 32 31 39 | 20 20 20 20 31 20 20 20 | 0.80219| 1 |
|00001cf0| 20 20 31 20 20 20 20 20 | 30 2e 37 37 0a 58 32 36 | 1 |0.77.X26|
|00001d00| 20 20 20 20 31 2e 30 31 | 30 34 36 20 20 20 20 30 | 1.01|046 0|
|00001d10| 20 20 20 20 20 31 20 20 | 20 20 20 33 2e 33 31 0a | 1 | 3.31.|
|00001d20| 58 32 37 20 20 20 20 30 | 2e 39 35 31 39 36 20 20 |X27 0|.95196 |
|00001d30| 20 20 31 20 20 20 20 20 | 31 20 20 20 20 20 34 2e | 1 |1 4.|
|00001d40| 35 31 0a 58 32 38 20 20 | 20 20 30 2e 39 37 36 35 |51.X28 | 0.9765|
|00001d50| 38 20 20 20 20 30 20 20 | 20 20 20 31 20 20 20 20 |8 0 | 1 |
|00001d60| 20 32 2e 36 35 0a 58 32 | 39 20 20 20 20 30 2e 35 | 2.65.X2|9 0.5|
|00001d70| 30 38 31 31 20 20 20 20 | 31 20 20 20 20 20 31 20 |0811 |1 1 |
|00001d80| 20 20 20 20 30 2e 30 38 | 0a 58 33 30 20 20 20 20 | 0.08|.X30 |
|00001d90| 30 2e 39 31 38 34 30 20 | 20 20 20 30 20 20 20 20 |0.91840 | 0 |
|00001da0| 20 31 20 20 20 20 20 36 | 2e 31 31 0a 45 4e 44 5f | 1 6|.11.END_|
|00001db0| 4f 46 5f 46 49 4c 45 0a | 69 66 20 74 65 73 74 20 |OF_FILE.|if test |
|00001dc0| 31 30 32 30 20 2d 6e 65 | 20 60 77 63 20 2d 63 20 |1020 -ne| `wc -c |
|00001dd0| 3c 27 63 68 31 65 67 31 | 2f 61 6d 73 74 61 74 2e |<'ch1eg1|/amstat.|
|00001de0| 64 61 74 27 60 3b 20 74 | 68 65 6e 0a 20 20 20 20 |dat'`; t|hen. |
|00001df0| 65 63 68 6f 20 73 68 61 | 72 3a 20 5c 22 27 63 68 |echo sha|r: \"'ch|
|00001e00| 31 65 67 31 2f 61 6d 73 | 74 61 74 2e 64 61 74 27 |1eg1/ams|tat.dat'|
|00001e10| 5c 22 20 75 6e 70 61 63 | 6b 65 64 20 77 69 74 68 |\" unpac|ked with|
|00001e20| 20 77 72 6f 6e 67 20 73 | 69 7a 65 21 0a 66 69 0a | wrong s|ize!.fi.|
|00001e30| 23 20 65 6e 64 20 6f 66 | 20 27 63 68 31 65 67 31 |# end of| 'ch1eg1|
|00001e40| 2f 61 6d 73 74 61 74 2e | 64 61 74 27 0a 66 69 0a |/amstat.|dat'.fi.|
|00001e50| 69 66 20 74 65 73 74 20 | 2d 66 20 27 63 68 31 65 |if test |-f 'ch1e|
|00001e60| 67 31 2f 65 6e 64 69 6e | 67 2e 62 61 6b 27 20 2d |g1/endin|g.bak' -|
|00001e70| 61 20 22 24 7b 31 7d 22 | 20 21 3d 20 22 2d 63 22 |a "${1}"| != "-c"|
|00001e80| 20 3b 20 74 68 65 6e 20 | 0a 20 20 65 63 68 6f 20 | ; then |. echo |
|00001e90| 73 68 61 72 3a 20 57 69 | 6c 6c 20 6e 6f 74 20 63 |shar: Wi|ll not c|
|00001ea0| 6c 6f 62 62 65 72 20 65 | 78 69 73 74 69 6e 67 20 |lobber e|xisting |
|00001eb0| 66 69 6c 65 20 5c 22 27 | 63 68 31 65 67 31 2f 65 |file \"'|ch1eg1/e|
|00001ec0| 6e 64 69 6e 67 2e 62 61 | 6b 27 5c 22 0a 65 6c 73 |nding.ba|k'\".els|
|00001ed0| 65 0a 65 63 68 6f 20 73 | 68 61 72 3a 20 45 78 74 |e.echo s|har: Ext|
|00001ee0| 72 61 63 74 69 6e 67 20 | 5c 22 27 63 68 31 65 67 |racting |\"'ch1eg|
|00001ef0| 31 2f 65 6e 64 69 6e 67 | 2e 62 61 6b 27 5c 22 20 |1/ending|.bak'\" |
|00001f00| 5c 28 32 35 32 33 20 63 | 68 61 72 61 63 74 65 72 |\(2523 c|haracter|
|00001f10| 73 5c 29 0a 73 65 64 20 | 22 73 2f 5e 58 2f 2f 22 |s\).sed |"s/^X//"|
|00001f20| 20 3e 27 63 68 31 65 67 | 31 2f 65 6e 64 69 6e 67 | >'ch1eg|1/ending|
|00001f30| 2e 62 61 6b 27 20 3c 3c | 27 45 4e 44 5f 4f 46 5f |.bak' <<|'END_OF_|
|00001f40| 46 49 4c 45 27 0a 58 47 | 61 6c 6c 61 6e 74 2c 20 |FILE'.XG|allant, |
|00001f50| 22 4e 6f 6e 6c 69 6e 65 | 61 72 20 53 74 61 74 69 |"Nonline|ar Stati|
|00001f60| 73 74 69 63 61 6c 20 4d | 6f 64 65 6c 73 2c 22 20 |stical M|odels," |
|00001f70| 43 68 61 70 74 65 72 20 | 31 2c 20 46 69 67 75 72 |Chapter |1, Figur|
|00001f80| 65 20 35 2c 20 70 2e 20 | 33 35 2d 33 36 2e 0a 58 |e 5, p. |35-36..X|
|00001f90| 53 55 52 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |SUR | |
|00001fa0| 20 20 20 57 68 61 74 20 | 65 73 74 69 6d 61 74 69 | What |estimati|
|00001fb0| 6f 6e 20 6d 65 74 68 6f | 64 3f 20 20 43 6f 64 65 |on metho|d? Code|
|00001fc0| 20 53 55 52 2c 20 54 53 | 4c 53 2c 20 6f 72 20 47 | SUR, TS|LS, or G|
|00001fd0| 4d 4d 2e 0a 58 33 30 20 | 20 20 20 20 20 20 20 20 |MM..X30 | |
|00001fe0| 20 20 20 20 20 20 20 20 | 4e 75 6d 62 65 72 20 6f | |Number o|
|00001ff0| 66 20 6f 62 73 65 72 76 | 61 74 69 6f 6e 73 2c 20 |f observ|ations, |
|00002000| 74 20 3d 20 31 2c 20 2e | 2e 2e 2c 20 6e 2e 0a 58 |t = 1, .|.., n..X|
|00002010| 31 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |1 | |
|00002020| 20 20 20 4e 75 6d 62 65 | 72 20 6f 66 20 65 71 75 | Numbe|r of equ|
|00002030| 61 74 69 6f 6e 73 2c 20 | 69 2e 65 2e 20 64 69 6d |ations, |i.e. dim|
|00002040| 65 6e 73 69 6f 6e 20 6f | 66 20 65 2e 0a 58 30 20 |ension o|f e..X0 |
|00002050| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002060| 20 4e 75 6d 62 65 72 20 | 6f 66 20 69 6e 73 74 72 | Number |of instr|
|00002070| 75 6d 65 6e 74 73 2c 20 | 69 2e 65 2e 20 64 69 6d |uments, |i.e. dim|
|00002080| 65 6e 73 69 6f 6e 20 6f | 66 20 5a 2e 0a 58 34 20 |ension o|f Z..X4 |
|00002090| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000020a0| 20 4e 75 6d 62 65 72 20 | 6f 66 20 70 61 72 61 6d | Number |of param|
|000020b0| 65 74 65 72 73 2c 20 69 | 2e 65 2e 20 64 69 6d 65 |eters, i|.e. dime|
|000020c0| 6e 73 69 6f 6e 20 6f 66 | 20 74 68 65 74 61 2e 0a |nsion of| theta..|
|000020d0| 58 35 30 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X50 | |
|000020e0| 20 20 20 20 55 70 70 65 | 72 20 6c 69 6d 69 74 20 | Uppe|r limit |
|000020f0| 6f 6e 20 47 61 75 73 73 | 2d 4e 65 77 74 6f 6e 20 |on Gauss|-Newton |
|00002100| 69 74 65 72 61 74 69 6f | 6e 73 2e 0a 58 31 20 20 |iteratio|ns..X1 |
|00002110| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002120| 4e 75 6d 62 65 72 20 76 | 61 72 20 69 74 65 72 61 |Number v|ar itera|
|00002130| 74 65 73 2c 20 69 76 61 | 72 3d 30 20 6d 65 61 6e |tes, iva|r=0 mean|
|00002140| 73 20 6e 6f 6e 65 2e 0a | 58 68 6f 6d 6f 73 6b 65 |s none..|Xhomoske|
|00002150| 64 61 73 74 69 63 20 20 | 20 20 20 20 43 6f 64 65 |dastic | Code|
|00002160| 20 68 6f 6d 6f 73 6b 65 | 64 61 73 74 69 63 20 6f | homoske|dastic o|
|00002170| 72 20 68 65 74 65 72 6f | 73 6b 65 64 61 73 74 69 |r hetero|skedasti|
|00002180| 63 2e 0a 58 30 20 20 20 | 20 20 20 20 20 20 20 20 |c..X0 | |
|00002190| 20 20 20 20 20 20 20 4e | 75 6d 62 65 72 20 6f 66 | N|umber of|
|000021a0| 20 6d 6f 76 69 6e 67 20 | 61 76 65 72 61 67 65 20 | moving |average |
|000021b0| 74 65 72 6d 73 20 4d 41 | 20 66 6f 72 20 76 61 72 |terms MA| for var|
|000021c0| 20 65 73 74 69 6d 61 74 | 65 2e 0a 58 6e 6f 6e 65 | estimat|e..Xnone|
|000021d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 43 | | C|
|000021e0| 6f 64 65 20 6e 6f 6e 65 | 20 6f 72 20 50 61 72 7a |ode none| or Parz|
|000021f0| 65 6e 2c 20 6e 6f 6e 65 | 20 77 68 65 6e 20 4d 41 |en, none| when MA|
|00002200| 3e 30 20 69 73 20 75 6e | 77 69 73 65 2e 0a 58 31 |>0 is un|wise..X1|
|00002210| 2e 30 30 30 30 30 30 65 | 2d 31 33 20 20 20 20 20 |.000000e|-13 |
|00002220| 20 20 43 6f 6e 76 65 72 | 67 65 6e 63 65 20 74 6f | Conver|gence to|
|00002230| 6c 65 72 61 6e 63 65 2c | 20 74 6f 6c 3d 31 2e 30 |lerance,| tol=1.0|
|00002240| 65 2d 38 20 69 73 20 72 | 65 61 73 6f 6e 61 62 6c |e-8 is r|easonabl|
|00002250| 65 2e 0a 58 31 2e 30 30 | 30 30 30 30 65 2d 31 30 |e..X1.00|0000e-10|
|00002260| 20 20 20 20 20 20 20 49 | 6e 76 65 72 73 69 6f 6e | I|nversion|
|00002270| 20 74 6f 6c 65 72 61 6e | 63 65 2c 20 65 70 73 3d | toleran|ce, eps=|
|00002280| 31 2e 30 65 2d 31 33 20 | 69 73 20 72 65 61 73 6f |1.0e-13 |is reaso|
|00002290| 6e 61 62 6c 65 0a 58 66 | 75 6c 6c 20 20 20 20 20 |nable.Xf|ull |
|000022a0| 20 20 20 20 20 20 20 20 | 20 20 48 6f 77 20 6d 75 | | How mu|
|000022b0| 63 68 20 6f 75 74 70 75 | 74 3f 20 20 43 6f 64 65 |ch outpu|t? Code|
|000022c0| 20 6e 6f 6e 65 2c 20 6d | 69 6e 69 6d 61 6c 2c 20 | none, m|inimal, |
|000022d0| 6f 72 20 66 75 6c 6c 2e | 0a 58 34 20 20 20 20 20 |or full.|.X4 |
|000022e0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 43 6f 6d | | Com|
|000022f0| 70 75 74 65 64 20 72 61 | 6e 6b 20 6f 66 20 56 2e |puted ra|nk of V.|
|00002300| 0a 58 63 6f 72 72 65 63 | 74 65 64 20 20 20 20 20 |.Xcorrec|ted |
|00002310| 20 20 20 20 20 44 69 76 | 69 73 6f 72 20 6f 66 20 | Div|isor of |
|00002320| 76 61 72 2c 20 63 6f 72 | 72 65 63 74 65 64 20 6f |var, cor|rected o|
|00002330| 72 20 75 6e 63 6f 72 72 | 65 63 74 65 64 2e 0a 58 |r uncorr|ected..X|
|00002340| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002350| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|00002360| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002370| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|00002380| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00002390| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|000023a0| 20 20 2d 32 2e 35 38 38 | 39 36 39 38 30 38 33 32 | -2.588|96980832|
|000023b0| 37 38 32 33 31 65 2d 30 | 32 20 20 20 20 20 20 20 |78231e-0|2 |
|000023c0| 20 20 20 2d 30 2e 30 32 | 35 38 38 39 37 30 20 20 | -0.02|588970 |
|000023d0| 74 68 65 74 61 28 31 29 | 0a 58 20 20 20 31 2e 30 |theta(1)|.X 1.0|
|000023e0| 31 35 36 37 39 36 36 39 | 33 31 39 36 38 35 39 65 |15679669|3196859e|
|000023f0| 2b 30 30 20 20 20 20 20 | 20 20 20 20 20 20 31 2e |+00 | 1.|
|00002400| 30 31 35 36 37 39 36 37 | 20 20 74 68 65 74 61 28 |01567967| theta(|
|00002410| 32 29 0a 58 20 20 2d 31 | 2e 31 31 35 36 39 37 31 |2).X -1|.1156971|
|00002420| 33 39 38 39 38 31 33 31 | 30 65 2b 30 30 20 20 20 |39898131|0e+00 |
|00002430| 20 20 20 20 20 20 20 2d | 31 2e 31 31 35 36 39 37 | -|1.115697|
|00002440| 31 34 20 20 74 68 65 74 | 61 28 33 29 0a 58 20 20 |14 thet|a(3).X |
|00002450| 2d 35 2e 30 34 39 30 32 | 38 36 31 35 31 31 35 34 |-5.04902|86151154|
|00002460| 34 39 35 65 2d 30 31 20 | 20 20 20 20 20 20 20 20 |495e-01 | |
|00002470| 20 2d 30 2e 35 30 34 39 | 30 32 38 36 20 20 74 68 | -0.5049|0286 th|
|00002480| 65 74 61 28 34 29 0a 58 | 20 20 20 31 2e 31 37 32 |eta(4).X| 1.172|
|00002490| 39 30 35 32 37 36 36 35 | 35 33 37 39 35 65 2d 30 |90527665|53795e-0|
|000024a0| 33 20 20 20 20 20 20 20 | 20 20 20 20 30 2e 30 30 |3 | 0.00|
|000024b0| 31 31 37 32 39 31 20 20 | 76 61 72 28 31 2c 31 29 |117291 |var(1,1)|
|000024c0| 0a 58 20 20 20 31 2e 35 | 39 33 36 31 33 30 31 31 |.X 1.5|93613011|
|000024d0| 38 36 33 38 32 36 35 65 | 2d 30 34 20 20 20 20 20 |8638265e|-04 |
|000024e0| 20 20 20 20 20 20 30 2e | 30 30 30 31 35 39 33 36 | 0.|00015936|
|000024f0| 20 20 56 28 31 2c 31 29 | 0a 58 20 20 2d 37 2e 38 | V(1,1)|.X -7.8|
|00002500| 37 31 35 37 32 35 31 31 | 32 32 36 32 31 38 30 65 |71572511|2262180e|
|00002510| 2d 30 35 20 20 20 20 20 | 20 20 20 20 20 2d 30 2e |-05 | -0.|
|00002520| 30 30 30 30 37 38 37 32 | 20 20 56 28 32 2c 31 29 |00007872| V(2,1)|
|00002530| 0a 58 20 20 2d 31 2e 37 | 37 31 30 37 30 34 37 37 |.X -1.7|71070477|
|00002540| 33 32 37 37 38 32 33 65 | 2d 30 34 20 20 20 20 20 |3277823e|-04 |
|00002550| 20 20 20 20 20 2d 30 2e | 30 30 30 31 37 37 31 31 | -0.|00017711|
|00002560| 20 20 56 28 33 2c 31 29 | 0a 58 20 20 2d 34 2e 34 | V(3,1)|.X -4.4|
|00002570| 30 39 34 38 33 36 39 36 | 33 39 34 30 39 31 38 65 |09483696|3940918e|
|00002580| 2d 30 35 20 20 20 20 20 | 20 20 20 20 20 2d 30 2e |-05 | -0.|
|00002590| 30 30 30 30 34 34 30 39 | 20 20 56 28 34 2c 31 29 |00004409| V(4,1)|
|000025a0| 0a 58 20 20 2d 37 2e 38 | 37 31 35 37 32 35 31 31 |.X -7.8|71572511|
|000025b0| 32 32 36 32 31 38 30 65 | 2d 30 35 20 20 20 20 20 |2262180e|-05 |
|000025c0| 20 20 20 20 20 2d 30 2e | 30 30 30 30 37 38 37 32 | -0.|00007872|
|000025d0| 20 20 56 28 31 2c 32 29 | 0a 58 20 20 20 39 2e 38 | V(1,2)|.X 9.8|
|000025e0| 37 36 32 34 30 31 36 36 | 34 35 38 31 30 30 34 65 |76240166|4581004e|
|000025f0| 2d 30 35 20 20 20 20 20 | 20 20 20 20 20 20 30 2e |-05 | 0.|
|00002600| 30 30 30 30 39 38 37 36 | 20 20 56 28 32 2c 32 29 |00009876| V(2,2)|
|00002610| 0a 58 20 20 20 36 2e 30 | 37 30 32 34 38 36 30 35 |.X 6.0|70248605|
|00002620| 37 39 37 33 35 31 34 65 | 2d 30 34 20 20 20 20 20 |7973514e|-04 |
|00002630| 20 20 20 20 20 20 30 2e | 30 30 30 36 30 37 30 32 | 0.|00060702|
|00002640| 20 20 56 28 33 2c 32 29 | 0a 58 20 20 2d 31 2e 38 | V(3,2)|.X -1.8|
|00002650| 35 31 33 38 31 30 35 33 | 38 37 38 35 38 30 32 65 |51381053|8785802e|
|00002660| 2d 30 36 20 20 20 20 20 | 20 20 20 20 20 2d 30 2e |-06 | -0.|
|00002670| 30 30 30 30 30 31 38 35 | 20 20 56 28 34 2c 32 29 |00000185| V(4,2)|
|00002680| 0a 58 20 20 2d 31 2e 37 | 37 31 30 37 30 34 37 37 |.X -1.7|71070477|
|00002690| 33 32 37 37 38 32 33 65 | 2d 30 34 20 20 20 20 20 |3277823e|-04 |
|000026a0| 20 20 20 20 20 2d 30 2e | 30 30 30 31 37 37 31 31 | -0.|00017711|
|000026b0| 20 20 56 28 31 2c 33 29 | 0a 58 20 20 20 36 2e 30 | V(1,3)|.X 6.0|
|000026c0| 37 30 32 34 38 36 30 35 | 37 39 37 33 35 31 34 65 |70248605|7973514e|
|000026d0| 2d 30 34 20 20 20 20 20 | 20 20 20 20 20 20 30 2e |-04 | 0.|
|000026e0| 30 30 30 36 30 37 30 32 | 20 20 56 28 32 2c 33 29 |00060702| V(2,3)|
|000026f0| 0a 58 20 20 20 32 2e 36 | 37 34 35 39 38 33 32 33 |.X 2.6|74598323|
|00002700| 37 37 31 31 34 32 38 65 | 2d 30 32 20 20 20 20 20 |7711428e|-02 |
|00002710| 20 20 20 20 20 20 30 2e | 30 32 36 37 34 35 39 38 | 0.|02674598|
|00002720| 20 20 56 28 33 2c 33 29 | 0a 58 20 20 20 32 2e 33 | V(3,3)|.X 2.3|
|00002730| 35 36 32 31 31 39 35 30 | 30 30 39 39 39 37 37 65 |56211950|0099977e|
|00002740| 2d 30 33 20 20 20 20 20 | 20 20 20 20 20 20 30 2e |-03 | 0.|
|00002750| 30 30 32 33 35 36 32 31 | 20 20 56 28 34 2c 33 29 |00235621| V(4,3)|
|00002760| 0a 58 20 20 2d 34 2e 34 | 30 39 34 38 33 36 39 36 |.X -4.4|09483696|
|00002770| 33 39 34 30 39 31 38 65 | 2d 30 35 20 20 20 20 20 |3940918e|-05 |
|00002780| 20 20 20 20 20 2d 30 2e | 30 30 30 30 34 34 30 39 | -0.|00004409|
|00002790| 20 20 56 28 31 2c 34 29 | 0a 58 20 20 2d 31 2e 38 | V(1,4)|.X -1.8|
|000027a0| 35 31 33 38 31 30 35 33 | 38 37 38 35 38 30 32 65 |51381053|8785802e|
|000027b0| 2d 30 36 20 20 20 20 20 | 20 20 20 20 20 2d 30 2e |-06 | -0.|
|000027c0| 30 30 30 30 30 31 38 35 | 20 20 56 28 32 2c 34 29 |00000185| V(2,4)|
|000027d0| 0a 58 20 20 20 32 2e 33 | 35 36 32 31 31 39 35 30 |.X 2.3|56211950|
|000027e0| 30 30 39 39 39 37 37 65 | 2d 30 33 20 20 20 20 20 |0099977e|-03 |
|000027f0| 20 20 20 20 20 20 30 2e | 30 30 32 33 35 36 32 31 | 0.|00235621|
|00002800| 20 20 56 28 33 2c 34 29 | 0a 58 20 20 20 36 2e 35 | V(3,4)|.X 6.5|
|00002810| 38 32 39 32 34 31 35 38 | 38 37 32 31 30 30 38 65 |82924158|8721008e|
|00002820| 2d 30 34 20 20 20 20 20 | 20 20 20 20 20 20 30 2e |-04 | 0.|
|00002830| 30 30 30 36 35 38 32 39 | 20 20 56 28 34 2c 34 29 |00065829| V(4,4)|
|00002840| 0a 58 20 20 2d 33 2e 31 | 31 32 32 36 35 33 31 33 |.X -3.1|12265313|
|00002850| 30 31 38 31 32 35 33 65 | 2d 31 33 20 20 20 20 20 |0181253e|-13 |
|00002860| 20 20 20 20 20 2d 30 2e | 30 30 30 30 30 30 30 30 | -0.|00000000|
|00002870| 20 20 44 28 31 29 0a 58 | 20 20 20 31 2e 30 36 36 | D(1).X| 1.066|
|00002880| 36 30 36 33 31 30 37 34 | 38 39 30 34 38 65 2d 31 |60631074|89048e-1|
|00002890| 32 20 20 20 20 20 20 20 | 20 20 20 20 30 2e 30 30 |2 | 0.00|
|000028a0| 30 30 30 30 30 30 20 20 | 44 28 32 29 0a 58 20 20 |000000 |D(2).X |
|000028b0| 20 34 2e 36 39 39 39 33 | 39 36 38 30 31 34 39 38 | 4.69993|96801498|
|000028c0| 33 30 32 65 2d 31 31 20 | 20 20 20 20 20 20 20 20 |302e-11 | |
|000028d0| 20 20 30 2e 30 30 30 30 | 30 30 30 30 20 20 44 28 | 0.0000|0000 D(|
|000028e0| 33 29 0a 58 20 20 20 34 | 2e 31 34 30 35 30 38 35 |3).X 4|.1405085|
|000028f0| 39 36 37 39 37 32 33 34 | 33 65 2d 31 32 20 20 20 |96797234|3e-12 |
|00002900| 20 20 20 20 20 20 20 20 | 30 2e 30 30 30 30 30 30 | |0.000000|
|00002910| 30 30 20 20 44 28 34 29 | 0a 58 20 20 20 32 2e 35 |00 D(4)|.X 2.5|
|00002920| 39 39 39 39 39 39 39 39 | 39 39 39 39 39 39 36 65 |99999999|9999996e|
|00002930| 2b 30 31 20 20 20 20 20 | 20 20 20 20 20 32 36 2e |+01 | 26.|
|00002940| 30 30 30 30 30 30 30 30 | 20 20 6f 62 6a 0a 45 4e |00000000| obj.EN|
|00002950| 44 5f 4f 46 5f 46 49 4c | 45 0a 69 66 20 74 65 73 |D_OF_FIL|E.if tes|
|00002960| 74 20 32 35 32 33 20 2d | 6e 65 20 60 77 63 20 2d |t 2523 -|ne `wc -|
|00002970| 63 20 3c 27 63 68 31 65 | 67 31 2f 65 6e 64 69 6e |c <'ch1e|g1/endin|
|00002980| 67 2e 62 61 6b 27 60 3b | 20 74 68 65 6e 0a 20 20 |g.bak'`;| then. |
|00002990| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 27 | echo s|har: \"'|
|000029a0| 63 68 31 65 67 31 2f 65 | 6e 64 69 6e 67 2e 62 61 |ch1eg1/e|nding.ba|
|000029b0| 6b 27 5c 22 20 75 6e 70 | 61 63 6b 65 64 20 77 69 |k'\" unp|acked wi|
|000029c0| 74 68 20 77 72 6f 6e 67 | 20 73 69 7a 65 21 0a 66 |th wrong| size!.f|
|000029d0| 69 0a 23 20 65 6e 64 20 | 6f 66 20 27 63 68 31 65 |i.# end |of 'ch1e|
|000029e0| 67 31 2f 65 6e 64 69 6e | 67 2e 62 61 6b 27 0a 66 |g1/endin|g.bak'.f|
|000029f0| 69 0a 69 66 20 74 65 73 | 74 20 2d 66 20 27 63 68 |i.if tes|t -f 'ch|
|00002a00| 31 65 67 31 2f 6d 6f 64 | 65 6c 2e 63 63 27 20 2d |1eg1/mod|el.cc' -|
|00002a10| 61 20 22 24 7b 31 7d 22 | 20 21 3d 20 22 2d 63 22 |a "${1}"| != "-c"|
|00002a20| 20 3b 20 74 68 65 6e 20 | 0a 20 20 65 63 68 6f 20 | ; then |. echo |
|00002a30| 73 68 61 72 3a 20 57 69 | 6c 6c 20 6e 6f 74 20 63 |shar: Wi|ll not c|
|00002a40| 6c 6f 62 62 65 72 20 65 | 78 69 73 74 69 6e 67 20 |lobber e|xisting |
|00002a50| 66 69 6c 65 20 5c 22 27 | 63 68 31 65 67 31 2f 6d |file \"'|ch1eg1/m|
|00002a60| 6f 64 65 6c 2e 63 63 27 | 5c 22 0a 65 6c 73 65 0a |odel.cc'|\".else.|
|00002a70| 65 63 68 6f 20 73 68 61 | 72 3a 20 45 78 74 72 61 |echo sha|r: Extra|
|00002a80| 63 74 69 6e 67 20 5c 22 | 27 63 68 31 65 67 31 2f |cting \"|'ch1eg1/|
|00002a90| 6d 6f 64 65 6c 2e 63 63 | 27 5c 22 20 5c 28 31 37 |model.cc|'\" \(17|
|00002aa0| 39 34 20 63 68 61 72 61 | 63 74 65 72 73 5c 29 0a |94 chara|cters\).|
|00002ab0| 73 65 64 20 22 73 2f 5e | 58 2f 2f 22 20 3e 27 63 |sed "s/^|X//" >'c|
|00002ac0| 68 31 65 67 31 2f 6d 6f | 64 65 6c 2e 63 63 27 20 |h1eg1/mo|del.cc' |
|00002ad0| 3c 3c 27 45 4e 44 5f 4f | 46 5f 46 49 4c 45 27 0a |<<'END_O|F_FILE'.|
|00002ae0| 58 23 69 6e 63 6c 75 64 | 65 20 22 6d 6f 64 65 6c |X#includ|e "model|
|00002af0| 2e 68 22 0a 58 0a 58 6d | 6f 64 65 6c 3a 3a 6d 6f |.h".X.Xm|odel::mo|
|00002b00| 64 65 6c 28 29 20 7b 20 | 7d 0a 58 0a 58 6d 6f 64 |del() { |}.X.Xmod|
|00002b10| 65 6c 3a 3a 7e 6d 6f 64 | 65 6c 28 29 20 7b 20 7d |el::~mod|el() { }|
|00002b20| 0a 58 0a 58 0a 58 72 65 | 61 6c 6d 61 74 20 6d 6f |.X.X.Xre|almat mo|
|00002b30| 64 65 6c 3a 3a 65 28 49 | 4e 54 45 47 45 52 20 74 |del::e(I|NTEGER t|
|00002b40| 29 0a 58 7b 0a 58 20 20 | 72 65 61 6c 6d 61 74 20 |).X{.X |realmat |
|00002b50| 65 72 72 28 31 2c 31 29 | 3b 0a 58 0a 58 20 20 65 |err(1,1)|;.X.X e|
|00002b60| 72 72 2e 65 6c 65 6d 28 | 31 2c 31 29 0a 58 20 20 |rr.elem(|1,1).X |
|00002b70| 20 20 3d 20 64 61 74 61 | 2e 65 6c 65 6d 28 74 2c | = data|.elem(t,|
|00002b80| 32 29 0a 58 20 20 20 20 | 20 20 2d 20 64 61 74 61 |2).X | - data|
|00002b90| 2e 65 6c 65 6d 28 74 2c | 33 29 2a 73 2e 74 68 65 |.elem(t,|3)*s.the|
|00002ba0| 74 61 2e 65 6c 65 6d 28 | 31 2c 31 29 0a 58 20 20 |ta.elem(|1,1).X |
|00002bb0| 20 20 20 20 2d 20 64 61 | 74 61 2e 65 6c 65 6d 28 | - da|ta.elem(|
|00002bc0| 74 2c 34 29 2a 73 2e 74 | 68 65 74 61 2e 65 6c 65 |t,4)*s.t|heta.ele|
|00002bd0| 6d 28 32 2c 31 29 0a 58 | 20 20 20 20 20 20 2d 20 |m(2,1).X| - |
|00002be0| 73 2e 74 68 65 74 61 2e | 65 6c 65 6d 28 34 2c 31 |s.theta.|elem(4,1|
|00002bf0| 29 2a 65 78 70 28 64 61 | 74 61 2e 65 6c 65 6d 28 |)*exp(da|ta.elem(|
|00002c00| 74 2c 35 29 2a 73 2e 74 | 68 65 74 61 2e 65 6c 65 |t,5)*s.t|heta.ele|
|00002c10| 6d 28 33 2c 31 29 29 3b | 0a 58 0a 58 20 20 72 65 |m(3,1));|.X.X re|
|00002c20| 74 75 72 6e 20 65 72 72 | 3b 0a 58 7d 0a 58 0a 58 |turn err|;.X}.X.X|
|00002c30| 72 65 61 6c 6d 61 74 20 | 6d 6f 64 65 6c 3a 3a 64 |realmat |model::d|
|00002c40| 65 6c 65 28 49 4e 54 45 | 47 45 52 20 74 29 20 0a |ele(INTE|GER t) .|
|00002c50| 58 7b 0a 58 20 20 72 65 | 61 6c 6d 61 74 20 64 65 |X{.X re|almat de|
|00002c60| 6c 65 72 72 28 31 2c 34 | 29 3b 0a 58 20 20 52 45 |lerr(1,4|);.X RE|
|00002c70| 41 4c 20 20 20 20 74 6d | 70 3b 0a 58 0a 58 20 20 |AL tm|p;.X.X |
|00002c80| 74 6d 70 20 3d 20 65 78 | 70 28 64 61 74 61 2e 65 |tmp = ex|p(data.e|
|00002c90| 6c 65 6d 28 74 2c 35 29 | 2a 73 2e 74 68 65 74 61 |lem(t,5)|*s.theta|
|00002ca0| 2e 65 6c 65 6d 28 33 2c | 31 29 29 3b 0a 58 0a 58 |.elem(3,|1));.X.X|
|00002cb0| 20 20 64 65 6c 65 72 72 | 2e 65 6c 65 6d 28 31 2c | delerr|.elem(1,|
|00002cc0| 31 29 20 3d 20 2d 20 64 | 61 74 61 2e 65 6c 65 6d |1) = - d|ata.elem|
|00002cd0| 28 74 2c 33 29 3b 0a 58 | 20 20 64 65 6c 65 72 72 |(t,3);.X| delerr|
|00002ce0| 2e 65 6c 65 6d 28 31 2c | 32 29 20 3d 20 2d 20 64 |.elem(1,|2) = - d|
|00002cf0| 61 74 61 2e 65 6c 65 6d | 28 74 2c 34 29 3b 0a 58 |ata.elem|(t,4);.X|
|00002d00| 20 20 64 65 6c 65 72 72 | 2e 65 6c 65 6d 28 31 2c | delerr|.elem(1,|
|00002d10| 33 29 20 3d 20 2d 20 73 | 2e 74 68 65 74 61 2e 65 |3) = - s|.theta.e|
|00002d20| 6c 65 6d 28 34 2c 31 29 | 2a 64 61 74 61 2e 65 6c |lem(4,1)|*data.el|
|00002d30| 65 6d 28 74 2c 35 29 2a | 74 6d 70 3b 0a 58 20 20 |em(t,5)*|tmp;.X |
|00002d40| 64 65 6c 65 72 72 2e 65 | 6c 65 6d 28 31 2c 34 29 |delerr.e|lem(1,4)|
|00002d50| 20 3d 20 2d 20 74 6d 70 | 3b 0a 58 0a 58 20 20 72 | = - tmp|;.X.X r|
|00002d60| 65 74 75 72 6e 20 64 65 | 6c 65 72 72 3b 0a 58 0a |eturn de|lerr;.X.|
|00002d70| 58 7d 0a 58 0a 58 0a 58 | 69 6e 74 20 6d 6f 64 65 |X}.X.X.X|int mode|
|00002d80| 6c 3a 3a 69 6e 69 74 69 | 61 6c 69 7a 65 28 29 0a |l::initi|alize().|
|00002d90| 58 7b 0a 58 23 69 66 64 | 65 66 20 47 4e 55 5f 47 |X{.X#ifd|ef GNU_G|
|00002da0| 50 50 5f 43 4f 4d 50 49 | 4c 45 52 0a 58 0a 58 23 |PP_COMPI|LER.X.X#|
|00002db0| 69 66 64 65 66 20 55 53 | 45 5f 41 54 54 5f 53 54 |ifdef US|E_ATT_ST|
|00002dc0| 59 4c 45 5f 49 4f 5f 57 | 49 54 48 5f 47 4e 55 0a |YLE_IO_W|ITH_GNU.|
|00002dd0| 58 20 20 66 69 6c 65 62 | 75 66 20 64 61 74 61 62 |X fileb|uf datab|
|00002de0| 75 66 3b 0a 58 20 20 69 | 66 28 20 64 61 74 61 62 |uf;.X i|f( datab|
|00002df0| 75 66 2e 6f 70 65 6e 28 | 22 61 6d 73 74 61 74 2e |uf.open(|"amstat.|
|00002e00| 64 61 74 22 2c 20 69 6e | 70 75 74 29 20 3d 3d 20 |dat", in|put) == |
|00002e10| 30 20 29 7b 0a 58 20 20 | 20 20 63 65 72 72 20 3c |0 ){.X | cerr <|
|00002e20| 3c 20 22 45 72 72 6f 72 | 2c 20 6d 6f 64 65 6c 3a |< "Error|, model:|
|00002e30| 3a 6d 6f 64 65 6c 2c 20 | 43 61 6e 6e 6f 74 20 6f |:model, |Cannot o|
|00002e40| 70 65 6e 20 61 6d 73 74 | 61 74 2e 64 61 74 5c 6e |pen amst|at.dat\n|
|00002e50| 22 3b 0a 58 20 20 20 20 | 65 78 69 74 28 31 29 3b |";.X |exit(1);|
|00002e60| 0a 58 20 20 7d 0a 58 20 | 20 69 73 74 72 65 61 6d |.X }.X | istream|
|00002e70| 20 64 61 74 61 67 65 74 | 28 26 64 61 74 61 62 75 | dataget|(&databu|
|00002e80| 66 29 3b 0a 58 23 65 6e | 64 69 66 20 0a 58 0a 58 |f);.X#en|dif .X.X|
|00002e90| 23 69 66 64 65 66 20 55 | 53 45 5f 47 4e 55 5f 53 |#ifdef U|SE_GNU_S|
|00002ea0| 54 59 4c 45 5f 49 4f 5f | 57 49 54 48 5f 47 4e 55 |TYLE_IO_|WITH_GNU|
|00002eb0| 0a 58 20 20 69 73 74 72 | 65 61 6d 20 64 61 74 61 |.X istr|eam data|
|00002ec0| 67 65 74 28 22 61 6d 73 | 74 61 74 2e 64 61 74 22 |get("ams|tat.dat"|
|00002ed0| 2c 69 6f 5f 72 65 61 64 | 6f 6e 6c 79 2c 61 5f 75 |,io_read|only,a_u|
|00002ee0| 73 65 6f 6e 6c 79 29 3b | 0a 58 20 20 69 66 20 28 |seonly);|.X if (|
|00002ef0| 21 64 61 74 61 67 65 74 | 29 20 7b 0a 58 20 20 20 |!dataget|) {.X |
|00002f00| 20 63 65 72 72 20 3c 3c | 20 22 45 72 72 6f 72 2c | cerr <<| "Error,|
|00002f10| 20 6d 6f 64 65 6c 3a 3a | 6d 6f 64 65 6c 5c 6e 22 | model::|model\n"|
|00002f20| 3b 0a 58 20 20 20 20 65 | 78 69 74 28 31 29 3b 0a |;.X e|xit(1);.|
|00002f30| 58 20 20 7d 20 0a 58 23 | 65 6e 64 69 66 20 0a 58 |X } .X#|endif .X|
|00002f40| 0a 58 23 65 6e 64 69 66 | 20 0a 58 0a 58 23 69 66 |.X#endif| .X.X#if|
|00002f50| 64 65 66 20 54 55 52 42 | 4f 5f 43 50 50 5f 43 4f |def TURB|O_CPP_CO|
|00002f60| 4d 50 49 4c 45 52 0a 58 | 20 20 69 66 73 74 72 65 |MPILER.X| ifstre|
|00002f70| 61 6d 20 64 61 74 61 67 | 65 74 28 22 61 6d 73 74 |am datag|et("amst|
|00002f80| 61 74 2e 64 61 74 22 29 | 3b 0a 58 20 20 69 66 20 |at.dat")|;.X if |
|00002f90| 28 21 64 61 74 61 67 65 | 74 29 20 7b 0a 58 20 20 |(!datage|t) {.X |
|00002fa0| 20 20 63 65 72 72 20 3c | 3c 20 22 45 72 72 6f 72 | cerr <|< "Error|
|00002fb0| 2c 20 6d 6f 64 65 6c 3a | 3a 6d 6f 64 65 6c 5c 6e |, model:|:model\n|
|00002fc0| 22 3b 0a 58 20 20 20 20 | 65 78 69 74 28 31 29 3b |";.X |exit(1);|
|00002fd0| 0a 58 20 20 7d 20 0a 58 | 23 65 6e 64 69 66 20 0a |.X } .X|#endif .|
|00002fe0| 58 0a 58 20 20 64 61 74 | 61 2e 72 65 73 69 7a 65 |X.X dat|a.resize|
|00002ff0| 28 33 30 2c 35 29 3b 0a | 58 20 20 49 4e 54 45 47 |(30,5);.|X INTEG|
|00003000| 45 52 20 69 2c 6a 3b 0a | 58 0a 58 20 20 64 6f 75 |ER i,j;.|X.X dou|
|00003010| 62 6c 65 20 78 3b 0a 58 | 20 20 66 6f 72 20 28 69 |ble x;.X| for (i|
|00003020| 3d 31 3b 20 69 3c 3d 33 | 30 3b 20 69 2b 2b 29 20 |=1; i<=3|0; i++) |
|00003030| 7b 20 20 0a 58 20 20 66 | 6f 72 20 28 6a 3d 31 3b |{ .X f|or (j=1;|
|00003040| 20 6a 3c 3d 35 3b 20 20 | 6a 2b 2b 29 20 7b 20 20 | j<=5; |j++) { |
|00003050| 20 20 0a 58 20 20 20 20 | 64 61 74 61 67 65 74 20 | .X |dataget |
|00003060| 3e 3e 20 78 3b 20 0a 58 | 20 20 20 20 64 61 74 61 |>> x; .X| data|
|00003070| 2e 65 6c 65 6d 28 69 2c | 6a 29 20 3d 20 28 52 45 |.elem(i,|j) = (RE|
|00003080| 41 4c 29 78 3b 0a 58 20 | 20 7d 0a 58 20 20 7d 0a |AL)x;.X | }.X }.|
|00003090| 58 0a 58 20 20 69 66 28 | 73 74 72 63 6d 70 28 73 |X.X if(|strcmp(s|
|000030a0| 2e 64 65 74 61 69 6c 2c | 22 66 75 6c 6c 22 29 3d |.detail,|"full")=|
|000030b0| 3d 30 29 20 63 6f 75 74 | 20 3c 3c 73 74 61 72 62 |=0) cout| <<starb|
|000030c0| 6f 78 28 22 2f 64 61 74 | 61 2f 2f 5f 22 29 3c 3c |ox("/dat|a//_")<<|
|000030d0| 20 64 61 74 61 3b 0a 58 | 0a 58 20 20 64 61 74 61 | data;.X|.X data|
|000030e0| 67 65 74 2e 63 6c 6f 73 | 65 28 29 3b 0a 58 0a 58 |get.clos|e();.X.X|
|000030f0| 20 20 72 65 74 75 72 6e | 20 31 3b 0a 58 7d 0a 58 | return| 1;.X}.X|
|00003100| 0a 58 72 65 61 6c 6d 61 | 74 20 6d 6f 64 65 6c 3a |.Xrealma|t model:|
|00003110| 3a 5a 28 49 4e 54 45 47 | 45 52 20 74 29 0a 58 7b |:Z(INTEG|ER t).X{|
|00003120| 0a 58 20 20 72 65 61 6c | 6d 61 74 20 7a 28 73 2e |.X real|mat z(s.|
|00003130| 4b 2c 31 29 3b 0a 58 20 | 20 72 65 61 6c 6d 61 74 |K,1);.X | realmat|
|00003140| 20 74 6d 70 28 73 2e 70 | 2c 31 29 3b 0a 58 20 20 | tmp(s.p|,1);.X |
|00003150| 74 6d 70 20 3d 20 64 65 | 6c 65 28 74 29 3b 0a 58 |tmp = de|le(t);.X|
|00003160| 20 20 66 6f 72 20 28 49 | 4e 54 45 47 45 52 20 69 | for (I|NTEGER i|
|00003170| 20 3d 20 31 3b 20 69 3c | 3d 73 2e 4b 3b 20 69 2b | = 1; i<|=s.K; i+|
|00003180| 2b 29 0a 58 20 20 20 20 | 69 66 20 28 69 20 3c 3d |+).X |if (i <=|
|00003190| 20 73 2e 70 29 20 0a 58 | 20 20 20 20 20 20 7a 5b | s.p) .X| z[|
|000031a0| 69 5d 20 3d 20 74 6d 70 | 5b 69 5d 3b 20 20 2f 2f |i] = tmp|[i]; //|
|000031b0| 20 77 69 6c 6c 20 70 72 | 6f 64 75 63 65 20 53 55 | will pr|oduce SU|
|000031c0| 52 20 65 73 74 69 6d 61 | 74 65 73 20 77 68 65 6e |R estima|tes when|
|000031d0| 20 73 2e 4b 20 3d 20 73 | 2e 70 20 0a 58 20 20 20 | s.K = s|.p .X |
|000031e0| 20 65 6c 73 65 20 0a 58 | 20 20 20 20 20 20 7a 5b | else .X| z[|
|000031f0| 69 5d 3d 70 6f 77 28 64 | 61 74 61 2e 65 6c 65 6d |i]=pow(d|ata.elem|
|00003200| 28 74 2c 35 29 2c 69 2d | 73 2e 70 29 3b 0a 58 20 |(t,5),i-|s.p);.X |
|00003210| 20 72 65 74 75 72 6e 20 | 7a 3b 0a 58 7d 0a 58 0a | return |z;.X}.X.|
|00003220| 58 69 6e 74 20 6d 6f 64 | 65 6c 3a 3a 74 65 72 6d |Xint mod|el::term|
|00003230| 69 6e 61 74 65 28 29 20 | 7b 20 72 65 74 75 72 6e |inate() |{ return|
|00003240| 20 30 3b 20 7d 0a 45 4e | 44 5f 4f 46 5f 46 49 4c | 0; }.EN|D_OF_FIL|
|00003250| 45 0a 69 66 20 74 65 73 | 74 20 31 37 39 34 20 2d |E.if tes|t 1794 -|
|00003260| 6e 65 20 60 77 63 20 2d | 63 20 3c 27 63 68 31 65 |ne `wc -|c <'ch1e|
|00003270| 67 31 2f 6d 6f 64 65 6c | 2e 63 63 27 60 3b 20 74 |g1/model|.cc'`; t|
|00003280| 68 65 6e 0a 20 20 20 20 | 65 63 68 6f 20 73 68 61 |hen. |echo sha|
|00003290| 72 3a 20 5c 22 27 63 68 | 31 65 67 31 2f 6d 6f 64 |r: \"'ch|1eg1/mod|
|000032a0| 65 6c 2e 63 63 27 5c 22 | 20 75 6e 70 61 63 6b 65 |el.cc'\"| unpacke|
|000032b0| 64 20 77 69 74 68 20 77 | 72 6f 6e 67 20 73 69 7a |d with w|rong siz|
|000032c0| 65 21 0a 66 69 0a 23 20 | 65 6e 64 20 6f 66 20 27 |e!.fi.# |end of '|
|000032d0| 63 68 31 65 67 31 2f 6d | 6f 64 65 6c 2e 63 63 27 |ch1eg1/m|odel.cc'|
|000032e0| 0a 66 69 0a 69 66 20 74 | 65 73 74 20 2d 66 20 27 |.fi.if t|est -f '|
|000032f0| 63 68 31 65 67 31 2f 6d | 6f 64 65 6c 2e 68 27 20 |ch1eg1/m|odel.h' |
|00003300| 2d 61 20 22 24 7b 31 7d | 22 20 21 3d 20 22 2d 63 |-a "${1}|" != "-c|
|00003310| 22 20 3b 20 74 68 65 6e | 20 0a 20 20 65 63 68 6f |" ; then| . echo|
|00003320| 20 73 68 61 72 3a 20 57 | 69 6c 6c 20 6e 6f 74 20 | shar: W|ill not |
|00003330| 63 6c 6f 62 62 65 72 20 | 65 78 69 73 74 69 6e 67 |clobber |existing|
|00003340| 20 66 69 6c 65 20 5c 22 | 27 63 68 31 65 67 31 2f | file \"|'ch1eg1/|
|00003350| 6d 6f 64 65 6c 2e 68 27 | 5c 22 0a 65 6c 73 65 0a |model.h'|\".else.|
|00003360| 65 63 68 6f 20 73 68 61 | 72 3a 20 45 78 74 72 61 |echo sha|r: Extra|
|00003370| 63 74 69 6e 67 20 5c 22 | 27 63 68 31 65 67 31 2f |cting \"|'ch1eg1/|
|00003380| 6d 6f 64 65 6c 2e 68 27 | 5c 22 20 5c 28 33 36 33 |model.h'|\" \(363|
|00003390| 20 63 68 61 72 61 63 74 | 65 72 73 5c 29 0a 73 65 | charact|ers\).se|
|000033a0| 64 20 22 73 2f 5e 58 2f | 2f 22 20 3e 27 63 68 31 |d "s/^X/|/" >'ch1|
|000033b0| 65 67 31 2f 6d 6f 64 65 | 6c 2e 68 27 20 3c 3c 27 |eg1/mode|l.h' <<'|
|000033c0| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 27 0a 58 23 69 |END_OF_F|ILE'.X#i|
|000033d0| 66 6e 64 65 66 20 5f 5f | 46 49 4c 45 5f 4d 4f 44 |fndef __|FILE_MOD|
|000033e0| 45 4c 5f 48 5f 53 45 45 | 4e 5f 5f 0a 58 23 70 72 |EL_H_SEE|N__.X#pr|
|000033f0| 61 67 6d 61 20 6f 6e 63 | 65 0a 58 23 64 65 66 69 |agma onc|e.X#defi|
|00003400| 6e 65 20 5f 5f 46 49 4c | 45 5f 4d 4f 44 45 4c 5f |ne __FIL|E_MODEL_|
|00003410| 48 5f 53 45 45 4e 5f 5f | 20 31 0a 58 0a 58 23 69 |H_SEEN__| 1.X.X#i|
|00003420| 6e 63 6c 75 64 65 20 22 | 73 74 61 74 75 73 2e 68 |nclude "|status.h|
|00003430| 22 0a 58 0a 58 65 78 74 | 65 72 6e 20 73 74 61 74 |".X.Xext|ern stat|
|00003440| 75 73 20 73 3b 0a 58 0a | 58 63 6c 61 73 73 20 6d |us s;.X.|Xclass m|
|00003450| 6f 64 65 6c 0a 58 7b 0a | 58 0a 58 70 72 6f 74 65 |odel.X{.|X.Xprote|
|00003460| 63 74 65 64 3a 0a 58 20 | 20 72 65 61 6c 6d 61 74 |cted:.X | realmat|
|00003470| 20 20 20 64 61 74 61 3b | 0a 58 0a 58 70 75 62 6c | data;|.X.Xpubl|
|00003480| 69 63 3a 0a 58 20 20 20 | 20 20 20 20 20 20 20 20 |ic:.X | |
|00003490| 20 6d 6f 64 65 6c 28 29 | 3b 0a 58 20 20 20 20 20 | model()|;.X |
|000034a0| 20 20 20 20 20 20 20 7e | 6d 6f 64 65 6c 28 29 3b | ~|model();|
|000034b0| 0a 58 0a 58 20 20 69 6e | 74 20 20 20 20 20 20 20 |.X.X in|t |
|000034c0| 69 6e 69 74 69 61 6c 69 | 7a 65 28 29 3b 0a 58 20 |initiali|ze();.X |
|000034d0| 20 69 6e 74 20 20 20 20 | 20 20 20 74 65 72 6d 69 | int | termi|
|000034e0| 6e 61 74 65 28 29 3b 0a | 58 0a 58 20 20 72 65 61 |nate();.|X.X rea|
|000034f0| 6c 6d 61 74 20 20 20 65 | 28 49 4e 54 45 47 45 52 |lmat e|(INTEGER|
|00003500| 20 74 29 3b 20 20 20 20 | 0a 58 20 20 72 65 61 6c | t); |.X real|
|00003510| 6d 61 74 20 20 20 64 65 | 6c 65 28 49 4e 54 45 47 |mat de|le(INTEG|
|00003520| 45 52 20 74 29 3b 0a 58 | 20 20 72 65 61 6c 6d 61 |ER t);.X| realma|
|00003530| 74 20 20 20 5a 28 49 4e | 54 45 47 45 52 20 74 29 |t Z(IN|TEGER t)|
|00003540| 3b 20 20 0a 58 0a 58 7d | 3b 0a 58 0a 58 23 65 6e |; .X.X}|;.X.X#en|
|00003550| 64 69 66 0a 45 4e 44 5f | 4f 46 5f 46 49 4c 45 0a |dif.END_|OF_FILE.|
|00003560| 69 66 20 74 65 73 74 20 | 33 36 33 20 2d 6e 65 20 |if test |363 -ne |
|00003570| 60 77 63 20 2d 63 20 3c | 27 63 68 31 65 67 31 2f |`wc -c <|'ch1eg1/|
|00003580| 6d 6f 64 65 6c 2e 68 27 | 60 3b 20 74 68 65 6e 0a |model.h'|`; then.|
|00003590| 20 20 20 20 65 63 68 6f | 20 73 68 61 72 3a 20 5c | echo| shar: \|
|000035a0| 22 27 63 68 31 65 67 31 | 2f 6d 6f 64 65 6c 2e 68 |"'ch1eg1|/model.h|
|000035b0| 27 5c 22 20 75 6e 70 61 | 63 6b 65 64 20 77 69 74 |'\" unpa|cked wit|
|000035c0| 68 20 77 72 6f 6e 67 20 | 73 69 7a 65 21 0a 66 69 |h wrong |size!.fi|
|000035d0| 0a 23 20 65 6e 64 20 6f | 66 20 27 63 68 31 65 67 |.# end o|f 'ch1eg|
|000035e0| 31 2f 6d 6f 64 65 6c 2e | 68 27 0a 66 69 0a 69 66 |1/model.|h'.fi.if|
|000035f0| 20 74 65 73 74 20 2d 66 | 20 27 63 68 31 65 67 31 | test -f| 'ch1eg1|
|00003600| 2f 73 74 61 72 74 69 6e | 67 2e 64 61 74 27 20 2d |/startin|g.dat' -|
|00003610| 61 20 22 24 7b 31 7d 22 | 20 21 3d 20 22 2d 63 22 |a "${1}"| != "-c"|
|00003620| 20 3b 20 74 68 65 6e 20 | 0a 20 20 65 63 68 6f 20 | ; then |. echo |
|00003630| 73 68 61 72 3a 20 57 69 | 6c 6c 20 6e 6f 74 20 63 |shar: Wi|ll not c|
|00003640| 6c 6f 62 62 65 72 20 65 | 78 69 73 74 69 6e 67 20 |lobber e|xisting |
|00003650| 66 69 6c 65 20 5c 22 27 | 63 68 31 65 67 31 2f 73 |file \"'|ch1eg1/s|
|00003660| 74 61 72 74 69 6e 67 2e | 64 61 74 27 5c 22 0a 65 |tarting.|dat'\".e|
|00003670| 6c 73 65 0a 65 63 68 6f | 20 73 68 61 72 3a 20 45 |lse.echo| shar: E|
|00003680| 78 74 72 61 63 74 69 6e | 67 20 5c 22 27 63 68 31 |xtractin|g \"'ch1|
|00003690| 65 67 31 2f 73 74 61 72 | 74 69 6e 67 2e 64 61 74 |eg1/star|ting.dat|
|000036a0| 27 5c 22 20 5c 28 31 30 | 39 37 20 63 68 61 72 61 |'\" \(10|97 chara|
|000036b0| 63 74 65 72 73 5c 29 0a | 73 65 64 20 22 73 2f 5e |cters\).|sed "s/^|
|000036c0| 58 2f 2f 22 20 3e 27 63 | 68 31 65 67 31 2f 73 74 |X//" >'c|h1eg1/st|
|000036d0| 61 72 74 69 6e 67 2e 64 | 61 74 27 20 3c 3c 27 45 |arting.d|at' <<'E|
|000036e0| 4e 44 5f 4f 46 5f 46 49 | 4c 45 27 0a 58 47 61 6c |ND_OF_FI|LE'.XGal|
|000036f0| 6c 61 6e 74 2c 20 22 4e | 6f 6e 6c 69 6e 65 61 72 |lant, "N|onlinear|
|00003700| 20 53 74 61 74 69 73 74 | 69 63 61 6c 20 4d 6f 64 | Statist|ical Mod|
|00003710| 65 6c 73 2c 22 20 43 68 | 61 70 74 65 72 20 31 2c |els," Ch|apter 1,|
|00003720| 20 46 69 67 75 72 65 20 | 35 2c 20 70 2e 20 33 35 | Figure |5, p. 35|
|00003730| 2d 33 36 2e 0a 58 53 55 | 52 20 20 20 20 20 20 20 |-36..XSU|R |
|00003740| 20 20 20 20 20 20 20 20 | 20 57 68 61 74 20 65 73 | | What es|
|00003750| 74 69 6d 61 74 69 6f 6e | 20 6d 65 74 68 6f 64 3f |timation| method?|
|00003760| 20 20 43 6f 64 65 20 53 | 55 52 2c 20 54 53 4c 53 | Code S|UR, TSLS|
|00003770| 2c 20 6f 72 20 47 4d 4d | 2e 0a 58 33 30 20 20 20 |, or GMM|..X30 |
|00003780| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 4e 75 | | Nu|
|00003790| 6d 62 65 72 20 6f 66 20 | 6f 62 73 65 72 76 61 74 |mber of |observat|
|000037a0| 69 6f 6e 73 2c 20 74 20 | 3d 20 31 2c 20 2e 2e 2e |ions, t |= 1, ...|
|000037b0| 2c 20 6e 2e 0a 58 31 20 | 20 20 20 20 20 20 20 20 |, n..X1 | |
|000037c0| 20 20 20 20 20 20 20 20 | 20 4e 75 6d 62 65 72 20 | | Number |
|000037d0| 6f 66 20 65 71 75 61 74 | 69 6f 6e 73 2c 20 69 2e |of equat|ions, i.|
|000037e0| 65 2e 20 64 69 6d 65 6e | 73 69 6f 6e 20 6f 66 20 |e. dimen|sion of |
|000037f0| 65 2e 0a 58 30 20 20 20 | 20 20 20 20 20 20 20 20 |e..X0 | |
|00003800| 20 20 20 20 20 20 20 4e | 75 6d 62 65 72 20 6f 66 | N|umber of|
|00003810| 20 69 6e 73 74 72 75 6d | 65 6e 74 73 2c 20 69 2e | instrum|ents, i.|
|00003820| 65 2e 20 64 69 6d 65 6e | 73 69 6f 6e 20 6f 66 20 |e. dimen|sion of |
|00003830| 5a 2e 0a 58 34 20 20 20 | 20 20 20 20 20 20 20 20 |Z..X4 | |
|00003840| 20 20 20 20 20 20 20 4e | 75 6d 62 65 72 20 6f 66 | N|umber of|
|00003850| 20 70 61 72 61 6d 65 74 | 65 72 73 2c 20 69 2e 65 | paramet|ers, i.e|
|00003860| 2e 20 64 69 6d 65 6e 73 | 69 6f 6e 20 6f 66 20 74 |. dimens|ion of t|
|00003870| 68 65 74 61 2e 0a 58 35 | 30 20 20 20 20 20 20 20 |heta..X5|0 |
|00003880| 20 20 20 20 20 20 20 20 | 20 20 55 70 70 65 72 20 | | Upper |
|00003890| 6c 69 6d 69 74 20 6f 6e | 20 47 61 75 73 73 2d 4e |limit on| Gauss-N|
|000038a0| 65 77 74 6f 6e 20 69 74 | 65 72 61 74 69 6f 6e 73 |ewton it|erations|
|000038b0| 2e 0a 58 31 20 20 20 20 | 20 20 20 20 20 20 20 20 |..X1 | |
|000038c0| 20 20 20 20 20 20 4e 75 | 6d 62 65 72 20 76 61 72 | Nu|mber var|
|000038d0| 20 69 74 65 72 61 74 65 | 73 2c 20 69 76 61 72 3d | iterate|s, ivar=|
|000038e0| 30 20 6d 65 61 6e 73 20 | 6e 6f 6e 65 2e 0a 58 68 |0 means |none..Xh|
|000038f0| 6f 6d 6f 73 6b 65 64 61 | 73 74 69 63 20 20 20 20 |omoskeda|stic |
|00003900| 20 20 43 6f 64 65 20 68 | 6f 6d 6f 73 6b 65 64 61 | Code h|omoskeda|
|00003910| 73 74 69 63 20 6f 72 20 | 68 65 74 65 72 6f 73 6b |stic or |heterosk|
|00003920| 65 64 61 73 74 69 63 2e | 0a 58 30 20 20 20 20 20 |edastic.|.X0 |
|00003930| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 4e 75 6d | | Num|
|00003940| 62 65 72 20 6f 66 20 6d | 6f 76 69 6e 67 20 61 76 |ber of m|oving av|
|00003950| 65 72 61 67 65 20 74 65 | 72 6d 73 20 4d 41 20 66 |erage te|rms MA f|
|00003960| 6f 72 20 76 61 72 20 65 | 73 74 69 6d 61 74 65 2e |or var e|stimate.|
|00003970| 0a 58 6e 6f 6e 65 20 20 | 20 20 20 20 20 20 20 20 |.Xnone | |
|00003980| 20 20 20 20 20 43 6f 64 | 65 20 6e 6f 6e 65 20 6f | Cod|e none o|
|00003990| 72 20 50 61 72 7a 65 6e | 2c 20 6e 6f 6e 65 20 77 |r Parzen|, none w|
|000039a0| 68 65 6e 20 4d 41 3e 30 | 20 69 73 20 75 6e 77 69 |hen MA>0| is unwi|
|000039b0| 73 65 2e 0a 58 31 2e 30 | 30 30 30 30 30 65 2d 31 |se..X1.0|00000e-1|
|000039c0| 33 20 20 20 20 20 20 20 | 43 6f 6e 76 65 72 67 65 |3 |Converge|
|000039d0| 6e 63 65 20 74 6f 6c 65 | 72 61 6e 63 65 2c 20 74 |nce tole|rance, t|
|000039e0| 6f 6c 3d 31 2e 30 65 2d | 38 20 69 73 20 72 65 61 |ol=1.0e-|8 is rea|
|000039f0| 73 6f 6e 61 62 6c 65 2e | 0a 58 31 2e 30 30 30 30 |sonable.|.X1.0000|
|00003a00| 30 30 65 2d 31 30 20 20 | 20 20 20 20 20 49 6e 76 |00e-10 | Inv|
|00003a10| 65 72 73 69 6f 6e 20 74 | 6f 6c 65 72 61 6e 63 65 |ersion t|olerance|
|00003a20| 2c 20 65 70 73 3d 31 2e | 30 65 2d 31 33 20 69 73 |, eps=1.|0e-13 is|
|00003a30| 20 72 65 61 73 6f 6e 61 | 62 6c 65 0a 58 66 75 6c | reasona|ble.Xful|
|00003a40| 6c 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |l | |
|00003a50| 48 6f 77 20 6d 75 63 68 | 20 6f 75 74 70 75 74 3f |How much| output?|
|00003a60| 20 20 43 6f 64 65 20 6e | 6f 6e 65 2c 20 6d 69 6e | Code n|one, min|
|00003a70| 69 6d 61 6c 2c 20 6f 72 | 20 66 75 6c 6c 2e 0a 58 |imal, or| full..X|
|00003a80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003a90| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|00003aa0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003ab0| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|00003ac0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003ad0| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|00003ae0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003af0| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|00003b00| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003b10| 20 20 20 42 6c 61 6e 6b | 20 6c 69 6e 65 2e 0a 58 | Blank| line..X|
|00003b20| 20 2d 30 2e 30 34 38 36 | 36 0a 58 20 20 31 2e 30 | -0.0486|6.X 1.0|
|00003b30| 33 38 38 34 0a 58 20 2d | 30 2e 37 33 37 39 32 0a |3884.X -|0.73792.|
|00003b40| 58 20 2d 30 2e 35 31 33 | 36 32 20 0a 45 4e 44 5f |X -0.513|62 .END_|
|00003b50| 4f 46 5f 46 49 4c 45 0a | 69 66 20 74 65 73 74 20 |OF_FILE.|if test |
|00003b60| 31 30 39 37 20 2d 6e 65 | 20 60 77 63 20 2d 63 20 |1097 -ne| `wc -c |
|00003b70| 3c 27 63 68 31 65 67 31 | 2f 73 74 61 72 74 69 6e |<'ch1eg1|/startin|
|00003b80| 67 2e 64 61 74 27 60 3b | 20 74 68 65 6e 0a 20 20 |g.dat'`;| then. |
|00003b90| 20 20 65 63 68 6f 20 73 | 68 61 72 3a 20 5c 22 27 | echo s|har: \"'|
|00003ba0| 63 68 31 65 67 31 2f 73 | 74 61 72 74 69 6e 67 2e |ch1eg1/s|tarting.|
|00003bb0| 64 61 74 27 5c 22 20 75 | 6e 70 61 63 6b 65 64 20 |dat'\" u|npacked |
|00003bc0| 77 69 74 68 20 77 72 6f | 6e 67 20 73 69 7a 65 21 |with wro|ng size!|
|00003bd0| 0a 66 69 0a 23 20 65 6e | 64 20 6f 66 20 27 63 68 |.fi.# en|d of 'ch|
|00003be0| 31 65 67 31 2f 73 74 61 | 72 74 69 6e 67 2e 64 61 |1eg1/sta|rting.da|
|00003bf0| 74 27 0a 66 69 0a 69 66 | 20 74 65 73 74 20 21 20 |t'.fi.if| test ! |
|00003c00| 2d 64 20 27 63 68 35 65 | 67 31 27 20 3b 20 74 68 |-d 'ch5e|g1' ; th|
|00003c10| 65 6e 0a 20 20 20 20 65 | 63 68 6f 20 73 68 61 72 |en. e|cho shar|
|00003c20| 3a 20 43 72 65 61 74 69 | 6e 67 20 64 69 72 65 63 |: Creati|ng direc|
|00003c30| 74 6f 72 79 20 5c 22 27 | 63 68 35 65 67 31 27 5c |tory \"'|ch5eg1'\|
|00003c40| 22 0a 20 20 20 20 6d 6b | 64 69 72 20 27 63 68 35 |". mk|dir 'ch5|
|00003c50| 65 67 31 27 0a 66 69 0a | 69 66 20 74 65 73 74 20 |eg1'.fi.|if test |
|00003c60| 2d 66 20 27 63 68 35 65 | 67 31 2f 65 6c 65 63 74 |-f 'ch5e|g1/elect|
|00003c70| 62 2e 64 61 74 27 20 2d | 61 20 22 24 7b 31 7d 22 |b.dat' -|a "${1}"|
|00003c80| 20 21 3d 20 22 2d 63 22 | 20 3b 20 74 68 65 6e 20 | != "-c"| ; then |
|00003c90| 0a 20 20 65 63 68 6f 20 | 73 68 61 72 3a 20 57 69 |. echo |shar: Wi|
|00003ca0| 6c 6c 20 6e 6f 74 20 63 | 6c 6f 62 62 65 72 20 65 |ll not c|lobber e|
|00003cb0| 78 69 73 74 69 6e 67 20 | 66 69 6c 65 20 5c 22 27 |xisting |file \"'|
|00003cc0| 63 68 35 65 67 31 2f 65 | 6c 65 63 74 62 2e 64 61 |ch5eg1/e|lectb.da|
|00003cd0| 74 27 5c 22 0a 65 6c 73 | 65 0a 65 63 68 6f 20 73 |t'\".els|e.echo s|
|00003ce0| 68 61 72 3a 20 45 78 74 | 72 61 63 74 69 6e 67 20 |har: Ext|racting |
|00003cf0| 5c 22 27 63 68 35 65 67 | 31 2f 65 6c 65 63 74 62 |\"'ch5eg|1/electb|
|00003d00| 2e 64 61 74 27 5c 22 20 | 5c 28 36 31 32 20 63 68 |.dat'\" |\(612 ch|
|00003d10| 61 72 61 63 74 65 72 73 | 5c 29 0a 73 65 64 20 22 |aracters|\).sed "|
|00003d20| 73 2f 5e 58 2f 2f 22 20 | 3e 27 63 68 35 65 67 31 |s/^X//" |>'ch5eg1|
|00003d30| 2f 65 6c 65 63 74 62 2e | 64 61 74 27 20 3c 3c 27 |/electb.|dat' <<'|
|00003d40| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 27 0a 58 20 20 |END_OF_F|ILE'.X |
|00003d50| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003d60| 20 20 20 20 31 20 20 20 | 20 20 20 20 20 20 20 20 | 1 | |
|00003d70| 20 20 20 20 20 20 20 20 | 20 20 31 2e 30 36 20 20 | | 1.06 |
|00003d80| 20 20 20 32 2e 38 36 20 | 20 20 20 20 20 33 2e 39 | 2.86 | 3.9|
|00003d90| 30 0a 58 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |0.X | |
|00003da0| 20 20 20 20 20 20 20 20 | 20 32 20 20 20 20 20 20 | | 2 |
|00003db0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 31 | | 1|
|00003dc0| 2e 37 38 20 20 20 20 20 | 32 2e 38 36 20 20 20 20 |.78 |2.86 |
|00003dd0| 20 20 33 2e 39 30 0a 58 | 20 20 20 20 20 20 20 20 | 3.90.X| |
|00003de0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 33 20 | | 3 |
|00003df0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003e00| 20 20 20 20 31 2e 30 36 | 20 20 20 20 20 33 2e 39 | 1.06| 3.9|
|00003e10| 30 20 20 20 20 20 20 33 | 2e 39 30 0a 58 20 20 20 |0 3|.90.X |
|00003e20| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003e30| 20 20 20 34 20 20 20 20 | 20 20 20 20 20 20 20 20 | 4 | |
|00003e40| 20 20 20 20 20 20 20 20 | 20 31 2e 37 38 20 20 20 | | 1.78 |
|00003e50| 20 20 33 2e 39 30 20 20 | 20 20 20 20 33 2e 39 30 | 3.90 | 3.90|
|00003e60| 0a 58 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |.X | |
|00003e70| 20 20 20 20 20 20 20 20 | 35 20 20 20 20 20 20 20 | |5 |
|00003e80| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 31 2e | | 1.|
|00003e90| 33 37 20 20 20 20 20 33 | 2e 33 34 20 20 20 20 20 |37 3|.34 |
|00003ea0| 20 35 2e 30 36 0a 58 20 | 20 20 20 20 20 20 20 20 | 5.06.X | |
|00003eb0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 36 20 20 | | 6 |
|00003ec0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003ed0| 20 20 20 31 2e 30 36 20 | 20 20 20 20 32 2e 38 36 | 1.06 | 2.86|
|00003ee0| 20 20 20 20 20 20 36 2e | 35 36 0a 58 20 20 20 20 | 6.|56.X |
|00003ef0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003f00| 20 20 37 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | 7 | |
|00003f10| 20 20 20 20 20 20 20 20 | 31 2e 37 38 20 20 20 20 | |1.78 |
|00003f20| 20 32 2e 38 36 20 20 20 | 20 20 20 36 2e 35 36 0a | 2.86 | 6.56.|
|00003f30| 58 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X | |
|00003f40| 20 20 20 20 20 20 20 38 | 20 20 20 20 20 20 20 20 | 8| |
|00003f50| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 31 2e 30 | | 1.0|
|00003f60| 36 20 20 20 20 20 33 2e | 39 30 20 20 20 20 20 20 |6 3.|90 |
|00003f70| 36 2e 35 36 0a 58 20 20 | 20 20 20 20 20 20 20 20 |6.56.X | |
|00003f80| 20 20 20 20 20 20 20 20 | 20 20 20 20 39 20 20 20 | | 9 |
|00003f90| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00003fa0| 20 20 31 2e 37 38 20 20 | 20 20 20 33 2e 39 30 20 | 1.78 | 3.90 |
|00003fb0| 20 20 20 20 20 36 2e 35 | 36 0a 45 4e 44 5f 4f 46 | 6.5|6.END_OF|
|00003fc0| 5f 46 49 4c 45 0a 69 66 | 20 74 65 73 74 20 36 31 |_FILE.if| test 61|
|00003fd0| 32 20 2d 6e 65 20 60 77 | 63 20 2d 63 20 3c 27 63 |2 -ne `w|c -c <'c|
|00003fe0| 68 35 65 67 31 2f 65 6c | 65 63 74 62 2e 64 61 74 |h5eg1/el|ectb.dat|
|00003ff0| 27 60 3b 20 74 68 65 6e | 0a 20 20 20 20 65 63 68 |'`; then|. ech|
|00004000| 6f 20 73 68 61 72 3a 20 | 5c 22 27 63 68 35 65 67 |o shar: |\"'ch5eg|
|00004010| 31 2f 65 6c 65 63 74 62 | 2e 64 61 74 27 5c 22 20 |1/electb|.dat'\" |
|00004020| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|00004030| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|00004040| 6e 64 20 6f 66 20 27 63 | 68 35 65 67 31 2f 65 6c |nd of 'c|h5eg1/el|
|00004050| 65 63 74 62 2e 64 61 74 | 27 0a 66 69 0a 69 66 20 |ectb.dat|'.fi.if |
|00004060| 74 65 73 74 20 2d 66 20 | 27 63 68 35 65 67 31 2f |test -f |'ch5eg1/|
|00004070| 6d 6f 64 65 6c 2e 68 27 | 20 2d 61 20 22 24 7b 31 |model.h'| -a "${1|
|00004080| 7d 22 20 21 3d 20 22 2d | 63 22 20 3b 20 74 68 65 |}" != "-|c" ; the|
|00004090| 6e 20 0a 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |n . ech|o shar: |
|000040a0| 57 69 6c 6c 20 6e 6f 74 | 20 63 6c 6f 62 62 65 72 |Will not| clobber|
|000040b0| 20 65 78 69 73 74 69 6e | 67 20 66 69 6c 65 20 5c | existin|g file \|
|000040c0| 22 27 63 68 35 65 67 31 | 2f 6d 6f 64 65 6c 2e 68 |"'ch5eg1|/model.h|
|000040d0| 27 5c 22 0a 65 6c 73 65 | 0a 65 63 68 6f 20 73 68 |'\".else|.echo sh|
|000040e0| 61 72 3a 20 45 78 74 72 | 61 63 74 69 6e 67 20 5c |ar: Extr|acting \|
|000040f0| 22 27 63 68 35 65 67 31 | 2f 6d 6f 64 65 6c 2e 68 |"'ch5eg1|/model.h|
|00004100| 27 5c 22 20 5c 28 34 32 | 33 20 63 68 61 72 61 63 |'\" \(42|3 charac|
|00004110| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00004120| 2f 2f 22 20 3e 27 63 68 | 35 65 67 31 2f 6d 6f 64 |//" >'ch|5eg1/mod|
|00004130| 65 6c 2e 68 27 20 3c 3c | 27 45 4e 44 5f 4f 46 5f |el.h' <<|'END_OF_|
|00004140| 46 49 4c 45 27 0a 58 23 | 69 66 6e 64 65 66 20 5f |FILE'.X#|ifndef _|
|00004150| 5f 46 49 4c 45 5f 4d 4f | 44 45 4c 5f 48 5f 53 45 |_FILE_MO|DEL_H_SE|
|00004160| 45 4e 5f 5f 0a 58 23 70 | 72 61 67 6d 61 20 6f 6e |EN__.X#p|ragma on|
|00004170| 63 65 0a 58 23 64 65 66 | 69 6e 65 20 5f 5f 46 49 |ce.X#def|ine __FI|
|00004180| 4c 45 5f 4d 4f 44 45 4c | 5f 48 5f 53 45 45 4e 5f |LE_MODEL|_H_SEEN_|
|00004190| 5f 20 31 0a 58 0a 58 23 | 69 6e 63 6c 75 64 65 20 |_ 1.X.X#|include |
|000041a0| 22 73 74 61 74 75 73 2e | 68 22 0a 58 0a 58 65 78 |"status.|h".X.Xex|
|000041b0| 74 65 72 6e 20 73 74 61 | 74 75 73 20 73 3b 0a 58 |tern sta|tus s;.X|
|000041c0| 0a 58 63 6c 61 73 73 20 | 6d 6f 64 65 6c 0a 58 7b |.Xclass |model.X{|
|000041d0| 0a 58 0a 58 70 72 6f 74 | 65 63 74 65 64 3a 0a 58 |.X.Xprot|ected:.X|
|000041e0| 20 20 72 65 61 6c 6d 61 | 74 20 20 20 64 61 74 61 | realma|t data|
|000041f0| 3b 0a 58 0a 58 70 75 62 | 6c 69 63 3a 0a 58 20 20 |;.X.Xpub|lic:.X |
|00004200| 20 20 20 20 20 20 20 20 | 20 20 6d 6f 64 65 6c 28 | | model(|
|00004210| 29 3b 0a 58 20 20 20 20 | 20 20 20 20 20 20 20 20 |);.X | |
|00004220| 7e 6d 6f 64 65 6c 28 29 | 3b 0a 58 0a 58 20 20 69 |~model()|;.X.X i|
|00004230| 6e 74 20 20 20 20 20 20 | 20 69 6e 69 74 69 61 6c |nt | initial|
|00004240| 69 7a 65 28 29 3b 0a 58 | 20 20 69 6e 74 20 20 20 |ize();.X| int |
|00004250| 20 20 20 20 74 65 72 6d | 69 6e 61 74 65 28 29 3b | term|inate();|
|00004260| 0a 58 0a 58 20 20 72 65 | 61 6c 6d 61 74 20 20 20 |.X.X re|almat |
|00004270| 65 28 49 4e 54 45 47 45 | 52 20 74 29 3b 20 20 20 |e(INTEGE|R t); |
|00004280| 20 20 20 2f 2f 65 20 20 | 20 20 69 73 20 20 32 20 | //e | is 2 |
|00004290| 62 79 20 31 0a 58 20 20 | 72 65 61 6c 6d 61 74 20 |by 1.X |realmat |
|000042a0| 20 20 64 65 6c 65 28 49 | 4e 54 45 47 45 52 20 74 | dele(I|NTEGER t|
|000042b0| 29 3b 20 20 20 2f 2f 64 | 65 6c 65 20 69 73 20 20 |); //d|ele is |
|000042c0| 32 20 62 79 20 38 0a 58 | 20 20 72 65 61 6c 6d 61 |2 by 8.X| realma|
|000042d0| 74 20 20 20 5a 28 49 4e | 54 45 47 45 52 20 74 29 |t Z(IN|TEGER t)|
|000042e0| 3b 20 20 20 20 20 20 2f | 2f 5a 20 20 20 20 69 73 |; /|/Z is|
|000042f0| 20 31 37 20 62 79 20 31 | 0a 58 0a 58 7d 3b 0a 58 | 17 by 1|.X.X};.X|
|00004300| 0a 58 23 65 6e 64 69 66 | 0a 45 4e 44 5f 4f 46 5f |.X#endif|.END_OF_|
|00004310| 46 49 4c 45 0a 69 66 20 | 74 65 73 74 20 34 32 33 |FILE.if |test 423|
|00004320| 20 2d 6e 65 20 60 77 63 | 20 2d 63 20 3c 27 63 68 | -ne `wc| -c <'ch|
|00004330| 35 65 67 31 2f 6d 6f 64 | 65 6c 2e 68 27 60 3b 20 |5eg1/mod|el.h'`; |
|00004340| 74 68 65 6e 0a 20 20 20 | 20 65 63 68 6f 20 73 68 |then. | echo sh|
|00004350| 61 72 3a 20 5c 22 27 63 | 68 35 65 67 31 2f 6d 6f |ar: \"'c|h5eg1/mo|
|00004360| 64 65 6c 2e 68 27 5c 22 | 20 75 6e 70 61 63 6b 65 |del.h'\"| unpacke|
|00004370| 64 20 77 69 74 68 20 77 | 72 6f 6e 67 20 73 69 7a |d with w|rong siz|
|00004380| 65 21 0a 66 69 0a 23 20 | 65 6e 64 20 6f 66 20 27 |e!.fi.# |end of '|
|00004390| 63 68 35 65 67 31 2f 6d | 6f 64 65 6c 2e 68 27 0a |ch5eg1/m|odel.h'.|
|000043a0| 66 69 0a 69 66 20 74 65 | 73 74 20 2d 66 20 27 63 |fi.if te|st -f 'c|
|000043b0| 68 35 65 67 31 2f 73 74 | 61 72 74 69 6e 67 2e 64 |h5eg1/st|arting.d|
|000043c0| 61 74 27 20 2d 61 20 22 | 24 7b 31 7d 22 20 21 3d |at' -a "|${1}" !=|
|000043d0| 20 22 2d 63 22 20 3b 20 | 74 68 65 6e 20 0a 20 20 | "-c" ; |then . |
|000043e0| 65 63 68 6f 20 73 68 61 | 72 3a 20 57 69 6c 6c 20 |echo sha|r: Will |
|000043f0| 6e 6f 74 20 63 6c 6f 62 | 62 65 72 20 65 78 69 73 |not clob|ber exis|
|00004400| 74 69 6e 67 20 66 69 6c | 65 20 5c 22 27 63 68 35 |ting fil|e \"'ch5|
|00004410| 65 67 31 2f 73 74 61 72 | 74 69 6e 67 2e 64 61 74 |eg1/star|ting.dat|
|00004420| 27 5c 22 0a 65 6c 73 65 | 0a 65 63 68 6f 20 73 68 |'\".else|.echo sh|
|00004430| 61 72 3a 20 45 78 74 72 | 61 63 74 69 6e 67 20 5c |ar: Extr|acting \|
|00004440| 22 27 63 68 35 65 67 31 | 2f 73 74 61 72 74 69 6e |"'ch5eg1|/startin|
|00004450| 67 2e 64 61 74 27 5c 22 | 20 5c 28 31 31 31 33 20 |g.dat'\"| \(1113 |
|00004460| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|00004470| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 27 63 68 35 65 | "s/^X//|" >'ch5e|
|00004480| 67 31 2f 73 74 61 72 74 | 69 6e 67 2e 64 61 74 27 |g1/start|ing.dat'|
|00004490| 20 3c 3c 27 45 4e 44 5f | 4f 46 5f 46 49 4c 45 27 | <<'END_|OF_FILE'|
|000044a0| 0a 58 47 61 6c 6c 61 6e | 74 2c 20 22 4e 6f 6e 6c |.XGallan|t, "Nonl|
|000044b0| 69 6e 65 61 72 20 53 74 | 61 74 69 73 74 69 63 61 |inear St|atistica|
|000044c0| 6c 20 4d 6f 64 65 6c 73 | 2c 22 20 43 68 61 70 74 |l Models|," Chapt|
|000044d0| 65 72 20 35 2c 20 46 69 | 67 75 72 65 20 33 2c 20 |er 5, Fi|gure 3, |
|000044e0| 70 2e 20 33 32 34 2d 33 | 32 36 0a 58 53 55 52 20 |p. 324-3|26.XSUR |
|000044f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 57 | | W|
|00004500| 68 61 74 20 65 73 74 69 | 6d 61 74 69 6f 6e 20 6d |hat esti|mation m|
|00004510| 65 74 68 6f 64 3f 20 20 | 43 6f 64 65 20 53 55 52 |ethod? |Code SUR|
|00004520| 2c 20 54 53 4c 53 2c 20 | 6f 72 20 47 4d 4d 2e 0a |, TSLS, |or GMM..|
|00004530| 58 32 32 34 20 20 20 20 | 20 20 20 20 20 20 20 20 |X224 | |
|00004540| 20 20 20 20 4e 75 6d 62 | 65 72 20 6f 66 20 6f 62 | Numb|er of ob|
|00004550| 73 65 72 76 61 74 69 6f | 6e 73 2c 20 74 20 3d 20 |servatio|ns, t = |
|00004560| 31 2c 20 2e 2e 2e 2c 20 | 6e 2e 0a 58 32 20 20 20 |1, ..., |n..X2 |
|00004570| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 4e | | N|
|00004580| 75 6d 62 65 72 20 6f 66 | 20 65 71 75 61 74 69 6f |umber of| equatio|
|00004590| 6e 73 2c 20 69 2e 65 2e | 20 64 69 6d 65 6e 73 69 |ns, i.e.| dimensi|
|000045a0| 6f 6e 20 6f 66 20 65 2e | 0a 58 30 20 20 20 20 20 |on of e.|.X0 |
|000045b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 4e 75 6d | | Num|
|000045c0| 62 65 72 20 6f 66 20 69 | 6e 73 74 72 75 6d 65 6e |ber of i|nstrumen|
|000045d0| 74 73 2c 20 69 2e 65 2e | 20 64 69 6d 65 6e 73 69 |ts, i.e.| dimensi|
|000045e0| 6f 6e 20 6f 66 20 5a 2e | 0a 58 38 20 20 20 20 20 |on of Z.|.X8 |
|000045f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 4e 75 6d | | Num|
|00004600| 62 65 72 20 6f 66 20 70 | 61 72 61 6d 65 74 65 72 |ber of p|arameter|
|00004610| 73 2c 20 69 2e 65 2e 20 | 64 69 6d 65 6e 73 69 6f |s, i.e. |dimensio|
|00004620| 6e 20 6f 66 20 74 68 65 | 74 61 2e 0a 58 35 30 20 |n of the|ta..X50 |
|00004630| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00004640| 55 70 70 65 72 20 6c 69 | 6d 69 74 20 6f 6e 20 47 |Upper li|mit on G|
|00004650| 61 75 73 73 2d 4e 65 77 | 74 6f 6e 20 69 74 65 72 |auss-New|ton iter|
|00004660| 61 74 69 6f 6e 73 2e 0a | 58 31 20 20 20 20 20 20 |ations..|X1 |
|00004670| 20 20 20 20 20 20 20 20 | 20 20 20 20 4e 75 6d 62 | | Numb|
|00004680| 65 72 20 76 61 72 20 69 | 74 65 72 61 74 65 73 2c |er var i|terates,|
|00004690| 20 69 76 61 72 3d 30 20 | 6d 65 61 6e 73 20 6e 6f | ivar=0 |means no|
|000046a0| 6e 65 2e 0a 58 68 6f 6d | 6f 73 6b 65 64 61 73 74 |ne..Xhom|oskedast|
|000046b0| 69 63 20 20 20 20 20 20 | 43 6f 64 65 20 68 6f 6d |ic |Code hom|
|000046c0| 6f 73 6b 65 64 61 73 74 | 69 63 20 6f 72 20 68 65 |oskedast|ic or he|
|000046d0| 74 65 72 6f 73 6b 65 64 | 61 73 74 69 63 2e 0a 58 |terosked|astic..X|
|000046e0| 30 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |0 | |
|000046f0| 20 20 20 4e 75 6d 62 65 | 72 20 6f 66 20 6d 6f 76 | Numbe|r of mov|
|00004700| 69 6e 67 20 61 76 65 72 | 61 67 65 20 74 65 72 6d |ing aver|age term|
|00004710| 73 20 4d 41 20 66 6f 72 | 20 76 61 72 20 65 73 74 |s MA for| var est|
|00004720| 69 6d 61 74 65 2e 0a 58 | 6e 6f 6e 65 20 20 20 20 |imate..X|none |
|00004730| 20 20 20 20 20 20 20 20 | 20 20 20 43 6f 64 65 20 | | Code |
|00004740| 6e 6f 6e 65 20 6f 72 20 | 50 61 72 7a 65 6e 2c 20 |none or |Parzen, |
|00004750| 6e 6f 6e 65 20 77 68 65 | 6e 20 4d 41 3e 30 20 69 |none whe|n MA>0 i|
|00004760| 73 20 75 6e 77 69 73 65 | 2e 0a 58 31 2e 30 30 30 |s unwise|..X1.000|
|00004770| 30 30 30 65 2d 31 33 20 | 20 20 20 20 20 20 43 6f |000e-13 | Co|
|00004780| 6e 76 65 72 67 65 6e 63 | 65 20 74 6f 6c 65 72 61 |nvergenc|e tolera|
|00004790| 6e 63 65 2c 20 74 6f 6c | 3d 31 2e 30 65 2d 38 20 |nce, tol|=1.0e-8 |
|000047a0| 69 73 20 72 65 61 73 6f | 6e 61 62 6c 65 2e 0a 58 |is reaso|nable..X|
|000047b0| 31 2e 30 30 30 30 30 30 | 65 2d 31 30 20 20 20 20 |1.000000|e-10 |
|000047c0| 20 20 20 49 6e 76 65 72 | 73 69 6f 6e 20 74 6f 6c | Inver|sion tol|
|000047d0| 65 72 61 6e 63 65 2c 20 | 65 70 73 3d 31 2e 30 65 |erance, |eps=1.0e|
|000047e0| 2d 31 33 20 69 73 20 72 | 65 61 73 6f 6e 61 62 6c |-13 is r|easonabl|
|000047f0| 65 0a 58 6d 69 6e 69 6d | 61 6c 20 20 20 20 20 20 |e.Xminim|al |
|00004800| 20 20 20 20 20 20 48 6f | 77 20 6d 75 63 68 20 6f | Ho|w much o|
|00004810| 75 74 70 75 74 3f 20 20 | 43 6f 64 65 20 6e 6f 6e |utput? |Code non|
|00004820| 65 2c 20 6d 69 6e 69 6d | 61 6c 2c 20 6f 72 20 66 |e, minim|al, or f|
|00004830| 75 6c 6c 2e 0a 58 20 20 | 20 20 20 20 20 20 20 20 |ull..X | |
|00004840| 20 20 20 20 20 20 20 20 | 20 42 6c 61 6e 6b 20 6c | | Blank l|
|00004850| 69 6e 65 2e 0a 58 20 20 | 20 20 20 20 20 20 20 20 |ine..X | |
|00004860| 20 20 20 20 20 20 20 20 | 20 42 6c 61 6e 6b 20 6c | | Blank l|
|00004870| 69 6e 65 2e 0a 58 20 20 | 20 20 20 20 20 20 20 20 |ine..X | |
|00004880| 20 20 20 20 20 20 20 20 | 20 42 6c 61 6e 6b 20 6c | | Blank l|
|00004890| 69 6e 65 2e 0a 58 20 20 | 20 20 20 20 20 20 20 20 |ine..X | |
|000048a0| 20 20 20 20 20 20 20 20 | 20 42 6c 61 6e 6b 20 6c | | Blank l|
|000048b0| 69 6e 65 2e 0a 58 20 20 | 20 20 20 20 20 20 20 20 |ine..X | |
|000048c0| 20 20 20 20 20 20 20 20 | 20 42 6c 61 6e 6b 20 6c | | Blank l|
|000048d0| 69 6e 65 2e 0a 58 20 2d | 32 2e 39 30 0a 58 20 2d |ine..X -|2.90.X -|
|000048e0| 31 2e 33 30 0a 58 20 20 | 30 2e 38 32 0a 58 20 20 |1.30.X |0.82.X |
|000048f0| 30 2e 33 36 0a 58 20 2d | 31 2e 35 30 0a 58 20 2d |0.36.X -|1.50.X -|
|00004900| 31 2e 30 30 0a 58 20 2d | 30 2e 30 33 0a 58 20 2d |1.00.X -|0.03.X -|
|00004910| 30 2e 34 37 0a 45 4e 44 | 5f 4f 46 5f 46 49 4c 45 |0.47.END|_OF_FILE|
|00004920| 0a 69 66 20 74 65 73 74 | 20 31 31 31 33 20 2d 6e |.if test| 1113 -n|
|00004930| 65 20 60 77 63 20 2d 63 | 20 3c 27 63 68 35 65 67 |e `wc -c| <'ch5eg|
|00004940| 31 2f 73 74 61 72 74 69 | 6e 67 2e 64 61 74 27 60 |1/starti|ng.dat'`|
|00004950| 3b 20 74 68 65 6e 0a 20 | 20 20 20 65 63 68 6f 20 |; then. | echo |
|00004960| 73 68 61 72 3a 20 5c 22 | 27 63 68 35 65 67 31 2f |shar: \"|'ch5eg1/|
|00004970| 73 74 61 72 74 69 6e 67 | 2e 64 61 74 27 5c 22 20 |starting|.dat'\" |
|00004980| 75 6e 70 61 63 6b 65 64 | 20 77 69 74 68 20 77 72 |unpacked| with wr|
|00004990| 6f 6e 67 20 73 69 7a 65 | 21 0a 66 69 0a 23 20 65 |ong size|!.fi.# e|
|000049a0| 6e 64 20 6f 66 20 27 63 | 68 35 65 67 31 2f 73 74 |nd of 'c|h5eg1/st|
|000049b0| 61 72 74 69 6e 67 2e 64 | 61 74 27 0a 66 69 0a 69 |arting.d|at'.fi.i|
|000049c0| 66 20 74 65 73 74 20 21 | 20 2d 64 20 27 63 68 36 |f test !| -d 'ch6|
|000049d0| 65 67 31 27 20 3b 20 74 | 68 65 6e 0a 20 20 20 20 |eg1' ; t|hen. |
|000049e0| 65 63 68 6f 20 73 68 61 | 72 3a 20 43 72 65 61 74 |echo sha|r: Creat|
|000049f0| 69 6e 67 20 64 69 72 65 | 63 74 6f 72 79 20 5c 22 |ing dire|ctory \"|
|00004a00| 27 63 68 36 65 67 31 27 | 5c 22 0a 20 20 20 20 6d |'ch6eg1'|\". m|
|00004a10| 6b 64 69 72 20 27 63 68 | 36 65 67 31 27 0a 66 69 |kdir 'ch|6eg1'.fi|
|00004a20| 0a 69 66 20 74 65 73 74 | 20 2d 66 20 27 63 68 36 |.if test| -f 'ch6|
|00004a30| 65 67 31 2f 52 45 41 44 | 4d 45 27 20 2d 61 20 22 |eg1/READ|ME' -a "|
|00004a40| 24 7b 31 7d 22 20 21 3d | 20 22 2d 63 22 20 3b 20 |${1}" !=| "-c" ; |
|00004a50| 74 68 65 6e 20 0a 20 20 | 65 63 68 6f 20 73 68 61 |then . |echo sha|
|00004a60| 72 3a 20 57 69 6c 6c 20 | 6e 6f 74 20 63 6c 6f 62 |r: Will |not clob|
|00004a70| 62 65 72 20 65 78 69 73 | 74 69 6e 67 20 66 69 6c |ber exis|ting fil|
|00004a80| 65 20 5c 22 27 63 68 36 | 65 67 31 2f 52 45 41 44 |e \"'ch6|eg1/READ|
|00004a90| 4d 45 27 5c 22 0a 65 6c | 73 65 0a 65 63 68 6f 20 |ME'\".el|se.echo |
|00004aa0| 73 68 61 72 3a 20 45 78 | 74 72 61 63 74 69 6e 67 |shar: Ex|tracting|
|00004ab0| 20 5c 22 27 63 68 36 65 | 67 31 2f 52 45 41 44 4d | \"'ch6e|g1/READM|
|00004ac0| 45 27 5c 22 20 5c 28 39 | 32 20 63 68 61 72 61 63 |E'\" \(9|2 charac|
|00004ad0| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00004ae0| 2f 2f 22 20 3e 27 63 68 | 36 65 67 31 2f 52 45 41 |//" >'ch|6eg1/REA|
|00004af0| 44 4d 45 27 20 3c 3c 27 | 45 4e 44 5f 4f 46 5f 46 |DME' <<'|END_OF_F|
|00004b00| 49 4c 45 27 0a 58 0a 58 | 44 6f 20 74 68 69 73 20 |ILE'.X.X|Do this |
|00004b10| 74 6f 20 67 65 74 20 74 | 68 65 20 6d 69 73 73 69 |to get t|he missi|
|00004b20| 6e 67 20 66 69 6c 65 73 | 3a 0a 58 0a 58 20 20 63 |ng files|:.X.X c|
|00004b30| 70 20 2e 2e 2f 63 68 35 | 65 67 31 2f 65 6c 2a 2e |p ../ch5|eg1/el*.|
|00004b40| 64 61 74 20 2e 2f 0a 58 | 20 20 63 70 20 2e 2e 2f |dat ./.X| cp ../|
|00004b50| 63 68 35 65 67 31 2f 6d | 6f 64 65 6c 2e 2a 20 2e |ch5eg1/m|odel.* .|
|00004b60| 2f 20 0a 58 20 20 0a 45 | 4e 44 5f 4f 46 5f 46 49 |/ .X .E|ND_OF_FI|
|00004b70| 4c 45 0a 69 66 20 74 65 | 73 74 20 39 32 20 2d 6e |LE.if te|st 92 -n|
|00004b80| 65 20 60 77 63 20 2d 63 | 20 3c 27 63 68 36 65 67 |e `wc -c| <'ch6eg|
|00004b90| 31 2f 52 45 41 44 4d 45 | 27 60 3b 20 74 68 65 6e |1/README|'`; then|
|00004ba0| 0a 20 20 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |. ech|o shar: |
|00004bb0| 5c 22 27 63 68 36 65 67 | 31 2f 52 45 41 44 4d 45 |\"'ch6eg|1/README|
|00004bc0| 27 5c 22 20 75 6e 70 61 | 63 6b 65 64 20 77 69 74 |'\" unpa|cked wit|
|00004bd0| 68 20 77 72 6f 6e 67 20 | 73 69 7a 65 21 0a 66 69 |h wrong |size!.fi|
|00004be0| 0a 23 20 65 6e 64 20 6f | 66 20 27 63 68 36 65 67 |.# end o|f 'ch6eg|
|00004bf0| 31 2f 52 45 41 44 4d 45 | 27 0a 66 69 0a 69 66 20 |1/README|'.fi.if |
|00004c00| 74 65 73 74 20 2d 66 20 | 27 63 68 36 65 67 31 2f |test -f |'ch6eg1/|
|00004c10| 73 74 61 72 74 69 6e 67 | 2e 64 61 74 27 20 2d 61 |starting|.dat' -a|
|00004c20| 20 22 24 7b 31 7d 22 20 | 21 3d 20 22 2d 63 22 20 | "${1}" |!= "-c" |
|00004c30| 3b 20 74 68 65 6e 20 0a | 20 20 65 63 68 6f 20 73 |; then .| echo s|
|00004c40| 68 61 72 3a 20 57 69 6c | 6c 20 6e 6f 74 20 63 6c |har: Wil|l not cl|
|00004c50| 6f 62 62 65 72 20 65 78 | 69 73 74 69 6e 67 20 66 |obber ex|isting f|
|00004c60| 69 6c 65 20 5c 22 27 63 | 68 36 65 67 31 2f 73 74 |ile \"'c|h6eg1/st|
|00004c70| 61 72 74 69 6e 67 2e 64 | 61 74 27 5c 22 0a 65 6c |arting.d|at'\".el|
|00004c80| 73 65 0a 65 63 68 6f 20 | 73 68 61 72 3a 20 45 78 |se.echo |shar: Ex|
|00004c90| 74 72 61 63 74 69 6e 67 | 20 5c 22 27 63 68 36 65 |tracting| \"'ch6e|
|00004ca0| 67 31 2f 73 74 61 72 74 | 69 6e 67 2e 64 61 74 27 |g1/start|ing.dat'|
|00004cb0| 5c 22 20 5c 28 31 31 31 | 33 20 63 68 61 72 61 63 |\" \(111|3 charac|
|00004cc0| 74 65 72 73 5c 29 0a 73 | 65 64 20 22 73 2f 5e 58 |ters\).s|ed "s/^X|
|00004cd0| 2f 2f 22 20 3e 27 63 68 | 36 65 67 31 2f 73 74 61 |//" >'ch|6eg1/sta|
|00004ce0| 72 74 69 6e 67 2e 64 61 | 74 27 20 3c 3c 27 45 4e |rting.da|t' <<'EN|
|00004cf0| 44 5f 4f 46 5f 46 49 4c | 45 27 0a 58 47 61 6c 6c |D_OF_FIL|E'.XGall|
|00004d00| 61 6e 74 2c 20 22 4e 6f | 6e 6c 69 6e 65 61 72 20 |ant, "No|nlinear |
|00004d10| 53 74 61 74 69 73 74 69 | 63 61 6c 20 4d 6f 64 65 |Statisti|cal Mode|
|00004d20| 6c 73 2c 22 20 43 68 61 | 70 74 65 72 20 36 2c 20 |ls," Cha|pter 6, |
|00004d30| 46 69 67 75 72 65 20 31 | 2c 20 70 2e 20 34 33 36 |Figure 1|, p. 436|
|00004d40| 2e 0a 58 54 53 4c 53 20 | 20 20 20 20 20 20 20 20 |..XTSLS | |
|00004d50| 20 20 20 20 20 20 57 68 | 61 74 20 65 73 74 69 6d | Wh|at estim|
|00004d60| 61 74 69 6f 6e 20 6d 65 | 74 68 6f 64 3f 20 20 43 |ation me|thod? C|
|00004d70| 6f 64 65 20 53 55 52 2c | 20 54 53 4c 53 2c 20 6f |ode SUR,| TSLS, o|
|00004d80| 72 20 47 4d 4d 2e 0a 58 | 32 32 34 20 20 20 20 20 |r GMM..X|224 |
|00004d90| 20 20 20 20 20 20 20 20 | 20 20 20 4e 75 6d 62 65 | | Numbe|
|00004da0| 72 20 6f 66 20 6f 62 73 | 65 72 76 61 74 69 6f 6e |r of obs|ervation|
|00004db0| 73 2c 20 74 20 3d 20 31 | 2c 20 2e 2e 2e 2c 20 6e |s, t = 1|, ..., n|
|00004dc0| 2e 0a 58 32 20 20 20 20 | 20 20 20 20 20 20 20 20 |..X2 | |
|00004dd0| 20 20 20 20 20 20 4e 75 | 6d 62 65 72 20 6f 66 20 | Nu|mber of |
|00004de0| 65 71 75 61 74 69 6f 6e | 73 2c 20 69 2e 65 2e 20 |equation|s, i.e. |
|00004df0| 64 69 6d 65 6e 73 69 6f | 6e 20 6f 66 20 65 2e 0a |dimensio|n of e..|
|00004e00| 58 31 37 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X17 | |
|00004e10| 20 20 20 20 4e 75 6d 62 | 65 72 20 6f 66 20 69 6e | Numb|er of in|
|00004e20| 73 74 72 75 6d 65 6e 74 | 73 2c 20 69 2e 65 2e 20 |strument|s, i.e. |
|00004e30| 64 69 6d 65 6e 73 69 6f | 6e 20 6f 66 20 5a 2e 0a |dimensio|n of Z..|
|00004e40| 58 38 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X8 | |
|00004e50| 20 20 20 20 4e 75 6d 62 | 65 72 20 6f 66 20 70 61 | Numb|er of pa|
|00004e60| 72 61 6d 65 74 65 72 73 | 2c 20 69 2e 65 2e 20 64 |rameters|, i.e. d|
|00004e70| 69 6d 65 6e 73 69 6f 6e | 20 6f 66 20 74 68 65 74 |imension| of thet|
|00004e80| 61 2e 0a 58 32 30 20 20 | 20 20 20 20 20 20 20 20 |a..X20 | |
|00004e90| 20 20 20 20 20 20 20 55 | 70 70 65 72 20 6c 69 6d | U|pper lim|
|00004ea0| 69 74 20 6f 6e 20 47 61 | 75 73 73 2d 4e 65 77 74 |it on Ga|uss-Newt|
|00004eb0| 6f 6e 20 69 74 65 72 61 | 74 69 6f 6e 73 2e 0a 58 |on itera|tions..X|
|00004ec0| 31 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |1 | |
|00004ed0| 20 20 20 4e 75 6d 62 65 | 72 20 76 61 72 20 69 74 | Numbe|r var it|
|00004ee0| 65 72 61 74 65 73 2c 20 | 69 76 61 72 3d 30 20 6d |erates, |ivar=0 m|
|00004ef0| 65 61 6e 73 20 6e 6f 6e | 65 2e 0a 58 68 6f 6d 6f |eans non|e..Xhomo|
|00004f00| 73 6b 65 64 61 73 74 69 | 63 20 20 20 20 20 20 43 |skedasti|c C|
|00004f10| 6f 64 65 20 68 6f 6d 6f | 73 6b 65 64 61 73 74 69 |ode homo|skedasti|
|00004f20| 63 20 6f 72 20 68 65 74 | 65 72 6f 73 6b 65 64 61 |c or het|eroskeda|
|00004f30| 73 74 69 63 2e 0a 58 30 | 20 20 20 20 20 20 20 20 |stic..X0| |
|00004f40| 20 20 20 20 20 20 20 20 | 20 20 4e 75 6d 62 65 72 | | Number|
|00004f50| 20 6f 66 20 6d 6f 76 69 | 6e 67 20 61 76 65 72 61 | of movi|ng avera|
|00004f60| 67 65 20 74 65 72 6d 73 | 20 4d 41 20 66 6f 72 20 |ge terms| MA for |
|00004f70| 76 61 72 20 65 73 74 69 | 6d 61 74 65 2e 0a 58 6e |var esti|mate..Xn|
|00004f80| 6f 6e 65 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |one | |
|00004f90| 20 20 43 6f 64 65 20 6e | 6f 6e 65 20 6f 72 20 50 | Code n|one or P|
|00004fa0| 61 72 7a 65 6e 2c 20 6e | 6f 6e 65 20 77 68 65 6e |arzen, n|one when|
|00004fb0| 20 4d 41 3e 30 20 69 73 | 20 75 6e 77 69 73 65 2e | MA>0 is| unwise.|
|00004fc0| 0a 58 31 2e 30 30 30 30 | 30 30 65 2d 30 38 20 20 |.X1.0000|00e-08 |
|00004fd0| 20 20 20 20 20 43 6f 6e | 76 65 72 67 65 6e 63 65 | Con|vergence|
|00004fe0| 20 74 6f 6c 65 72 61 6e | 63 65 2c 20 74 6f 6c 3d | toleran|ce, tol=|
|00004ff0| 31 2e 30 65 2d 38 20 69 | 73 20 72 65 61 73 6f 6e |1.0e-8 i|s reason|
|00005000| 61 62 6c 65 2e 0a 58 31 | 2e 30 30 30 30 30 30 65 |able..X1|.000000e|
|00005010| 2d 31 30 20 20 20 20 20 | 20 20 49 6e 76 65 72 73 |-10 | Invers|
|00005020| 69 6f 6e 20 74 6f 6c 65 | 72 61 6e 63 65 2c 20 65 |ion tole|rance, e|
|00005030| 70 73 3d 31 2e 30 65 2d | 31 33 20 69 73 20 72 65 |ps=1.0e-|13 is re|
|00005040| 61 73 6f 6e 61 62 6c 65 | 0a 58 66 75 6c 6c 20 20 |asonable|.Xfull |
|00005050| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 48 6f 77 | | How|
|00005060| 20 6d 75 63 68 20 6f 75 | 74 70 75 74 3f 20 20 43 | much ou|tput? C|
|00005070| 6f 64 65 20 6e 6f 6e 65 | 2c 20 6d 69 6e 69 6d 61 |ode none|, minima|
|00005080| 6c 2c 20 6f 72 20 66 75 | 6c 6c 2e 0a 58 20 20 20 |l, or fu|ll..X |
|00005090| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000050a0| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 20 20 |Blank li|ne..X |
|000050b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000050c0| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 20 20 |Blank li|ne..X |
|000050d0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000050e0| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 20 20 |Blank li|ne..X |
|000050f0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005100| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 20 20 |Blank li|ne..X |
|00005110| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005120| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 2d 32 |Blank li|ne..X -2|
|00005130| 2e 39 38 0a 58 20 2d 31 | 2e 31 36 0a 58 20 20 30 |.98.X -1|.16.X 0|
|00005140| 2e 37 38 37 0a 58 20 20 | 30 2e 33 35 33 0a 58 20 |.787.X |0.353.X |
|00005150| 2d 31 2e 35 31 0a 58 20 | 2d 31 2e 30 0a 58 20 20 |-1.51.X |-1.0.X |
|00005160| 30 2e 30 35 34 0a 58 20 | 2d 30 2e 34 37 34 0a 45 |0.054.X |-0.474.E|
|00005170| 4e 44 5f 4f 46 5f 46 49 | 4c 45 0a 69 66 20 74 65 |ND_OF_FI|LE.if te|
|00005180| 73 74 20 31 31 31 33 20 | 2d 6e 65 20 60 77 63 20 |st 1113 |-ne `wc |
|00005190| 2d 63 20 3c 27 63 68 36 | 65 67 31 2f 73 74 61 72 |-c <'ch6|eg1/star|
|000051a0| 74 69 6e 67 2e 64 61 74 | 27 60 3b 20 74 68 65 6e |ting.dat|'`; then|
|000051b0| 0a 20 20 20 20 65 63 68 | 6f 20 73 68 61 72 3a 20 |. ech|o shar: |
|000051c0| 5c 22 27 63 68 36 65 67 | 31 2f 73 74 61 72 74 69 |\"'ch6eg|1/starti|
|000051d0| 6e 67 2e 64 61 74 27 5c | 22 20 75 6e 70 61 63 6b |ng.dat'\|" unpack|
|000051e0| 65 64 20 77 69 74 68 20 | 77 72 6f 6e 67 20 73 69 |ed with |wrong si|
|000051f0| 7a 65 21 0a 66 69 0a 23 | 20 65 6e 64 20 6f 66 20 |ze!.fi.#| end of |
|00005200| 27 63 68 36 65 67 31 2f | 73 74 61 72 74 69 6e 67 |'ch6eg1/|starting|
|00005210| 2e 64 61 74 27 0a 66 69 | 0a 69 66 20 74 65 73 74 |.dat'.fi|.if test|
|00005220| 20 21 20 2d 64 20 27 63 | 68 36 65 67 32 27 20 3b | ! -d 'c|h6eg2' ;|
|00005230| 20 74 68 65 6e 0a 20 20 | 20 20 65 63 68 6f 20 73 | then. | echo s|
|00005240| 68 61 72 3a 20 43 72 65 | 61 74 69 6e 67 20 64 69 |har: Cre|ating di|
|00005250| 72 65 63 74 6f 72 79 20 | 5c 22 27 63 68 36 65 67 |rectory |\"'ch6eg|
|00005260| 32 27 5c 22 0a 20 20 20 | 20 6d 6b 64 69 72 20 27 |2'\". | mkdir '|
|00005270| 63 68 36 65 67 32 27 0a | 66 69 0a 69 66 20 74 65 |ch6eg2'.|fi.if te|
|00005280| 73 74 20 2d 66 20 27 63 | 68 36 65 67 32 2f 65 6e |st -f 'c|h6eg2/en|
|00005290| 64 69 6e 67 2e 62 61 6b | 27 20 2d 61 20 22 24 7b |ding.bak|' -a "${|
|000052a0| 31 7d 22 20 21 3d 20 22 | 2d 63 22 20 3b 20 74 68 |1}" != "|-c" ; th|
|000052b0| 65 6e 20 0a 20 20 65 63 | 68 6f 20 73 68 61 72 3a |en . ec|ho shar:|
|000052c0| 20 57 69 6c 6c 20 6e 6f | 74 20 63 6c 6f 62 62 65 | Will no|t clobbe|
|000052d0| 72 20 65 78 69 73 74 69 | 6e 67 20 66 69 6c 65 20 |r existi|ng file |
|000052e0| 5c 22 27 63 68 36 65 67 | 32 2f 65 6e 64 69 6e 67 |\"'ch6eg|2/ending|
|000052f0| 2e 62 61 6b 27 5c 22 0a | 65 6c 73 65 0a 65 63 68 |.bak'\".|else.ech|
|00005300| 6f 20 73 68 61 72 3a 20 | 45 78 74 72 61 63 74 69 |o shar: |Extracti|
|00005310| 6e 67 20 5c 22 27 63 68 | 36 65 67 32 2f 65 6e 64 |ng \"'ch|6eg2/end|
|00005320| 69 6e 67 2e 62 61 6b 27 | 5c 22 20 5c 28 32 31 30 |ing.bak'|\" \(210|
|00005330| 31 20 63 68 61 72 61 63 | 74 65 72 73 5c 29 0a 73 |1 charac|ters\).s|
|00005340| 65 64 20 22 73 2f 5e 58 | 2f 2f 22 20 3e 27 63 68 |ed "s/^X|//" >'ch|
|00005350| 36 65 67 32 2f 65 6e 64 | 69 6e 67 2e 62 61 6b 27 |6eg2/end|ing.bak'|
|00005360| 20 3c 3c 27 45 4e 44 5f | 4f 46 5f 46 49 4c 45 27 | <<'END_|OF_FILE'|
|00005370| 0a 58 47 61 6c 6c 61 6e | 74 2c 20 22 4e 6f 6e 6c |.XGallan|t, "Nonl|
|00005380| 69 6e 65 61 72 20 53 74 | 61 74 69 73 74 69 63 61 |inear St|atistica|
|00005390| 6c 20 4d 6f 64 65 6c 73 | 2c 22 20 43 68 61 70 74 |l Models|," Chapt|
|000053a0| 65 72 20 36 2c 20 46 69 | 67 75 72 65 20 32 2c 20 |er 6, Fi|gure 2, |
|000053b0| 70 2e 20 34 34 39 2d 34 | 35 30 2e 0a 58 47 4d 4d |p. 449-4|50..XGMM|
|000053c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000053d0| 57 68 61 74 20 65 73 74 | 69 6d 61 74 69 6f 6e 20 |What est|imation |
|000053e0| 6d 65 74 68 6f 64 3f 20 | 20 43 6f 64 65 20 53 55 |method? | Code SU|
|000053f0| 52 2c 20 54 53 4c 53 2c | 20 6f 72 20 47 4d 4d 2e |R, TSLS,| or GMM.|
|00005400| 0a 58 32 33 38 20 20 20 | 20 20 20 20 20 20 20 20 |.X238 | |
|00005410| 20 20 20 20 20 4e 75 6d | 62 65 72 20 6f 66 20 6f | Num|ber of o|
|00005420| 62 73 65 72 76 61 74 69 | 6f 6e 73 2c 20 74 20 3d |bservati|ons, t =|
|00005430| 20 31 2c 20 2e 2e 2e 2c | 20 6e 2e 0a 58 31 20 20 | 1, ...,| n..X1 |
|00005440| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005450| 4e 75 6d 62 65 72 20 6f | 66 20 65 71 75 61 74 69 |Number o|f equati|
|00005460| 6f 6e 73 2c 20 69 2e 65 | 2e 20 64 69 6d 65 6e 73 |ons, i.e|. dimens|
|00005470| 69 6f 6e 20 6f 66 20 65 | 2e 0a 58 33 20 20 20 20 |ion of e|..X3 |
|00005480| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 4e 75 | | Nu|
|00005490| 6d 62 65 72 20 6f 66 20 | 69 6e 73 74 72 75 6d 65 |mber of |instrume|
|000054a0| 6e 74 73 2c 20 69 2e 65 | 2e 20 64 69 6d 65 6e 73 |nts, i.e|. dimens|
|000054b0| 69 6f 6e 20 6f 66 20 5a | 2e 0a 58 32 20 20 20 20 |ion of Z|..X2 |
|000054c0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 4e 75 | | Nu|
|000054d0| 6d 62 65 72 20 6f 66 20 | 70 61 72 61 6d 65 74 65 |mber of |paramete|
|000054e0| 72 73 2c 20 69 2e 65 2e | 20 64 69 6d 65 6e 73 69 |rs, i.e.| dimensi|
|000054f0| 6f 6e 20 6f 66 20 74 68 | 65 74 61 2e 0a 58 32 30 |on of th|eta..X20|
|00005500| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005510| 20 55 70 70 65 72 20 6c | 69 6d 69 74 20 6f 6e 20 | Upper l|imit on |
|00005520| 47 61 75 73 73 2d 4e 65 | 77 74 6f 6e 20 69 74 65 |Gauss-Ne|wton ite|
|00005530| 72 61 74 69 6f 6e 73 2e | 0a 58 31 20 20 20 20 20 |rations.|.X1 |
|00005540| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 4e 75 6d | | Num|
|00005550| 62 65 72 20 76 61 72 20 | 69 74 65 72 61 74 65 73 |ber var |iterates|
|00005560| 2c 20 69 76 61 72 3d 30 | 20 6d 65 61 6e 73 20 6e |, ivar=0| means n|
|00005570| 6f 6e 65 2e 0a 58 68 6f | 6d 6f 73 6b 65 64 61 73 |one..Xho|moskedas|
|00005580| 74 69 63 20 20 20 20 20 | 20 43 6f 64 65 20 68 6f |tic | Code ho|
|00005590| 6d 6f 73 6b 65 64 61 73 | 74 69 63 20 6f 72 20 68 |moskedas|tic or h|
|000055a0| 65 74 65 72 6f 73 6b 65 | 64 61 73 74 69 63 2e 0a |eteroske|dastic..|
|000055b0| 58 30 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 |X0 | |
|000055c0| 20 20 20 20 4e 75 6d 62 | 65 72 20 6f 66 20 6d 6f | Numb|er of mo|
|000055d0| 76 69 6e 67 20 61 76 65 | 72 61 67 65 20 74 65 72 |ving ave|rage ter|
|000055e0| 6d 73 20 4d 41 20 66 6f | 72 20 76 61 72 20 65 73 |ms MA fo|r var es|
|000055f0| 74 69 6d 61 74 65 2e 0a | 58 50 61 72 7a 65 6e 20 |timate..|XParzen |
|00005600| 20 20 20 20 20 20 20 20 | 20 20 20 20 43 6f 64 65 | | Code|
|00005610| 20 6e 6f 6e 65 20 6f 72 | 20 50 61 72 7a 65 6e 2c | none or| Parzen,|
|00005620| 20 6e 6f 6e 65 20 77 68 | 65 6e 20 4d 41 3e 30 20 | none wh|en MA>0 |
|00005630| 69 73 20 75 6e 77 69 73 | 65 2e 0a 58 31 2e 30 30 |is unwis|e..X1.00|
|00005640| 30 30 30 30 65 2d 30 35 | 20 20 20 20 20 20 20 43 |0000e-05| C|
|00005650| 6f 6e 76 65 72 67 65 6e | 63 65 20 74 6f 6c 65 72 |onvergen|ce toler|
|00005660| 61 6e 63 65 2c 20 74 6f | 6c 3d 31 2e 30 65 2d 38 |ance, to|l=1.0e-8|
|00005670| 20 69 73 20 72 65 61 73 | 6f 6e 61 62 6c 65 2e 0a | is reas|onable..|
|00005680| 58 31 2e 30 30 30 30 30 | 30 65 2d 31 33 20 20 20 |X1.00000|0e-13 |
|00005690| 20 20 20 20 49 6e 76 65 | 72 73 69 6f 6e 20 74 6f | Inve|rsion to|
|000056a0| 6c 65 72 61 6e 63 65 2c | 20 65 70 73 3d 31 2e 30 |lerance,| eps=1.0|
|000056b0| 65 2d 31 33 20 69 73 20 | 72 65 61 73 6f 6e 61 62 |e-13 is |reasonab|
|000056c0| 6c 65 0a 58 66 75 6c 6c | 20 20 20 20 20 20 20 20 |le.Xfull| |
|000056d0| 20 20 20 20 20 20 20 48 | 6f 77 20 6d 75 63 68 20 | H|ow much |
|000056e0| 6f 75 74 70 75 74 3f 20 | 20 43 6f 64 65 20 6e 6f |output? | Code no|
|000056f0| 6e 65 2c 20 6d 69 6e 69 | 6d 61 6c 2c 20 6f 72 20 |ne, mini|mal, or |
|00005700| 66 75 6c 6c 2e 0a 58 32 | 20 20 20 20 20 20 20 20 |full..X2| |
|00005710| 20 20 20 20 20 20 20 20 | 20 20 43 6f 6d 70 75 74 | | Comput|
|00005720| 65 64 20 72 61 6e 6b 20 | 6f 66 20 56 2e 0a 58 75 |ed rank |of V..Xu|
|00005730| 6e 63 6f 72 72 65 63 74 | 65 64 20 20 20 20 20 20 |ncorrect|ed |
|00005740| 20 20 44 69 76 69 73 6f | 72 20 6f 66 20 76 61 72 | Diviso|r of var|
|00005750| 2c 20 63 6f 72 72 65 63 | 74 65 64 20 6f 72 20 75 |, correc|ted or u|
|00005760| 6e 63 6f 72 72 65 63 74 | 65 64 2e 0a 58 20 20 20 |ncorrect|ed..X |
|00005770| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|00005780| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 20 20 |Blank li|ne..X |
|00005790| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000057a0| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 20 20 |Blank li|ne..X |
|000057b0| 20 20 20 20 20 20 20 20 | 20 20 20 20 20 20 20 20 | | |
|000057c0| 42 6c 61 6e 6b 20 6c 69 | 6e 65 2e 0a 58 20 20 2d |Blank li|ne..X -|
|000057d0| 31 2e 30 33 33 35 32 31 | 37 33 39 32 30 33 35 33 |1.033521|73920353|
|000057e0| 39 30 65 2b 30 30 20 20 | 20 20 20 20 20 20 20 20 |90e+00 | |
|000057f0| 2d 31 2e 30 33 33 35 32 | 31 37 34 20 20 74 68 65 |-1.03352|174 the|
|00005800| 74 61 28 31 29 0a 58 20 | 20 20 39 2e 39 38 32 35 |ta(1).X | 9.9825|
|00005810| 35 37 32 33 36 34 32 34 | 30 32 30 30 65 2d 30 31 |57236424|0200e-01|
|00005820| 20 20 20 20 20 20 20 20 | 20 20 20 30 2e 39 39 38 | | 0.998|
|00005830| 32 35 35 37 32 20 20 74 | 68 65 74 61 28 32 29 0a |25572 t|heta(2).|
|00005840| 58 20 20 20 34 2e 30 35 | 38 32 31 36 38 33 34 36 |X 4.05|82168346|
|00005850| 36 36 35 38 32 36 65 2d | 30 31 20 20 20 20 20 20 |665826e-|01 |
|00005860| 20 20 20 20 20 30 2e 34 | 30 35 38 32 31 36 38 20 | 0.4|0582168 |
|00005870| 20 76 61 72 28 31 2c 31 | 29 0a 58 20 20 20 34 2e | var(1,1|).X 4.|
|00005880| 30 36 34 33 34 31 30 39 | 36 35 39 36 30 32 33 31 |06434109|65960231|
|00005890| 65 2d 30 31 20 20 20 20 | 20 20 20 20 20 20 20 30 |e-01 | 0|
|000058a0| 2e 34 30 36 34 33 34 31 | 31 20 20 76 61 72 28 32 |.4064341|1 var(2|
|000058b0| 2c 31 29 0a 58 20 20 20 | 33 2e 39 38 37 33 36 35 |,1).X |3.987365|
|000058c0| 38 32 36 37 36 39 32 34 | 38 30 65 2d 30 31 20 20 |82676924|80e-01 |
|000058d0| 20 20 20 20 20 20 20 20 | 20 30 2e 33 39 38 37 33 | | 0.39873|
|000058e0| 36 35 38 20 20 76 61 72 | 28 33 2c 31 29 0a 58 20 |658 var|(3,1).X |
|000058f0| 20 20 34 2e 30 36 34 33 | 34 31 30 39 36 35 39 36 | 4.0643|41096596|
|00005900| 30 32 33 31 65 2d 30 31 | 20 20 20 20 20 20 20 20 |0231e-01| |
|00005910| 20 20 20 30 2e 34 30 36 | 34 33 34 31 31 20 20 76 | 0.406|43411 v|
|00005920| 61 72 28 31 2c 32 29 0a | 58 20 20 20 34 2e 30 37 |ar(1,2).|X 4.07|
|00005930| 30 35 34 38 37 34 31 34 | 38 36 30 36 35 38 65 2d |05487414|860658e-|
|00005940| 30 31 20 20 20 20 20 20 | 20 20 20 20 20 30 2e 34 |01 | 0.4|
|00005950| 30 37 30 35 34 38 37 20 | 20 76 61 72 28 32 2c 32 |0705487 | var(2,2|
|00005960| 29 0a 58 20 20 20 33 2e | 39 39 33 36 33 32 36 32 |).X 3.|99363262|
|00005970| 33 31 39 33 32 32 32 33 | 65 2d 30 31 20 20 20 20 |31932223|e-01 |
|00005980| 20 20 20 20 20 20 20 30 | 2e 33 39 39 33 36 33 32 | 0|.3993632|
|00005990| 36 20 20 76 61 72 28 33 | 2c 32 29 0a 58 20 20 20 |6 var(3|,2).X |
|000059a0| 33 2e 39 38 37 33 36 35 | 38 32 36 37 36 39 32 34 |3.987365|82676924|
|000059b0| 38 30 65 2d 30 31 20 20 | 20 20 20 20 20 20 20 20 |80e-01 | |
|000059c0| 20 30 2e 33 39 38 37 33 | 36 35 38 20 20 76 61 72 | 0.39873|658 var|
|000059d0| 28 31 2c 33 29 0a 58 20 | 20 20 33 2e 39 39 33 36 |(1,3).X | 3.9936|
|000059e0| 33 32 36 32 33 31 39 33 | 32 32 32 33 65 2d 30 31 |32623193|2223e-01|
|000059f0| 20 20 20 20 20 20 20 20 | 20 20 20 30 2e 33 39 39 | | 0.399|
|00005a00| 33 36 33 32 36 20 20 76 | 61 72 28 32 2c 33 29 0a |36326 v|ar(2,3).|
|00005a10| 58 20 20 20 33 2e 39 32 | 37 32 33 32 37 31 33 36 |X 3.92|72327136|
|00005a20| 33 33 36 30 32 37 65 2d | 30 31 20 20 20 20 20 20 |336027e-|01 |
|00005a30| 20 20 20 20 20 30 2e 33 | 39 32 37 32 33 32 37 20 | 0.3|9272327 |
|00005a40| 20 76 61 72 28 33 2c 33 | 29 0a 58 20 20 20 33 2e | var(3,3|).X 3.|
|00005a50| 35 38 30 30 39 33 32 30 | 37 32 34 37 37 31 32 31 |58009320|72477121|
|00005a60| 65 2b 30 30 20 20 20 20 | 20 20 20 20 20 20 20 33 |e+00 | 3|
|00005a70| 2e 35 38 30 30 39 33 32 | 31 20 20 56 28 31 2c 31 |.5800932|1 V(1,1|
|00005a80| 29 0a 58 20 20 2d 37 2e | 32 31 32 36 38 31 37 36 |).X -7.|21268176|
|00005a90| 39 31 33 39 30 39 31 33 | 65 2d 30 33 20 20 20 20 |91390913|e-03 |
|00005aa0| 20 20 20 20 20 20 2d 30 | 2e 30 30 37 32 31 32 36 | -0|.0072126|
|00005ab0| 38 20 20 56 28 32 2c 31 | 29 0a 58 20 20 2d 37 2e |8 V(2,1|).X -7.|
|00005ac0| 32 31 32 36 38 31 37 36 | 39 31 33 39 30 39 31 33 |21268176|91390913|
|00005ad0| 65 2d 30 33 20 20 20 20 | 20 20 20 20 20 20 2d 30 |e-03 | -0|
|00005ae0| 2e 30 30 37 32 31 32 36 | 38 20 20 56 28 31 2c 32 |.0072126|8 V(1,2|
|00005af0| 29 0a 58 20 20 20 32 2e | 30 36 30 33 31 38 30 39 |).X 2.|06031809|
|00005b00| 30 30 34 35 37 38 33 31 | 65 2d 30 35 20 20 20 20 |00457831|e-05 |
|00005b10| 20 20 20 20 20 20 20 30 | 2e 30 30 30 30 32 30 36 | 0|.0000206|
|00005b20| 30 20 20 56 28 32 2c 32 | 29 0a 58 20 20 2d 39 2e |0 V(2,2|).X -9.|
|00005b30| 31 39 33 31 38 35 39 31 | 37 39 35 34 34 34 34 36 |19318591|79544446|
|00005b40| 65 2d 30 38 20 20 20 20 | 20 20 20 20 20 20 2d 30 |e-08 | -0|
|00005b50| 2e 30 30 30 30 30 30 30 | 39 20 20 44 28 31 29 0a |.0000000|9 D(1).|
|00005b60| 58 20 20 20 31 2e 38 35 | 32 35 39 39 39 37 35 33 |X 1.85|25999753|
|00005b70| 36 30 38 36 30 37 65 2d | 31 30 20 20 20 20 20 20 |608607e-|10 |
|00005b80| 20 20 20 20 20 30 2e 30 | 30 30 30 30 30 30 30 20 | 0.0|0000000 |
|00005b90| 20 44 28 32 29 0a 58 20 | 20 20 31 2e 30 35 36 39 | D(2).X | 1.0569|
|00005ba0| 32 31 37 31 34 38 35 30 | 35 30 34 31 65 2b 30 30 |21714850|5041e+00|
|00005bb0| 20 20 20 20 20 20 20 20 | 20 20 20 31 2e 30 35 36 | | 1.056|
|00005bc0| 39 32 31 37 31 20 20 6f | 62 6a 0a 45 4e 44 5f 4f |92171 o|bj.END_O|
|00005bd0| 46 5f 46 49 4c 45 0a 69 | 66 20 74 65 73 74 20 32 |F_FILE.i|f test 2|
|00005be0| 31 30 31 20 2d 6e 65 20 | 60 77 63 20 2d 63 20 3c |101 -ne |`wc -c <|
|00005bf0| 27 63 68 36 65 67 32 2f | 65 6e 64 69 6e 67 2e 62 |'ch6eg2/|ending.b|
|00005c00| 61 6b 27 60 3b 20 74 68 | 65 6e 0a 20 20 20 20 65 |ak'`; th|en. e|
|00005c10| 63 68 6f 20 73 68 61 72 | 3a 20 5c 22 27 63 68 36 |cho shar|: \"'ch6|
|00005c20| 65 67 32 2f 65 6e 64 69 | 6e 67 2e 62 61 6b 27 5c |eg2/endi|ng.bak'\|
|00005c30| 22 20 75 6e 70 61 63 6b | 65 64 20 77 69 74 68 20 |" unpack|ed with |
|00005c40| 77 72 6f 6e 67 20 73 69 | 7a 65 21 0a 66 69 0a 23 |wrong si|ze!.fi.#|
|00005c50| 20 65 6e 64 20 6f 66 20 | 27 63 68 36 65 67 32 2f | end of |'ch6eg2/|
|00005c60| 65 6e 64 69 6e 67 2e 62 | 61 6b 27 0a 66 69 0a 69 |ending.b|ak'.fi.i|
|00005c70| 66 20 74 65 73 74 20 2d | 66 20 27 63 68 36 65 67 |f test -|f 'ch6eg|
|00005c80| 32 2f 6d 6f 64 65 6c 2e | 63 63 27 20 2d 61 20 22 |2/model.|cc' -a "|
|00005c90| 24 7b 31 7d 22 20 21 3d | 20 22 2d 63 22 20 3b 20 |${1}" !=| "-c" ; |
|00005ca0| 74 68 65 6e 20 0a 20 20 | 65 63 68 6f 20 73 68 61 |then . |echo sha|
|00005cb0| 72 3a 20 57 69 6c 6c 20 | 6e 6f 74 20 63 6c 6f 62 |r: Will |not clob|
|00005cc0| 62 65 72 20 65 78 69 73 | 74 69 6e 67 20 66 69 6c |ber exis|ting fil|
|00005cd0| 65 20 5c 22 27 63 68 36 | 65 67 32 2f 6d 6f 64 65 |e \"'ch6|eg2/mode|
|00005ce0| 6c 2e 63 63 27 5c 22 0a | 65 6c 73 65 0a 65 63 68 |l.cc'\".|else.ech|
|00005cf0| 6f 20 73 68 61 72 3a 20 | 45 78 74 72 61 63 74 69 |o shar: |Extracti|
|00005d00| 6e 67 20 5c 22 27 63 68 | 36 65 67 32 2f 6d 6f 64 |ng \"'ch|6eg2/mod|
|00005d10| 65 6c 2e 63 63 27 5c 22 | 20 5c 28 31 37 39 36 20 |el.cc'\"| \(1796 |
|00005d20| 63 68 61 72 61 63 74 65 | 72 73 5c 29 0a 73 65 64 |characte|rs\).sed|
|00005d30| 20 22 73 2f 5e 58 2f 2f | 22 20 3e 27 63 68 36 65 | "s/^X//|" >'ch6e|
|00005d40| 67 32 2f 6d 6f 64 65 6c | 2e 63 63 27 20 3c 3c 27 |g2/model|.cc' <<'|
|00005d50| 45 4e 44 5f 4f 46 5f 46 | 49 4c 45 27 0a 58 23 69 |END_OF_F|ILE'.X#i|
|00005d60| 6e 63 6c 75 64 65 20 22 | 6d 6f 64 65 6c 2e 68 22 |nclude "|model.h"|
|00005d70| 0a 58 0a 58 6d 6f 64 65 | 6c 3a 3a 6d 6f 64 65 6c |.X.Xmode|l::model|
|00005d80| 28 29 20 7b 20 7d 0a 58 | 0a 58 6d 6f 64 65 6c 3a |() { }.X|.Xmodel:|
|00005d90| 3a 7e 6d 6f 64 65 6c 28 | 29 20 7b 20 7d 0a 58 0a |:~model(|) { }.X.|
|00005da0| 58 0a 58 72 65 61 6c 6d | 61 74 20 6d 6f 64 65 6c |X.Xrealm|at model|
|00005db0| 3a 3a 65 28 49 4e 54 45 | 47 45 52 20 74 29 0a 58 |::e(INTE|GER t).X|
|00005dc0| 7b 0a 58 20 20 72 65 61 | 6c 6d 61 74 20 65 74 28 |{.X rea|lmat et(|
|00005dd0| 31 2c 31 29 3b 0a 58 20 | 20 52 45 41 4c 20 20 20 |1,1);.X | REAL |
|00005de0| 20 61 2c 62 2c 79 74 61 | 2c 79 74 2c 78 74 3b 0a | a,b,yta|,yt,xt;.|
|00005df0| 58 20 20 49 4e 54 45 47 | 45 52 20 69 3b 0a 58 0a |X INTEG|ER i;.X.|
|00005e00| 58 20 20 69 3d 74 2b 32 | 3b 0a 58 20 20 0a 58 20 |X i=t+2|;.X .X |
|00005e10| 20 79 74 20 3d 20 28 6e | 64 73 5b 69 5d 2f 70 6f | yt = (n|ds[i]/po|
|00005e20| 70 5b 69 5d 29 2f 28 6e | 64 73 5b 69 2d 31 5d 2f |p[i])/(n|ds[i-1]/|
|00005e30| 70 6f 70 5b 69 2d 31 5d | 29 3b 0a 58 20 20 78 74 |pop[i-1]|);.X xt|
|00005e40| 20 3d 20 28 31 2e 30 2b | 6e 79 73 65 5b 69 5d 29 | = (1.0+|nyse[i])|
|00005e50| 2a 28 64 66 6c 61 74 5b | 69 2d 31 5d 2f 64 66 6c |*(dflat[|i-1]/dfl|
|00005e60| 61 74 5b 69 5d 29 3b 0a | 58 0a 58 20 20 61 3d 73 |at[i]);.|X.X a=s|
|00005e70| 2e 74 68 65 74 61 5b 31 | 5d 3b 0a 58 20 20 62 3d |.theta[1|];.X b=|
|00005e80| 73 2e 74 68 65 74 61 5b | 32 5d 3b 0a 58 0a 58 20 |s.theta[|2];.X.X |
|00005e90| 20 79 74 61 20 3d 20 70 | 6f 77 28 79 74 2c 61 29 | yta = p|ow(yt,a)|
|00005ea0| 3b 0a 58 0a 58 20 20 65 | 74 5b 31 5d 20 3d 20 62 |;.X.X e|t[1] = b|
|00005eb0| 20 2a 20 79 74 61 20 2a | 20 78 74 20 2d 20 31 2e | * yta *| xt - 1.|
|00005ec0| 30 3b 0a 58 0a 58 20 20 | 72 65 74 75 72 6e 20 65 |0;.X.X |return e|
|00005ed0| 74 3b 0a 58 7d 0a 58 0a | 58 72 65 61 6c 6d 61 74 |t;.X}.X.|Xrealmat|
|00005ee0| 20 6d 6f 64 65 6c 3a 3a | 64 65 6c 65 28 49 4e 54 | model::|dele(INT|
|00005ef0| 45 47 45 52 20 74 29 20 | 0a 58 7b 0a 58 20 20 72 |EGER t) |.X{.X r|
|00005f00| 65 61 6c 6d 61 74 20 64 | 65 72 71 74 28 31 2c 32 |ealmat d|erqt(1,2|
|00005f10| 29 3b 0a 58 0a 58 20 20 | 52 45 41 4c 20 20 20 20 |);.X.X |REAL |
|00005f20| 61 2c 62 2c 79 74 61 2c | 79 74 2c 78 74 3b 0a 58 |a,b,yta,|yt,xt;.X|
|00005f30| 20 20 49 4e 54 45 47 45 | 52 20 69 3b 0a 58 0a 58 | INTEGE|R i;.X.X|
|00005f40| 20 20 69 3d 74 2b 32 3b | 0a 58 20 20 0a 58 20 20 | i=t+2;|.X .X |
|00005f50| 79 74 20 3d 20 28 6e 64 | 73 5b 69 5d 2f 70 6f 70 |yt = (nd|s[i]/pop|
|00005f60| 5b 69 5d 29 2f 28 6e 64 | 73 5b 69 2d 31 5d 2f 70 |[i])/(nd|s[i-1]/p|
|00005f70| 6f 70 5b 69 2d 31 5d 29 | 3b 0a 58 20 20 78 74 20 |op[i-1])|;.X xt |
|00005f80| 3d 20 28 31 2e 30 2b 6e | 79 73 65 5b 69 5d 29 2a |= (1.0+n|yse[i])*|
|00005f90| 28 64 66 6c 61 74 5b 69 | 2d 31 5d 2f 64 66 6c 61 |(dflat[i|-1]/dfla|
|00005fa0| 74 5b 69 5d 29 3b 0a 58 | 0a 58 20 20 61 3d 73 2e |t[i]);.X|.X a=s.|
|00005fb0| 74 68 65 74 61 5b 31 5d | 3b 0a 58 20 20 62 3d 73 |theta[1]|;.X b=s|
|00005fc0| 2e 74 68 65 74 61 5b 32 | 5d 3b 0a 58 0a 58 20 20 |.theta[2|];.X.X |
|00005fd0| 79 74 61 20 3d 20 70 6f | 77 28 79 74 2c 61 29 3b |yta = po|w(yt,a);|
|00005fe0| 0a 58 0a 58 20 20 64 65 | 72 71 74 2e 65 6c 65 6d |.X.X de|rqt.elem|
|00005ff0| 28 31 2c 31 29 20 3d 20 | 62 20 2a 20 6c 6f 67 28 |(1,1) = |b * log(|
|00006000| 79 74 29 20 2a 20 79 74 | 61 20 2a 20 78 74 3b 0a |yt) * yt|a * xt;.|
|00006010| 58 0a 58 20 20 64 65 72 | 71 74 2e 65 6c 65 6d 28 |X.X der|qt.elem(|
|00006020| 31 2c 32 29 20 3d 20 79 | 74 61 2a 78 74 3b 0a 58 |1,2) = y|ta*xt;.X|
|00006030| 0a 58 20 20 72 65 74 75 | 72 6e 20 64 65 72 71 74 |.X retu|rn derqt|
|00006040| 3b 0a 58 0a 58 7d 0a 58 | 0a 58 0a 58 69 6e 74 20 |;.X.X}.X|.X.Xint |
|00006050| 6d 6f 64 65 6c 3a 3a 69 | 6e 69 74 69 61 6c 69 7a |model::i|nitializ|
|00006060| 65 28 29 0a 58 7b 0a 58 | 20 20 49 4e 54 45 47 45 |e().X{.X| INTEGE|
|00006070| 52 20 69 3b 0a 58 20 20 | 64 6f 75 62 6c 65 20 74 |R i;.X |double t|
|00006080| 6d 70 3b 0a 58 0a 58 23 | 69 66 64 65 66 20 20 47 |mp;.X.X#|ifdef G|
|00006090| 4e 55 5f 47 50 50 5f 43 | 4f 4d 50 49 4c 45 52 0a |NU_GPP_C|OMPILER.|
|000060a0| 58 23 69 66 64 65 66 20 | 20 55 53 45 5f 41 54 54 |X#ifdef | USE_ATT|
|000060b0| 5f 53 54 59 4c 45 5f 49 | 4f 5f 57 49 54 48 5f 47 |_STYLE_I|O_WITH_G|
|000060c0| 4e 55 0a 58 0a 58 20 20 | 66 69 6c 65 62 75 66 20 |NU.X.X |filebuf |
|000060d0| 63 64 61 74 61 5f 62 75 | 66 3b 0a 58 20 20 69 66 |cdata_bu|f;.X if|
|000060e0| 20 28 63 64 61 74 61 5f | 62 75 66 2e 6f 70 65 6e | (cdata_|buf.open|
|000060f0| 28 22 68 61 6e 73 65 6e | 2e 64 61 74 22 2c 69 6e |("hansen|.dat",in|
|00006100| 70 75 74 29 20 3d 3d 30 | 29 20 7b 0a 58 20 20 20 |put) ==0|) {.X |
|00006110| 20 63 65 72 72 20 3c 3c | 20 22 45 72 72 6f 72 2c | cerr <<| "Error,|
|00006120| 20 6d 6f 64 65 6c 3a 3a | 6d 6f 64 65 6c 2c 20 43 | model::|model, C|
|00006130| 61 6e 6e 6f 74 20 6f 70 | 65 6e 20 68 61 6e 73 65 |annot op|en hanse|
|00006140| 6e 2e 64 61 74 2e 5c 6e | 22 3b 0a 58 20 20 20 20 |n.dat.\n|";.X |
|00006150| 65 78 69 74 28 31 29 3b | 0a 58 20 20 7d 0a 58 20 |exit(1);|.X }.X |
|00006160| 20 69 73 74 72 65 61 6d | 20 63 64 61 74 61 28 26 | istream| cdata(&|
|00006170| 63 64 61 74 61 5f 62 75 | 66 29 3b 0a 58 0a 58 23 |cdata_bu|f);.X.X#|
|00006180| 65 6e 64 69 66 0a 58 0a | 58 23 69 66 64 65 66 20 |endif.X.|X#ifdef |
|00006190| 20 55 53 45 5f 47 4e 55 | 5f 53 54 59 4c 45 5f 49 | USE_GNU|_STYLE_I|
|000061a0| 4f 5f 57 49 54 48 5f 47 | 4e 55 0a 58 0a 58 20 20 |O_WITH_G|NU.X.X |
|000061b0| 69 73 74 72 65 61 6d 20 | 63 64 61 74 61 28 22 68 |istream |cdata("h|
|000061c0| 61 6e 73 65 6e 2e 64 61 | 74 22 2c 69 6f 5f 72 65 |ansen.da|t",io_re|
|000061d0| 61 64 6f 6e 6c 79 2c 61 | 5f 75 73 65 6f 6e 6c 79 |adonly,a|_useonly|
|000061e0| 29 3b 0a 58 20 20 69 66 | 20 28 20 21 63 64 61 74 |);.X if| ( !cdat|
|000061f0| 61 20 29 20 7b 0a 58 20 | 20 20 20 63 65 72 72 20 |a ) {.X | cerr |
|00006200| 3c 3c 20 22 45 72 72 6f | 72 2c 20 6d 6f 64 65 6c |<< "Erro|r, model|
|00006210| 3a 3a 6d 6f 64 65 6c 2c | 20 43 61 6e 6e 6f 74 20 |::model,| Cannot |
|00006220| 6f 70 65 6e 20 68 61 6e | 73 65 6e 2e 64 61 74 2e |open han|sen.dat.|
|00006230| 5c 6e 22 3b 0a 58 20 20 | 20 20 65 78 69 74 28 31 |\n";.X | exit(1|
|00006240| 29 3b 0a 58 20 20 7d 0a | 58 0a 58 23 65 6e 64 69 |);.X }.|X.X#endi|
|00006250| 66 0a 58 0a 58 23 65 6e | 64 69 66 0a 58 0a 58 23 |f.X.X#en|dif.X.X#|
|00006260| 69 66 64 65 66 20 20 54 | 55 52 42 4f 5f 43 50 50 |ifdef T|URBO_CPP|
|00006270| 5f 43 4f 4d 50 49 4c 45 | 52 0a 58 0a 58 20 20 69 |_COMPILE|R.X.X i|
|00006280| 66 73 74 72 65 61 6d 20 | 63 64 61 74 61 28 22 68 |fstream |cdata("h|
|00006290| 61 6e 73 65 6e 2e 64 61 | 74 22 29 3b 0a 58 20 20 |ansen.da|t");.X |
|000062a0| 69 66 20 28 21 63 64 61 | 74 61 29 20 7b 0a 58 20 |if (!cda|ta) {.X |
|000062b0| 20 20 20 63 65 72 72 20 | 3c 3c 20 22 45 72 72 6f | cerr |<< "Erro|
|000062c0| 72 2c 20 6d 6f 64 65 6c | 3a 3a 6d 6f 64 65 6c 2c |r, model|::model,|
|000062d0| 20 43 61 6e 6e 6f 74 20 | 6f 70 65 6e 20 68 61 6e | Cannot |open han|
|000062e0| 73 65 6e 2e 64 61 74 2e | 5c 6e 22 3b 0a 58 20 20 |sen.dat.|\n";.X |
|000062f0| 20 20 65 78 69 74 28 31 | 29 3b 0a 58 20 20 7d 0a | exit(1|);.X }.|
|00006300| 58 0a 58 23 65 6e 64 69 | 66 0a 58 0a 58 0a 58 20 |X.X#endi|f.X.X.X |
|00006310| 20 66 6f 72 20 28 69 3d | 31 3b 20 69 3c 3d 32 34 | for (i=|1; i<=24|
|00006320| 30 3b 20 69 2b 2b 29 20 | 7b 0a 58 20 20 20 20 63 |0; i++) |{.X c|
|00006330| 64 61 74 61 20 3e 3e 20 | 74 6d 70 3b 20 6e 64 73 |data >> |tmp; nds|
|00006340| 5b 69 5d 3d 28 52 45 41 | 4c 29 74 6d 70 3b 0a 58 |[i]=(REA|L)tmp;.X|
|00006350| 20 20 20 20 63 64 61 74 | 61 20 3e 3e 20 74 6d 70 | cdat|a >> tmp|
|00006360| 3b 20 70 6f 70 5b 69 5d | 3d 28 52 45 41 4c 29 74 |; pop[i]|=(REAL)t|
|00006370| 6d 70 3b 0a 58 20 20 20 | 20 63 64 61 74 61 20 3e |mp;.X | cdata >|
|00006380| 3e 20 74 6d 70 3b 20 6e | 79 73 65 5b 69 5d 3d 28 |> tmp; n|yse[i]=(|
|00006390| 52 45 41 4c 29 74 6d 70 | 3b 0a 58 20 20 20 20 63 |REAL)tmp|;.X c|
|000063a0| 64 61 74 61 20 3e 3e 20 | 74 6d 70 3b 20 64 66 6c |data >> |tmp; dfl|
|000063b0| 61 74 5b 69 5d 3d 28 52 | 45 41 4c 29 74 6d 70 3b |at[i]=(R|EAL)tmp;|
|000063c0| 0a 58 20 20 7d 0a 58 0a | 58 20 20 72 65 74 75 72 |.X }.X.|X retur|
|000063d0| 6e 20 31 3b 0a 58 7d 0a | 58 0a 58 72 65 61 6c 6d |n 1;.X}.|X.Xrealm|
|000063e0| 61 74 20 6d 6f 64 65 6c | 3a 3a 5a 28 49 4e 54 45 |at model|::Z(INTE|
|000063f0| 47 45 52 20 74 29 0a 58 | 7b 0a 58 20 20 49 4e 54 |GER t).X|{.X INT|
+--------+-------------------------+-------------------------+--------+--------+
Only 25.0 KB of data is shown above.
