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

< prev

next >

view JSON data

view as text

open on a Mac

open on a PC

+--------+-------------------------+-------------------------+--------+--------+
|00000000| 44 45 52 49 56 45 20 66 | 6f 72 20 57 69 6e 64 6f |DERIVE f|or Windo|
|00000010| 77 73 20 76 65 72 73 69 | 6f 6e 20 35 2e 78 78 20 |ws versi|on 5.xx |
|00000020| 44 66 57 20 66 69 6c 65 | 20 73 61 76 65 64 20 6f |DfW file| saved o|
|00000030| 6e 20 30 32 20 4e 6f 76 | 20 32 30 30 31 0d 0a 1a |n 02 Nov| 2001...|
|00000040| 1b 00 00 00 41 50 50 45 | 4e 44 5f 43 4f 4c 55 4d |....APPE|ND_COLUM|
|00000050| 4e 53 28 61 2c 20 62 29 | 3a 3d 41 50 50 45 4e 44 |NS(a, b)|:=APPEND|
|00000060| 28 61 60 2c 20 62 60 29 | 60 0d 0a 43 48 49 5f 53 |(a`, b`)|`..CHI_S|
|00000070| 51 55 41 52 45 28 78 32 | 2c 20 76 29 3a 3d 49 4e |QUARE(x2|, v):=IN|
|00000080| 43 4f 4d 50 4c 45 54 45 | 5f 47 41 4d 4d 41 28 76 |COMPLETE|_GAMMA(v|
|00000090| 2f 32 2c 20 78 32 2f 32 | 29 0d 0a 49 4e 43 4f 4d |/2, x2/2|)..INCOM|
|000000a0| 50 4c 45 54 45 5f 47 41 | 4d 4d 41 28 7a 2c 20 77 |PLETE_GA|MMA(z, w|
|000000b0| 29 3a 3d 28 77 5e 7a 2f | 7a 20 2b 20 49 4e 54 28 |):=(w^z/|z + INT(|
|000000c0| 74 5f 5e 28 7a 20 2d 20 | 31 29 b7 28 ea 5e 28 2d |t_^(z - |1).(.^(-|
|000000d0| 74 5f 29 20 2d 20 31 29 | 2c 20 74 5f 2c 20 30 2c |t_) - 1)|, t_, 0,|
|000000e0| 20 77 29 29 2f 9b 28 7a | 29 0d 0a 50 72 65 64 69 | w))/.(z|)..Predi|
|000000f0| 63 74 65 64 5f 56 61 6c | 75 65 73 28 76 61 6c 2c |cted_Val|ues(val,|
|00000100| 20 63 69 20 3a 3d 20 30 | 2e 39 35 29 3a 3d 50 52 | ci := 0|.95):=PR|
|00000110| 4f 47 28 6c 20 3a 3d 20 | 53 55 42 53 54 28 66 65 |OG(l := |SUBST(fe|
|00000120| 71 2c 20 69 76 61 72 73 | 5f 2c 20 76 61 6c 29 2c |q, ivars|_, val),|
|00000130| 20 78 31 20 3a 3d 20 5b | 41 44 4a 4f 49 4e 28 31 | x1 := [|ADJOIN(1|
|00000140| 2c 20 76 61 6c 29 5d 2c | 20 73 20 3a 3d 20 8b 28 |, val)],| s := .(|
|00000150| 28 78 31 20 95 20 78 70 | 78 69 20 95 20 78 31 60 |(x1 . xp|xi . x1`|
|00000160| 29 99 31 99 31 29 2c 20 | 69 6e 76 6e 20 3a 3d 20 |).1.1), |invn := |
|00000170| 52 48 53 28 4e 53 4f 4c | 56 45 28 4e 4f 52 4d 41 |RHS(NSOL|VE(NORMA|
|00000180| 4c 28 7a 29 20 2d 20 63 | 69 2f 32 20 2d 20 30 2e |L(z) - c|i/2 - 0.|
|00000190| 35 2c 20 7a 2c 20 30 2c | 20 96 29 29 2c 20 6c 6c |5, z, 0,| .)), ll|
|000001a0| 20 3a 3d 20 6c 20 2d 20 | 69 6e 76 6e b7 73 2c 20 | := l - |invn.s, |
|000001b0| 75 6c 20 3a 3d 20 6c 20 | 2b 20 69 6e 76 6e b7 73 |ul := l |+ invn.s|
|000001c0| 2c 20 49 46 28 6c 69 6e | 6b 5f 20 3d 20 31 2c 20 |, IF(lin|k_ = 1, |
|000001d0| 50 52 4f 47 28 6c 5f 20 | 3a 3d 20 22 6c 6f 67 69 |PROG(l_ |:= "logi|
|000001e0| 74 28 70 72 6f 62 29 22 | 20 3d 20 6c 2c 20 70 20 |t(prob)"| = l, p |
|000001f0| 3a 3d 20 31 2f 28 31 20 | 2b 20 ea 5e 28 2d 6c 29 |:= 1/(1 |+ .^(-l)|
|00000200| 29 2c 20 6c 70 20 3a 3d | 20 31 2f 28 31 20 2b 20 |), lp :=| 1/(1 + |
|00000210| ea 5e 28 2d 6c 6c 29 29 | 2c 20 75 70 20 3a 3d 20 |.^(-ll))|, up := |
|00000220| 31 2f 28 31 20 2b 20 ea | 5e 28 2d 75 6c 29 29 29 |1/(1 + .|^(-ul)))|
|00000230| 2c 20 49 46 28 6c 69 6e | 6b 5f 20 3d 20 32 2c 20 |, IF(lin|k_ = 2, |
|00000240| 50 52 4f 47 28 6c 5f 20 | 3a 3d 20 22 70 72 6f 62 |PROG(l_ |:= "prob|
|00000250| 69 74 28 70 72 6f 62 29 | 22 20 3d 20 6c 2c 20 70 |it(prob)|" = l, p|
|00000260| 20 3a 3d 20 4e 4f 52 4d | 41 4c 28 6c 29 2c 20 6c | := NORM|AL(l), l|
|00000270| 70 20 3a 3d 20 4e 4f 52 | 4d 41 4c 28 6c 6c 29 2c |p := NOR|MAL(ll),|
|00000280| 20 75 70 20 3a 3d 20 4e | 4f 52 4d 41 4c 28 75 6c | up := N|ORMAL(ul|
|00000290| 29 29 2c 20 50 52 4f 47 | 28 6c 5f 20 3a 3d 20 22 |)), PROG|(l_ := "|
|000002a0| 63 6c 6f 67 6c 6f 67 28 | 70 72 6f 62 29 22 20 3d |cloglog(|prob)" =|
|000002b0| 20 6c 2c 20 70 20 3a 3d | 20 31 20 2d 20 ea 5e 28 | l, p :=| 1 - .^(|
|000002c0| 2d 20 ea 5e 6c 29 2c 20 | 6c 70 20 3a 3d 20 31 20 |- .^l), |lp := 1 |
|000002d0| 2d 20 ea 5e 28 2d 20 ea | 5e 6c 6c 29 2c 20 75 70 |- .^(- .|^ll), up|
|000002e0| 20 3a 3d 20 31 20 2d 20 | ea 5e 28 2d 20 ea 5e 75 | := 1 - |.^(- .^u|
|000002f0| 6c 29 29 29 29 2c 20 5b | 6c 5f 2c 20 22 73 65 22 |l)))), [|l_, "se"|
|00000300| 20 3d 20 73 2c 20 22 70 | 72 6f 62 22 20 3d 20 70 | = s, "p|rob" = p|
|00000310| 3b 20 41 50 50 45 4e 44 | 28 53 54 52 49 4e 47 28 |; APPEND|(STRING(|
|00000320| 31 30 30 b7 63 69 29 2c | 20 22 25 20 43 2e 49 2e |100.ci),| "% C.I.|
|00000330| 20 70 72 6f 62 22 29 2c | 20 6c 70 2c 20 75 70 5d | prob"),| lp, up]|
|00000340| 29 0d 0a 6c 6f 67 69 73 | 74 69 63 28 64 76 61 72 |)..logis|tic(dvar|
|00000350| 2c 20 69 76 61 72 73 2c | 20 64 61 74 61 2c 20 6c |, ivars,| data, l|
|00000360| 69 6e 6b 20 3a 3d 20 31 | 2c 20 66 72 65 71 20 3a |ink := 1|, freq :|
|00000370| 3d 20 5b 5d 2c 20 6e 6e | 20 3a 3d 20 33 30 2c 20 |= [], nn| := 30, |
|00000380| 63 63 20 3a 3d 20 31 30 | 5e 28 2d 38 29 29 3a 3d |cc := 10|^(-8)):=|
|00000390| 50 52 4f 47 28 76 61 72 | 73 20 3a 3d 20 64 61 74 |PROG(var|s := dat|
|000003a0| 61 99 31 2c 20 6f 62 73 | 5f 20 3a 3d 20 44 49 4d |a.1, obs|_ := DIM|
|000003b0| 28 64 61 74 61 29 20 2d | 20 31 2c 20 69 76 61 72 |(data) -| 1, ivar|
|000003c0| 73 5f 20 3a 3d 20 69 76 | 61 72 73 2c 20 6c 69 6e |s_ := iv|ars, lin|
|000003d0| 6b 5f 20 3a 3d 20 6c 69 | 6e 6b 2c 20 64 61 74 61 |k_ := li|nk, data|
|000003e0| 31 20 3a 3d 20 64 61 74 | 61 99 5b 32 2c 20 2e 2e |1 := dat|a.[2, ..|
|000003f0| 2e 2c 20 6f 62 73 5f 20 | 2b 20 31 5d 2c 20 6b 5f |., obs_ |+ 1], k_|
|00000400| 20 3a 3d 20 44 49 4d 28 | 69 76 61 72 73 29 2c 20 | := DIM(|ivars), |
|00000410| 66 72 65 71 5f 20 3a 3d | 20 44 49 4d 28 66 72 65 |freq_ :=| DIM(fre|
|00000420| 71 29 2c 20 49 46 28 66 | 72 65 71 5f 20 3d 20 30 |q), IF(f|req_ = 0|
|00000430| 2c 20 77 67 74 20 3a 3d | 20 56 45 43 54 4f 52 28 |, wgt :=| VECTOR(|
|00000440| 31 2c 20 69 5f 2c 20 31 | 2c 20 6f 62 73 5f 29 2c |1, i_, 1|, obs_),|
|00000450| 20 49 46 28 66 72 65 71 | 5f 20 3d 20 31 2c 20 50 | IF(freq|_ = 1, P|
|00000460| 52 4f 47 28 70 6f 73 5f | 20 3a 3d 20 70 6f 73 69 |ROG(pos_| := posi|
|00000470| 74 69 6f 6e 28 66 72 65 | 71 99 31 2c 20 76 61 72 |tion(fre|q.1, var|
|00000480| 73 29 2c 20 77 67 74 20 | 3a 3d 20 64 61 74 61 99 |s), wgt |:= data.|
|00000490| 99 70 6f 73 5f 29 2c 20 | 49 46 28 66 72 65 71 5f |.pos_), |IF(freq_|
|000004a0| 20 3d 20 32 2c 20 50 52 | 4f 47 28 70 6f 73 5f 20 | = 2, PR|OG(pos_ |
|000004b0| 3a 3d 20 70 6f 73 69 74 | 69 6f 6e 28 66 72 65 71 |:= posit|ion(freq|
|000004c0| 99 31 2c 20 76 61 72 73 | 29 2c 20 65 76 74 20 3a |.1, vars|), evt :|
|000004d0| 3d 20 64 61 74 61 31 99 | 99 70 6f 73 5f 2c 20 70 |= data1.|.pos_, p|
|000004e0| 6f 73 5f 20 3a 3d 20 70 | 6f 73 69 74 69 6f 6e 28 |os_ := p|osition(|
|000004f0| 66 72 65 71 99 32 2c 20 | 76 61 72 73 29 2c 20 74 |freq.2, |vars), t|
|00000500| 72 69 20 3a 3d 20 64 61 | 74 61 31 99 99 70 6f 73 |ri := da|ta1..pos|
|00000510| 5f 2c 20 77 67 74 20 3a | 3d 20 41 50 50 45 4e 44 |_, wgt :|= APPEND|
|00000520| 28 65 76 74 2c 20 74 72 | 69 20 2d 20 65 76 74 29 |(evt, tr|i - evt)|
|00000530| 29 2c 20 52 45 54 55 52 | 4e 20 22 49 6e 76 61 6c |), RETUR|N "Inval|
|00000540| 69 64 20 66 72 65 71 75 | 65 6e 63 79 20 76 61 72 |id frequ|ency var|
|00000550| 69 61 62 6c 65 22 29 29 | 29 2c 20 49 46 28 66 72 |iable"))|), IF(fr|
|00000560| 65 71 5f 20 3d 20 30 20 | 90 20 66 72 65 71 5f 20 |eq_ = 0 |. freq_ |
|00000570| 3d 20 31 2c 20 50 52 4f | 47 28 70 6f 73 5f 20 3a |= 1, PRO|G(pos_ :|
|00000580| 3d 20 70 6f 73 69 74 69 | 6f 6e 28 64 76 61 72 2c |= positi|on(dvar,|
|00000590| 20 76 61 72 73 29 2c 20 | 79 5f 20 3a 3d 20 64 61 | vars), |y_ := da|
|000005a0| 74 61 31 99 99 70 6f 73 | 5f 29 2c 20 49 46 28 66 |ta1..pos|_), IF(f|
|000005b0| 72 65 71 5f 20 3d 20 32 | 2c 20 79 5f 20 3a 3d 20 |req_ = 2|, y_ := |
|000005c0| 41 50 50 45 4e 44 28 56 | 45 43 54 4f 52 28 31 2c |APPEND(V|ECTOR(1,|
|000005d0| 20 69 5f 2c 20 31 2c 20 | 6f 62 73 5f 29 2c 20 56 | i_, 1, |obs_), V|
|000005e0| 45 43 54 4f 52 28 30 2c | 20 69 5f 2c 20 31 2c 20 |ECTOR(0,| i_, 1, |
|000005f0| 6f 62 73 5f 29 29 29 29 | 2c 20 70 6f 73 5f 20 3a |obs_))))|, pos_ :|
|00000600| 3d 20 56 45 43 54 4f 52 | 28 70 6f 73 69 74 69 6f |= VECTOR|(positio|
|00000610| 6e 28 69 76 61 72 73 5f | 99 69 5f 2c 20 76 61 72 |n(ivars_|.i_, var|
|00000620| 73 29 2c 20 69 5f 2c 20 | 31 2c 20 6b 5f 29 2c 20 |s), i_, |1, k_), |
|00000630| 49 46 28 66 72 65 71 5f | 20 3d 20 32 2c 20 50 52 |IF(freq_| = 2, PR|
|00000640| 4f 47 28 78 5f 20 3a 3d | 20 41 50 50 45 4e 44 28 |OG(x_ :=| APPEND(|
|00000650| 64 61 74 61 31 99 99 70 | 6f 73 5f 2c 20 64 61 74 |data1..p|os_, dat|
|00000660| 61 31 99 99 70 6f 73 5f | 29 2c 20 6f 62 73 5f 20 |a1..pos_|), obs_ |
|00000670| 3a 3d 20 32 b7 6f 62 73 | 5f 29 2c 20 78 5f 20 3a |:= 2.obs|_), x_ :|
|00000680| 3d 20 64 61 74 61 31 99 | 99 70 6f 73 5f 29 2c 20 |= data1.|.pos_), |
|00000690| 78 5f 20 3a 3d 20 56 45 | 43 54 4f 52 28 53 55 42 |x_ := VE|CTOR(SUB|
|000006a0| 53 54 28 69 76 61 72 73 | 2c 20 69 76 61 72 73 5f |ST(ivars|, ivars_|
|000006b0| 2c 20 78 5f 99 69 5f 29 | 2c 20 69 5f 2c 20 31 2c |, x_.i_)|, i_, 1,|
|000006c0| 20 6f 62 73 5f 29 2c 20 | 78 31 20 3a 3d 20 56 45 | obs_), |x1 := VE|
|000006d0| 43 54 4f 52 28 56 45 43 | 54 4f 52 28 31 2c 20 69 |CTOR(VEC|TOR(1, i|
|000006e0| 5f 2c 20 31 2c 20 6f 62 | 73 5f 29 2c 20 6a 2c 20 |_, 1, ob|s_), j, |
|000006f0| 31 2c 20 31 29 60 2c 20 | 78 78 5f 20 3a 3d 20 41 |1, 1)`, |xx_ := A|
|00000700| 50 50 45 4e 44 5f 43 4f | 4c 55 4d 4e 53 28 78 31 |PPEND_CO|LUMNS(x1|
|00000710| 2c 20 78 5f 29 2c 20 df | 5f 20 3a 3d 20 56 45 43 |, x_), .|_ := VEC|
|00000720| 54 4f 52 28 30 2c 20 69 | 5f 2c 20 31 2c 20 6b 5f |TOR(0, i|_, 1, k_|
|00000730| 20 2b 20 31 29 2c 20 6f | 6c 64 df 20 3a 3d 20 56 | + 1), o|ld. := V|
|00000740| 45 43 54 4f 52 28 31 2c | 20 69 5f 2c 20 31 2c 20 |ECTOR(1,| i_, 1, |
|00000750| 6b 5f 20 2b 20 31 29 2c | 20 69 74 65 72 5f 20 3a |k_ + 1),| iter_ :|
|00000760| 3d 20 30 2c 20 4c 4f 4f | 50 28 63 63 5f 20 3a 3d |= 0, LOO|P(cc_ :=|
|00000770| 20 4d 41 58 28 56 45 43 | 54 4f 52 28 41 42 53 28 | MAX(VEC|TOR(ABS(|
|00000780| df 5f 99 69 5f 20 2d 20 | 6f 6c 64 df 99 69 5f 29 |._.i_ - |old..i_)|
|00000790| 2c 20 69 5f 2c 20 31 2c | 20 6b 5f 20 2b 20 31 29 |, i_, 1,| k_ + 1)|
|000007a0| 29 2c 20 49 46 28 69 74 | 65 72 5f 20 3e 20 6e 6e |), IF(it|er_ > nn|
|000007b0| 20 90 20 63 63 5f 20 3c | 20 63 63 2c 20 65 78 69 | . cc_ <| cc, exi|
|000007c0| 74 29 2c 20 6f 6c 64 df | 20 3a 3d 20 df 5f 2c 20 |t), old.| := ._, |
|000007d0| 7a 5f 20 3a 3d 20 78 78 | 5f 20 95 20 df 5f 2c 20 |z_ := xx|_ . ._, |
|000007e0| 49 46 28 6c 69 6e 6b 20 | 3d 20 31 2c 20 50 52 4f |IF(link |= 1, PRO|
|000007f0| 47 28 70 5f 20 3a 3d 20 | 56 45 43 54 4f 52 28 31 |G(p_ := |VECTOR(1|
|00000800| 2f 28 31 20 2b 20 ea 5e | 28 2d 20 7a 5f 99 69 5f |/(1 + .^|(- z_.i_|
|00000810| 29 29 2c 20 69 5f 2c 20 | 31 2c 20 6f 62 73 5f 29 |)), i_, |1, obs_)|
|00000820| 2c 20 66 5f 20 3a 3d 20 | 56 45 43 54 4f 52 28 ea |, f_ := |VECTOR(.|
|00000830| 5e 28 2d 20 7a 5f 99 69 | 5f 29 b7 70 5f 99 69 5f |^(- z_.i|_).p_.i_|
|00000840| 5e 32 2c 20 69 5f 2c 20 | 31 2c 20 6f 62 73 5f 29 |^2, i_, |1, obs_)|
|00000850| 29 2c 20 49 46 28 6c 69 | 6e 6b 20 3d 20 32 2c 20 |), IF(li|nk = 2, |
|00000860| 50 52 4f 47 28 70 5f 20 | 3a 3d 20 56 45 43 54 4f |PROG(p_ |:= VECTO|
|00000870| 52 28 4e 4f 52 4d 41 4c | 28 7a 5f 99 69 5f 29 2c |R(NORMAL|(z_.i_),|
|00000880| 20 69 5f 2c 20 31 2c 20 | 6f 62 73 5f 29 2c 20 66 | i_, 1, |obs_), f|
|00000890| 5f 20 3a 3d 20 56 45 43 | 54 4f 52 28 6e 6f 72 6d |_ := VEC|TOR(norm|
|000008a0| 70 64 66 28 7a 5f 99 69 | 5f 29 2c 20 69 5f 2c 20 |pdf(z_.i|_), i_, |
|000008b0| 31 2c 20 6f 62 73 5f 29 | 29 2c 20 49 46 28 6c 69 |1, obs_)|), IF(li|
|000008c0| 6e 6b 20 3d 20 33 2c 20 | 50 52 4f 47 28 70 5f 20 |nk = 3, |PROG(p_ |
|000008d0| 3a 3d 20 56 45 43 54 4f | 52 28 31 20 2d 20 ea 5e |:= VECTO|R(1 - .^|
|000008e0| 28 2d 20 ea 5e 7a 5f 99 | 69 5f 29 2c 20 69 5f 2c |(- .^z_.|i_), i_,|
|000008f0| 20 31 2c 20 6f 62 73 5f | 29 2c 20 66 5f 20 3a 3d | 1, obs_|), f_ :=|
|00000900| 20 56 45 43 54 4f 52 28 | ea 5e 28 7a 5f 99 69 5f | VECTOR(|.^(z_.i_|
|00000910| 20 2d 20 ea 5e 7a 5f 99 | 69 5f 29 2c 20 69 5f 2c | - .^z_.|i_), i_,|
|00000920| 20 31 2c 20 6f 62 73 5f | 29 29 2c 20 52 45 54 55 | 1, obs_|)), RETU|
|00000930| 52 4e 20 22 55 6e 6b 6e | 6f 77 6e 20 6c 69 6e 6b |RN "Unkn|own link|
|00000940| 20 66 75 6e 63 74 69 6f | 6e 22 2c 20 52 45 54 55 | functio|n", RETU|
|00000950| 52 4e 20 22 55 6e 6b 6e | 6f 77 6e 20 6c 69 6e 6b |RN "Unkn|own link|
|00000960| 20 66 75 6e 63 74 69 6f | 6e 22 29 29 29 2c 20 6c | functio|n"))), l|
|00000970| 6f 67 6c 69 6b 20 3a 3d | 20 a4 28 56 45 43 54 4f |oglik :=| .(VECTO|
|00000980| 52 28 28 79 5f 99 69 5f | b7 4c 4e 28 70 5f 99 69 |R((y_.i_|.LN(p_.i|
|00000990| 5f 29 20 2b 20 41 42 53 | 28 79 5f 99 69 5f 20 2d |_) + ABS|(y_.i_ -|
|000009a0| 20 31 29 b7 4c 4e 28 31 | 20 2d 20 70 5f 99 69 5f | 1).LN(1| - p_.i_|
|000009b0| 29 29 b7 77 67 74 99 69 | 5f 2c 20 69 5f 2c 20 31 |)).wgt.i|_, i_, 1|
|000009c0| 2c 20 6f 62 73 5f 29 29 | 2c 20 77 5f 20 3a 3d 20 |, obs_))|, w_ := |
|000009d0| 56 45 43 54 4f 52 28 77 | 67 74 99 69 5f 2f 28 70 |VECTOR(w|gt.i_/(p|
|000009e0| 5f 99 69 5f b7 28 31 20 | 2d 20 70 5f 99 69 5f 29 |_.i_.(1 |- p_.i_)|
|000009f0| 29 2c 20 69 5f 2c 20 31 | 2c 20 6f 62 73 5f 29 2c |), i_, 1|, obs_),|
|00000a00| 20 78 78 20 3a 3d 20 56 | 45 43 54 4f 52 28 78 78 | xx := V|ECTOR(xx|
|00000a10| 5f 99 69 5f b7 66 5f 99 | 69 5f 2c 20 69 5f 2c 20 |_.i_.f_.|i_, i_, |
|00000a20| 31 2c 20 6f 62 73 5f 29 | 2c 20 77 78 78 20 3a 3d |1, obs_)|, wxx :=|
|00000a30| 20 56 45 43 54 4f 52 28 | 78 78 99 69 5f b7 77 5f | VECTOR(|xx.i_.w_|
|00000a40| 99 69 5f 2c 20 69 5f 2c | 20 31 2c 20 6f 62 73 5f |.i_, i_,| 1, obs_|
|00000a50| 29 2c 20 78 70 78 69 20 | 3a 3d 20 31 2f 28 78 78 |), xpxi |:= 1/(xx|
|00000a60| 60 20 95 20 77 78 78 29 | 2c 20 77 31 20 3a 3d 20 |` . wxx)|, w1 := |
|00000a70| 56 45 43 54 4f 52 28 77 | 5f 99 69 5f b7 28 79 5f |VECTOR(w|_.i_.(y_|
|00000a80| 99 69 5f 20 2d 20 70 5f | 99 69 5f 29 2c 20 69 5f |.i_ - p_|.i_), i_|
|00000a90| 2c 20 31 2c 20 6f 62 73 | 5f 29 2c 20 df 5f 20 3a |, 1, obs|_), ._ :|
|00000aa0| 2b 20 78 70 78 69 20 95 | 20 78 78 60 20 95 20 77 |+ xpxi .| xx` . w|
|00000ab0| 31 2c 20 69 74 65 72 5f | 20 3a 2b 20 31 29 2c 20 |1, iter_| :+ 1), |
|00000ac0| 49 46 28 69 74 65 72 5f | 20 93 20 6e 6e 2c 20 6d |IF(iter_| . nn, m|
|00000ad0| 73 67 20 3a 3d 20 22 43 | 6f 6e 76 65 72 67 65 6e |sg := "C|onvergen|
|00000ae0| 63 65 20 63 72 69 74 65 | 72 69 61 20 6d 65 74 21 |ce crite|ria met!|
|00000af0| 22 2c 20 6d 73 67 20 3a | 3d 20 22 43 6f 6e 76 65 |", msg :|= "Conve|
|00000b00| 72 67 65 6e 63 65 20 63 | 72 69 74 65 72 69 61 20 |rgence c|riteria |
|00000b10| 6e 6f 74 20 6d 65 74 21 | 22 29 2c 20 70 30 20 3a |not met!|"), p0 :|
|00000b20| 3d 20 28 79 5f 20 95 20 | 77 67 74 29 2f a4 28 77 |= (y_ . |wgt)/.(w|
|00000b30| 67 74 29 2c 20 6c 6f 67 | 6c 69 6b 30 20 3a 3d 20 |gt), log|lik0 := |
|00000b40| 28 79 5f b7 4c 4e 28 70 | 30 29 20 2b 20 56 45 43 |(y_.LN(p|0) + VEC|
|00000b50| 54 4f 52 28 41 42 53 28 | 79 5f 99 69 5f 20 2d 20 |TOR(ABS(|y_.i_ - |
|00000b60| 31 29 2c 20 69 5f 2c 20 | 31 2c 20 6f 62 73 5f 29 |1), i_, |1, obs_)|
|00000b70| b7 4c 4e 28 31 20 2d 20 | 70 30 29 29 20 95 20 77 |.LN(1 - |p0)) . w|
|00000b80| 67 74 2c 20 63 68 69 73 | 71 20 3a 3d 20 32 b7 28 |gt, chis|q := 2.(|
|00000b90| 6c 6f 67 6c 69 6b 20 2d | 20 6c 6f 67 6c 69 6b 30 |loglik -| loglik0|
|00000ba0| 29 2c 20 70 72 6f 62 20 | 3a 3d 20 31 20 2d 20 43 |), prob |:= 1 - C|
|00000bb0| 48 49 5f 53 51 55 41 52 | 45 28 63 68 69 73 71 2c |HI_SQUAR|E(chisq,|
|00000bc0| 20 6b 5f 29 2c 20 6c 6c | 72 61 74 69 6f 20 3a 3d | k_), ll|ratio :=|
|00000bd0| 20 5b 22 43 68 69 5e 32 | 22 2c 20 22 44 46 22 2c | ["Chi^2|", "DF",|
|00000be0| 20 22 50 72 6f 62 61 62 | 69 6c 69 74 79 22 3b 20 | "Probab|ility"; |
|00000bf0| 63 68 69 73 71 2c 20 6b | 5f 2c 20 70 72 6f 62 5d |chisq, k|_, prob]|
|00000c00| 2c 20 73 74 64 65 72 72 | 20 3a 3d 20 56 45 43 54 |, stderr| := VECT|
|00000c10| 4f 52 28 8b 28 78 70 78 | 69 99 69 5f 99 69 5f 29 |OR(.(xpx|i.i_.i_)|
|00000c20| 2c 20 69 5f 2c 20 31 2c | 20 6b 5f 20 2b 20 31 29 |, i_, 1,| k_ + 1)|
|00000c30| 2c 20 77 61 6c 64 20 3a | 3d 20 56 45 43 54 4f 52 |, wald :|= VECTOR|
|00000c40| 28 28 df 5f 99 69 5f 2f | 73 74 64 65 72 72 99 69 |((._.i_/|stderr.i|
|00000c50| 5f 29 5e 32 2c 20 69 5f | 2c 20 31 2c 20 6b 5f 20 |_)^2, i_|, 1, k_ |
|00000c60| 2b 20 31 29 2c 20 70 77 | 20 3a 3d 20 56 45 43 54 |+ 1), pw| := VECT|
|00000c70| 4f 52 28 31 20 2d 20 43 | 48 49 5f 53 51 55 41 52 |OR(1 - C|HI_SQUAR|
|00000c80| 45 28 77 61 6c 64 99 69 | 5f 2c 20 31 29 2c 20 69 |E(wald.i|_, 1), i|
|00000c90| 5f 2c 20 31 2c 20 6b 5f | 20 2b 20 31 29 2c 20 74 |_, 1, k_| + 1), t|
|00000ca0| 69 74 6c 20 3a 3d 20 5b | 22 50 61 72 6d 22 3b 20 |itl := [|"Parm"; |
|00000cb0| 22 56 61 6c 75 65 22 3b | 20 22 53 74 64 45 72 72 |"Value";| "StdErr|
|00000cc0| 22 3b 20 22 57 43 68 69 | 5e 32 22 3b 20 22 50 72 |"; "WChi|^2"; "Pr|
|00000cd0| 6f 62 28 57 29 22 5d 2c | 20 6f 75 74 20 3a 3d 20 |ob(W)"],| out := |
|00000ce0| 41 50 50 45 4e 44 5f 43 | 4f 4c 55 4d 4e 53 28 74 |APPEND_C|OLUMNS(t|
|00000cf0| 69 74 6c 2c 20 41 50 50 | 45 4e 44 28 5b 41 44 4a |itl, APP|END([ADJ|
|00000d00| 4f 49 4e 28 22 49 6e 74 | 65 72 63 65 70 74 22 2c |OIN("Int|ercept",|
|00000d10| 20 69 76 61 72 73 29 5d | 2c 20 5b df 5f 5d 2c 20 | ivars)]|, [._], |
|00000d20| 5b 73 74 64 65 72 72 5d | 2c 20 5b 77 61 6c 64 5d |[stderr]|, [wald]|
|00000d30| 2c 20 5b 70 77 5d 29 29 | 60 2c 20 66 65 71 20 3a |, [pw]))|`, feq :|
|00000d40| 3d 20 df 5f 20 95 20 41 | 44 4a 4f 49 4e 28 31 2c |= ._ . A|DJOIN(1,|
|00000d50| 20 69 76 61 72 73 29 2c | 20 49 46 28 6c 69 6e 6b | ivars),| IF(link|
|00000d60| 20 3d 20 31 2c 20 66 75 | 6e 63 20 3a 3d 20 22 4c | = 1, fu|nc := "L|
|00000d70| 6f 67 69 74 20 6c 69 6e | 6b 20 66 75 6e 63 74 69 |ogit lin|k functi|
|00000d80| 6f 6e 22 2c 20 49 46 28 | 6c 69 6e 6b 20 3d 20 32 |on", IF(|link = 2|
|00000d90| 2c 20 66 75 6e 63 20 3a | 3d 20 22 50 72 6f 62 69 |, func :|= "Probi|
|00000da0| 74 20 6c 69 6e 6b 20 66 | 75 6e 63 74 69 6f 6e 22 |t link f|unction"|
|00000db0| 2c 20 66 75 6e 63 20 3a | 3d 20 22 43 6f 6d 70 6c |, func :|= "Compl|
|00000dc0| 65 6d 65 6e 74 61 72 79 | 20 4c 6f 67 2d 4c 6f 67 |ementary| Log-Log|
|00000dd0| 20 6c 69 6e 6b 20 66 75 | 6e 63 74 69 6f 6e 22 29 | link fu|nction")|
|00000de0| 29 2c 20 5b 22 4c 6f 67 | 69 73 74 69 63 20 52 65 |), ["Log|istic Re|
|00000df0| 67 72 65 73 73 69 6f 6e | 22 3b 20 22 20 22 3b 20 |gression|"; " "; |
|00000e00| 66 75 6e 63 3b 20 22 20 | 22 3b 20 6d 73 67 3b 20 |func; " |"; msg; |
|00000e10| 22 20 22 3b 20 6f 75 74 | 3b 20 22 20 22 3b 20 6c |" "; out|; " "; l|
|00000e20| 6c 72 61 74 69 6f 3b 20 | 22 20 22 3b 20 66 65 71 |lratio; |" "; feq|
|00000e30| 5d 29 0d 0a 6e 6f 72 6d | 70 64 66 28 7a 29 3a 3d |])..norm|pdf(z):=|
|00000e40| ea 5e 28 2d 20 7a 5e 32 | 2f 32 29 b7 8b 32 2f 28 |.^(- z^2|/2)..2/(|
|00000e50| 32 b7 8b b9 29 0d 0a 70 | 6f 73 69 74 69 6f 6e 28 |2...)..p|osition(|
|00000e60| 65 2c 20 75 29 3a 3d 50 | 52 4f 47 28 64 20 3a 3d |e, u):=P|ROG(d :=|
|00000e70| 20 44 49 4d 28 75 29 2c | 20 69 5f 20 3a 3d 20 31 | DIM(u),| i_ := 1|
|00000e80| 2c 20 4c 4f 4f 50 28 49 | 46 28 65 20 3d 20 75 99 |, LOOP(I|F(e = u.|
|00000e90| 69 5f 2c 20 52 45 54 55 | 52 4e 20 69 5f 29 2c 20 |i_, RETU|RN i_), |
|00000ea0| 49 46 28 69 5f 20 3e 20 | 64 2c 20 52 45 54 55 52 |IF(i_ > |d, RETUR|
|00000eb0| 4e 20 30 29 2c 20 69 5f | 20 3a 2b 20 31 29 29 0d |N 0), i_| :+ 1)).|
|00000ec0| 0a 76 61 72 69 28 76 29 | 3a 3d 50 52 4f 47 28 6c |.vari(v)|:=PROG(l|
|00000ed0| 5f 20 3a 3d 20 44 49 4d | 28 76 29 2c 20 77 5f 20 |_ := DIM|(v), w_ |
|00000ee0| 3a 3d 20 5b 5d 2c 20 63 | 5f 20 3a 3d 20 31 2c 20 |:= [], c|_ := 1, |
|00000ef0| 4c 4f 4f 50 28 49 46 28 | 63 5f 20 3e 20 6c 5f 2c |LOOP(IF(|c_ > l_,|
|00000f00| 20 52 45 54 55 52 4e 20 | 77 5f 29 2c 20 65 5f 20 | RETURN |w_), e_ |
|00000f10| 3a 3d 20 76 99 63 5f 2c | 20 77 5f 20 3a 3d 20 41 |:= v.c_,| w_ := A|
|00000f20| 50 50 45 4e 44 28 77 5f | 2c 20 56 41 52 49 41 42 |PPEND(w_|, VARIAB|
|00000f30| 4c 45 53 28 65 5f 29 29 | 2c 20 63 5f 20 3a 2b 20 |LES(e_))|, c_ :+ |
|00000f40| 31 29 29 0d 0a 63 63 3a | 3d 0d 0a 63 69 3a 3d 0d |1))..cc:|=..ci:=.|
|00000f50| 0a 64 61 74 61 3a 3d 0d | 0a 64 76 61 72 3a 3d 0d |.data:=.|.dvar:=.|
|00000f60| 0a 66 65 71 3a 3d 30 2e | 30 38 32 30 33 30 38 30 |.feq:=0.|08203080|
|00000f70| 32 b7 68 20 2b 20 30 2e | 30 35 36 37 37 31 33 31 |2.h + 0.|05677131|
|00000f80| 33 b7 73 20 2d 20 35 2e | 35 35 39 31 36 36 34 0d |3.s - 5.|5591664.|
|00000f90| 0a 66 72 65 71 3a 3d 0d | 0a 66 75 6e 63 3a 3d 22 |.freq:=.|.func:="|
|00000fa0| 4c 6f 67 69 74 20 6c 69 | 6e 6b 20 66 75 6e 63 74 |Logit li|nk funct|
|00000fb0| 69 6f 6e 22 0d 0a 69 6e | 67 6f 74 3a 3d 5b 65 2c |ion"..in|got:=[e,|
|00000fc0| 20 74 2c 20 68 2c 20 73 | 3b 20 30 2c 20 31 30 2c | t, h, s|; 0, 10,|
|00000fd0| 20 37 2c 20 31 3b 20 30 | 2c 20 31 37 2c 20 37 2c | 7, 1; 0|, 17, 7,|
|00000fe0| 20 31 2e 37 3b 20 30 2c | 20 37 2c 20 37 2c 20 32 | 1.7; 0,| 7, 7, 2|
|00000ff0| 2e 32 3b 20 30 2c 20 31 | 32 2c 20 37 2c 20 32 2e |.2; 0, 1|2, 7, 2.|
|00001000| 38 3b 20 30 2c 20 39 2c | 20 37 2c 20 34 3b 20 30 |8; 0, 9,| 7, 4; 0|
|00001010| 2c 20 33 31 2c 20 31 34 | 2c 20 31 3b 20 30 2c 20 |, 31, 14|, 1; 0, |
|00001020| 34 33 2c 20 31 34 2c 20 | 31 2e 37 3b 20 32 2c 20 |43, 14, |1.7; 2, |
|00001030| 33 33 2c 20 31 34 2c 20 | 32 2e 32 3b 20 30 2c 20 |33, 14, |2.2; 0, |
|00001040| 33 31 2c 20 31 34 2c 20 | 32 2e 38 3b 20 30 2c 20 |31, 14, |2.8; 0, |
|00001050| 31 39 2c 20 31 34 2c 20 | 34 3b 20 31 2c 20 35 36 |19, 14, |4; 1, 56|
|00001060| 2c 20 32 37 2c 20 31 3b | 20 34 2c 20 34 34 2c 20 |, 27, 1;| 4, 44, |
|00001070| 32 37 2c 20 31 2e 37 3b | 20 30 2c 20 32 31 2c 20 |27, 1.7;| 0, 21, |
|00001080| 32 37 2c 20 32 2e 32 3b | 20 31 2c 20 32 32 2c 20 |27, 2.2;| 1, 22, |
|00001090| 32 37 2c 20 32 2e 38 3b | 20 31 2c 20 31 36 2c 20 |27, 2.8;| 1, 16, |
|000010a0| 32 37 2c 20 34 3b 20 33 | 2c 20 31 33 2c 20 35 31 |27, 4; 3|, 13, 51|
|000010b0| 2c 20 31 3b 20 30 2c 20 | 31 2c 20 35 31 2c 20 31 |, 1; 0, |1, 51, 1|
|000010c0| 2e 37 3b 20 30 2c 20 31 | 2c 20 35 31 2c 20 32 2e |.7; 0, 1|, 51, 2.|
|000010d0| 32 3b 20 30 2c 20 31 2c | 20 35 31 2c 20 34 5d 0d |2; 0, 1,| 51, 4].|
|000010e0| 0a 69 76 61 72 73 3a 3d | 0d 0a 69 76 61 72 73 5f |.ivars:=|..ivars_|
|000010f0| 3a 3d 5b 68 2c 20 73 5d | 0d 0a 6c 69 6e 6b 3a 3d |:=[h, s]|..link:=|
|00001100| 0d 0a 6c 69 6e 6b 5f 3a | 3d 31 0d 0a 6c 6c 72 61 |..link_:|=1..llra|
|00001110| 74 69 6f 3a 3d 5b 22 43 | 68 69 5e 32 22 2c 20 22 |tio:=["C|hi^2", "|
|00001120| 44 46 22 2c 20 22 50 72 | 6f 62 61 62 69 6c 69 74 |DF", "Pr|obabilit|
|00001130| 79 22 3b 20 31 31 2e 36 | 34 32 38 32 30 2c 20 32 |y"; 11.6|42820, 2|
|00001140| 2c 20 30 2e 30 30 32 39 | 36 33 34 32 33 36 5d 0d |, 0.0029|634236].|
|00001150| 0a 6d 73 67 3a 3d 22 43 | 6f 6e 76 65 72 67 65 6e |.msg:="C|onvergen|
|00001160| 63 65 20 63 72 69 74 65 | 72 69 61 20 6d 65 74 21 |ce crite|ria met!|
|00001170| 22 0d 0a 6e 6e 3a 3d 0d | 0a 6f 75 74 3a 3d 5b 22 |"..nn:=.|.out:=["|
|00001180| 50 61 72 6d 22 2c 20 22 | 56 61 6c 75 65 22 2c 20 |Parm", "|Value", |
|00001190| 22 53 74 64 45 72 72 22 | 2c 20 22 57 43 68 69 5e |"StdErr"|, "WChi^|
|000011a0| 32 22 2c 20 22 50 72 6f | 62 28 57 29 22 3b 20 22 |2", "Pro|b(W)"; "|
|000011b0| 49 6e 74 65 72 63 65 70 | 74 22 2c 20 2d 35 2e 35 |Intercep|t", -5.5|
|000011c0| 35 39 31 36 36 34 2c 20 | 31 2e 31 31 39 36 39 34 |591664, |1.119694|
|000011d0| 36 2c 20 32 34 2e 36 35 | 30 31 38 31 2c 20 37 2e |6, 24.65|0181, 7.|
|000011e0| 30 32 39 36 39 34 31 b7 | 31 30 5e 28 2d 37 29 3b |0296941.|10^(-7);|
|000011f0| 20 68 2c 20 30 2e 30 38 | 32 30 33 30 38 30 32 2c | h, 0.08|2030802,|
|00001200| 20 30 2e 30 32 33 37 33 | 34 34 38 33 2c 20 31 31 | 0.02373|4483, 11|
|00001210| 2e 39 34 35 32 32 36 2c | 20 30 2e 30 30 30 35 34 |.945226,| 0.00054|
|00001220| 37 38 39 30 37 39 3b 20 | 73 2c 20 30 2e 30 35 36 |789079; |s, 0.056|
|00001230| 37 37 31 33 31 33 2c 20 | 30 2e 33 33 31 32 31 33 |771313, |0.331213|
|00001240| 31 34 2c 20 30 2e 30 32 | 39 33 37 39 33 38 39 2c |14, 0.02|9379389,|
|00001250| 20 30 2e 38 36 33 39 30 | 35 39 37 5d 0d 0a 76 61 | 0.86390|597]..va|
|00001260| 6c 3a 3d 0d 0a 78 32 3a | 3d 0d 0a 78 70 78 69 3a |l:=..x2:|=..xpxi:|
|00001270| 3d 5b 31 2e 32 35 33 37 | 31 36 32 2c 20 2d 30 2e |=[1.2537|162, -0.|
|00001280| 30 32 31 35 36 36 34 38 | 39 2c 20 2d 30 2e 32 38 |02156648|9, -0.28|
|00001290| 31 37 36 35 33 34 3b 20 | 2d 30 2e 30 32 31 35 36 |176534; |-0.02156|
|000012a0| 36 34 38 39 2c 20 30 2e | 30 30 30 35 36 33 33 32 |6489, 0.|00056332|
|000012b0| 35 36 38 2c 20 30 2e 30 | 30 32 36 32 34 32 36 36 |568, 0.0|02624266|
|000012c0| 31 3b 20 2d 30 2e 32 38 | 31 37 36 35 33 34 2c 20 |1; -0.28|176534, |
|000012d0| 30 2e 30 30 32 36 32 34 | 32 36 36 31 2c 20 30 2e |0.002624|2661, 0.|
|000012e0| 31 30 39 37 30 32 31 34 | 5d 0d 0a 12 00 ff ff 00 |10970214|].......|
|000012f0| 00 08 00 43 54 65 78 74 | 4f 62 6a 08 00 00 00 0c |...CText|Obj.....|
|00001300| 00 00 00 cd 03 00 00 70 | 00 00 00 00 ff 04 02 7b |.......p|.......{|
|00001310| 5c 72 74 66 31 5c 61 6e | 73 69 5c 64 65 66 66 30 |\rtf1\an|si\deff0|
|00001320| 5c 64 65 66 74 61 62 37 | 32 30 7b 5c 66 6f 6e 74 |\deftab7|20{\font|
|00001330| 74 62 6c 7b 5c 66 30 5c | 66 73 77 69 73 73 20 4d |tbl{\f0\|fswiss M|
|00001340| 53 20 53 61 6e 73 20 53 | 65 72 69 66 3b 7d 7b 5c |S Sans S|erif;}{\|
|00001350| 66 31 5c 66 64 65 63 6f | 72 5c 66 63 68 61 72 73 |f1\fdeco|r\fchars|
|00001360| 65 74 32 20 53 79 6d 62 | 6f 6c 3b 7d 7b 5c 66 32 |et2 Symb|ol;}{\f2|
|00001370| 5c 66 73 77 69 73 73 5c | 66 70 72 71 32 20 53 79 |\fswiss\|fprq2 Sy|
|00001380| 73 74 65 6d 3b 7d 7b 5c | 66 33 5c 66 73 77 69 73 |stem;}{\|f3\fswis|
|00001390| 73 5c 66 70 72 71 32 20 | 41 72 69 61 6c 3b 7d 7b |s\fprq2 |Arial;}{|
|000013a0| 5c 66 34 5c 66 73 77 69 | 73 73 5c 66 63 68 61 72 |\f4\fswi|ss\fchar|
|000013b0| 73 65 74 31 20 41 72 69 | 61 6c 3b 7d 7b 5c 66 35 |set1 Ari|al;}{\f5|
|000013c0| 5c 66 6d 6f 64 65 72 6e | 5c 66 63 68 61 72 73 65 |\fmodern|\fcharse|
|000013d0| 74 32 20 44 66 57 35 20 | 50 72 69 6e 74 65 72 3b |t2 DfW5 |Printer;|
|000013e0| 7d 7d 0d 0a 7b 5c 63 6f | 6c 6f 72 74 62 6c 5c 72 |}}..{\co|lortbl\r|
|000013f0| 65 64 30 5c 67 72 65 65 | 6e 30 5c 62 6c 75 65 30 |ed0\gree|n0\blue0|
|00001400| 3b 5c 72 65 64 32 35 35 | 5c 67 72 65 65 6e 30 5c |;\red255|\green0\|
|00001410| 62 6c 75 65 30 3b 7d 0d | 0a 5c 64 65 66 6c 61 6e |blue0;}.|.\deflan|
|00001420| 67 31 30 33 33 5c 70 61 | 72 64 5c 71 63 5c 70 6c |g1033\pa|rd\qc\pl|
|00001430| 61 69 6e 5c 66 33 5c 66 | 73 32 34 5c 63 66 31 20 |ain\f3\f|s24\cf1 |
|00001440| 4c 6f 67 69 73 74 69 63 | 20 52 65 67 72 65 73 73 |Logistic| Regress|
|00001450| 69 6f 6e 0d 0a 5c 70 61 | 72 20 46 6f 72 20 42 69 |ion..\pa|r For Bi|
|00001460| 6e 61 72 79 20 44 65 70 | 65 6e 64 65 6e 74 20 44 |nary Dep|endent D|
|00001470| 61 74 61 0d 0a 5c 70 61 | 72 20 5c 70 61 72 64 5c |ata..\pa|r \pard\|
|00001480| 71 72 5c 70 6c 61 69 6e | 5c 66 33 5c 66 73 32 30 |qr\plain|\f3\fs20|
|00001490| 5c 63 66 31 20 4d 61 63 | 44 6f 6e 61 6c 64 20 52 |\cf1 Mac|Donald R|
|000014a0| 2e 20 50 68 69 6c 6c 69 | 70 73 0d 0a 5c 70 61 72 |. Philli|ps..\par|
|000014b0| 20 50 68 69 6c 6c 69 70 | 73 4d 40 67 61 6f 2e 67 | Phillip|sM@gao.g|
|000014c0| 6f 76 0d 0a 5c 70 61 72 | 20 64 6f 6e 70 68 69 6c |ov..\par| donphil|
|000014d0| 6c 69 70 73 40 73 74 61 | 72 70 6f 77 65 72 2e 6e |lips@sta|rpower.n|
|000014e0| 65 74 0d 0a 5c 70 61 72 | 20 4e 6f 76 65 6d 62 65 |et..\par| Novembe|
|000014f0| 72 20 32 30 30 31 5c 70 | 6c 61 69 6e 5c 66 33 5c |r 2001\p|lain\f3\|
|00001500| 66 73 32 30 5c 63 66 30 | 20 0d 0a 5c 70 61 72 20 |fs20\cf0| ..\par |
|00001510| 7d 0d 0a ff ff 00 00 08 | 00 43 45 78 70 6e 4f 62 |}.......|.CExpnOb|
|00001520| 6a 38 00 00 00 7c 00 00 | 00 c0 03 00 00 88 00 00 |j8...|..|........|
|00001530| 00 00 00 00 00 00 00 00 | 00 f0 bf 01 00 00 00 01 |........|........|
|00001540| 00 00 00 63 5b 49 6e 70 | 75 74 4d 6f 64 65 3a 3d |...c[Inp|utMode:=|
|00001550| 57 6f 72 64 2c 50 72 65 | 63 69 73 69 6f 6e 3a 3d |Word,Pre|cision:=|
|00001560| 41 70 70 72 6f 78 69 6d | 61 74 65 2c 50 72 65 63 |Approxim|ate,Prec|
|00001570| 69 73 69 6f 6e 44 69 67 | 69 74 73 3a 3d 31 30 2c |isionDig|its:=10,|
|00001580| 4e 6f 74 61 74 69 6f 6e | 3a 3d 53 63 69 65 6e 74 |Notation|:=Scient|
|00001590| 69 66 69 63 2c 4e 6f 74 | 61 74 69 6f 6e 44 69 67 |ific,Not|ationDig|
|000015a0| 69 74 73 3a 3d 38 5d 03 | 80 38 00 00 00 94 00 00 |its:=8].|.8......|
|000015b0| 00 90 02 00 00 a0 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000015c0| 00 f0 bf 02 00 00 00 01 | 00 00 00 3c 5b 78 70 78 |........|...<[xpx|
|000015d0| 69 3a 3d 2c 66 75 6e 63 | 3a 3d 2c 6d 73 67 3a 3d |i:=,func|:=,msg:=|
|000015e0| 2c 6f 75 74 3a 3d 2c 6c | 6c 72 61 74 69 6f 3a 3d |,out:=,l|lratio:=|
|000015f0| 2c 66 65 71 3a 3d 2c 6c | 69 6e 6b 5f 3a 3d 2c 69 |,feq:=,l|ink_:=,i|
|00001600| 76 61 72 73 5f 3a 3d 5d | 01 80 08 00 00 00 ac 00 |vars_:=]|........|
|00001610| 00 00 cd 03 00 00 8c 02 | 00 00 00 ff 5b 10 7b 5c |........|....[.{\|
|00001620| 72 74 66 31 5c 61 6e 73 | 69 5c 64 65 66 66 30 5c |rtf1\ans|i\deff0\|
|00001630| 64 65 66 74 61 62 37 32 | 30 7b 5c 66 6f 6e 74 74 |deftab72|0{\fontt|
|00001640| 62 6c 7b 5c 66 30 5c 66 | 73 77 69 73 73 20 4d 53 |bl{\f0\f|swiss MS|
|00001650| 20 53 61 6e 73 20 53 65 | 72 69 66 3b 7d 7b 5c 66 | Sans Se|rif;}{\f|
|00001660| 31 5c 66 64 65 63 6f 72 | 5c 66 63 68 61 72 73 65 |1\fdecor|\fcharse|
|00001670| 74 32 20 53 79 6d 62 6f | 6c 3b 7d 7b 5c 66 32 5c |t2 Symbo|l;}{\f2\|
|00001680| 66 73 77 69 73 73 5c 66 | 70 72 71 32 20 53 79 73 |fswiss\f|prq2 Sys|
|00001690| 74 65 6d 3b 7d 7b 5c 66 | 33 5c 66 73 77 69 73 73 |tem;}{\f|3\fswiss|
|000016a0| 5c 66 70 72 71 32 20 41 | 72 69 61 6c 3b 7d 7b 5c |\fprq2 A|rial;}{\|
|000016b0| 66 34 5c 66 73 77 69 73 | 73 5c 66 63 68 61 72 73 |f4\fswis|s\fchars|
|000016c0| 65 74 31 20 41 72 69 61 | 6c 3b 7d 7b 5c 66 35 5c |et1 Aria|l;}{\f5\|
|000016d0| 66 6d 6f 64 65 72 6e 5c | 66 63 68 61 72 73 65 74 |fmodern\|fcharset|
|000016e0| 32 20 44 66 57 35 20 50 | 72 69 6e 74 65 72 3b 7d |2 DfW5 P|rinter;}|
|000016f0| 7d 0d 0a 7b 5c 63 6f 6c | 6f 72 74 62 6c 5c 72 65 |}..{\col|ortbl\re|
|00001700| 64 30 5c 67 72 65 65 6e | 30 5c 62 6c 75 65 30 3b |d0\green|0\blue0;|
|00001710| 5c 72 65 64 32 35 35 5c | 67 72 65 65 6e 30 5c 62 |\red255\|green0\b|
|00001720| 6c 75 65 30 3b 7d 0d 0a | 5c 64 65 66 6c 61 6e 67 |lue0;}..|\deflang|
|00001730| 31 30 33 33 5c 70 61 72 | 64 5c 70 6c 61 69 6e 5c |1033\par|d\plain\|
|00001740| 66 34 5c 66 73 32 30 20 | 4c 6f 67 69 73 74 69 63 |f4\fs20 |Logistic|
|00001750| 20 72 65 67 72 65 73 73 | 69 6f 6e 20 66 6f 72 20 | regress|ion for |
|00001760| 62 69 6e 61 72 79 20 64 | 65 70 65 6e 64 65 6e 74 |binary d|ependent|
|00001770| 20 64 61 74 61 20 69 73 | 20 6f 66 66 65 72 65 64 | data is| offered|
|00001780| 20 22 61 73 20 69 73 2e | 22 20 20 49 20 6d 61 6b | "as is.|"  I mak|
|00001790| 65 20 6e 6f 20 63 6c 61 | 69 6d 20 74 68 61 74 20 |e no cla|im that |
|000017a0| 69 74 20 69 73 20 62 75 | 67 20 66 72 65 65 2c 20 |it is bu|g free, |
|000017b0| 61 6c 74 68 6f 75 67 68 | 20 49 20 62 65 6c 69 65 |although| I belie|
|000017c0| 76 65 20 69 74 20 69 73 | 2e 20 20 49 66 20 79 6f |ve it is|.  If yo|
|000017d0| 75 20 66 69 6e 64 20 61 | 6e 79 20 62 75 67 73 2c |u find a|ny bugs,|
|000017e0| 20 70 6c 65 61 73 65 20 | 6c 65 74 20 6d 65 20 6b | please |let me k|
|000017f0| 6e 6f 77 20 73 6f 20 49 | 20 63 61 6e 20 63 6f 72 |now so I| can cor|
|00001800| 72 65 63 74 20 74 68 65 | 6d 2e 0d 0a 5c 70 61 72 |rect the|m...\par|
|00001810| 20 0d 0a 5c 70 61 72 20 | 54 68 65 20 64 65 70 65 | ..\par |The depe|
|00001820| 6e 64 65 6e 74 20 76 61 | 72 69 61 62 6c 65 20 66 |ndent va|riable f|
|00001830| 6f 72 20 65 61 63 68 20 | 6f 62 73 65 72 76 61 74 |or each |observat|
|00001840| 69 6f 6e 20 69 73 20 65 | 69 74 68 65 72 20 61 20 |ion is e|ither a |
|00001850| 31 20 6f 72 20 61 20 30 | 2e 20 20 41 20 31 20 73 |1 or a 0|.  A 1 s|
|00001860| 69 67 6e 69 66 69 65 73 | 20 74 68 61 74 20 61 6e |ignifies| that an|
|00001870| 20 65 76 65 6e 74 20 68 | 61 70 70 65 6e 65 64 20 | event h|appened |
|00001880| 67 69 76 65 6e 20 74 68 | 65 20 73 74 61 74 65 20 |given th|e state |
|00001890| 6f 66 20 74 68 65 20 69 | 6e 64 65 70 65 6e 64 65 |of the i|ndepende|
|000018a0| 6e 74 20 76 61 72 69 61 | 62 6c 65 73 3b 20 61 20 |nt varia|bles; a |
|000018b0| 30 20 73 69 67 6e 69 66 | 69 65 73 20 74 68 65 20 |0 signif|ies the |
|000018c0| 65 76 65 6e 74 20 64 69 | 64 20 6e 6f 74 20 68 61 |event di|d not ha|
|000018d0| 70 70 65 6e 2e 20 20 54 | 68 65 20 6d 61 69 6e 20 |ppen.  T|he main |
|000018e0| 70 72 6f 67 72 61 6d 2c | 20 5c 70 6c 61 69 6e 5c |program,| \plain\|
|000018f0| 66 34 5c 66 73 32 30 5c | 63 66 31 20 6c 6f 67 69 |f4\fs20\|cf1 logi|
|00001900| 73 74 69 63 28 64 76 61 | 72 2c 69 76 61 72 73 2c |stic(dva|r,ivars,|
|00001910| 64 61 74 61 2c 6c 69 6e | 6b 2c 66 72 65 71 2c 6e |data,lin|k,freq,n|
|00001920| 6e 2c 63 63 29 5c 70 6c | 61 69 6e 5c 66 34 5c 66 |n,cc)\pl|ain\f4\f|
|00001930| 73 32 30 5c 63 66 30 20 | 2c 20 75 73 65 73 20 6f |s20\cf0 |, uses o|
|00001940| 6e 65 20 6f 66 20 74 68 | 72 65 65 20 6c 69 6e 6b |ne of th|ree link|
|00001950| 20 66 75 6e 63 74 69 6f | 6e 73 3a 20 6c 6f 67 69 | functio|ns: logi|
|00001960| 74 2c 20 70 72 6f 62 69 | 74 20 28 6f 72 20 6e 6f |t, probi|t (or no|
|00001970| 72 6d 69 74 29 2c 20 61 | 6e 64 20 63 6f 6d 70 6c |rmit), a|nd compl|
|00001980| 65 6d 65 6e 74 61 72 79 | 20 6c 6f 67 2d 6c 6f 67 |ementary| log-log|
|00001990| 2e 20 20 49 66 20 5c 70 | 6c 61 69 6e 5c 66 34 5c |.  If \p|lain\f4\|
|000019a0| 66 73 32 30 5c 63 66 31 | 20 6c 69 6e 6b 5c 70 6c |fs20\cf1| link\pl|
|000019b0| 61 69 6e 5c 66 34 5c 66 | 73 32 30 5c 63 66 30 20 |ain\f4\f|s20\cf0 |
|000019c0| 20 69 73 20 67 69 76 65 | 6e 20 74 68 65 20 76 61 | is give|n the va|
|000019d0| 6c 75 65 20 6f 66 20 31 | 2c 20 74 68 65 20 6c 6f |lue of 1|, the lo|
|000019e0| 67 69 74 20 6c 69 6e 6b | 20 66 75 6e 63 74 69 6f |git link| functio|
|000019f0| 6e 20 69 73 20 75 73 65 | 64 3b 20 61 20 76 61 6c |n is use|d; a val|
|00001a00| 75 65 20 6f 66 20 32 2c | 20 74 68 65 20 70 72 6f |ue of 2,| the pro|
|00001a10| 62 69 74 20 6c 69 6e 6b | 20 66 75 6e 63 74 69 6f |bit link| functio|
|00001a20| 6e 3b 20 61 6e 64 20 61 | 20 76 61 6c 75 65 20 6f |n; and a| value o|
|00001a30| 66 20 33 2c 20 74 68 65 | 20 63 6f 6d 70 6c 65 6d |f 3, the| complem|
|00001a40| 65 6e 74 61 72 79 20 6c | 6f 67 2d 6c 6f 67 20 6c |entary l|og-log l|
|00001a50| 69 6e 6b 20 66 75 6e 63 | 74 69 6f 6e 2e 20 20 54 |ink func|tion.  T|
|00001a60| 68 65 20 64 65 66 61 75 | 6c 74 20 76 61 6c 75 65 |he defau|lt value|
|00001a70| 20 6f 66 20 5c 70 6c 61 | 69 6e 5c 66 34 5c 66 73 | of \pla|in\f4\fs|
|00001a80| 32 30 5c 63 66 31 20 6c | 69 6e 6b 5c 70 6c 61 69 |20\cf1 l|ink\plai|
|00001a90| 6e 5c 66 34 5c 66 73 32 | 30 5c 63 66 30 20 20 69 |n\f4\fs2|0\cf0  i|
|00001aa0| 73 20 31 2e 0d 0a 5c 70 | 61 72 20 0d 0a 5c 70 61 |s 1...\p|ar ..\pa|
|00001ab0| 72 20 54 68 65 20 64 61 | 74 61 20 6d 61 79 20 6f |r The da|ta may o|
|00001ac0| 72 20 6d 61 79 20 6e 6f | 74 20 75 73 65 20 61 20 |r may no|t use a |
|00001ad0| 66 72 65 71 75 65 6e 63 | 79 20 76 61 72 69 61 62 |frequenc|y variab|
|00001ae0| 6c 65 3b 20 74 68 65 20 | 64 65 66 61 75 6c 74 20 |le; the |default |
|00001af0| 76 61 6c 75 65 20 6f 66 | 20 5c 70 6c 61 69 6e 5c |value of| \plain\|
|00001b00| 66 34 5c 66 73 32 30 5c | 63 66 31 20 66 72 65 71 |f4\fs20\|cf1 freq|
|00001b10| 5c 70 6c 61 69 6e 5c 66 | 34 5c 66 73 32 30 5c 63 |\plain\f|4\fs20\c|
|00001b20| 66 30 20 20 69 73 20 5b | 5d 2c 20 73 69 67 6e 69 |f0  is [|], signi|
|00001b30| 66 79 69 6e 67 20 6e 6f | 20 66 72 65 71 75 65 6e |fying no| frequen|
|00001b40| 63 79 20 76 61 72 69 61 | 62 6c 65 2e 20 20 54 68 |cy varia|ble.  Th|
|00001b50| 65 72 65 20 61 72 65 20 | 74 77 6f 20 74 79 70 65 |ere are |two type|
|00001b60| 73 20 6f 66 20 66 72 65 | 71 75 65 6e 63 69 65 73 |s of fre|quencies|
|00001b70| 20 75 73 65 64 20 69 6e | 20 6c 6f 67 69 73 74 69 | used in| logisti|
|00001b80| 63 20 72 65 67 72 65 73 | 73 69 6f 6e 2c 20 61 20 |c regres|sion, a |
|00001b90| 73 74 72 61 69 67 68 74 | 20 66 72 65 71 75 65 6e |straight| frequen|
|00001ba0| 63 79 20 61 6e 64 20 63 | 6f 75 6e 74 20 64 61 74 |cy and c|ount dat|
|00001bb0| 61 20 66 72 6f 6d 20 61 | 20 62 69 6e 6f 6d 69 61 |a from a| binomia|
|00001bc0| 6c 20 65 78 70 65 72 69 | 6d 65 6e 74 2e 20 20 54 |l experi|ment.  T|
|00001bd0| 68 65 20 63 6f 75 6e 74 | 20 64 61 74 61 20 69 73 |he count| data is|
|00001be0| 20 69 6e 20 74 68 65 20 | 66 6f 72 6d 20 6f 66 20 | in the |form of |
|00001bf0| 65 2c 20 74 68 65 20 6e | 75 6d 62 65 72 20 6f 66 |e, the n|umber of|
|00001c00| 20 73 75 63 65 73 73 66 | 75 6c 20 65 76 65 6e 74 | sucessf|ul event|
|00001c10| 73 2c 20 61 6e 64 20 74 | 2c 20 74 68 65 20 6e 75 |s, and t|, the nu|
|00001c20| 6d 62 65 72 20 6f 66 20 | 74 72 69 61 6c 73 2e 20 |mber of |trials. |
|00001c30| 20 49 66 20 63 6f 75 6e | 74 20 64 61 74 61 20 69 | If coun|t data i|
|00001c40| 73 20 75 73 65 64 2c 20 | 74 68 65 72 65 20 69 73 |s used, |there is|
|00001c50| 20 6e 6f 20 6e 65 65 64 | 20 74 6f 20 68 61 76 65 | no need| to have|
|00001c60| 20 61 20 64 65 70 65 6e | 64 65 6e 74 20 76 61 72 | a depen|dent var|
|00001c70| 69 61 62 6c 65 3b 20 69 | 74 20 77 69 6c 6c 20 62 |iable; i|t will b|
|00001c80| 65 20 63 72 65 61 74 65 | 64 20 62 79 20 74 68 65 |e create|d by the|
|00001c90| 20 70 72 6f 67 72 61 6d | 2e 20 20 49 66 20 61 20 | program|.  If a |
|00001ca0| 73 74 72 61 69 67 68 74 | 20 66 72 65 71 75 65 6e |straight| frequen|
|00001cb0| 63 79 20 69 73 20 75 73 | 65 64 2c 20 65 6e 74 65 |cy is us|ed, ente|
|00001cc0| 72 20 69 74 20 61 73 20 | 5b 76 61 72 69 61 62 6c |r it as |[variabl|
|00001cd0| 65 5d 2c 20 77 68 65 72 | 65 20 76 61 72 69 61 62 |e], wher|e variab|
|00001ce0| 6c 65 20 69 73 20 74 68 | 65 20 6e 61 6d 65 20 6f |le is th|e name o|
|00001cf0| 66 20 74 68 65 20 66 72 | 65 71 75 65 6e 63 79 20 |f the fr|equency |
|00001d00| 76 61 72 69 61 62 6c 65 | 2e 20 20 49 66 20 63 6f |variable|.  If co|
|00001d10| 75 6e 74 20 64 61 74 61 | 20 69 73 20 75 73 65 64 |unt data| is used|
|00001d20| 2c 20 65 6e 74 65 72 20 | 69 74 20 61 73 20 5b 65 |, enter |it as [e|
|00001d30| 76 65 6e 74 2c 74 72 69 | 61 6c 73 5d 2c 20 77 68 |vent,tri|als], wh|
|00001d40| 65 72 65 20 65 76 65 6e | 74 20 69 73 20 74 68 65 |ere even|t is the|
|00001d50| 20 6e 61 6d 65 20 6f 66 | 20 74 68 65 20 65 76 65 | name of| the eve|
|00001d60| 6e 74 20 76 61 72 69 61 | 62 6c 65 20 61 6e 64 20 |nt varia|ble and |
|00001d70| 74 72 69 61 6c 73 20 69 | 73 20 74 68 65 20 6e 61 |trials i|s the na|
|00001d80| 6d 65 20 6f 66 20 74 68 | 65 20 74 72 69 61 6c 73 |me of th|e trials|
|00001d90| 20 76 61 72 69 61 62 6c | 65 2e 0d 0a 5c 70 61 72 | variabl|e...\par|
|00001da0| 20 0d 0a 5c 70 61 72 20 | 54 68 65 20 64 61 74 61 | ..\par |The data|
|00001db0| 20 73 68 6f 75 6c 64 20 | 62 65 20 65 6e 74 65 72 | should |be enter|
|00001dc0| 65 64 20 69 6e 20 6d 61 | 74 72 69 78 20 66 6f 72 |ed in ma|trix for|
|00001dd0| 6d 61 74 20 77 69 74 68 | 20 74 68 65 20 66 69 72 |mat with| the fir|
|00001de0| 73 74 20 72 6f 77 20 6f | 66 20 74 68 65 20 6d 61 |st row o|f the ma|
|00001df0| 74 72 69 78 20 63 6f 6e | 74 61 69 6e 69 6e 67 20 |trix con|taining |
|00001e00| 74 68 65 20 6e 61 6d 65 | 73 20 6f 66 20 74 68 65 |the name|s of the|
|00001e10| 20 76 61 72 69 61 62 6c | 65 73 2e 20 20 49 20 73 | variabl|es.  I s|
|00001e20| 75 67 67 65 73 74 20 6f | 6e 65 2d 6c 65 74 74 65 |uggest o|ne-lette|
|00001e30| 72 20 76 61 72 69 61 62 | 6c 65 73 20 73 6f 20 61 |r variab|les so a|
|00001e40| 73 20 6e 6f 74 20 74 6f | 20 63 6f 6e 66 6c 69 63 |s not to| conflic|
|00001e50| 74 20 77 69 74 68 20 61 | 6e 79 20 72 65 73 65 72 |t with a|ny reser|
|00001e60| 76 65 64 20 6f 72 20 66 | 75 6e 63 74 69 6f 6e 20 |ved or f|unction |
|00001e70| 6e 61 6d 65 73 2e 20 20 | 41 73 73 69 67 6e 20 74 |names.  |Assign t|
|00001e80| 68 65 20 6d 61 74 72 69 | 78 20 74 6f 20 61 20 76 |he matri|x to a v|
|00001e90| 61 72 69 61 62 6c 65 2e | 0d 0a 5c 70 61 72 20 0d |ariable.|..\par .|
|00001ea0| 0a 5c 70 61 72 20 54 6f | 20 72 75 6e 20 74 68 65 |.\par To| run the|
|00001eb0| 20 70 72 6f 67 72 61 6d | 2c 20 65 6e 74 65 72 20 | program|, enter |
|00001ec0| 5c 70 6c 61 69 6e 5c 66 | 34 5c 66 73 32 30 5c 63 |\plain\f|4\fs20\c|
|00001ed0| 66 31 20 64 76 61 72 5c | 70 6c 61 69 6e 5c 66 34 |f1 dvar\|plain\f4|
|00001ee0| 5c 66 73 32 30 5c 63 66 | 30 20 2c 20 74 68 65 20 |\fs20\cf|0 , the |
|00001ef0| 6e 61 6d 65 20 6f 66 20 | 74 68 65 20 64 65 70 65 |name of |the depe|
|00001f00| 6e 64 65 6e 74 20 76 61 | 72 69 61 62 6c 65 3b 20 |ndent va|riable; |
|00001f10| 69 66 20 63 6f 75 6e 74 | 20 64 61 74 61 20 69 73 |if count| data is|
|00001f20| 20 75 73 65 64 20 5c 70 | 6c 61 69 6e 5c 66 34 5c | used \p|lain\f4\|
|00001f30| 66 73 32 30 5c 63 66 31 | 20 64 76 61 72 5c 70 6c |fs20\cf1| dvar\pl|
|00001f40| 61 69 6e 5c 66 34 5c 66 | 73 32 30 5c 63 66 30 20 |ain\f4\f|s20\cf0 |
|00001f50| 20 69 73 20 6e 6f 74 20 | 75 73 65 64 2e 20 20 49 | is not |used.  I|
|00001f60| 6e 20 74 68 69 73 20 63 | 61 73 65 20 61 6e 79 20 |n this c|ase any |
|00001f70| 6e 61 6d 65 20 6d 61 79 | 20 62 65 20 75 73 65 64 |name may| be used|
|00001f80| 3b 20 49 20 6a 75 73 74 | 20 65 6e 74 65 72 20 61 |; I just| enter a|
|00001f90| 20 30 20 28 7a 65 72 6f | 29 2e 20 20 45 6e 74 65 | 0 (zero|).  Ente|
|00001fa0| 72 20 5c 70 6c 61 69 6e | 5c 66 34 5c 66 73 32 30 |r \plain|\f4\fs20|
|00001fb0| 5c 63 66 31 20 69 76 61 | 72 73 5c 70 6c 61 69 6e |\cf1 iva|rs\plain|
|00001fc0| 5c 66 34 5c 66 73 32 30 | 5c 63 66 30 20 2c 20 61 |\f4\fs20|\cf0 , a|
|00001fd0| 20 76 65 63 74 6f 72 20 | 6f 66 20 74 68 65 20 69 | vector |of the i|
|00001fe0| 6e 64 65 70 65 6e 64 65 | 6e 74 20 76 61 72 69 61 |ndepende|nt varia|
|00001ff0| 62 6c 65 73 3b 20 65 6e | 74 65 72 20 69 74 20 61 |bles; en|ter it a|
|00002000| 73 20 61 20 76 65 63 74 | 6f 72 20 65 76 65 6e 20 |s a vect|or even |
|00002010| 69 66 20 74 68 65 72 65 | 20 69 73 20 6f 6e 6c 79 |if there| is only|
|00002020| 20 6f 6e 65 20 76 61 72 | 69 61 62 6c 65 2e 20 20 | one var|iable.  |
|00002030| 45 6e 74 65 72 20 5c 70 | 6c 61 69 6e 5c 66 34 5c |Enter \p|lain\f4\|
|00002040| 66 73 32 30 5c 63 66 31 | 20 64 61 74 61 5c 70 6c |fs20\cf1| data\pl|
|00002050| 61 69 6e 5c 66 34 5c 66 | 73 32 30 5c 63 66 30 20 |ain\f4\f|s20\cf0 |
|00002060| 2c 20 74 68 65 20 6e 61 | 6d 65 20 6f 66 20 74 68 |, the na|me of th|
|00002070| 65 20 64 61 74 61 73 65 | 74 2e 20 20 45 6e 74 65 |e datase|t.  Ente|
|00002080| 72 20 5c 70 6c 61 69 6e | 5c 66 34 5c 66 73 32 30 |r \plain|\f4\fs20|
|00002090| 5c 63 66 31 20 6c 69 6e | 6b 5c 70 6c 61 69 6e 5c |\cf1 lin|k\plain\|
|000020a0| 66 34 5c 66 73 32 30 5c | 63 66 30 20 20 61 73 20 |f4\fs20\|cf0  as |
|000020b0| 32 20 69 66 20 70 72 6f | 62 69 74 20 69 73 20 75 |2 if pro|bit is u|
|000020c0| 73 65 64 20 6f 72 20 33 | 20 69 66 20 63 6f 6d 70 |sed or 3| if comp|
|000020d0| 6c 65 6d 65 6e 74 61 72 | 79 20 6c 6f 67 2d 6c 6f |lementar|y log-lo|
|000020e0| 67 20 69 73 20 75 73 65 | 64 3b 20 74 68 65 20 64 |g is use|d; the d|
|000020f0| 65 66 61 75 6c 74 20 69 | 73 20 31 2c 20 6c 6f 67 |efault i|s 1, log|
|00002100| 69 74 2e 20 20 45 6e 74 | 65 72 20 5c 70 6c 61 69 |it.  Ent|er \plai|
|00002110| 6e 5c 66 34 5c 66 73 32 | 30 5c 63 66 31 20 66 72 |n\f4\fs2|0\cf1 fr|
|00002120| 65 71 5c 70 6c 61 69 6e | 5c 66 34 5c 66 73 32 30 |eq\plain|\f4\fs20|
|00002130| 5c 63 66 30 20 20 61 73 | 20 61 20 76 65 63 74 6f |\cf0  as| a vecto|
|00002140| 72 20 6f 66 20 74 68 65 | 20 6e 61 6d 65 20 6f 66 |r of the| name of|
|00002150| 20 74 68 65 20 66 72 65 | 71 75 65 6e 63 79 20 76 | the fre|quency v|
|00002160| 61 72 69 61 62 6c 65 20 | 6f 72 20 65 76 65 6e 74 |ariable |or event|
|00002170| 73 2d 74 72 69 61 6c 73 | 20 76 61 72 69 61 62 6c |s-trials| variabl|
|00002180| 65 73 2e 20 20 5c 70 6c | 61 69 6e 5c 66 34 5c 66 |es.  \pl|ain\f4\f|
|00002190| 73 32 30 5c 63 66 31 20 | 6e 6e 5c 70 6c 61 69 6e |s20\cf1 |nn\plain|
|000021a0| 5c 66 34 5c 66 73 32 30 | 5c 63 66 30 20 20 61 6e |\f4\fs20|\cf0  an|
|000021b0| 64 20 5c 70 6c 61 69 6e | 5c 66 34 5c 66 73 32 30 |d \plain|\f4\fs20|
|000021c0| 5c 63 66 31 20 63 63 5c | 70 6c 61 69 6e 5c 66 34 |\cf1 cc\|plain\f4|
|000021d0| 5c 66 73 32 30 5c 63 66 | 30 20 20 63 6f 6e 74 72 |\fs20\cf|0  contr|
|000021e0| 6f 6c 20 74 68 65 20 63 | 6f 6e 76 65 72 67 65 6e |ol the c|onvergen|
|000021f0| 63 65 20 63 72 69 74 65 | 72 69 61 2e 20 20 5c 70 |ce crite|ria.  \p|
|00002200| 6c 61 69 6e 5c 66 34 5c | 66 73 32 30 5c 63 66 31 |lain\f4\|fs20\cf1|
|00002210| 20 6e 6e 5c 70 6c 61 69 | 6e 5c 66 34 5c 66 73 32 | nn\plai|n\f4\fs2|
|00002220| 30 5c 63 66 30 20 20 73 | 70 65 63 69 66 69 65 73 |0\cf0  s|pecifies|
|00002230| 20 74 68 65 20 6d 61 78 | 69 6d 75 6d 20 6e 75 6d | the max|imum num|
|00002240| 62 65 72 20 6f 66 20 69 | 74 65 72 61 74 69 6f 6e |ber of i|teration|
|00002250| 73 20 74 6f 20 70 65 72 | 66 6f 72 6d 2e 20 20 5c |s to per|form.  \|
|00002260| 70 6c 61 69 6e 5c 66 34 | 5c 66 73 32 30 5c 63 66 |plain\f4|\fs20\cf|
|00002270| 31 20 63 63 5c 70 6c 61 | 69 6e 5c 66 34 5c 66 73 |1 cc\pla|in\f4\fs|
|00002280| 32 30 5c 63 66 30 20 20 | 73 70 65 63 69 66 69 65 |20\cf0  |specifie|
|00002290| 73 20 74 68 65 20 6d 61 | 78 69 6d 75 6d 20 61 62 |s the ma|ximum ab|
|000022a0| 73 6f 6c 75 74 65 20 64 | 69 66 66 65 72 65 6e 63 |solute d|ifferenc|
|000022b0| 65 20 62 65 74 77 65 65 | 6e 20 73 75 63 65 73 73 |e betwee|n sucess|
|000022c0| 69 76 65 20 69 74 65 72 | 61 74 69 6f 6e 73 20 6f |ive iter|ations o|
|000022d0| 66 20 74 68 65 20 69 6e | 64 65 70 65 6e 64 65 6e |f the in|dependen|
|000022e0| 74 20 76 61 72 69 61 62 | 6c 65 73 2e 20 20 54 68 |t variab|les.  Th|
|000022f0| 61 74 20 69 73 2c 20 74 | 68 65 20 63 68 61 6e 67 |at is, t|he chang|
|00002300| 65 20 69 6e 20 65 61 63 | 68 20 69 6e 64 65 70 65 |e in eac|h indepe|
|00002310| 6e 64 65 6e 74 20 76 61 | 72 69 61 62 6c 65 20 6d |ndent va|riable m|
|00002320| 75 73 74 20 62 65 20 6c | 65 73 73 20 74 68 61 6e |ust be l|ess than|
|00002330| 20 6f 72 20 65 71 75 61 | 6c 20 74 6f 20 5c 70 6c | or equa|l to \pl|
|00002340| 61 69 6e 5c 66 34 5c 66 | 73 32 30 5c 63 66 31 20 |ain\f4\f|s20\cf1 |
|00002350| 63 63 5c 70 6c 61 69 6e | 5c 66 34 5c 66 73 32 30 |cc\plain|\f4\fs20|
|00002360| 5c 63 66 30 20 2e 0d 0a | 5c 70 61 72 20 0d 0a 5c |\cf0 ...|\par ..\|
|00002370| 70 61 72 20 5c 70 6c 61 | 69 6e 5c 66 34 5c 66 73 |par \pla|in\f4\fs|
|00002380| 32 30 5c 63 66 31 20 50 | 72 65 64 69 63 74 65 64 |20\cf1 P|redicted|
|00002390| 5f 56 61 6c 75 65 73 28 | 76 61 6c 2c 63 69 29 5c |_Values(|val,ci)\|
|000023a0| 70 6c 61 69 6e 5c 66 34 | 5c 66 73 32 30 5c 63 66 |plain\f4|\fs20\cf|
|000023b0| 30 20 20 63 6f 6d 70 75 | 74 65 73 20 74 68 65 20 |0  compu|tes the |
|000023c0| 6c 6f 67 69 74 2c 20 70 | 72 6f 62 69 74 2c 20 6f |logit, p|robit, o|
|000023d0| 72 20 63 6c 6f 67 6c 6f | 67 20 6f 66 20 74 68 65 |r cloglo|g of the|
|000023e0| 20 70 72 6f 62 61 62 69 | 6c 74 79 20 61 73 73 6f | probabi|lty asso|
|000023f0| 63 69 61 74 65 64 20 77 | 69 74 68 20 61 20 67 69 |ciated w|ith a gi|
|00002400| 76 65 6e 20 73 65 74 20 | 6f 66 20 69 6e 64 65 70 |ven set |of indep|
|00002410| 65 6e 64 65 6e 74 20 76 | 61 72 69 61 62 6c 65 73 |endent v|ariables|
|00002420| 2e 20 20 49 74 20 61 6c | 73 6f 20 63 6f 6d 70 75 |.  It al|so compu|
|00002430| 74 65 73 20 74 68 65 20 | 73 74 61 6e 64 61 72 64 |tes the |standard|
|00002440| 20 65 72 72 6f 72 2c 20 | 73 65 2c 20 6f 66 20 74 | error, |se, of t|
|00002450| 68 65 20 6c 6f 67 69 74 | 2c 20 65 74 63 2c 20 69 |he logit|, etc, i|
|00002460| 74 73 20 70 72 6f 62 61 | 62 69 6c 69 74 79 2c 20 |ts proba|bility, |
|00002470| 70 2c 20 61 6e 64 20 61 | 20 63 6f 6e 66 69 64 65 |p, and a| confide|
|00002480| 6e 63 65 20 69 6e 74 65 | 72 76 61 6c 20 6f 66 20 |nce inte|rval of |
|00002490| 74 68 65 20 61 6d 6f 75 | 6e 74 20 5c 70 6c 61 69 |the amou|nt \plai|
|000024a0| 6e 5c 66 34 5c 66 73 32 | 30 5c 63 66 31 20 63 69 |n\f4\fs2|0\cf1 ci|
|000024b0| 5c 70 6c 61 69 6e 5c 66 | 34 5c 66 73 32 30 5c 63 |\plain\f|4\fs20\c|
|000024c0| 66 30 20 20 61 72 6f 75 | 6e 64 20 70 2e 0d 0a 5c |f0  arou|nd p...\|
|000024d0| 70 61 72 20 0d 0a 5c 70 | 61 72 20 54 68 65 20 70 |par ..\p|ar The p|
|000024e0| 72 6f 67 72 61 6d 20 61 | 6c 6c 6f 77 73 20 73 69 |rogram a|llows si|
|000024f0| 6d 70 6c 65 20 74 72 61 | 6e 73 66 6f 72 6d 61 74 |mple tra|nsformat|
|00002500| 69 6f 6e 73 20 6f 66 20 | 74 68 65 20 69 6e 64 65 |ions of |the inde|
|00002510| 70 65 6e 64 65 6e 74 20 | 76 61 72 69 61 62 6c 65 |pendent |variable|
|00002520| 73 2e 20 20 46 6f 72 20 | 69 6e 73 74 61 6e 63 65 |s.  For |instance|
|00002530| 2c 20 61 20 72 65 67 72 | 65 73 73 69 6f 6e 20 6d |, a regr|ession m|
|00002540| 61 79 20 64 65 70 65 6e | 64 20 6f 6e 20 74 68 65 |ay depen|d on the|
|00002550| 20 6e 61 74 75 72 61 6c | 20 6c 6f 67 20 6f 66 20 | natural| log of |
|00002560| 6f 6e 65 20 6f 72 20 6d | 6f 72 65 20 6f 66 20 74 |one or m|ore of t|
|00002570| 68 65 20 69 6e 64 65 70 | 65 6e 64 65 6e 74 20 76 |he indep|endent v|
|00002580| 61 72 69 61 62 6c 65 73 | 2e 20 20 57 68 65 6e 20 |ariables|.  When |
|00002590| 65 6e 74 65 72 69 6e 67 | 20 74 68 65 20 76 65 63 |entering| the vec|
|000025a0| 74 6f 72 20 6f 66 20 5c | 70 6c 61 69 6e 5c 66 34 |tor of \|plain\f4|
|000025b0| 5c 66 73 32 30 5c 63 66 | 31 20 69 76 61 72 73 5c |\fs20\cf|1 ivars\|
|000025c0| 70 6c 61 69 6e 5c 66 34 | 5c 66 73 32 30 5c 63 66 |plain\f4|\fs20\cf|
|000025d0| 30 20 20 65 6e 74 65 72 | 20 74 68 65 20 76 61 72 |0  enter| the var|
|000025e0| 69 61 62 6c 65 20 6e 61 | 6d 65 28 73 29 20 61 73 |iable na|me(s) as|
|000025f0| 20 4c 4e 28 6e 61 6d 65 | 29 2e 0d 0a 5c 70 61 72 | LN(name|)...\par|
|00002600| 20 0d 0a 5c 70 61 72 20 | 41 6e 20 65 78 61 6d 70 | ..\par |An examp|
|00002610| 6c 65 20 69 73 20 67 69 | 76 65 6e 20 77 69 74 68 |le is gi|ven with|
|00002620| 20 61 20 64 61 74 61 73 | 65 74 20 75 73 65 64 20 | a datas|et used |
|00002630| 62 79 20 53 41 53 20 61 | 73 20 61 6e 20 65 78 61 |by SAS a|s an exa|
|00002640| 6d 70 6c 65 20 66 6f 72 | 20 69 74 73 20 4c 4f 47 |mple for| its LOG|
|00002650| 49 53 54 49 43 20 70 72 | 6f 63 65 64 75 72 65 2e |ISTIC pr|ocedure.|
|00002660| 5c 70 6c 61 69 6e 5c 66 | 35 5c 66 73 32 34 20 0d |\plain\f|5\fs24 .|
|00002670| 0a 5c 70 61 72 20 7d 0d | 0a 03 80 38 00 00 00 98 |.\par }.|...8....|
|00002680| 02 00 00 50 02 00 00 c8 | 02 00 00 00 04 55 73 65 |...P....|.....Use|
|00002690| 72 00 00 00 00 00 00 f0 | bf 03 00 00 00 01 00 00 |r.......|........|
|000026a0| 00 4b 49 4e 43 4f 4d 50 | 4c 45 54 45 5f 47 41 4d |.KINCOMP|LETE_GAM|
|000026b0| 4d 41 28 7a 2c 77 29 3a | 3d 31 2f 47 41 4d 4d 41 |MA(z,w):|=1/GAMMA|
|000026c0| 28 7a 29 2a 28 77 5e 7a | 2f 7a 2b 49 4e 54 28 28 |(z)*(w^z|/z+INT((|
|000026d0| 23 65 5e 28 2d 74 5f 29 | 2d 31 29 2a 74 5f 5e 28 |#e^(-t_)|-1)*t_^(|
|000026e0| 7a 2d 31 29 2c 74 5f 2c | 30 2c 77 29 29 03 80 38 |z-1),t_,|0,w))..8|
|000026f0| 00 00 00 d4 02 00 00 b8 | 01 00 00 f8 02 00 00 00 |........|........|
|00002700| 04 55 73 65 72 00 00 00 | 00 00 00 f0 bf 04 00 00 |.User...|........|
|00002710| 00 01 00 00 00 2c 43 48 | 49 5f 53 51 55 41 52 45 |.....,CH|I_SQUARE|
|00002720| 28 78 32 2c 76 29 3a 3d | 49 4e 43 4f 4d 50 4c 45 |(x2,v):=|INCOMPLE|
|00002730| 54 45 5f 47 41 4d 4d 41 | 28 76 2f 32 2c 78 32 2f |TE_GAMMA|(v/2,x2/|
|00002740| 32 29 03 80 38 00 00 00 | 04 03 00 00 70 01 00 00 |2)..8...|....p...|
|00002750| 10 03 00 00 00 04 55 73 | 65 72 00 00 00 00 00 00 |......Us|er......|
|00002760| f0 bf 05 00 00 00 01 00 | 00 00 23 41 50 50 45 4e |........|..#APPEN|
|00002770| 44 5f 43 4f 4c 55 4d 4e | 53 28 61 2c 62 29 3a 3d |D_COLUMN|S(a,b):=|
|00002780| 41 50 50 45 4e 44 28 61 | 60 2c 62 60 29 60 03 80 |APPEND(a|`,b`)`..|
|00002790| 38 00 00 00 1c 03 00 00 | c8 00 00 00 94 03 00 00 |8.......|........|
|000027a0| 00 00 00 00 00 00 00 00 | f0 bf 06 00 00 00 01 00 |........|........|
|000027b0| 00 00 5d 70 6f 73 69 74 | 69 6f 6e 28 65 2c 75 29 |..]posit|ion(e,u)|
|000027c0| 3a 3d 50 52 4f 47 28 64 | 3a 3d 44 49 4d 28 75 29 |:=PROG(d|:=DIM(u)|
|000027d0| 2c 69 5f 3a 3d 31 2c 4c | 4f 4f 50 28 49 46 28 65 |,i_:=1,L|OOP(IF(e|
|000027e0| 3d 75 20 53 55 42 20 69 | 5f 2c 52 45 54 55 52 4e |=u SUB i|_,RETURN|
|000027f0| 28 69 5f 29 29 2c 49 46 | 28 69 5f 3e 64 2c 52 45 |(i_)),IF|(i_>d,RE|
|00002800| 54 55 52 4e 28 30 29 29 | 2c 69 5f 3a 2b 31 29 29 |TURN(0))|,i_:+1))|
|00002810| 03 80 38 00 00 00 a0 03 | 00 00 08 01 00 00 dc 03 |..8.....|........|
|00002820| 00 00 00 00 00 00 00 00 | 00 00 f0 bf 07 00 00 00 |........|........|
|00002830| 01 00 00 00 2c 6e 6f 72 | 6d 70 64 66 28 7a 29 3a |....,nor|mpdf(z):|
|00002840| 3d 53 51 52 54 28 32 29 | 2a 23 65 5e 28 2d 7a 5e |=SQRT(2)|*#e^(-z^|
|00002850| 32 2f 32 29 2f 28 32 2a | 53 51 52 54 28 70 69 29 |2/2)/(2*|SQRT(pi)|
|00002860| 29 03 80 38 00 00 00 e8 | 03 00 00 58 03 00 00 08 |)..8....|...X....|
|00002870| 08 00 00 00 04 55 73 65 | 72 00 00 00 00 00 00 f0 |.....Use|r.......|
|00002880| bf 08 00 00 00 01 00 00 | 00 ff 76 0a 6c 6f 67 69 |........|..v.logi|
|00002890| 73 74 69 63 28 64 76 61 | 72 2c 69 76 61 72 73 2c |stic(dva|r,ivars,|
|000028a0| 64 61 74 61 2c 6c 69 6e | 6b 3a 3d 31 2c 66 72 65 |data,lin|k:=1,fre|
|000028b0| 71 3a 3d 5b 5d 2c 6e 6e | 3a 3d 33 30 2c 63 63 3a |q:=[],nn|:=30,cc:|
|000028c0| 3d 31 30 5e 28 2d 38 29 | 29 3a 3d 50 52 4f 47 28 |=10^(-8)|):=PROG(|
|000028d0| 76 61 72 73 3a 3d 64 61 | 74 61 20 53 55 42 20 31 |vars:=da|ta SUB 1|
|000028e0| 2c 6f 62 73 5f 3a 3d 44 | 49 4d 28 64 61 74 61 29 |,obs_:=D|IM(data)|
|000028f0| 2d 31 2c 69 76 61 72 73 | 5f 3a 3d 69 76 61 72 73 |-1,ivars|_:=ivars|
|00002900| 2c 6c 69 6e 6b 5f 3a 3d | 6c 69 6e 6b 2c 64 61 74 |,link_:=|link,dat|
|00002910| 61 31 3a 3d 64 61 74 61 | 20 53 55 42 20 5b 32 2c |a1:=data| SUB [2,|
|00002920| 2e 2e 2e 2c 6f 62 73 5f | 2b 31 5d 2c 6b 5f 3a 3d |...,obs_|+1],k_:=|
|00002930| 44 49 4d 28 69 76 61 72 | 73 29 2c 66 72 65 71 5f |DIM(ivar|s),freq_|
|00002940| 3a 3d 44 49 4d 28 66 72 | 65 71 29 2c 49 46 28 66 |:=DIM(fr|eq),IF(f|
|00002950| 72 65 71 5f 3d 30 2c 77 | 67 74 3a 3d 56 45 43 54 |req_=0,w|gt:=VECT|
|00002960| 4f 52 28 31 2c 69 5f 2c | 31 2c 6f 62 73 5f 29 2c |OR(1,i_,|1,obs_),|
|00002970| 49 46 28 66 72 65 71 5f | 3d 31 2c 50 52 4f 47 28 |IF(freq_|=1,PROG(|
|00002980| 70 6f 73 5f 3a 3d 70 6f | 73 69 74 69 6f 6e 28 66 |pos_:=po|sition(f|
|00002990| 72 65 71 20 53 55 42 20 | 31 2c 76 61 72 73 29 2c |req SUB |1,vars),|
|000029a0| 77 67 74 3a 3d 64 61 74 | 61 20 53 55 42 20 20 53 |wgt:=dat|a SUB  S|
|000029b0| 55 42 20 70 6f 73 5f 29 | 2c 49 46 28 66 72 65 71 |UB pos_)|,IF(freq|
|000029c0| 5f 3d 32 2c 50 52 4f 47 | 28 70 6f 73 5f 3a 3d 70 |_=2,PROG|(pos_:=p|
|000029d0| 6f 73 69 74 69 6f 6e 28 | 66 72 65 71 20 53 55 42 |osition(|freq SUB|
|000029e0| 20 31 2c 76 61 72 73 29 | 2c 65 76 74 3a 3d 64 61 | 1,vars)|,evt:=da|
|000029f0| 74 61 31 20 53 55 42 20 | 20 53 55 42 20 70 6f 73 |ta1 SUB | SUB pos|
|00002a00| 5f 2c 70 6f 73 5f 3a 3d | 70 6f 73 69 74 69 6f 6e |_,pos_:=|position|
|00002a10| 28 66 72 65 71 20 53 55 | 42 20 32 2c 76 61 72 73 |(freq SU|B 2,vars|
|00002a20| 29 2c 74 72 69 3a 3d 64 | 61 74 61 31 20 53 55 42 |),tri:=d|ata1 SUB|
|00002a30| 20 20 53 55 42 20 70 6f | 73 5f 2c 77 67 74 3a 3d |  SUB po|s_,wgt:=|
|00002a40| 41 50 50 45 4e 44 28 65 | 76 74 2c 74 72 69 2d 65 |APPEND(e|vt,tri-e|
|00002a50| 76 74 29 29 2c 52 45 54 | 55 52 4e 28 22 49 6e 76 |vt)),RET|URN("Inv|
|00002a60| 61 6c 69 64 20 66 72 65 | 71 75 65 6e 63 79 20 76 |alid fre|quency v|
|00002a70| 61 72 69 61 62 6c 65 22 | 29 29 29 29 2c 49 46 28 |ariable"|)))),IF(|
|00002a80| 66 72 65 71 5f 3d 30 20 | 4f 52 20 66 72 65 71 5f |freq_=0 |OR freq_|
|00002a90| 3d 31 2c 50 52 4f 47 28 | 70 6f 73 5f 3a 3d 70 6f |=1,PROG(|pos_:=po|
|00002aa0| 73 69 74 69 6f 6e 28 64 | 76 61 72 2c 76 61 72 73 |sition(d|var,vars|
|00002ab0| 29 2c 79 5f 3a 3d 64 61 | 74 61 31 20 53 55 42 20 |),y_:=da|ta1 SUB |
|00002ac0| 20 53 55 42 20 70 6f 73 | 5f 29 2c 49 46 28 66 72 | SUB pos|_),IF(fr|
|00002ad0| 65 71 5f 3d 32 2c 79 5f | 3a 3d 41 50 50 45 4e 44 |eq_=2,y_|:=APPEND|
|00002ae0| 28 56 45 43 54 4f 52 28 | 31 2c 69 5f 2c 31 2c 6f |(VECTOR(|1,i_,1,o|
|00002af0| 62 73 5f 29 2c 56 45 43 | 54 4f 52 28 30 2c 69 5f |bs_),VEC|TOR(0,i_|
|00002b00| 2c 31 2c 6f 62 73 5f 29 | 29 29 29 2c 70 6f 73 5f |,1,obs_)|))),pos_|
|00002b10| 3a 3d 56 45 43 54 4f 52 | 28 70 6f 73 69 74 69 6f |:=VECTOR|(positio|
|00002b20| 6e 28 69 76 61 72 73 5f | 20 53 55 42 20 69 5f 2c |n(ivars_| SUB i_,|
|00002b30| 76 61 72 73 29 2c 69 5f | 2c 31 2c 6b 5f 29 2c 49 |vars),i_|,1,k_),I|
|00002b40| 46 28 66 72 65 71 5f 3d | 32 2c 50 52 4f 47 28 78 |F(freq_=|2,PROG(x|
|00002b50| 5f 3a 3d 41 50 50 45 4e | 44 28 64 61 74 61 31 20 |_:=APPEN|D(data1 |
|00002b60| 53 55 42 20 20 53 55 42 | 20 70 6f 73 5f 2c 64 61 |SUB  SUB| pos_,da|
|00002b70| 74 61 31 20 53 55 42 20 | 20 53 55 42 20 70 6f 73 |ta1 SUB | SUB pos|
|00002b80| 5f 29 2c 6f 62 73 5f 3a | 3d 32 2a 6f 62 73 5f 29 |_),obs_:|=2*obs_)|
|00002b90| 2c 78 5f 3a 3d 64 61 74 | 61 31 20 53 55 42 20 20 |,x_:=dat|a1 SUB  |
|00002ba0| 53 55 42 20 70 6f 73 5f | 29 2c 78 5f 3a 3d 56 45 |SUB pos_|),x_:=VE|
|00002bb0| 43 54 4f 52 28 53 55 42 | 53 54 28 69 76 61 72 73 |CTOR(SUB|ST(ivars|
|00002bc0| 2c 69 76 61 72 73 5f 2c | 78 5f 20 53 55 42 20 69 |,ivars_,|x_ SUB i|
|00002bd0| 5f 29 2c 69 5f 2c 31 2c | 6f 62 73 5f 29 2c 78 31 |_),i_,1,|obs_),x1|
|00002be0| 3a 3d 56 45 43 54 4f 52 | 28 56 45 43 54 4f 52 28 |:=VECTOR|(VECTOR(|
|00002bf0| 31 2c 69 5f 2c 31 2c 6f | 62 73 5f 29 2c 6a 2c 31 |1,i_,1,o|bs_),j,1|
|00002c00| 2c 31 29 60 2c 78 78 5f | 3a 3d 41 50 50 45 4e 44 |,1)`,xx_|:=APPEND|
|00002c10| 5f 43 4f 4c 55 4d 4e 53 | 28 78 31 2c 78 5f 29 2c |_COLUMNS|(x1,x_),|
|00002c20| 62 65 74 61 5f 3a 3d 56 | 45 43 54 4f 52 28 30 2c |beta_:=V|ECTOR(0,|
|00002c30| 69 5f 2c 31 2c 6b 5f 2b | 31 29 2c 6f 6c 64 df 3a |i_,1,k_+|1),old.:|
|00002c40| 3d 56 45 43 54 4f 52 28 | 31 2c 69 5f 2c 31 2c 6b |=VECTOR(|1,i_,1,k|
|00002c50| 5f 2b 31 29 2c 69 74 65 | 72 5f 3a 3d 30 2c 4c 4f |_+1),ite|r_:=0,LO|
|00002c60| 4f 50 28 63 63 5f 3a 3d | 4d 41 58 28 56 45 43 54 |OP(cc_:=|MAX(VECT|
|00002c70| 4f 52 28 41 42 53 28 62 | 65 74 61 5f 20 53 55 42 |OR(ABS(b|eta_ SUB|
|00002c80| 20 69 5f 2d 6f 6c 64 df | 20 53 55 42 20 69 5f 29 | i_-old.| SUB i_)|
|00002c90| 2c 69 5f 2c 31 2c 6b 5f | 2b 31 29 29 2c 49 46 28 |,i_,1,k_|+1)),IF(|
|00002ca0| 69 74 65 72 5f 3e 6e 6e | 20 4f 52 20 63 63 5f 3c |iter_>nn| OR cc_<|
|00002cb0| 63 63 2c 65 78 69 74 29 | 2c 6f 6c 64 df 3a 3d 62 |cc,exit)|,old.:=b|
|00002cc0| 65 74 61 5f 2c 7a 5f 3a | 3d 78 78 5f 20 2e 20 62 |eta_,z_:|=xx_ . b|
|00002cd0| 65 74 61 5f 2c 49 46 28 | 6c 69 6e 6b 3d 31 2c 50 |eta_,IF(|link=1,P|
|00002ce0| 52 4f 47 28 70 5f 3a 3d | 56 45 43 54 4f 52 28 31 |ROG(p_:=|VECTOR(1|
|00002cf0| 2f 28 31 2b 23 65 5e 28 | 2d 7a 5f 20 53 55 42 20 |/(1+#e^(|-z_ SUB |
|00002d00| 69 5f 29 29 2c 69 5f 2c | 31 2c 6f 62 73 5f 29 2c |i_)),i_,|1,obs_),|
|00002d10| 66 5f 3a 3d 56 45 43 54 | 4f 52 28 70 5f 20 53 55 |f_:=VECT|OR(p_ SU|
|00002d20| 42 20 69 5f 5e 32 2a 23 | 65 5e 28 2d 7a 5f 20 53 |B i_^2*#|e^(-z_ S|
|00002d30| 55 42 20 69 5f 29 2c 69 | 5f 2c 31 2c 6f 62 73 5f |UB i_),i|_,1,obs_|
|00002d40| 29 29 2c 49 46 28 6c 69 | 6e 6b 3d 32 2c 50 52 4f |)),IF(li|nk=2,PRO|
|00002d50| 47 28 70 5f 3a 3d 56 45 | 43 54 4f 52 28 4e 4f 52 |G(p_:=VE|CTOR(NOR|
|00002d60| 4d 41 4c 28 7a 5f 20 53 | 55 42 20 69 5f 29 2c 69 |MAL(z_ S|UB i_),i|
|00002d70| 5f 2c 31 2c 6f 62 73 5f | 29 2c 66 5f 3a 3d 56 45 |_,1,obs_|),f_:=VE|
|00002d80| 43 54 4f 52 28 6e 6f 72 | 6d 70 64 66 28 7a 5f 20 |CTOR(nor|mpdf(z_ |
|00002d90| 53 55 42 20 69 5f 29 2c | 69 5f 2c 31 2c 6f 62 73 |SUB i_),|i_,1,obs|
|00002da0| 5f 29 29 2c 49 46 28 6c | 69 6e 6b 3d 33 2c 50 52 |_)),IF(l|ink=3,PR|
|00002db0| 4f 47 28 70 5f 3a 3d 56 | 45 43 54 4f 52 28 31 2d |OG(p_:=V|ECTOR(1-|
|00002dc0| 23 65 5e 28 2d 23 65 5e | 7a 5f 20 53 55 42 20 69 |#e^(-#e^|z_ SUB i|
|00002dd0| 5f 29 2c 69 5f 2c 31 2c | 6f 62 73 5f 29 2c 66 5f |_),i_,1,|obs_),f_|
|00002de0| 3a 3d 56 45 43 54 4f 52 | 28 23 65 5e 28 7a 5f 20 |:=VECTOR|(#e^(z_ |
|00002df0| 53 55 42 20 69 5f 2d 23 | 65 5e 7a 5f 20 53 55 42 |SUB i_-#|e^z_ SUB|
|00002e00| 20 69 5f 29 2c 69 5f 2c | 31 2c 6f 62 73 5f 29 29 | i_),i_,|1,obs_))|
|00002e10| 2c 52 45 54 55 52 4e 28 | 22 55 6e 6b 6e 6f 77 6e |,RETURN(|"Unknown|
|00002e20| 20 6c 69 6e 6b 20 66 75 | 6e 63 74 69 6f 6e 22 29 | link fu|nction")|
|00002e30| 2c 52 45 54 55 52 4e 28 | 22 55 6e 6b 6e 6f 77 6e |,RETURN(|"Unknown|
|00002e40| 20 6c 69 6e 6b 20 66 75 | 6e 63 74 69 6f 6e 22 29 | link fu|nction")|
|00002e50| 29 29 29 2c 6c 6f 67 6c | 69 6b 3a 3d 53 55 4d 28 |))),logl|ik:=SUM(|
|00002e60| 56 45 43 54 4f 52 28 28 | 79 5f 20 53 55 42 20 69 |VECTOR((|y_ SUB i|
|00002e70| 5f 2a 4c 4e 28 70 5f 20 | 53 55 42 20 69 5f 29 2b |_*LN(p_ |SUB i_)+|
|00002e80| 41 42 53 28 79 5f 20 53 | 55 42 20 69 5f 2d 31 29 |ABS(y_ S|UB i_-1)|
|00002e90| 2a 4c 4e 28 31 2d 70 5f | 20 53 55 42 20 69 5f 29 |*LN(1-p_| SUB i_)|
|00002ea0| 29 2a 77 67 74 20 53 55 | 42 20 69 5f 2c 69 5f 2c |)*wgt SU|B i_,i_,|
|00002eb0| 31 2c 6f 62 73 5f 29 29 | 2c 77 5f 3a 3d 56 45 43 |1,obs_))|,w_:=VEC|
|00002ec0| 54 4f 52 28 77 67 74 20 | 53 55 42 20 69 5f 2f 28 |TOR(wgt |SUB i_/(|
|00002ed0| 70 5f 20 53 55 42 20 69 | 5f 2a 28 31 2d 70 5f 20 |p_ SUB i|_*(1-p_ |
|00002ee0| 53 55 42 20 69 5f 29 29 | 2c 69 5f 2c 31 2c 6f 62 |SUB i_))|,i_,1,ob|
|00002ef0| 73 5f 29 2c 78 78 3a 3d | 56 45 43 54 4f 52 28 78 |s_),xx:=|VECTOR(x|
|00002f00| 78 5f 20 53 55 42 20 69 | 5f 2a 66 5f 20 53 55 42 |x_ SUB i|_*f_ SUB|
|00002f10| 20 69 5f 2c 69 5f 2c 31 | 2c 6f 62 73 5f 29 2c 77 | i_,i_,1|,obs_),w|
|00002f20| 78 78 3a 3d 56 45 43 54 | 4f 52 28 78 78 20 53 55 |xx:=VECT|OR(xx SU|
|00002f30| 42 20 69 5f 2a 77 5f 20 | 53 55 42 20 69 5f 2c 69 |B i_*w_ |SUB i_,i|
|00002f40| 5f 2c 31 2c 6f 62 73 5f | 29 2c 78 70 78 69 3a 3d |_,1,obs_|),xpxi:=|
|00002f50| 28 78 78 60 20 2e 20 77 | 78 78 29 5e 28 2d 31 29 |(xx` . w|xx)^(-1)|
|00002f60| 2c 77 31 3a 3d 56 45 43 | 54 4f 52 28 77 5f 20 53 |,w1:=VEC|TOR(w_ S|
|00002f70| 55 42 20 69 5f 2a 28 79 | 5f 20 53 55 42 20 69 5f |UB i_*(y|_ SUB i_|
|00002f80| 2d 70 5f 20 53 55 42 20 | 69 5f 29 2c 69 5f 2c 31 |-p_ SUB |i_),i_,1|
|00002f90| 2c 6f 62 73 5f 29 2c 62 | 65 74 61 5f 3a 2b 78 70 |,obs_),b|eta_:+xp|
|00002fa0| 78 69 20 2e 20 78 78 60 | 20 2e 20 77 31 2c 69 74 |xi . xx`| . w1,it|
|00002fb0| 65 72 5f 3a 2b 31 29 2c | 49 46 28 69 74 65 72 5f |er_:+1),|IF(iter_|
|00002fc0| 3c 3d 6e 6e 2c 6d 73 67 | 3a 3d 22 43 6f 6e 76 65 |<=nn,msg|:="Conve|
|00002fd0| 72 67 65 6e 63 65 20 63 | 72 69 74 65 72 69 61 20 |rgence c|riteria |
|00002fe0| 6d 65 74 21 22 2c 6d 73 | 67 3a 3d 22 43 6f 6e 76 |met!",ms|g:="Conv|
|00002ff0| 65 72 67 65 6e 63 65 20 | 63 72 69 74 65 72 69 61 |ergence |criteria|
|00003000| 20 6e 6f 74 20 6d 65 74 | 21 22 29 2c 70 30 3a 3d | not met|!"),p0:=|
|00003010| 28 79 5f 20 2e 20 77 67 | 74 29 2f 53 55 4d 28 77 |(y_ . wg|t)/SUM(w|
|00003020| 67 74 29 2c 6c 6f 67 6c | 69 6b 30 3a 3d 28 79 5f |gt),logl|ik0:=(y_|
|00003030| 2a 4c 4e 28 70 30 29 2b | 56 45 43 54 4f 52 28 41 |*LN(p0)+|VECTOR(A|
|00003040| 42 53 28 79 5f 20 53 55 | 42 20 69 5f 2d 31 29 2c |BS(y_ SU|B i_-1),|
|00003050| 69 5f 2c 31 2c 6f 62 73 | 5f 29 2a 4c 4e 28 31 2d |i_,1,obs|_)*LN(1-|
|00003060| 70 30 29 29 20 2e 20 77 | 67 74 2c 63 68 69 73 71 |p0)) . w|gt,chisq|
|00003070| 3a 3d 32 2a 28 6c 6f 67 | 6c 69 6b 2d 6c 6f 67 6c |:=2*(log|lik-logl|
|00003080| 69 6b 30 29 2c 70 72 6f | 62 3a 3d 31 2d 43 48 49 |ik0),pro|b:=1-CHI|
|00003090| 5f 53 51 55 41 52 45 28 | 63 68 69 73 71 2c 6b 5f |_SQUARE(|chisq,k_|
|000030a0| 29 2c 6c 6c 72 61 74 69 | 6f 3a 3d 5b 5b 22 43 68 |),llrati|o:=[["Ch|
|000030b0| 69 5e 32 22 2c 22 44 46 | 22 2c 22 50 72 6f 62 61 |i^2","DF|","Proba|
|000030c0| 62 69 6c 69 74 79 22 5d | 2c 5b 63 68 69 73 71 2c |bility"]|,[chisq,|
|000030d0| 6b 5f 2c 70 72 6f 62 5d | 5d 2c 73 74 64 65 72 72 |k_,prob]|],stderr|
|000030e0| 3a 3d 56 45 43 54 4f 52 | 28 53 51 52 54 28 78 70 |:=VECTOR|(SQRT(xp|
|000030f0| 78 69 20 53 55 42 20 69 | 5f 20 53 55 42 20 69 5f |xi SUB i|_ SUB i_|
|00003100| 29 2c 69 5f 2c 31 2c 6b | 5f 2b 31 29 2c 77 61 6c |),i_,1,k|_+1),wal|
|00003110| 64 3a 3d 56 45 43 54 4f | 52 28 28 62 65 74 61 5f |d:=VECTO|R((beta_|
|00003120| 20 53 55 42 20 69 5f 2f | 73 74 64 65 72 72 20 53 | SUB i_/|stderr S|
|00003130| 55 42 20 69 5f 29 5e 32 | 2c 69 5f 2c 31 2c 6b 5f |UB i_)^2|,i_,1,k_|
|00003140| 2b 31 29 2c 70 77 3a 3d | 56 45 43 54 4f 52 28 31 |+1),pw:=|VECTOR(1|
|00003150| 2d 43 48 49 5f 53 51 55 | 41 52 45 28 77 61 6c 64 |-CHI_SQU|ARE(wald|
|00003160| 20 53 55 42 20 69 5f 2c | 31 29 2c 69 5f 2c 31 2c | SUB i_,|1),i_,1,|
|00003170| 6b 5f 2b 31 29 2c 74 69 | 74 6c 3a 3d 5b 5b 22 50 |k_+1),ti|tl:=[["P|
|00003180| 61 72 6d 22 5d 2c 5b 22 | 56 61 6c 75 65 22 5d 2c |arm"],["|Value"],|
|00003190| 5b 22 53 74 64 45 72 72 | 22 5d 2c 5b 22 57 43 68 |["StdErr|"],["WCh|
|000031a0| 69 5e 32 22 5d 2c 5b 22 | 50 72 6f 62 28 57 29 22 |i^2"],["|Prob(W)"|
|000031b0| 5d 5d 2c 6f 75 74 3a 3d | 41 50 50 45 4e 44 5f 43 |]],out:=|APPEND_C|
|000031c0| 4f 4c 55 4d 4e 53 28 74 | 69 74 6c 2c 41 50 50 45 |OLUMNS(t|itl,APPE|
|000031d0| 4e 44 28 5b 41 44 4a 4f | 49 4e 28 22 49 6e 74 65 |ND([ADJO|IN("Inte|
|000031e0| 72 63 65 70 74 22 2c 69 | 76 61 72 73 29 5d 2c 5b |rcept",i|vars)],[|
|000031f0| 62 65 74 61 5f 5d 2c 5b | 73 74 64 65 72 72 5d 2c |beta_],[|stderr],|
|00003200| 5b 77 61 6c 64 5d 2c 5b | 70 77 5d 29 29 60 2c 66 |[wald],[|pw]))`,f|
|00003210| 65 71 3a 3d 62 65 74 61 | 5f 20 2e 20 41 44 4a 4f |eq:=beta|_ . ADJO|
|00003220| 49 4e 28 31 2c 69 76 61 | 72 73 29 2c 49 46 28 6c |IN(1,iva|rs),IF(l|
|00003230| 69 6e 6b 3d 31 2c 66 75 | 6e 63 3a 3d 22 4c 6f 67 |ink=1,fu|nc:="Log|
|00003240| 69 74 20 6c 69 6e 6b 20 | 66 75 6e 63 74 69 6f 6e |it link |function|
|00003250| 22 2c 49 46 28 6c 69 6e | 6b 3d 32 2c 66 75 6e 63 |",IF(lin|k=2,func|
|00003260| 3a 3d 22 50 72 6f 62 69 | 74 20 6c 69 6e 6b 20 66 |:="Probi|t link f|
|00003270| 75 6e 63 74 69 6f 6e 22 | 2c 66 75 6e 63 3a 3d 22 |unction"|,func:="|
|00003280| 43 6f 6d 70 6c 65 6d 65 | 6e 74 61 72 79 20 4c 6f |Compleme|ntary Lo|
|00003290| 67 2d 4c 6f 67 20 6c 69 | 6e 6b 20 66 75 6e 63 74 |g-Log li|nk funct|
|000032a0| 69 6f 6e 22 29 29 2c 5b | 5b 22 4c 6f 67 69 73 74 |ion")),[|["Logist|
|000032b0| 69 63 20 52 65 67 72 65 | 73 73 69 6f 6e 22 5d 2c |ic Regre|ssion"],|
|000032c0| 5b 22 20 22 5d 2c 5b 66 | 75 6e 63 5d 2c 5b 22 20 |[" "],[f|unc],[" |
|000032d0| 22 5d 2c 5b 6d 73 67 5d | 2c 5b 22 20 22 5d 2c 5b |"],[msg]|,[" "],[|
|000032e0| 6f 75 74 5d 2c 5b 22 20 | 22 5d 2c 5b 6c 6c 72 61 |out],[" |"],[llra|
|000032f0| 74 69 6f 5d 2c 5b 22 20 | 22 5d 2c 5b 66 65 71 5d |tio],[" |"],[feq]|
|00003300| 5d 29 03 80 38 00 00 00 | 14 08 00 00 a0 02 00 00 |])..8...|........|
|00003310| 4c 09 00 00 00 04 55 73 | 65 72 00 00 00 00 00 00 |L.....Us|er......|
|00003320| f0 bf 09 00 00 00 01 00 | 00 00 ff 06 02 50 72 65 |........|.....Pre|
|00003330| 64 69 63 74 65 64 5f 56 | 61 6c 75 65 73 28 76 61 |dicted_V|alues(va|
|00003340| 6c 2c 63 69 3a 3d 30 2e | 39 35 29 3a 3d 50 52 4f |l,ci:=0.|95):=PRO|
|00003350| 47 28 6c 3a 3d 53 55 42 | 53 54 28 66 65 71 2c 69 |G(l:=SUB|ST(feq,i|
|00003360| 76 61 72 73 5f 2c 76 61 | 6c 29 2c 78 31 3a 3d 5b |vars_,va|l),x1:=[|
|00003370| 41 44 4a 4f 49 4e 28 31 | 2c 76 61 6c 29 5d 2c 73 |ADJOIN(1|,val)],s|
|00003380| 3a 3d 53 51 52 54 28 28 | 78 31 20 2e 20 78 70 78 |:=SQRT((|x1 . xpx|
|00003390| 69 20 2e 20 78 31 60 29 | 20 53 55 42 20 31 20 53 |i . x1`)| SUB 1 S|
|000033a0| 55 42 20 31 29 2c 69 6e | 76 6e 3a 3d 52 48 53 28 |UB 1),in|vn:=RHS(|
|000033b0| 4e 53 4f 4c 56 45 28 4e | 4f 52 4d 41 4c 28 7a 29 |NSOLVE(N|ORMAL(z)|
|000033c0| 2d 28 63 69 2f 32 2b 30 | 2e 35 29 2c 7a 2c 30 2c |-(ci/2+0|.5),z,0,|
|000033d0| 69 6e 66 29 29 2c 6c 6c | 3a 3d 6c 2d 69 6e 76 6e |inf)),ll|:=l-invn|
|000033e0| 2a 73 2c 75 6c 3a 3d 6c | 2b 69 6e 76 6e 2a 73 2c |*s,ul:=l|+invn*s,|
|000033f0| 49 46 28 6c 69 6e 6b 5f | 3d 31 2c 50 52 4f 47 28 |IF(link_|=1,PROG(|
|00003400| 6c 5f 3a 3d 22 6c 6f 67 | 69 74 28 70 72 6f 62 29 |l_:="log|it(prob)|
|00003410| 22 3d 6c 2c 70 3a 3d 31 | 2f 28 31 2b 23 65 5e 28 |"=l,p:=1|/(1+#e^(|
|00003420| 2d 6c 29 29 2c 6c 70 3a | 3d 31 2f 28 31 2b 23 65 |-l)),lp:|=1/(1+#e|
|00003430| 5e 28 2d 6c 6c 29 29 2c | 75 70 3a 3d 31 2f 28 31 |^(-ll)),|up:=1/(1|
|00003440| 2b 23 65 5e 28 2d 75 6c | 29 29 29 2c 49 46 28 6c |+#e^(-ul|))),IF(l|
|00003450| 69 6e 6b 5f 3d 32 2c 50 | 52 4f 47 28 6c 5f 3a 3d |ink_=2,P|ROG(l_:=|
|00003460| 22 70 72 6f 62 69 74 28 | 70 72 6f 62 29 22 3d 6c |"probit(|prob)"=l|
|00003470| 2c 70 3a 3d 4e 4f 52 4d | 41 4c 28 6c 29 2c 6c 70 |,p:=NORM|AL(l),lp|
|00003480| 3a 3d 4e 4f 52 4d 41 4c | 28 6c 6c 29 2c 75 70 3a |:=NORMAL|(ll),up:|
|00003490| 3d 4e 4f 52 4d 41 4c 28 | 75 6c 29 29 2c 50 52 4f |=NORMAL(|ul)),PRO|
|000034a0| 47 28 6c 5f 3a 3d 22 63 | 6c 6f 67 6c 6f 67 28 70 |G(l_:="c|loglog(p|
|000034b0| 72 6f 62 29 22 3d 6c 2c | 70 3a 3d 31 2d 23 65 5e |rob)"=l,|p:=1-#e^|
|000034c0| 28 2d 23 65 5e 6c 29 2c | 6c 70 3a 3d 31 2d 23 65 |(-#e^l),|lp:=1-#e|
|000034d0| 5e 28 2d 23 65 5e 6c 6c | 29 2c 75 70 3a 3d 31 2d |^(-#e^ll|),up:=1-|
|000034e0| 23 65 5e 28 2d 23 65 5e | 75 6c 29 29 29 29 2c 5b |#e^(-#e^|ul)))),[|
|000034f0| 5b 6c 5f 2c 22 73 65 22 | 3d 73 2c 22 70 72 6f 62 |[l_,"se"|=s,"prob|
|00003500| 22 3d 70 5d 2c 5b 41 50 | 50 45 4e 44 28 53 54 52 |"=p],[AP|PEND(STR|
|00003510| 49 4e 47 28 63 69 2a 31 | 30 30 29 2c 22 25 20 43 |ING(ci*1|00),"% C|
|00003520| 2e 49 2e 20 70 72 6f 62 | 22 29 2c 6c 70 2c 75 70 |.I. prob|"),lp,up|
|00003530| 5d 5d 29 03 80 38 00 00 | 00 58 09 00 00 10 01 00 |]])..8..|.X......|
|00003540| 00 2c 0b 00 00 00 04 55 | 73 65 72 00 00 00 00 00 |.,.....U|ser.....|
|00003550| 00 f0 bf 0a 00 00 00 01 | 00 00 00 ff 02 01 69 6e |........|......in|
|00003560| 67 6f 74 3a 3d 5b 5b 65 | 2c 74 2c 68 2c 73 5d 2c |got:=[[e|,t,h,s],|
|00003570| 5b 30 2c 31 30 2c 37 2c | 31 5d 2c 5b 30 2c 31 37 |[0,10,7,|1],[0,17|
|00003580| 2c 37 2c 31 2e 37 5d 2c | 5b 30 2c 37 2c 37 2c 32 |,7,1.7],|[0,7,7,2|
|00003590| 2e 32 5d 2c 5b 30 2c 31 | 32 2c 37 2c 32 2e 38 5d |.2],[0,1|2,7,2.8]|
|000035a0| 2c 5b 30 2c 39 2c 37 2c | 34 5d 2c 5b 30 2c 33 31 |,[0,9,7,|4],[0,31|
|000035b0| 2c 31 34 2c 31 5d 2c 5b | 30 2c 34 33 2c 31 34 2c |,14,1],[|0,43,14,|
|000035c0| 31 2e 37 5d 2c 5b 32 2c | 33 33 2c 31 34 2c 32 2e |1.7],[2,|33,14,2.|
|000035d0| 32 5d 2c 5b 30 2c 33 31 | 2c 31 34 2c 32 2e 38 5d |2],[0,31|,14,2.8]|
|000035e0| 2c 5b 30 2c 31 39 2c 31 | 34 2c 34 5d 2c 5b 31 2c |,[0,19,1|4,4],[1,|
|000035f0| 35 36 2c 32 37 2c 31 5d | 2c 5b 34 2c 34 34 2c 32 |56,27,1]|,[4,44,2|
|00003600| 37 2c 31 2e 37 5d 2c 5b | 30 2c 32 31 2c 32 37 2c |7,1.7],[|0,21,27,|
|00003610| 32 2e 32 5d 2c 5b 31 2c | 32 32 2c 32 37 2c 32 2e |2.2],[1,|22,27,2.|
|00003620| 38 5d 2c 5b 31 2c 31 36 | 2c 32 37 2c 34 5d 2c 5b |8],[1,16|,27,4],[|
|00003630| 33 2c 31 33 2c 35 31 2c | 31 5d 2c 5b 30 2c 31 2c |3,13,51,|1],[0,1,|
|00003640| 35 31 2c 31 2e 37 5d 2c | 5b 30 2c 31 2c 35 31 2c |51,1.7],|[0,1,51,|
|00003650| 32 2e 32 5d 2c 5b 30 2c | 31 2c 35 31 2c 34 5d 5d |2.2],[0,|1,51,4]]|
|00003660| 01 80 08 00 00 00 38 0b | 00 00 cd 03 00 00 58 0b |......8.|......X.|
|00003670| 00 00 00 ff 66 02 7b 5c | 72 74 66 31 5c 61 6e 73 |....f.{\|rtf1\ans|
|00003680| 69 5c 64 65 66 66 30 5c | 64 65 66 74 61 62 37 32 |i\deff0\|deftab72|
|00003690| 30 7b 5c 66 6f 6e 74 74 | 62 6c 7b 5c 66 30 5c 66 |0{\fontt|bl{\f0\f|
|000036a0| 73 77 69 73 73 20 4d 53 | 20 53 61 6e 73 20 53 65 |swiss MS| Sans Se|
|000036b0| 72 69 66 3b 7d 7b 5c 66 | 31 5c 66 64 65 63 6f 72 |rif;}{\f|1\fdecor|
|000036c0| 5c 66 63 68 61 72 73 65 | 74 32 20 53 79 6d 62 6f |\fcharse|t2 Symbo|
|000036d0| 6c 3b 7d 7b 5c 66 32 5c | 66 73 77 69 73 73 5c 66 |l;}{\f2\|fswiss\f|
|000036e0| 70 72 71 32 20 53 79 73 | 74 65 6d 3b 7d 7b 5c 66 |prq2 Sys|tem;}{\f|
|000036f0| 33 5c 66 73 77 69 73 73 | 5c 66 70 72 71 32 20 41 |3\fswiss|\fprq2 A|
|00003700| 72 69 61 6c 3b 7d 7b 5c | 66 34 5c 66 73 77 69 73 |rial;}{\|f4\fswis|
|00003710| 73 5c 66 63 68 61 72 73 | 65 74 31 20 41 72 69 61 |s\fchars|et1 Aria|
|00003720| 6c 3b 7d 7b 5c 66 35 5c | 66 6d 6f 64 65 72 6e 5c |l;}{\f5\|fmodern\|
|00003730| 66 63 68 61 72 73 65 74 | 32 20 44 66 57 35 20 50 |fcharset|2 DfW5 P|
|00003740| 72 69 6e 74 65 72 3b 7d | 7d 0d 0a 7b 5c 63 6f 6c |rinter;}|}..{\col|
|00003750| 6f 72 74 62 6c 5c 72 65 | 64 30 5c 67 72 65 65 6e |ortbl\re|d0\green|
|00003760| 30 5c 62 6c 75 65 30 3b | 5c 72 65 64 32 35 35 5c |0\blue0;|\red255\|
|00003770| 67 72 65 65 6e 30 5c 62 | 6c 75 65 30 3b 7d 0d 0a |green0\b|lue0;}..|
|00003780| 5c 64 65 66 6c 61 6e 67 | 31 30 33 33 5c 70 61 72 |\deflang|1033\par|
|00003790| 64 5c 70 6c 61 69 6e 5c | 66 33 5c 66 73 32 30 5c |d\plain\|f3\fs20\|
|000037a0| 63 66 30 20 49 6e 20 74 | 68 69 73 20 64 61 74 61 |cf0 In t|his data|
|000037b0| 73 65 74 2c 20 65 20 69 | 73 20 74 68 65 20 65 76 |set, e i|s the ev|
|000037c0| 65 6e 74 73 20 76 61 72 | 69 61 62 6c 65 20 28 74 |ents var|iable (t|
|000037d0| 68 65 20 6e 75 6d 62 65 | 72 20 6f 66 20 73 75 63 |he numbe|r of suc|
|000037e0| 65 73 73 66 75 6c 20 74 | 72 69 61 6c 73 29 2c 20 |essful t|rials), |
|000037f0| 74 20 74 68 65 20 74 72 | 69 61 6c 73 20 76 61 72 |t the tr|ials var|
|00003800| 69 61 62 6c 65 2c 20 61 | 6e 64 20 68 20 61 6e 64 |iable, a|nd h and|
|00003810| 20 73 20 61 72 65 20 69 | 6e 64 65 70 65 6e 64 65 | s are i|ndepende|
|00003820| 6e 74 20 76 61 72 69 61 | 62 6c 65 73 2e 20 20 53 |nt varia|bles.  S|
|00003830| 69 6e 63 65 20 74 68 65 | 20 64 61 74 61 73 65 74 |ince the| dataset|
|00003840| 20 72 65 70 72 65 73 65 | 6e 74 73 20 74 68 65 20 | represe|nts the |
|00003850| 72 65 73 75 6c 74 73 20 | 6f 66 20 61 20 62 69 6e |results |of a bin|
|00003860| 6f 6d 69 61 6c 20 65 78 | 70 65 72 69 6d 65 6e 74 |omial ex|periment|
|00003870| 2c 20 74 68 65 72 65 20 | 69 73 20 6e 6f 20 64 65 |, there |is no de|
|00003880| 70 65 6e 64 65 6e 74 20 | 76 61 72 69 61 62 6c 65 |pendent |variable|
|00003890| 3b 20 74 68 61 74 20 69 | 73 20 63 72 65 61 74 65 |; that i|s create|
|000038a0| 64 20 62 79 20 74 68 65 | 20 70 72 6f 67 72 61 6d |d by the| program|
|000038b0| 20 66 72 6f 6d 20 74 68 | 65 20 65 76 65 6e 74 73 | from th|e events|
|000038c0| 2f 74 72 69 61 6c 73 20 | 76 61 72 69 61 62 6c 65 |/trials |variable|
|000038d0| 73 2e 0d 0a 5c 70 61 72 | 20 7d 0d 0a 03 80 38 00 |s...\par| }....8.|
|000038e0| 00 00 64 0b 00 00 60 01 | 00 00 70 0b 00 00 00 04 |..d...`.|..p.....|
|000038f0| 55 73 65 72 00 00 00 00 | 00 00 f0 bf 0b 00 00 00 |User....|........|
|00003900| 01 00 00 00 1f 6c 6f 67 | 69 73 74 69 63 28 30 2c |.....log|istic(0,|
|00003910| 5b 68 2c 73 5d 2c 69 6e | 67 6f 74 2c 31 2c 5b 65 |[h,s],in|got,1,[e|
|00003920| 2c 74 5d 29 03 80 e0 00 | 00 00 7c 0b 00 00 20 03 |,t])....|..|... .|
|00003930| 00 00 e4 0c 00 00 01 09 | 53 69 6d 70 28 23 31 31 |........|Simp(#11|
|00003940| 29 fa 7e 6a bc 74 93 00 | 40 0c 00 00 00 01 00 00 |).~j.t..|@.......|
|00003950| 00 ff a7 01 5b 5b 22 4c | 6f 67 69 73 74 69 63 20 |....[["L|ogistic |
|00003960| 52 65 67 72 65 73 73 69 | 6f 6e 22 5d 2c 5b 22 20 |Regressi|on"],[" |
|00003970| 22 5d 2c 5b 22 4c 6f 67 | 69 74 20 6c 69 6e 6b 20 |"],["Log|it link |
|00003980| 66 75 6e 63 74 69 6f 6e | 22 5d 2c 5b 22 20 22 5d |function|"],[" "]|
|00003990| 2c 5b 22 43 6f 6e 76 65 | 72 67 65 6e 63 65 20 63 |,["Conve|rgence c|
|000039a0| 72 69 74 65 72 69 61 20 | 6d 65 74 21 22 5d 2c 5b |riteria |met!"],[|
|000039b0| 22 20 22 5d 2c 5b 5b 5b | 22 50 61 72 6d 22 2c 22 |" "],[[[|"Parm","|
|000039c0| 56 61 6c 75 65 22 2c 22 | 53 74 64 45 72 72 22 2c |Value","|StdErr",|
|000039d0| 22 57 43 68 69 5e 32 22 | 2c 22 50 72 6f 62 28 57 |"WChi^2"|,"Prob(W|
|000039e0| 29 22 5d 2c 5b 22 49 6e | 74 65 72 63 65 70 74 22 |)"],["In|tercept"|
|000039f0| 2c 2d 35 2e 35 35 39 31 | 36 36 34 2c 31 2e 31 31 |,-5.5591|664,1.11|
|00003a00| 39 36 39 34 36 2c 32 34 | 2e 36 35 30 31 38 31 2c |96946,24|.650181,|
|00003a10| 37 2e 30 32 39 36 39 34 | 31 2a 31 30 5e 28 2d 37 |7.029694|1*10^(-7|
|00003a20| 29 5d 2c 5b 68 2c 30 2e | 30 38 32 30 33 30 38 30 |)],[h,0.|08203080|
|00003a30| 32 2c 30 2e 30 32 33 37 | 33 34 34 38 33 2c 31 31 |2,0.0237|34483,11|
|00003a40| 2e 39 34 35 32 32 36 2c | 30 2e 30 30 30 35 34 37 |.945226,|0.000547|
|00003a50| 38 39 30 37 39 5d 2c 5b | 73 2c 30 2e 30 35 36 37 |89079],[|s,0.0567|
|00003a60| 37 31 33 31 33 2c 30 2e | 33 33 31 32 31 33 31 34 |71313,0.|33121314|
|00003a70| 2c 30 2e 30 32 39 33 37 | 39 33 38 39 2c 30 2e 38 |,0.02937|9389,0.8|
|00003a80| 36 33 39 30 35 39 37 5d | 5d 5d 2c 5b 22 20 22 5d |6390597]|]],[" "]|
|00003a90| 2c 5b 5b 5b 22 43 68 69 | 5e 32 22 2c 22 44 46 22 |,[[["Chi|^2","DF"|
|00003aa0| 2c 22 50 72 6f 62 61 62 | 69 6c 69 74 79 22 5d 2c |,"Probab|ility"],|
|00003ab0| 5b 31 31 2e 36 34 32 38 | 32 30 2c 32 2c 30 2e 30 |[11.6428|20,2,0.0|
|00003ac0| 30 32 39 36 33 34 32 33 | 36 5d 5d 5d 2c 5b 22 20 |02963423|6]]],[" |
|00003ad0| 22 5d 2c 5b 30 2e 30 38 | 32 30 33 30 38 30 32 2a |"],[0.08|2030802*|
|00003ae0| 68 2b 30 2e 30 35 36 37 | 37 31 33 31 33 2a 73 2d |h+0.0567|71313*s-|
|00003af0| 35 2e 35 35 39 31 36 36 | 34 5d 5d 01 80 08 00 00 |5.559166|4]].....|
|00003b00| 00 f0 0c 00 00 cd 03 00 | 00 50 0d 00 00 00 ff cd |........|.P......|
|00003b10| 03 7b 5c 72 74 66 31 5c | 61 6e 73 69 5c 64 65 66 |.{\rtf1\|ansi\def|
|00003b20| 66 30 5c 64 65 66 74 61 | 62 37 32 30 7b 5c 66 6f |f0\defta|b720{\fo|
|00003b30| 6e 74 74 62 6c 7b 5c 66 | 30 5c 66 73 77 69 73 73 |nttbl{\f|0\fswiss|
|00003b40| 20 4d 53 20 53 61 6e 73 | 20 53 65 72 69 66 3b 7d | MS Sans| Serif;}|
|00003b50| 7b 5c 66 31 5c 66 64 65 | 63 6f 72 5c 66 63 68 61 |{\f1\fde|cor\fcha|
|00003b60| 72 73 65 74 32 20 53 79 | 6d 62 6f 6c 3b 7d 7b 5c |rset2 Sy|mbol;}{\|
|00003b70| 66 32 5c 66 73 77 69 73 | 73 5c 66 70 72 71 32 20 |f2\fswis|s\fprq2 |
|00003b80| 53 79 73 74 65 6d 3b 7d | 7b 5c 66 33 5c 66 73 77 |System;}|{\f3\fsw|
|00003b90| 69 73 73 5c 66 70 72 71 | 32 20 41 72 69 61 6c 3b |iss\fprq|2 Arial;|
|00003ba0| 7d 7b 5c 66 34 5c 66 73 | 77 69 73 73 5c 66 63 68 |}{\f4\fs|wiss\fch|
|00003bb0| 61 72 73 65 74 31 20 41 | 72 69 61 6c 3b 7d 7b 5c |arset1 A|rial;}{\|
|00003bc0| 66 35 5c 66 6d 6f 64 65 | 72 6e 5c 66 63 68 61 72 |f5\fmode|rn\fchar|
|00003bd0| 73 65 74 32 20 44 66 57 | 35 20 50 72 69 6e 74 65 |set2 DfW|5 Printe|
|00003be0| 72 3b 7d 7d 0d 0a 7b 5c | 63 6f 6c 6f 72 74 62 6c |r;}}..{\|colortbl|
|00003bf0| 5c 72 65 64 30 5c 67 72 | 65 65 6e 30 5c 62 6c 75 |\red0\gr|een0\blu|
|00003c00| 65 30 3b 5c 72 65 64 32 | 35 35 5c 67 72 65 65 6e |e0;\red2|55\green|
|00003c10| 30 5c 62 6c 75 65 30 3b | 7d 0d 0a 5c 64 65 66 6c |0\blue0;|}..\defl|
|00003c20| 61 6e 67 31 30 33 33 5c | 70 61 72 64 5c 70 6c 61 |ang1033\|pard\pla|
|00003c30| 69 6e 5c 66 33 5c 66 73 | 32 30 5c 63 66 30 20 54 |in\f3\fs|20\cf0 T|
|00003c40| 68 65 20 6f 75 74 70 75 | 74 20 6d 61 74 72 69 78 |he outpu|t matrix|
|00003c50| 20 74 65 6c 6c 73 20 79 | 6f 75 20 77 68 61 74 20 | tells y|ou what |
|00003c60| 6c 69 6e 6b 20 66 75 6e | 63 74 69 6f 6e 20 77 61 |link fun|ction wa|
|00003c70| 73 20 75 73 65 64 20 61 | 6e 64 20 77 68 65 74 68 |s used a|nd wheth|
|00003c80| 65 72 20 74 68 65 20 63 | 6f 6e 76 65 72 67 65 6e |er the c|onvergen|
|00003c90| 63 65 20 63 72 69 74 65 | 72 69 61 20 77 61 73 20 |ce crite|ria was |
|00003ca0| 6d 65 74 20 6f 72 20 6e | 6f 74 2e 20 20 49 74 20 |met or n|ot.  It |
|00003cb0| 74 68 65 6e 20 67 69 76 | 65 73 20 74 68 65 20 65 |then giv|es the e|
|00003cc0| 73 74 69 6d 61 74 65 64 | 20 70 61 72 61 6d 65 74 |stimated| paramet|
|00003cd0| 65 72 20 76 61 6c 75 65 | 73 2c 20 74 68 65 69 72 |er value|s, their|
|00003ce0| 20 73 74 61 6e 64 61 72 | 64 20 65 72 72 6f 72 73 | standar|d errors|
|00003cf0| 2c 20 57 61 6c 64 20 63 | 68 69 2d 73 71 75 61 72 |, Wald c|hi-squar|
|00003d00| 65 20 73 74 61 74 69 73 | 74 69 63 2c 20 61 6e 64 |e statis|tic, and|
|00003d10| 20 74 68 65 20 70 72 6f | 62 61 62 69 6c 69 74 79 | the pro|bability|
|00003d20| 20 6f 66 20 74 68 65 20 | 57 61 6c 64 20 63 68 69 | of the |Wald chi|
|00003d30| 2d 73 71 75 61 72 65 2e | 20 20 49 74 20 74 68 65 |-square.|  It the|
|00003d40| 6e 20 67 69 76 65 73 20 | 74 68 65 20 2d 32 20 6c |n gives |the -2 l|
|00003d50| 6f 67 20 6c 69 6b 6c 69 | 68 6f 6f 64 20 73 74 61 |og likli|hood sta|
|00003d60| 74 69 73 74 69 63 20 66 | 6f 72 20 74 68 65 20 63 |tistic f|or the c|
|00003d70| 6f 76 61 72 69 61 74 65 | 73 20 28 69 6e 64 65 70 |ovariate|s (indep|
|00003d80| 65 6e 64 65 6e 74 20 76 | 61 72 69 61 62 6c 65 73 |endent v|ariables|
|00003d90| 29 2c 20 74 68 65 20 64 | 65 67 72 65 65 73 20 66 |), the d|egrees f|
|00003da0| 6f 72 20 66 72 65 65 64 | 6f 6d 2c 20 61 6e 64 20 |or freed|om, and |
|00003db0| 70 72 6f 62 61 62 69 6c | 69 74 79 2e 20 20 46 69 |probabil|ity.  Fi|
|00003dc0| 6e 61 6c 6c 79 2c 20 74 | 68 65 20 66 69 74 74 65 |nally, t|he fitte|
|00003dd0| 64 20 65 71 75 61 74 69 | 6f 6e 20 69 73 20 70 72 |d equati|on is pr|
|00003de0| 65 73 65 6e 74 65 64 2e | 0d 0a 5c 70 61 72 20 0d |esented.|..\par .|
|00003df0| 0a 5c 70 61 72 20 49 66 | 20 68 20 3d 20 37 20 61 |.\par If| h = 7 a|
|00003e00| 6e 64 20 73 20 3d 20 31 | 20 74 68 65 20 6c 6f 67 |nd s = 1| the log|
|00003e10| 69 74 28 70 72 6f 62 29 | 20 61 6e 64 20 70 72 6f |it(prob)| and pro|
|00003e20| 62 20 63 61 6e 20 62 65 | 20 65 73 74 69 6d 61 74 |b can be| estimat|
|00003e30| 65 64 20 77 69 74 68 20 | 61 20 39 35 20 70 65 72 |ed with |a 95 per|
|00003e40| 63 65 6e 74 20 63 6f 6e | 66 69 64 65 6e 63 65 20 |cent con|fidence |
|00003e50| 69 6e 74 65 72 76 61 6c | 2e 20 20 54 68 65 20 76 |interval|.  The v|
|00003e60| 61 6c 75 65 73 20 6f 66 | 20 74 68 65 20 76 61 72 |alues of| the var|
|00003e70| 69 61 62 6c 65 73 20 6d | 75 73 74 20 62 65 20 65 |iables m|ust be e|
|00003e80| 6e 74 65 72 65 64 20 69 | 6e 20 74 68 65 20 73 61 |ntered i|n the sa|
|00003e90| 6d 65 20 6f 72 64 65 72 | 20 61 73 20 74 68 65 79 |me order| as they|
|00003ea0| 20 61 72 65 20 69 6e 20 | 5c 70 6c 61 69 6e 5c 66 | are in |\plain\f|
|00003eb0| 33 5c 66 73 32 30 5c 63 | 66 31 20 69 76 61 72 73 |3\fs20\c|f1 ivars|
|00003ec0| 5c 70 6c 61 69 6e 5c 66 | 33 5c 66 73 32 30 5c 63 |\plain\f|3\fs20\c|
|00003ed0| 66 30 20 2e 0d 0a 5c 70 | 61 72 20 7d 0d 0a 03 80 |f0 ...\p|ar }....|
|00003ee0| 38 00 00 00 5c 0d 00 00 | f8 00 00 00 68 0d 00 00 |8...\...|....h...|
|00003ef0| 00 04 55 73 65 72 00 00 | 00 00 00 00 f0 bf 0d 00 |..User..|........|
|00003f00| 00 00 01 00 00 00 17 50 | 72 65 64 69 63 74 65 64 |.......P|redicted|
|00003f10| 5f 56 61 6c 75 65 73 28 | 5b 37 2c 31 5d 29 03 80 |_Values(|[7,1])..|
|00003f20| f8 00 00 00 74 0d 00 00 | 08 03 00 00 98 0d 00 00 |....t...|........|
|00003f30| 01 09 53 69 6d 70 28 23 | 31 33 29 8b 6c e7 fb a9 |..Simp(#|13).l...|
|00003f40| f1 a2 3f 0e 00 00 00 01 | 00 00 00 6b 5b 5b 22 6c |..?.....|...k[["l|
|00003f50| 6f 67 69 74 28 70 72 6f | 62 29 22 3d 2d 34 2e 39 |ogit(pro|b)"=-4.9|
|00003f60| 32 38 31 37 39 35 2c 22 | 73 65 22 3d 30 2e 37 34 |281795,"|se"=0.74|
|00003f70| 39 38 36 36 33 32 2c 22 | 70 72 6f 62 22 3d 30 2e |986632,"|prob"=0.|
|00003f80| 30 30 37 31 38 37 36 33 | 34 37 5d 2c 5b 22 39 35 |00718763|47],["95|
|00003f90| 25 20 43 2e 49 2e 20 70 | 72 6f 62 22 2c 30 2e 30 |% C.I. p|rob",0.0|
|00003fa0| 30 31 36 36 32 32 39 38 | 31 2c 30 2e 30 33 30 35 |01662298|1,0.0305|
|00003fb0| 31 37 33 31 30 5d 5d 0f | 00 00 00 01 00 00 00 01 |17310]].|........|
|00003fc0| 00 00 00 ff ff ff 00 00 | 00 00 00 00 f0 00 00 00 |........|........|
|00003fd0| 00 00 00 00 00 00 00 00 | 00 00 00                |........|...     |
+--------+-------------------------+-------------------------+--------+--------+