data grey;
	input s Y;
	lnY=log(Y); lns=log(s);
	datalines;
0.8	80
1.6	65
2.4	47
3.2	43
4.0	45
4.8	38
5.6	36
6.4	29
;

data lateblight;
	input s t y; *--distance (s), time (t), and disease (y);
	logit=log(y/(1-y));
	datalines;
0	3	0.378
1	3	0.378
2	3	0.269
3	3	0.109
4	3	0.047
5	3	0.029
5	3	0.119
6	3	0.057
7	3	0.119
8	3	0.005
9	3	0.011
10	3	0.004
11	3	0.004
12	3	0.002
13	3	0.004
14	3	0.005
14	3	0.005
15	3	0.005
16	3	0.005
17	3	0.002
18	3	0.002
19	3	0.002
20	3	0.015
21	3	.
22	3	.
23	3	.
23	3	.
24	3	.
25	3	.
26	3	.
0	6	0.731
1	6	0.731
2	6	0.668
3	6	0.500
4	6	0.332
5	6	0.332
5	6	0.401
6	6	0.332
7	6	0.332
8	6	0.198
9	6	0.198
10	6	0.198
11	6	0.142
12	6	0.269
13	6	0.332
14	6	0.269
14	6	0.198
15	6	0.142
16	6	0.100
17	6	0.100
18	6	0.142
19	6	0.142
20	6	0.100
21	6	0.047
22	6	0.076
23	6	0.057
23	6	0.057
24	6	0.057
25	6	0.013
26	6	0.004
0	10	0.953
1	10	0.971
2	10	0.953
3	10	0.953
4	10	0.870
5	10	0.900
5	10	0.846
6	10	0.846
7	10	0.881
8	10	0.802
9	10	0.731
10	10	0.731
11	10	0.622
12	10	0.622
13	10	0.731
14	10	0.668
14	10	0.731
15	10	0.668
16	10	0.574
17	10	0.646
18	10	0.574
19	10	0.646
20	10	0.646
21	10	0.574
22	10	0.450
23	10	0.450
23	10	0.500
24	10	0.269
25	10	0.332
26	10	0.269
0	14	0.998
1	14	0.993
2	14	0.990
3	14	0.988
4	14	0.985
5	14	0.982
5	14	0.980
6	14	0.953
7	14	0.980
8	14	0.980
9	14	0.900
10	14	0.953
11	14	0.953
12	14	0.900
13	14	0.953
14	14	0.900
14	14	0.900
15	14	0.858
16	14	0.802
17	14	0.858
18	14	0.690
19	14	0.802
20	14	0.802
21	14	0.802
22	14	0.690
23	14	0.731
23	14	0.769
24	14	0.731
25	14	0.731
26	14	0.690
;

data strawberry;
	input Y;
	yy=Y/15;
	datalines;
6
4
1
7
2
5
4
0
2
1
9
10
7
8
11
7
8
3
15
7
5
11
5
5
3
10
2
5
5
8
6
5
4
6
1
0
0
2
3
4
6
0
4
8
2
1
2
1
0
3
4
6
1
6
3
1
1
8
0
1
1
0
6
2
3
0
3
3
0
;

proc reg data=grey;
	title 'Grey leaf spot of maize';
	model lnY = s;
	plot r.*s;
run;

proc reg data=lateblight;
	title 'Late blight data of Minogue and Fry; LOGISTIC each time';
	by t; /*separate for each time */
	model logit=s;
proc reg data=lateblight;
	title 'Late blight data of Minogue and Fry - LOGIST-LOGISTIC ';
	model logit= s t; /*combined analysis */
run;

proc means data=strawberry mean variance;
	title 'Summary statistics for Phomopsis leaf blight';
	var Y yy;
run;