%!PS
%%Version: 3.3.1
%%DocumentFonts: (atend)
%%Pages: (atend)
%%EndComments
%
% Version 3.3.1 prologue for troff files.
%

/#copies 1 store
/aspectratio 1 def
/formsperpage 1 def
/landscape false def
/linewidth .3 def
/magnification 1 def
/margin 0 def
/orientation 0 def
/resolution 720 def
/rotation 1 def
/xoffset 0 def
/yoffset 0 def

/roundpage true def
/useclippath true def
/pagebbox [0 0 612 792] def

/R  /Times-Roman def
/I  /Times-Italic def
/B  /Times-Bold def
/BI /Times-BoldItalic def
/H  /Helvetica def
/HI /Helvetica-Oblique def
/HB /Helvetica-Bold def
/HX /Helvetica-BoldOblique def
/CW /Courier def
/CO /Courier def
/CI /Courier-Oblique def
/CB /Courier-Bold def
/CX /Courier-BoldOblique def
/PA /Palatino-Roman def
/PI /Palatino-Italic def
/PB /Palatino-Bold def
/PX /Palatino-BoldItalic def
/Hr /Helvetica-Narrow def
/Hi /Helvetica-Narrow-Oblique def
/Hb /Helvetica-Narrow-Bold def
/Hx /Helvetica-Narrow-BoldOblique def
/KR /Bookman-Light def
/KI /Bookman-LightItalic def
/KB /Bookman-Demi def
/KX /Bookman-DemiItalic def
/AR /AvantGarde-Book def
/AI /AvantGarde-BookOblique def
/AB /AvantGarde-Demi def
/AX /AvantGarde-DemiOblique def
/NR /NewCenturySchlbk-Roman def
/NI /NewCenturySchlbk-Italic def
/NB /NewCenturySchlbk-Bold def
/NX /NewCenturySchlbk-BoldItalic def
/ZD /ZapfDingbats def
/ZI /ZapfChancery-MediumItalic def
/S  /S def
/S1 /S1 def
/GR /Symbol def

/inch {72 mul} bind def
/min {2 copy gt {exch} if pop} bind def

/setup {
	counttomark 2 idiv {def} repeat pop

	landscape {/orientation 90 orientation add def} if
	/scaling 72 resolution div def
	linewidth setlinewidth
	1 setlinecap

	pagedimensions
	xcenter ycenter translate
	orientation rotation mul rotate
	width 2 div neg height 2 div translate
	xoffset inch yoffset inch neg translate
	margin 2 div dup neg translate
	magnification dup aspectratio mul scale
	scaling scaling scale

	addmetrics
	0 0 moveto
} def

/pagedimensions {
	useclippath userdict /gotpagebbox known not and {
		/pagebbox [clippath pathbbox newpath] def
		roundpage currentdict /roundpagebbox known and {roundpagebbox} if
	} if
	pagebbox aload pop
	4 -1 roll exch 4 1 roll 4 copy
	landscape {4 2 roll} if
	sub /width exch def
	sub /height exch def
	add 2 div /xcenter exch def
	add 2 div /ycenter exch def
	userdict /gotpagebbox true put
} def

/addmetrics {
	/Symbol /S null Sdefs cf
	/Times-Roman /S1 StandardEncoding dup length array copy S1defs cf
} def

/pagesetup {
	/page exch def
	currentdict /pagedict known currentdict page known and {
		page load pagedict exch get cvx exec
	} if
} def

/decodingdefs [
	{counttomark 2 idiv {y moveto show} repeat}
	{neg /y exch def counttomark 2 idiv {y moveto show} repeat}
	{neg moveto {2 index stringwidth pop sub exch div 0 32 4 -1 roll widthshow} repeat}
	{neg moveto {spacewidth sub 0.0 32 4 -1 roll widthshow} repeat}
	{counttomark 2 idiv {y moveto show} repeat}
	{neg setfunnytext}
] def

/setdecoding {/t decodingdefs 3 -1 roll get bind def} bind def

/w {neg moveto show} bind def
/m {neg dup /y exch def moveto} bind def
/done {/lastpage where {pop lastpage} if} def

/f {
	dup /font exch def findfont exch
	dup /ptsize exch def scaling div dup /size exch def scalefont setfont
	linewidth ptsize mul scaling 10 mul div setlinewidth
	/spacewidth ( ) stringwidth pop def
} bind def

/changefont {
	/fontheight exch def
	/fontslant exch def
	currentfont [
		1 0
		fontheight ptsize div fontslant sin mul fontslant cos div
		fontheight ptsize div
		0 0
	] makefont setfont
} bind def

/sf {f} bind def

/cf {
	dup length 2 idiv
	/entries exch def
	/chtab exch def
	/newencoding exch def
	/newfont exch def

	findfont dup length 1 add dict
	/newdict exch def
	{1 index /FID ne {newdict 3 1 roll put}{pop pop} ifelse} forall

	newencoding type /arraytype eq {newdict /Encoding newencoding put} if

	newdict /Metrics entries dict put
	newdict /Metrics get
	begin
		chtab aload pop
		1 1 entries {pop def} for
		newfont newdict definefont pop
	end
} bind def

%
% A few arrays used to adjust reference points and character widths in some
% of the printer resident fonts. If square roots are too high try changing
% the lines describing /radical and /radicalex to,
%
%	/radical	[0 -75 550 0]
%	/radicalex	[-50 -75 500 0]
%
% Move braceleftbt a bit - default PostScript character is off a bit.
%

/Sdefs [
	/bracketlefttp		[201 500]
	/bracketleftbt		[201 500]
	/bracketrighttp		[-81 380]
	/bracketrightbt		[-83 380]
	/braceleftbt		[203 490]
	/bracketrightex		[220 -125 500 0]
	/radical		[0 0 550 0]
	/radicalex		[-50 0 500 0]
	/parenleftex		[-20 -170 0 0]
	/integral		[100 -50 500 0]
	/infinity		[10 -75 730 0]
] def

/S1defs [
	/underscore		[0 80 500 0]
	/endash			[7 90 650 0]
] def
%
% Tries to round clipping path dimensions, as stored in array pagebbox, so they
% match one of the known sizes in the papersizes array. Lower left coordinates
% are always set to 0.
%

/roundpagebbox {
    7 dict begin
	/papersizes [8.5 inch 11 inch 14 inch 17 inch] def

	/mappapersize {
		/val exch def
		/slop .5 inch def
		/diff slop def
		/j 0 def
		0 1 papersizes length 1 sub {
			/i exch def
			papersizes i get val sub abs
			dup diff le {/diff exch def /j i def} {pop} ifelse
		} for
		diff slop lt {papersizes j get} {val} ifelse
	} def

	pagebbox 0 0 put
	pagebbox 1 0 put
	pagebbox dup 2 get mappapersize 2 exch put
	pagebbox dup 3 get mappapersize 3 exch put
    end
} bind def

%%EndProlog
%%BeginSetup
mark
/xoffset 0 def
/yoffset 0 def
/#copies 1 store
/magnification 1 def
%%FormsPerPage: 1
/formsperpage 1 def

/landscape false def
/resolution 720 def
setup
2 setdecoding
%%EndSetup
%%Page: 1 1
/saveobj save def
mark
1 pagesetup
12 B f
(Driver PMAIN for DGL[FG][B ])4 1684 1 2038 1220 t
10 I f
(David M. Gay)2 568 1 2596 1480 t
10 R f
( variant DGLF\) address robust, maximum-quasi-likelihood, and)6 2697(DGLG \(and it's finite-difference)3 1373 2 970 2056 t
(maximum-likelihood regression problems that involve nuisance parameters: for)7 3182 1 720 2211 t
10 I f
(x)3927 2211 w
10 S f
(=)4020 2211 w
(\351)4132 2174 w
(\353)4132 2274 w
(g)4189 2271 w
(b)4182 2171 w
(\371)4237 2174 w
(\373)4237 2274 w
10 R f
(,)4275 2211 w
10 I f
(f)2120 2496 w
10 R f
(\()2164 2496 w
10 I f
(x)2205 2496 w
10 R f
(\) :)1 110 1 2257 2496 t
10 S f
(=)2383 2496 w
10 I f
(f)2495 2496 w
10 R f
(\()2539 2496 w
10 S f
(b)2580 2496 w
10 R f
(,)2643 2496 w
10 S f
(g)2676 2496 w
10 R f
(\))2725 2496 w
10 S f
(=)2815 2496 w
7 I f
(i)2919 2596 w
7 S f
(=)2955 2596 w
7 R f
(1)3005 2596 w
15 S f
(S)2935 2526 w
7 I f
(n)2962 2396 w
10 S f
(r)3081 2496 w
7 I f
(i)3147 2516 w
10 R f
(\()3183 2496 w
10 S f
(h)3224 2496 w
7 I f
(i)3295 2516 w
10 R f
(\()3331 2496 w
10 S f
(b)3372 2496 w
10 R f
(\) ,)1 74 1 3435 2496 t
10 S f
(g)3517 2496 w
10 R f
(\) ,)1 74 1 3566 2496 t
(where)720 2756 w
10 S f
(b \316)1 167 1 988 2756 t
10 R f
(I R)1 84 1 1196 2756 t
7 I f
(p)1285 2716 w
10 R f
(is the vector of parameters being estimated,)6 1744 1 1353 2756 t
10 S f
(h)3123 2756 w
7 I f
(i)3194 2776 w
10 R f
( R)1 51(: I)1 102 2 3254 2756 t
7 I f
(p)3412 2716 w
10 S f
(\256)3496 2756 w
10 R f
( is a model for the)5 729(I R)1 84 2 3603 2756 t
10 I f
(i)4442 2756 w
10 R f
(th observation)1 570 1 4470 2756 t
(\(typically)720 2876 w
10 S f
(h)1137 2876 w
7 I f
(i)1208 2896 w
10 R f
(has the form)2 517 1 1270 2876 t
10 S f
(h)1821 2876 w
7 I f
(i)1892 2896 w
10 R f
(\()1928 2876 w
10 S f
(b)1969 2876 w
10 R f
(\))2032 2876 w
10 S f
(= h)1 164 1 2122 2876 t
10 R f
(\()2294 2876 w
10 S f
(b)2335 2876 w
10 R f
(,)2398 2876 w
10 I f
(y)2431 2876 w
7 I f
(i)2486 2896 w
10 R f
(\), where)1 335 1 2522 2876 t
10 I f
(y)2891 2876 w
7 I f
(i)2946 2896 w
10 R f
(is the)1 223 1 3008 2876 t
10 I f
(i)3265 2876 w
10 R f
(th observation\), and)2 812 1 3293 2876 t
10 S f
(g \316)1 153 1 4138 2876 t
10 R f
(I R)1 84 1 4332 2876 t
7 S f
(n)4421 2836 w
10 R f
(is a vector of)3 542 1 4498 2876 t
( simply-bounded counterparts, DGLGB and)4 1812( of DGLG and DGLF and of their)7 1441( Use)1 218(nuisance parameters.)1 849 4 720 2996 t
(DGLFB, is summarized in the accompanying P)6 1974 1 720 3116 t
8 R f
(ORT)2694 3116 w
10 R f
( note gives details of my test)6 1216( This)1 240(-style description.)1 729 3 2855 3116 t
(program,)720 3236 w
10 CW f
(pmain.f)1154 3236 w
10 R f
(, and two subroutines,)3 913 1 1574 3236 t
10 CW f
(rhpoil.f)2558 3236 w
10 R f
(\(for RHO\) and)2 609 1 3074 3236 t
10 CW f
(brj.f)3719 3236 w
10 R f
( that this)2 369(\(for CALCRJ\),)1 616 2 4055 3236 t
( latter include several choices of the)6 1495( The)1 215(test program uses.)2 744 3 720 3356 t
10 I f
(r)3208 3356 w
10 R f
(and)3281 3356 w
10 S f
(r)3459 3356 w
10 R f
( pro-)1 200( The)1 214(, described in detail below.)4 1112 3 3514 3356 t
( the menu is:)3 514( Currently)1 433(gram is menu driven; menu item 1 causes the menu to be displayed.)12 2711 3 720 3476 t
10 CW f
(1 = LIST MENU)3 780 1 1500 3656 t
(2 = READ IV)3 660 1 1500 3776 t
(3 = READ V)3 600 1 1500 3896 t
(4 = READ ALG: 1 = GLG, 2 = GLF, 3 = GLGB, 4 = GLFB)15 3000 1 1500 4016 t
(5 = READ ALL OF X0)5 1080 1 1500 4136 t
(6 = COPY X TO X0)5 960 1 1500 4256 t
(7 = START)2 540 1 1500 4376 t
(8 = CONTINUE)2 720 1 1500 4496 t
(9 = READ COMMANDS FROM SPECIFIED FILE)6 2220 1 1500 4616 t
(10 = READ PROBLEM)3 1020 1 1500 4736 t
(11 = READ RHO)3 780 1 1500 4856 t
(12 = READ MODEL)3 900 1 1500 4976 t
(13 = CHECK RHO DERIVATIVES)4 1560 1 1500 5096 t
(14 = READ P)3 660 1 1500 5216 t
(15 = READ X0 COMPONENTWISE)4 1560 1 1500 5336 t
(16 = read new Y)4 900 1 1500 5456 t
(17 = negate RHO \(negative ==> use weights; see KW = 19\))11 3300 1 1500 5576 t
(18 = read KDIAG: 1 = from X*, 2 = from X0, 3 = both)14 3060 1 1500 5696 t
(19 = read KW: 1 = RHO", 2 = IRLS, 3 = score, 4 = deviance)15 3420 1 1500 5816 t
(20 = READ B \(format i, b\(1,i\), b\(2,i\)\))7 2280 1 1500 5936 t
(21,22 = Read,Show RHOI \(componentwise\))4 2280 1 1500 6056 t
( ")1 540(23,24 = Read,Show RHOR)3 1320 2 1500 6176 t
(25 = Show range of RHOR components)6 2040 1 1500 6296 t
(26,27 = Show IV, V components)5 1740 1 1500 6416 t
(28 = Read and echo comment)5 1560 1 1500 6536 t
(READ IV)1 422 1 970 6692 t
10 R f
(reads pairs of integers \()4 944 1 1419 6692 t
10 I f
(i)2363 6692 w
10 R f
(,)2391 6692 w
10 I f
(j)2443 6692 w
10 R f
(\) and sets IV\()3 546 1 2471 6692 t
10 I f
(i)3017 6692 w
10 R f
(\) :=)1 144 1 3045 6692 t
10 I f
(j)3216 6692 w
10 R f
(until)3271 6692 w
10 I f
(i)3482 6692 w
10 S f
(\243)3537 6692 w
10 R f
( is for overriding default IV)5 1116(0. This)1 305 2 3619 6692 t
(values.)720 6812 w
10 CW f
(READ V)1 370 1 970 6968 t
10 R f
(reads pairs of numbers \()4 1004 1 1375 6968 t
10 I f
(i)2379 6968 w
10 R f
(,)2407 6968 w
10 S f
(s)2467 6968 w
10 R f
(\) \()1 101 1 2527 6968 t
10 I f
(i)2628 6968 w
10 R f
(an integer,)1 431 1 2691 6968 t
10 S f
(s)3157 6968 w
10 R f
(a floating-point number\) and sets V\()5 1511 1 3252 6968 t
10 I f
(i)4763 6968 w
10 R f
(\) :=)1 153 1 4791 6968 t
10 S f
(s)4980 6968 w
10 R f
(until)720 7088 w
10 I f
(i)929 7088 w
10 S f
(\243)982 7088 w
10 R f
( is for overriding default V values.)6 1378(0. This)1 303 2 1062 7088 t
10 CW f
(READ ALG)1 480 1 970 7244 t
10 R f
(reads an integer that tells which solver to invoke.)8 1961 1 1475 7244 t
( 18:34:02)1 1772( 16, 1991)2 375( September)1 1307(DRAFT for comment)2 866 4 720 7680 t
cleartomark
showpage
saveobj restore
%%EndPage: 1 1
%%Page: 2 2
/saveobj save def
mark
2 pagesetup
10 R f
(- 2 -)2 166 1 2797 480 t
10 CW f
(READ ALL OF X0)3 840 1 970 875 t
10 R f
(lets you specify an alternate starting guess)6 1686 1 1835 875 t
10 I f
(x)3546 875 w
10 S f
(\260 =)1 144 1 3598 875 t
(\351)3799 838 w
(\353)3799 938 w
(g \260)1 89 1 3856 935 t
(b \260)1 103 1 3849 835 t
(\371)3952 838 w
(\373)3952 938 w
10 R f
(.)3990 875 w
10 CW f
(COPY X TO X0)3 720 1 970 1086 t
10 R f
(copies the current)2 709 1 1715 1086 t
10 I f
(x)2449 1086 w
10 R f
(vector to)1 352 1 2518 1086 t
10 I f
(x)2895 1086 w
10 S f
(\260)2947 1086 w
10 R f
(for use in a subsequent)4 915 1 3012 1086 t
10 CW f
(START)3952 1086 w
10 R f
(.)4252 1086 w
10 CW f
(START)970 1242 w
10 R f
(starts the algorithm from scratch \(from)5 1550 1 1295 1242 t
10 I f
(x)2870 1242 w
10 S f
(\260)2922 1242 w
10 R f
(\).)2962 1242 w
10 CW f
(CONTINUE)970 1398 w
10 R f
( the current)2 496(continues the iteration where it left off \(from)7 1928 2 1495 1398 t
10 I f
(x)3965 1398 w
10 R f
( you have)2 430(\). Presumably)1 601 2 4009 1398 t
(altered one or more IV or V components.)7 1646 1 720 1518 t
10 CW f
( FILE)1 308(READ COMMANDS FROM SPECIFIED)3 1701 2 970 1674 t
10 R f
(must be followed by the name of an input file on)10 2028 1 3012 1674 t
( exhausted, after which)3 955(the next input line; further commands are read from this new input file until it is)15 3365 2 720 1794 t
( input diversions can be nested to a maxi-)8 1747( These)1 299( file.)1 194(command reading resumes with the previous input)6 2080 4 720 1914 t
(mum depth of 5.)3 661 1 720 2034 t
10 CW f
(READ PROBLEM)1 732 1 970 2190 t
10 R f
( this)1 183( If)1 129( name of a problem input file on the next input line.)11 2208(expects to read the)3 781 4 1739 2190 t
(name is)1 311 1 720 2310 t
10 CW f
(-)1059 2310 w
10 R f
( problem file,)2 545( The)1 207( sign\), then the problem is read from the current input file.)11 2349(\(i.e., a single minus)3 792 4 1147 2310 t
(in turn, starts with)3 728 1 720 2430 t
( one-line description of the data set;)6 1431(1. a)1 169 2 936 2550 t
( line containing 8 integers \(format 8I4\):)6 1581(2. a)1 169 2 936 2670 t
(2.1)1080 2790 w
10 I f
(n)1255 2790 w
10 R f
(, the number of observations;)4 1168 1 1305 2790 t
(2.2)1080 2910 w
10 I f
(p)1255 2910 w
10 R f
(, the number of parameters, i.e.,)5 1269 1 1305 2910 t
10 I f
(x)2599 2910 w
10 R f
(components, including the nuisance parameters;)4 1922 1 2668 2910 t
(2.3)1080 3030 w
10 I f
(k)1255 3030 w
7 S f
(h)1310 3050 w
10 R f
(, the choice of)3 565 1 1360 3030 t
10 S f
(h)1950 3030 w
10 R f
(, described below;)2 729 1 2010 3030 t
(2.4)1080 3150 w
10 I f
(m)1255 3150 w
10 R f
(, the number of extra data per observation \(beyond)8 2020 1 1327 3150 t
10 I f
(y)3372 3150 w
7 I f
(i)3427 3170 w
10 R f
(and)3480 3150 w
10 I f
(c)3649 3150 w
7 I f
(i)3704 3170 w
10 R f
(\) to read;)2 360 1 3732 3150 t
(2.5)1080 3270 w
10 I f
(k)1255 3270 w
7 S f
(r)1310 3290 w
10 R f
(the choice of)2 515 1 1382 3270 t
10 S f
(r)1922 3270 w
10 R f
(, described below;)2 729 1 1977 3270 t
(2.6)1080 3390 w
10 I f
(dummy)1255 3390 w
10 R f
(\(ignored; for compatibility with Frome's data format\);)6 2170 1 1568 3390 t
(2.7)1080 3510 w
10 I f
(dummy)1255 3510 w
10 R f
(\(ignored; for compatibility with Frome's data format\);)6 2170 1 1568 3510 t
(2.8)1080 3630 w
10 I f
(p)1255 3630 w
7 I f
(S)1316 3650 w
10 R f
(\(the number of non-nuisance parameters; if)5 1726 1 1384 3630 t
10 I f
(p)3135 3630 w
7 I f
(S)3196 3650 w
10 R f
(is read as 0, then)4 668 1 3264 3630 t
10 I f
(p)3957 3630 w
7 I f
(S)4018 3650 w
10 R f
(is set to)2 306 1 4086 3630 t
10 I f
(p)4417 3630 w
10 R f
(\).)4467 3630 w
(3. For)1 269 1 936 3750 t
10 I f
(k)1235 3750 w
7 S f
(h)1290 3770 w
10 S f
(>)1389 3750 w
10 R f
( guess)1 253(2, the initial)2 491 2 1493 3750 t
10 I f
(x)2268 3750 w
10 S f
(\260)2320 3750 w
10 R f
(; for)1 175 1 2360 3750 t
10 I f
(k)2566 3750 w
7 S f
(h)2621 3770 w
10 S f
(\243)2712 3750 w
10 R f
(2 and)1 225 1 2808 3750 t
10 I f
(p)3064 3750 w
7 I f
(S)3125 3770 w
10 S f
(<)3217 3750 w
10 I f
(p)3321 3750 w
10 R f
(, the initial guess)3 696 1 3371 3750 t
10 I f
(x)4098 3750 w
10 S f
(\260)4150 3750 w
7 I f
(p)4161 3770 w
4 I f
(S)4202 3784 w
7 S f
(+)4238 3770 w
7 R f
(1 ,)1 58 1 4288 3770 t
(. . .)2 90 1 4369 3753 t
(,)4482 3770 w
7 I f
(p)4505 3770 w
10 R f
(for the nui-)2 461 1 4579 3750 t
(sance parameters.)1 708 1 936 3884 t
(4. A)1 197 1 936 4004 t
10 CW f
(FORMAT)1158 4004 w
10 R f
(for reading the next)3 784 1 1543 4004 t
10 I f
(n)2352 4004 w
10 R f
(lines.)2427 4004 w
10 S f
(*)936 4124 w
10 R f
(.)986 4124 w
10 I f
(n)1068 4124 w
10 R f
(data lines giving)2 675 1 1150 4124 t
10 I f
(y)1857 4124 w
7 I f
(i)1912 4144 w
10 R f
(,)1940 4124 w
10 I f
(c)1998 4124 w
7 I f
(i)2053 4144 w
10 R f
(, and)1 202 1 2081 4124 t
10 I f
(A)2316 4124 w
7 I f
(i)2388 4144 w
7 R f
(,)2413 4144 w
7 I f
(j)2442 4144 w
10 R f
(, 1)1 108 1 2470 4124 t
10 S f
(\243)2611 4124 w
10 I f
(j)2699 4124 w
10 S f
(\243)2760 4124 w
10 I f
(m)2848 4124 w
10 R f
( matrix)1 294(. The)1 238 2 2920 4124 t
10 I f
(A)3485 4124 w
7 I f
(T)3557 4084 w
10 R f
(thus defined is passed to BRJ \(the)6 1403 1 3637 4124 t
(CALCRJ subroutine\) as UR, while the \()6 1722 1 720 4244 t
10 I f
(y)2442 4244 w
7 I f
(i)2497 4264 w
10 R f
(,)2525 4244 w
10 I f
(c)2595 4244 w
7 I f
(i)2650 4264 w
10 R f
(\) pairs are passed to RHPOIL \(the RHO subroutine\) as)9 2362 1 2678 4244 t
(RHOR.)720 4364 w
10 CW f
(READ RHO)1 480 1 970 4520 t
10 R f
(causes a new)2 520 1 1475 4520 t
10 I f
(k)2020 4520 w
7 S f
(r)2075 4540 w
10 R f
(to be read from the next input line.)7 1387 1 2147 4520 t
10 CW f
(READ MODEL)1 600 1 970 4676 t
10 R f
(causes a new)2 520 1 1595 4676 t
10 I f
(k)2140 4676 w
7 S f
(h)2195 4696 w
10 R f
(to be read from the next input line.)7 1387 1 2270 4676 t
10 CW f
(CHECK RHO DERIVATIVES)2 1262 1 970 4832 t
10 R f
( forget whether it copes with nuisance parame-)7 1884( I)1 84(is a debugging item.)3 814 3 2258 4832 t
(ters.)720 4952 w
10 CW f
(READ P)1 360 1 970 5108 t
10 R f
(causes a new)2 520 1 1355 5108 t
10 I f
(p)1900 5108 w
10 R f
(to be read from the next input line.)7 1387 1 1975 5108 t
10 CW f
(READ X0 COMPONENTWISE)2 1260 1 970 5264 t
10 R f
(reads pairs \()2 487 1 2255 5264 t
10 I f
(i)2742 5264 w
10 R f
(,)2770 5264 w
10 S f
(s)2820 5264 w
10 R f
(\) and sets)2 377 1 2880 5264 t
10 I f
(x)3282 5264 w
10 S f
(\260)3334 5264 w
7 I f
(i)3385 5284 w
10 R f
(:=)3438 5264 w
10 S f
(s)3547 5264 w
10 R f
(until)3632 5264 w
10 I f
(i)3841 5264 w
10 S f
(\243)3894 5264 w
10 R f
(0.)3974 5264 w
10 CW f
(read new Y)2 614 1 970 5420 t
10 R f
(causes new \()2 523 1 1616 5420 t
10 I f
(y)2139 5420 w
7 I f
(i)2194 5440 w
10 R f
(,)2222 5420 w
10 I f
(c)2279 5420 w
7 I f
(i)2334 5440 w
10 R f
(\) values to be read \(under the previously read)8 1870 1 2362 5420 t
10 CW f
(FORMAT)4265 5420 w
10 R f
(\) from the)2 415 1 4625 5420 t
(next)720 5540 w
10 I f
(n)917 5540 w
10 R f
(input lines.)1 445 1 992 5540 t
10 CW f
(negate RHO)1 606 1 970 5696 t
10 R f
(causes)1607 5696 w
10 I f
(k)1898 5696 w
7 S f
(r)1953 5716 w
10 R f
(:)2041 5696 w
10 S f
(= -)1 167 1 2085 5696 t
10 I f
(k)2268 5696 w
7 S f
(r)2323 5716 w
10 R f
( values of)2 402(; for some)2 417 2 2370 5696 t
10 S f
(\357)3221 5713 w
10 I f
(k)3269 5696 w
7 S f
(r)3324 5716 w
10 S f
(\357)3371 5713 w
10 R f
(\(those for which I knew the appropriate)6 1621 1 3419 5696 t
10 I f
(w)720 5816 w
10 R f
(\),)787 5816 w
10 I f
(k)890 5816 w
7 S f
(r)945 5836 w
10 S f
(<)1037 5816 w
10 R f
(0 causes Roy's choice of the Gauss-Newton model, i.e.,)8 2389 1 1137 5816 t
10 I f
(J)3571 5816 w
7 I f
(T)3626 5776 w
10 R f
(diag \()1 213 1 3681 5816 t
10 I f
(w)3902 5816 w
10 R f
(\))3977 5816 w
10 I f
(J)4026 5816 w
10 R f
( of)1 127(, to be used instead)4 843 2 4070 5816 t
10 I f
(J)720 6036 w
7 I f
(T)775 5996 w
10 R f
(diag \()1 213 1 830 6036 t
10 S f
(\266 h)1 117 1 1076 6116 t
7 I f
(i)1198 6135 w
7 R f
(2)1198 6076 w
10 S f
(\266)1079 5976 w
7 R f
(2)1133 5936 w
10 S f
(r)1184 5976 w
10 S1 f
(_ ___)1 195 1 1062 6006 t
10 R f
(\))1275 6036 w
10 I f
(J)1324 6036 w
10 R f
(.)1368 6036 w
10 CW f
(read KDIAG)1 608 1 970 6272 t
10 R f
( iterated leave-one-out diagnostics with approxi-)5 1983(is for comparing one-step and fully)5 1446 2 1611 6272 t
( real testing rou-)3 657( [The)1 238( finite-difference Hessian approximation\).)3 1679(mations based on a quadratic model \(with a)7 1746 4 720 6392 t
(tine, louchk.c, is written in C and is not included in the Fortran-only version of this test program.])17 3906 1 720 6512 t
10 CW f
(read KW)1 420 1 970 6668 t
10 R f
(specifies the choice of diagonal matrix)5 1544 1 1415 6668 t
10 S f
(\341)2984 6668 w
10 I f
(w)3025 6668 w
10 S f
(\361)3100 6668 w
10 R f
(used in the Gauss-Newton Hessian approxima-)5 1882 1 3158 6668 t
(tion)720 6788 w
10 I f
(J)901 6788 w
7 I f
(T)956 6748 w
10 S f
(\341)1011 6788 w
10 I f
(w)1052 6788 w
10 S f
(\361)1127 6788 w
10 I f
(J)1168 6788 w
10 R f
(.)1212 6788 w
10 CW f
(READ B)1 368 1 970 6944 t
10 R f
( bounds passed to DGLFB and)5 1275(permits setting selected components of the array of simple)8 2394 2 1371 6944 t
( all lower bounds are)4 837(DGLGB. Initially)1 736 2 720 7064 t
10 S f
(- \245)1 144 1 2318 7064 t
10 R f
(and all upper bounds are)4 981 1 2487 7064 t
10 S f
(+ \245)1 144 1 3493 7064 t
10 R f
(.)3637 7064 w
10 CW f
(Read,Show RHOI/RHOR \(componentwise\))2 2108 1 970 7220 t
10 R f
(and)3107 7220 w
10 CW f
(Show range of RHOR components)4 1760 1 3280 7220 t
10 R f
( 18:34:02)1 1772( 16, 1991)2 375( September)1 1307(DRAFT for comment)2 866 4 720 7680 t
cleartomark
showpage
saveobj restore
%%EndPage: 2 2
%%Page: 3 3
/saveobj save def
mark
3 pagesetup
10 R f
(- 3 -)2 166 1 2797 480 t
(are occasionally useful for debugging or playing with block or fractional regression diagnostics.)12 3835 1 720 840 t
10 CW f
(Show IV, V components)3 1260 1 970 996 t
10 R f
(are mainly for debugging.)3 1037 1 2255 996 t
10 CW f
(Read and echo comment)3 1260 1 970 1152 t
10 R f
(may be useful in input scripts that solve several problems.)9 2320 1 2255 1152 t
10 B f
(Models)720 1392 w
10 S f
(h)1056 1392 w
10 BI f
(\()1124 1392 w
10 S f
(b)1165 1392 w
10 BI f
(\))1228 1392 w
10 R f
(Choices for)1 463 1 970 1548 t
10 I f
(k)1458 1548 w
7 S f
(h)1513 1568 w
10 R f
(include:)1588 1548 w
(\(1\))720 1792 w
10 I f
(Linear)1080 1792 w
10 R f
(:)1347 1792 w
10 S f
(h)1425 1792 w
7 I f
(i)1496 1812 w
10 R f
(\()1532 1792 w
10 S f
(b)1573 1792 w
10 R f
(\))1636 1792 w
10 S f
(=)1726 1792 w
10 I f
(A)1830 1792 w
7 I f
(i)1902 1812 w
10 S f
(b =)1 159 1 1938 1792 t
7 I f
(j)2148 1892 w
7 S f
(=)2179 1892 w
7 R f
(1)2229 1892 w
15 S f
(S)2162 1822 w
7 I f
(p)2174 1678 w
4 I f
(S)2215 1692 w
10 I f
(A)2267 1792 w
7 I f
(i)2339 1812 w
7 R f
(,)2364 1812 w
7 I f
(j)2393 1812 w
10 S f
(b)2429 1792 w
7 I f
(j)2495 1812 w
10 R f
(.)2523 1792 w
(\(2\))720 2116 w
10 I f
(Log-linear)1080 2116 w
10 R f
(:)1508 2116 w
10 S f
(h)1586 2116 w
7 I f
(i)1657 2136 w
10 R f
(\()1693 2116 w
10 S f
(b)1734 2116 w
10 R f
(\))1797 2116 w
10 S f
(=)1887 2116 w
10 R f
(exp \()1 185 1 1991 2116 t
10 I f
(A)2184 2116 w
7 I f
(i)2256 2136 w
10 S f
(b)2292 2116 w
10 R f
(\))2355 2116 w
10 S f
(=)2445 2116 w
10 R f
(exp \()1 185 1 2549 2116 t
7 I f
(j)2744 2216 w
7 S f
(=)2775 2216 w
7 R f
(1)2825 2216 w
15 S f
(S)2758 2146 w
7 I f
(p)2770 2002 w
4 I f
(S)2811 2016 w
10 I f
(A)2863 2116 w
7 I f
(i)2935 2136 w
7 R f
(,)2960 2136 w
7 I f
(j)2989 2136 w
10 S f
(b)3025 2116 w
7 I f
(j)3091 2136 w
10 R f
(\).)3127 2116 w
(\(3\))720 2376 w
10 I f
(Log of \(4\))2 400 1 1080 2376 t
10 R f
(.)1480 2376 w
(\(4\))720 2556 w
10 I f
(Nonlinear Poisson model from Frome's PREG manual)6 2199 1 1080 2556 t
10 R f
(.)3279 2556 w
(\(5\))720 2736 w
10 I f
(Caesium dose effect model \(2.6 of [Frome, 1984]\))7 2006 1 1080 2736 t
10 R f
(.)3086 2736 w
(\(6\))720 2916 w
10 I f
(Lung cancer model \(2.8 of [Frome, 1984]\))6 1709 1 1080 2916 t
10 R f
(.)2789 2916 w
(\(7\))720 3156 w
10 I f
(Logistic of linear)2 690 1 1080 3156 t
10 R f
(:)1770 3156 w
10 S f
(h)1848 3156 w
7 I f
(i)1919 3176 w
10 R f
(\()1955 3156 w
10 S f
(b)1996 3156 w
10 R f
(\))2059 3156 w
10 S f
(=)2149 3156 w
10 R f
(1)2278 3226 w
10 S f
(+)2377 3226 w
10 R f
(exp \()1 185 1 2481 3226 t
10 I f
(A)2674 3226 w
7 I f
(i)2746 3246 w
10 S f
(b)2782 3226 w
10 R f
(\))2845 3226 w
(exp \()1 185 1 2380 3076 t
10 I f
(A)2573 3076 w
7 I f
(i)2645 3096 w
10 S f
(b)2681 3076 w
10 R f
(\))2744 3076 w
10 S1 f
(_ ____________)1 630 1 2264 3126 t
10 R f
(.)2904 3156 w
(\(8\))720 3406 w
10 I f
(Exp-linear)1080 3406 w
10 R f
(:)1507 3406 w
10 S f
(h)1585 3406 w
7 I f
(i)1656 3426 w
10 R f
(\()1692 3406 w
10 I f
(x)1733 3406 w
10 R f
(\))1785 3406 w
10 S f
(=)1875 3406 w
10 R f
(log \()1 169 1 1979 3406 t
10 I f
(A)2156 3406 w
7 I f
(i)2228 3426 w
10 S f
(b)2264 3406 w
10 R f
(\).)2327 3406 w
(\(9\))720 3586 w
10 I f
(Model from p. 204\2615 of [McCullagh and Nelder])7 1972 1 1080 3586 t
10 R f
(:)3052 3586 w
10 S f
(h)1962 3826 w
7 I f
(i)2033 3846 w
10 R f
(\()2069 3826 w
10 S f
(b)2110 3826 w
10 R f
(\))2173 3826 w
10 S f
(= b)1 159 1 2263 3826 t
7 R f
(1)2433 3846 w
10 S f
(+ b)1 159 1 2525 3826 t
7 R f
(2)2695 3846 w
10 R f
(log \()1 169 1 2746 3826 t
10 I f
(A)2923 3826 w
7 I f
(i)2995 3846 w
7 R f
(, 1)1 58 1 3020 3846 t
10 R f
(\))3094 3826 w
10 S f
(+)3184 3826 w
(b)3313 3896 w
7 R f
(4)3379 3916 w
10 S f
(+)3471 3896 w
10 I f
(A)3575 3896 w
7 I f
(i)3647 3916 w
7 R f
(, 2)1 58 1 3672 3916 t
10 S f
(b)3386 3746 w
7 R f
(3)3452 3766 w
10 I f
(A)3503 3746 w
7 I f
(i)3575 3766 w
7 R f
(, 2)1 58 1 3600 3766 t
10 S1 f
(_ _________)1 455 1 3299 3796 t
10 I f
(.)3772 3826 w
10 R f
(In this case,)2 474 1 1080 4076 t
10 I f
(m)1579 4076 w
10 S f
(=)1700 4076 w
10 R f
(2.)1804 4076 w
(\(10\))720 4256 w
10 I f
(Model from p. 205\2616 of [McCullagh and Nelder])7 1972 1 1080 4256 t
10 R f
(:)3052 4256 w
10 S f
(h)1831 4496 w
7 I f
(i)1902 4516 w
10 R f
(\()1938 4496 w
10 S f
(b)1979 4496 w
10 R f
(\))2042 4496 w
10 S f
(= b)1 159 1 2132 4496 t
7 R f
(1)2302 4516 w
10 S f
(+ b)1 159 1 2394 4496 t
7 R f
(2)2564 4516 w
10 R f
(log \()1 169 1 2615 4496 t
10 I f
(A)2792 4496 w
7 I f
(i)2864 4516 w
7 R f
(, 1)1 58 1 2889 4516 t
10 S f
(- b)1 159 1 3004 4496 t
7 R f
(5)3174 4516 w
10 R f
(\))3225 4496 w
10 S f
(+)3315 4496 w
(b)3444 4566 w
7 R f
(4)3510 4586 w
10 S f
(+)3602 4566 w
10 I f
(A)3706 4566 w
7 I f
(i)3778 4586 w
7 R f
(, 2)1 58 1 3803 4586 t
10 S f
(b)3517 4416 w
7 R f
(3)3583 4436 w
10 I f
(A)3634 4416 w
7 I f
(i)3706 4436 w
7 R f
(, 2)1 58 1 3731 4436 t
10 S1 f
(_ _________)1 455 1 3430 4466 t
10 R f
(,)3903 4496 w
(where once again)2 697 1 1080 4746 t
10 I f
(m)1802 4746 w
10 S f
(=)1923 4746 w
10 R f
(2.)2027 4746 w
(\(11\))720 4926 w
10 I f
(Model from p. 202 of [McCullagh and Nelder])7 1872 1 1080 4926 t
10 R f
(:)2952 4926 w
10 S f
(h)2155 5166 w
7 I f
(i)2226 5186 w
10 R f
(\()2262 5166 w
10 S f
(b)2303 5166 w
10 R f
(\))2366 5166 w
10 S f
(= b)1 159 1 2456 5166 t
7 R f
(1)2626 5186 w
10 S f
(+)2718 5166 w
7 I f
(j)2824 5266 w
7 S f
(=)2855 5266 w
7 R f
(1)2905 5266 w
15 S f
(S)2838 5196 w
7 R f
(3)2865 5066 w
10 S f
(a)3000 5236 w
7 I f
(i)3074 5256 w
7 R f
(,)3099 5256 w
7 I f
(j)3128 5256 w
10 S f
(+ b)1 159 1 3205 5236 t
7 R f
(2)3375 5256 w
7 I f
(j)3421 5256 w
7 S f
(+)3452 5256 w
7 R f
(1)3502 5256 w
10 S f
(b)3203 5086 w
7 R f
(2)3269 5106 w
7 I f
(j)3315 5106 w
10 S1 f
(_ ___________)1 575 1 2986 5136 t
10 I f
(.)3579 5166 w
10 B f
(Loss Functions)1 643 1 1080 5546 t
10 S f
(r)1748 5546 w
7 BI f
(i)1814 5566 w
10 BI f
(\()1850 5546 w
10 S f
(h)1891 5546 w
7 BI f
(i)1962 5566 w
10 BI f
(,)1998 5546 w
10 S f
(g)2064 5546 w
10 BI f
(\))2113 5546 w
10 R f
(Choices for)1 463 1 1330 5702 t
10 I f
(k)1818 5702 w
7 S f
(r)1873 5722 w
10 R f
(include:)1945 5702 w
(\(1\))720 5952 w
10 I f
(IRLS Poisson)1 559 1 1080 5952 t
10 R f
(:)1639 5952 w
10 S f
(m = h)2 271 1 1734 5952 t
10 R f
(:)2005 5952 w
10 S f
(r)2100 5952 w
7 I f
(i)2166 5972 w
10 R f
(\()2202 5952 w
10 S f
(h)2243 5952 w
10 R f
(\))2311 5952 w
10 S f
(=)2401 5952 w
10 I f
(c)2505 5952 w
7 I f
(i)2560 5972 w
10 S f
(h -)1 164 1 2596 5952 t
10 I f
(y)2809 5952 w
7 I f
(i)2864 5972 w
10 R f
(log \()1 169 1 2900 5952 t
10 S f
(h)3077 5952 w
10 R f
(\), but we lie about)4 792 1 3145 5952 t
10 S f
(\266 h)1 117 1 4021 6032 t
7 R f
(2)4143 5992 w
10 S f
(\266)4004 5872 w
7 R f
(2)4058 5832 w
10 S f
(r)4109 5872 w
7 I f
(i)4175 5892 w
10 S1 f
(_ ____)1 229 1 3989 5922 t
10 R f
( use the IRLS)3 601(: we)1 211 2 4228 5952 t
(weights)1080 6202 w
10 I f
(w)1416 6202 w
7 I f
(i)1494 6222 w
10 S f
(=)1571 6202 w
10 I f
(c)1675 6202 w
7 I f
(i)1730 6222 w
10 I f
(/)1766 6202 w
10 S f
(h)1802 6202 w
10 R f
(instead of)1 391 1 1887 6202 t
10 S f
(\266 h)1 117 1 2345 6282 t
7 R f
(2)2467 6242 w
10 S f
(\266)2328 6122 w
7 R f
(2)2382 6082 w
10 S f
(r)2433 6122 w
7 I f
(i)2499 6142 w
10 S1 f
(_ ____)1 229 1 2313 6172 t
10 R f
(for both RD\()2 516 1 2577 6202 t
10 I f
(i)3093 6202 w
10 R f
( and W\()2 321(, 1\))1 116 2 3129 6202 t
10 I f
(i)3566 6202 w
10 R f
(\).)3594 6202 w
(\(2\))720 6512 w
10 I f
(Poisson)1080 6512 w
10 R f
(: Again)1 331 1 1397 6512 t
10 S f
(r)1762 6512 w
7 I f
(i)1828 6532 w
10 R f
(\()1864 6512 w
10 S f
(h)1905 6512 w
10 R f
(\))1973 6512 w
10 S f
(=)2063 6512 w
10 I f
(c)2167 6512 w
7 I f
(i)2222 6532 w
10 S f
(h -)1 164 1 2258 6512 t
10 I f
(y)2471 6512 w
7 I f
(i)2526 6532 w
10 R f
(log \()1 169 1 2562 6512 t
10 S f
(h)2739 6512 w
10 R f
( we supply the true)4 800(\), but now)2 426 2 2807 6512 t
10 S f
(\266 h)1 117 1 4110 6592 t
7 R f
(2)4232 6552 w
10 S f
(\266)4093 6432 w
7 R f
(2)4147 6392 w
10 S f
(r)4198 6432 w
7 I f
(i)4264 6452 w
10 S1 f
(_ ____)1 229 1 4078 6482 t
10 R f
(for both RD\()2 536 1 4352 6512 t
10 I f
(i)4888 6512 w
10 R f
(, 1\))1 116 1 4924 6512 t
(and W\()1 296 1 1080 6692 t
10 I f
(i)1376 6692 w
10 R f
(\).)1404 6692 w
(\(\2612\))720 6942 w
10 I f
(Poisson with IRLS GN model)4 1173 1 1080 6942 t
10 R f
( W\()1 177( is a mixture of \(1\) and \(2\):)7 1084(: This)1 256 3 2253 6942 t
10 I f
(i)3770 6942 w
10 R f
(\) =)1 114 1 3798 6942 t
10 I f
(w)3937 6942 w
7 I f
(i)4015 6962 w
10 R f
(and RD\()1 341 1 4068 6942 t
10 I f
(i)4409 6942 w
10 R f
( =)1 81(, 1\))1 116 2 4445 6942 t
10 S f
(\266 h)1 117 1 4709 7022 t
7 R f
(2)4831 6982 w
10 S f
(\266)4692 6862 w
7 R f
(2)4746 6822 w
10 S f
(r)4797 6862 w
7 I f
(i)4863 6882 w
10 S1 f
(_ ____)1 229 1 4677 6912 t
10 R f
(.)4916 6942 w
( 18:34:02)1 1772( 16, 1991)2 375( September)1 1307(DRAFT for comment)2 866 4 720 7680 t
cleartomark
showpage
saveobj restore
%%EndPage: 3 3
%%Page: 4 4
/saveobj save def
mark
4 pagesetup
10 R f
(- 4 -)2 166 1 2797 480 t
(\(3\))720 840 w
10 I f
(Log Linear Poisson)2 790 1 1080 840 t
10 R f
(:)1870 840 w
10 S f
(r)1948 840 w
7 I f
(i)2014 860 w
10 R f
(\()2050 840 w
10 S f
(h)2091 840 w
10 R f
(\))2159 840 w
10 S f
(=)2249 840 w
10 I f
(c)2353 840 w
7 I f
(i)2408 860 w
10 R f
(exp \()1 185 1 2444 840 t
10 S f
(h)2637 840 w
10 R f
(\))2705 840 w
10 S f
(-)2795 840 w
10 I f
(y)2899 840 w
7 I f
(i)2954 860 w
10 S f
(h)2990 840 w
10 R f
(.)3050 840 w
(\(\2613\))720 1020 w
10 I f
(IRLS GN = No change)4 912 1 1080 1020 t
10 R f
(:)1992 1020 w
10 S f
(m =)1 162 1 2070 1020 t
10 R f
(exp \()1 185 1 2281 1020 t
10 S f
(h)2474 1020 w
10 R f
(\) ==)1 170 1 2542 1020 t
10 S f
(>)2712 1020 w
10 I f
(w)2792 1020 w
7 I f
(i)2870 1040 w
10 S f
(=)2947 1020 w
10 I f
(c)3051 1020 w
7 I f
(i)3106 1040 w
10 R f
(exp \()1 185 1 3166 1020 t
10 S f
(h)3359 1020 w
10 R f
(\) ==)1 170 1 3427 1020 t
10 S f
(>)3597 1020 w
10 R f
(IRLS GN =)2 467 1 3677 1020 t
10 S f
(r \242 \242)2 121 1 4169 1020 t
10 R f
(GN.)4315 1020 w
(\(4\))720 1200 w
10 I f
(Square-root Linear for Poisson)3 1259 1 1080 1200 t
10 R f
(:)2339 1200 w
10 S f
(r)2417 1200 w
7 I f
(i)2483 1220 w
10 R f
(\()2519 1200 w
10 S f
(h)2560 1200 w
10 R f
(\))2628 1200 w
10 S f
(=)2718 1200 w
10 I f
(c)2822 1200 w
7 I f
(i)2877 1220 w
10 S f
(h)2913 1200 w
7 R f
(2)2978 1160 w
10 S f
(-)3070 1200 w
10 R f
(2)3174 1200 w
10 I f
(y)3232 1200 w
7 I f
(i)3287 1220 w
10 R f
(log)3323 1200 w
10 S f
(h)3459 1200 w
10 R f
(.)3519 1200 w
(\(\2614\))720 1380 w
10 I f
(IRLS)1080 1380 w
10 R f
(:)1280 1380 w
10 S f
(m = h)2 271 1 1358 1380 t
7 R f
(2)1634 1340 w
10 R f
(==)1702 1380 w
10 S f
(>)1814 1380 w
10 I f
(w)1894 1380 w
7 I f
(i)1972 1400 w
10 S f
(=)2049 1380 w
10 R f
(4)2153 1380 w
10 I f
(c)2211 1380 w
7 I f
(i)2266 1400 w
10 R f
(.)2294 1380 w
(\(5\))720 1560 w
10 I f
(Binomial IRLS)1 606 1 1080 1560 t
10 R f
(:)1686 1560 w
10 S f
(m = h)2 271 1 1778 1560 t
10 R f
(:)2049 1560 w
10 S f
(r)2141 1560 w
7 I f
(i)2207 1580 w
10 R f
(\()2243 1560 w
10 S f
(h)2284 1560 w
10 R f
(\))2352 1560 w
10 S f
(= -)1 167 1 2442 1560 t
10 I f
(y)2625 1560 w
7 I f
(i)2680 1580 w
10 R f
(log \()1 169 1 2716 1560 t
10 S f
(h)2893 1560 w
10 R f
(\))2961 1560 w
10 S f
(-)3051 1560 w
10 R f
(\()3155 1560 w
10 I f
(c)3196 1560 w
7 I f
(i)3251 1580 w
10 S f
(-)3319 1560 w
10 I f
(y)3414 1560 w
7 I f
(i)3469 1580 w
10 R f
( \( 1)2 99(\) log)1 177 2 3505 1560 t
10 S f
(- h)1 155 1 3821 1560 t
10 R f
(\), with RD\()2 486 1 3984 1560 t
10 I f
(i)4470 1560 w
10 R f
( = W\()2 261(, 1\))1 116 2 4506 1560 t
10 I f
(i)4883 1560 w
10 R f
(\) =)1 129 1 4911 1560 t
(IRLS weight =)2 595 1 1080 1680 t
10 S f
(h)1700 1680 w
7 S f
(-)1771 1640 w
7 R f
(1)1821 1640 w
10 R f
(\( 1)1 91 1 1872 1680 t
10 S f
(- h)1 155 1 2003 1680 t
10 R f
(\))2166 1680 w
7 S f
(-)2210 1640 w
7 R f
(1)2260 1640 w
10 R f
(.)2303 1680 w
(\(6\))720 1930 w
10 I f
(Binomial)1080 1930 w
10 R f
(:)1447 1930 w
10 S f
(r)1525 1930 w
7 I f
(i)1591 1950 w
10 R f
(as for \(5\), but RD\()4 740 1 1644 1930 t
10 I f
(i)2384 1930 w
10 R f
( = W\()2 233(, 1\))1 116 2 2420 1930 t
10 I f
(i)2769 1930 w
10 R f
(\) =)1 114 1 2797 1930 t
10 S f
(\266 h)1 117 1 2978 2010 t
7 R f
(2)3100 1970 w
10 S f
(\266)2961 1850 w
7 R f
(2)3015 1810 w
10 S f
(r)3066 1850 w
7 I f
(i)3132 1870 w
10 S1 f
(_ ____)1 229 1 2946 1900 t
10 R f
(.)3185 1930 w
(\(\2616\))720 2170 w
10 I f
(Binomial with IRLS GN model)4 1223 1 1080 2170 t
10 R f
( is a mixture of \(5\) and \(6\), analogous to \(\2612\).)10 1830(: This)1 256 2 2303 2170 t
(\(7\))720 2350 w
10 I f
(Binomial Logistic)1 715 1 1080 2350 t
10 R f
(:)1795 2350 w
10 S f
(r)1873 2350 w
7 I f
(i)1939 2370 w
10 R f
(\()1975 2350 w
10 S f
(h)2016 2350 w
10 R f
(\))2084 2350 w
10 S f
(=)2174 2350 w
10 I f
(c)2278 2350 w
7 I f
(i)2333 2370 w
10 R f
(log \( 1)2 227 1 2369 2350 t
10 S f
(+)2645 2350 w
10 I f
(e)2749 2350 w
7 S f
(h)2804 2310 w
10 R f
(\))2862 2350 w
10 S f
(-)2952 2350 w
10 I f
(y)3056 2350 w
7 I f
(i)3111 2370 w
10 S f
(h)3147 2350 w
10 R f
(.)3207 2350 w
(\(\2617\))720 2580 w
10 I f
(IRLS GN = No change)4 912 1 1080 2580 t
10 R f
(:)1992 2580 w
10 S f
(m =)1 162 1 2070 2580 t
10 R f
(1)2306 2650 w
10 S f
(+)2405 2650 w
10 R f
(exp \()1 185 1 2509 2650 t
10 S f
(h)2702 2650 w
10 R f
(\))2770 2650 w
(exp \()1 185 1 2408 2520 t
10 S f
(h)2601 2520 w
10 R f
(\))2669 2520 w
10 S1 f
(_ __________)1 527 1 2292 2550 t
10 R f
(==)2854 2580 w
10 S f
(>)2966 2580 w
10 I f
(w)3046 2580 w
7 I f
(i)3124 2600 w
10 S f
(= m)1 162 1 3201 2580 t
10 R f
(\( 1)1 91 1 3371 2580 t
10 S f
(- m)1 153 1 3502 2580 t
10 R f
(\))3663 2580 w
10 S f
(=)3753 2580 w
10 R f
(\( 1)1 91 1 3882 2660 t
10 S f
(+)4022 2660 w
10 I f
(e)4126 2660 w
7 S f
(h)4181 2620 w
10 R f
(\))4239 2660 w
7 R f
(2)4277 2620 w
10 I f
(e)4049 2520 w
7 S f
(h)4104 2480 w
10 S1 f
(_ _________)1 468 1 3868 2550 t
10 R f
(.)4346 2580 w
(\(8\))720 2820 w
10 I f
(Probit)1080 2820 w
10 R f
( with)1 203(: \(binomial)1 467 2 1336 2820 t
10 S f
(m = F)2 287 1 2031 2820 t
10 R f
(\()2326 2820 w
10 S f
(h)2367 2820 w
10 R f
(\),)2435 2820 w
10 S f
(F)2518 2820 w
10 R f
(= cumulative normal distribution\):)3 1381 1 2619 2820 t
10 S f
(r)1746 3000 w
7 I f
(i)1812 3020 w
10 R f
(\()1848 3000 w
10 S f
(h)1889 3000 w
10 R f
(\))1957 3000 w
10 S f
(= -)1 167 1 2047 3000 t
10 I f
(y)2230 3000 w
7 I f
(i)2285 3020 w
10 R f
(log [)1 169 1 2321 3000 t
10 S f
(F)2498 3000 w
10 R f
(\()2582 3000 w
10 S f
(h)2623 3000 w
10 R f
(\) ])1 74 1 2691 3000 t
10 S f
(-)2822 3000 w
10 R f
(\()2926 3000 w
10 I f
(c)2967 3000 w
7 I f
(i)3022 3020 w
10 S f
(-)3090 3000 w
10 I f
(y)3185 3000 w
7 I f
(i)3240 3020 w
10 R f
( [ 1)2 99(\) log)1 177 2 3276 3000 t
10 S f
(- F)1 180 1 3601 3000 t
10 R f
(\()3789 3000 w
10 S f
(h)3830 3000 w
10 R f
(\) ])1 74 1 3898 3000 t
10 I f
(.)3988 3000 w
10 R f
(\(\2618\))720 3290 w
10 I f
(IRLS Probit)1 481 1 1080 3290 t
10 R f
(:)1561 3290 w
10 S f
(m = F)2 287 1 1639 3290 t
10 R f
(\()1934 3290 w
10 S f
(h)1975 3290 w
10 R f
(\) ==)1 170 1 2043 3290 t
10 S f
(>)2213 3290 w
10 I f
(w)2293 3290 w
7 I f
(i)2371 3310 w
10 S f
(=)2448 3290 w
10 R f
(2)2577 3360 w
10 S f
(p m)1 121 1 2635 3360 t
10 R f
(\( 1)1 91 1 2764 3360 t
10 S f
(- m)1 153 1 2895 3360 t
10 R f
(\))3056 3360 w
10 I f
(e)2742 3230 w
7 S f
(- h)1 92 1 2797 3190 t
4 R f
(2)2892 3153 w
10 S1 f
(_ __________)1 542 1 2563 3260 t
10 R f
(.)3115 3290 w
(\(9\))720 3520 w
10 I f
(Weibull)1080 3520 w
10 R f
(:)1391 3520 w
10 S f
(r)1824 3700 w
7 I f
(i)1890 3720 w
10 R f
(\()1926 3700 w
10 S f
(h)1967 3700 w
10 R f
(\))2035 3700 w
10 S f
(= -)1 167 1 2125 3700 t
10 I f
(y)2308 3700 w
7 I f
(i)2363 3720 w
10 R f
(log [ 1)2 227 1 2399 3700 t
10 S f
(-)2675 3700 w
10 R f
(exp \()1 185 1 2779 3700 t
10 S f
(-)2980 3700 w
10 I f
(e)3051 3700 w
7 S f
(h)3106 3660 w
10 R f
(\) ])1 74 1 3164 3700 t
10 S f
(+)3295 3700 w
10 R f
(\()3399 3700 w
10 I f
(c)3440 3700 w
7 I f
(i)3495 3720 w
10 S f
(-)3563 3700 w
10 I f
(y)3658 3700 w
7 I f
(i)3713 3720 w
10 R f
(\))3749 3700 w
10 I f
(e)3798 3700 w
7 S f
(h)3853 3660 w
10 I f
(.)3911 3700 w
10 R f
(\(\2619\))720 3990 w
10 I f
(IRLS Weibull)1 536 1 1080 3990 t
10 R f
(:)1616 3990 w
10 S f
(m =)1 162 1 1694 3990 t
10 R f
(1)1905 3990 w
10 S f
(-)2004 3990 w
10 I f
(e)2108 3990 w
7 S f
(-)2163 3950 w
7 I f
(e)2213 3950 w
4 S f
(h)2250 3913 w
10 R f
(==)2312 3990 w
10 S f
(>)2424 3990 w
10 I f
(w)2504 3990 w
7 I f
(i)2582 4010 w
10 S f
(=)2659 3990 w
10 I f
(e)2788 4070 w
7 I f
(e)2843 4030 w
4 S f
(h)2880 3993 w
10 S f
(-)2966 4070 w
10 R f
(1)3070 4070 w
10 I f
(e)2882 3930 w
7 R f
(2)2937 3890 w
7 S f
(h)2977 3890 w
10 S1 f
(_ _______)1 362 1 2774 3960 t
10 R f
(.)3146 3990 w
(\(10\))720 4230 w
10 I f
(Gamma)1080 4230 w
10 R f
(:)1396 4230 w
10 S f
(r)1474 4230 w
7 I f
(i)1540 4250 w
10 R f
(\()1576 4230 w
10 S f
(h)1617 4230 w
10 R f
(\))1685 4230 w
10 S f
(=)1775 4230 w
10 I f
(c)1879 4230 w
7 I f
(i)1934 4250 w
10 R f
(\()1970 4230 w
10 I f
(y)2011 4230 w
7 I f
(i)2066 4250 w
10 S f
(h -)1 164 1 2102 4230 t
10 R f
(log \()1 169 1 2315 4230 t
10 S f
(h)2492 4230 w
10 R f
(\) \).)1 99 1 2560 4230 t
(\(\26110\))720 4410 w
10 I f
(IRLS GN = No change)4 912 1 1080 4410 t
10 R f
(:)1992 4410 w
10 S f
(m =)1 162 1 2070 4410 t
10 R f
(1)2281 4410 w
10 I f
(/)2363 4410 w
10 S f
(h)2423 4410 w
10 R f
(==)2508 4410 w
10 S f
(>)2620 4410 w
10 I f
(w)2700 4410 w
7 I f
(i)2778 4430 w
10 S f
(= m)1 162 1 2855 4410 t
7 R f
(2)3022 4370 w
10 S f
(= h)1 164 1 3114 4410 t
7 S f
(-)3289 4370 w
7 R f
(2)3339 4370 w
10 R f
(==)3407 4410 w
10 S f
(>)3519 4410 w
10 R f
(IRLS GN =)2 467 1 3599 4410 t
10 S f
(r \242 \242)2 121 1 4091 4410 t
10 R f
(GN.)4237 4410 w
(\(11\))720 4590 w
10 I f
(Nelder and Pregibon Quasi-Likelihood, replicate weights)5 2305 1 1080 4590 t
10 R f
(:)3385 4590 w
10 S f
(r)1287 4905 w
7 I f
(i)1353 4925 w
10 R f
(\()1389 4905 w
10 S f
(h)1430 4905 w
10 R f
(,)1498 4905 w
10 S f
(g)1531 4905 w
10 R f
(\))1580 4905 w
10 S f
(=)1670 4905 w
10 I f
(c)1774 4905 w
7 I f
(i)1829 4925 w
10 S f
(\351)1873 4768 w
(\357)1873 4868 w
(\357)1873 4968 w
(\353)1873 5068 w
10 S1 f
()1923 4905 w
cleartomark
saveobj restore
%%BeginGlobal
/build_12 {
    pop
    /optsize ptsize def
    /osize size def
    /ofont font def

    optsize 2 div dup R exch R f
    0 size 2 mul 3 div dup neg exch 0 exch rmoveto

    (1) show
    rmoveto
    optsize R f
    (\244) show
    f
    (2) show

    optsize ofont f
} def
%%EndGlobal
/saveobj save def
mark
10 S1 f
1923 4905 m
75 build_12
1998 4905 m
10 R f
(log [ 2)2 227 1 2006 4905 t
10 S f
(p f)1 115 1 2241 4905 t
10 I f
(y)2364 4905 w
7 I f
(i)2413 4924 w
7 S f
(q)2413 4865 w
10 R f
(])2465 4905 w
10 S f
(+)2555 4905 w
(\354)2667 4818 w
(\357)2667 4918 w
(\356)2667 5018 w
(f)2870 4995 w
10 R f
(\( 1)1 91 1 2930 4995 t
10 S f
(- q)1 147 1 3061 4995 t
10 R f
(\))3216 4995 w
10 I f
(y)2741 4845 w
7 I f
(i)2790 4864 w
7 R f
(1)2790 4805 w
7 S f
(- q)1 86 1 2836 4805 t
10 S f
(- h)1 155 1 2970 4845 t
7 S f
(y)3130 4805 w
7 R f
(\( 1)1 63 1 3183 4805 t
7 S f
(- q)1 86 1 3257 4805 t
7 R f
(\))3348 4805 w
10 S1 f
(_ _____________)1 668 1 2726 4895 t
10 S f
(\374)3404 4818 w
(\357)3404 4918 w
(\376)3404 5018 w
10 I f
(y)3469 4905 w
7 I f
(i)3524 4925 w
10 S f
(+)3601 4905 w
(f)3863 4975 w
10 R f
(\( 2)1 91 1 3923 4975 t
10 S f
(- q)1 147 1 4054 4975 t
10 R f
(\))4209 4975 w
10 S f
(h)3730 4825 w
7 S f
(y)3795 4785 w
7 R f
(\( 2)1 63 1 3848 4785 t
7 S f
(- q)1 86 1 3922 4785 t
7 R f
(\))4013 4785 w
10 S f
(-)4092 4825 w
10 I f
(y)4187 4825 w
7 I f
(i)4236 4844 w
7 R f
(2)4236 4785 w
7 S f
(- q)1 86 1 4282 4785 t
10 S1 f
(_ _____________)1 676 1 3715 4875 t
10 S f
(\371)4401 4768 w
(\357)4401 4868 w
(\357)4401 4968 w
(\373)4401 5068 w
10 R f
(,)4447 4905 w
(where)1080 5245 w
10 S f
(f = g)2 246 1 1355 5245 t
7 R f
(1)1612 5265 w
10 R f
(,)1655 5245 w
10 S f
(q = g)2 246 1 1712 5245 t
7 R f
(2)1969 5265 w
10 R f
(, and)1 201 1 2012 5245 t
10 S f
(y = g)2 263 1 2245 5245 t
7 R f
(3)2519 5265 w
10 R f
(. \(L)1 176 1 2562 5245 t
10 S f
(\242)2738 5245 w
10 R f
(Ho)2763 5245 w
11 R f
(\303)2842 5240 w
10 R f
( gives special forms for)4 969(pital's rule)1 437 2 2885 5245 t
10 S f
(r)4324 5245 w
7 I f
(i)4390 5265 w
10 R f
(when)4451 5245 w
10 S f
(q)4700 5245 w
10 R f
(= 1 or)2 255 1 4785 5245 t
(2.\))1080 5365 w
(\(\26111\))720 5545 w
10 I f
(IRLS)1080 5545 w
10 R f
(:)1280 5545 w
10 S f
(m = h)2 271 1 1358 5545 t
7 S f
(y)1634 5505 w
10 R f
(==)1715 5545 w
10 S f
(>)1827 5545 w
10 I f
(w)1907 5545 w
7 I f
(i)1985 5565 w
10 S f
(=)2062 5545 w
10 I f
(c)2166 5545 w
7 I f
(i)2221 5565 w
10 S f
(y)2257 5545 w
7 R f
(2)2331 5505 w
10 S f
(h)2382 5545 w
7 R f
(\( 2)1 63 1 2447 5505 t
7 S f
(- q)1 102 1 2537 5505 t
7 R f
(\))2644 5505 w
7 S f
(y -)1 121 1 2678 5505 t
7 R f
(2)2833 5505 w
10 R f
(.)2876 5545 w
(\(12\))720 5725 w
10 I f
(Least Squares)1 564 1 1080 5725 t
10 R f
(:)1644 5725 w
10 S f
(r)1722 5725 w
7 I f
(i)1788 5745 w
10 R f
(\()1824 5725 w
10 S f
(h)1865 5725 w
10 R f
(\))1933 5725 w
10 S f
(=)2023 5725 w
10 S1 f
()2127 5725 w
2127 5725 m
75 build_12
2202 5725 m
10 R f
(\()2210 5725 w
10 S f
(h -)1 164 1 2251 5725 t
10 I f
(y)2464 5725 w
7 I f
(i)2519 5745 w
10 R f
(\))2555 5725 w
7 R f
(2)2593 5685 w
10 R f
(.)2636 5725 w
(\(13\))720 5905 w
10 I f
(Nelder and Pregibon Quasi-Likelihood, variance weights)5 2299 1 1080 5905 t
10 R f
(:)3379 5905 w
10 S f
(r)1224 6220 w
7 I f
(i)1290 6240 w
10 R f
(\()1326 6220 w
10 S f
(h)1367 6220 w
10 R f
(,)1435 6220 w
10 S f
(g)1468 6220 w
10 R f
(\))1517 6220 w
10 S f
(=)1607 6220 w
10 S1 f
()1711 6220 w
1711 6220 m
75 build_12
1786 6220 m
10 R f
(log [ 2)2 227 1 1794 6220 t
10 S f
(p f)1 115 1 2029 6220 t
10 I f
(y)2152 6220 w
7 I f
(i)2201 6239 w
7 S f
(q)2201 6180 w
10 I f
(/ c)1 80 1 2253 6220 t
7 I f
(i)2344 6240 w
10 R f
(])2380 6220 w
10 S f
(+)2470 6220 w
10 I f
(c)2574 6220 w
7 I f
(i)2629 6240 w
10 S f
(\351)2673 6083 w
(\357)2673 6183 w
(\357)2673 6283 w
(\353)2673 6383 w
(\354)2731 6133 w
(\357)2731 6233 w
(\356)2731 6333 w
(f)2934 6310 w
10 R f
(\( 1)1 91 1 2994 6310 t
10 S f
(- q)1 147 1 3125 6310 t
10 R f
(\))3280 6310 w
10 I f
(y)2805 6160 w
7 I f
(i)2854 6179 w
7 R f
(1)2854 6120 w
7 S f
(- q)1 86 1 2900 6120 t
10 S f
(- h)1 155 1 3034 6160 t
7 S f
(y)3194 6120 w
7 R f
(\( 1)1 63 1 3247 6120 t
7 S f
(- q)1 86 1 3321 6120 t
7 R f
(\))3412 6120 w
10 S1 f
(_ _____________)1 668 1 2790 6210 t
10 S f
(\374)3468 6133 w
(\357)3468 6233 w
(\376)3468 6333 w
10 I f
(y)3533 6220 w
7 I f
(i)3588 6240 w
10 S f
(+)3665 6220 w
(f)3927 6290 w
10 R f
(\( 2)1 91 1 3987 6290 t
10 S f
(- q)1 147 1 4118 6290 t
10 R f
(\))4273 6290 w
10 S f
(h)3794 6140 w
7 S f
(y)3859 6100 w
7 R f
(\( 2)1 63 1 3912 6100 t
7 S f
(- q)1 86 1 3986 6100 t
7 R f
(\))4077 6100 w
10 S f
(-)4156 6140 w
10 I f
(y)4251 6140 w
7 I f
(i)4300 6159 w
7 R f
(2)4300 6100 w
7 S f
(- q)1 86 1 4346 6100 t
10 S1 f
(_ _____________)1 676 1 3779 6190 t
10 S f
(\371)4465 6083 w
(\357)4465 6183 w
(\357)4465 6283 w
(\373)4465 6383 w
10 R f
(,)4511 6220 w
(where)1080 6555 w
10 S f
(f)1348 6555 w
10 R f
(,)1400 6555 w
10 S f
(q)1450 6555 w
10 R f
(, and)1 194 1 1502 6555 t
10 S f
(y)1721 6555 w
10 R f
(are as for \(11\).)3 586 1 1815 6555 t
( \(\26111\).)1 266( see)1 177(\(\26113\) :)1 388 3 720 6735 t
( 18:34:02)1 1772( 16, 1991)2 375( September)1 1307(DRAFT for comment)2 866 4 720 7680 t
cleartomark
showpage
saveobj restore
%%EndPage: 4 4
%%Trailer
done
%%Pages: 4
%%DocumentFonts: Times-Roman Symbol Times-BoldItalic Courier Times-Roman Times-Bold Times-Italic

.