#
#  Matched case-control study data from
#  Trichopoulos et al. (1976) Br J Obs Gynaec 83: 645-650.
# 
#              coef exp(coef) se(coef)    z       p
#  spontaneous 1.99      7.29    0.352 5.63 1.8e-08
#  induced     1.41      4.09    0.361 3.91 9.4e-05
# 
#  Likelihood ratio test=53.1  on 2 df, p=2.87e-12  n= 248 
#
set locus education    qua . 1 = 0-5 2= 6-11 3= 12+
set locus age          qua . age in years
set locus parity       qua . count
set locus induced      qua . TOP 0=0 1=1 2=2+
set locus case         qua . 
set locus spontaneous  qua . miscarriage 0=0 1=1 2=2+
set locus stratum      qua . matching stratum
read pedigree inline
#
# strat id      education age parity induced case spontaneous stratum
1       1    f1 m2 f 1      26  6      1       1    2           1   
1       84   f1 m2 f 1      26  6      2       0    0           1   
1       166  f1 m2 f 1      26  6      2       0    0           1   
2       2    f1 m2 f 1      42  1      1       1    0           2   
2       85   f1 m2 f 1      42  1      0       0    0           2   
2       167  f1 m2 f 1      42  1      0       0    0           2   
3       3    f1 m2 f 1      39  6      2       1    0           3   
3       86   f1 m2 f 1      39  6      2       0    0           3   
3       168  f1 m2 f 1      39  6      2       0    0           3   
4       4    f1 m2 f 1      34  4      2       1    0           4   
4       87   f1 m2 f 1      34  4      0       0    1           4   
4       169  f1 m2 f 1      34  4      0       0    2           4   
5       5    f1 m2 f 2      35  3      1       1    1           5   
5       88   f1 m2 f 2      35  3      2       0    0           5   
5       170  f1 m2 f 2      35  3      0       0    0           5   
6       6    f1 m2 f 2      36  4      2       1    1           6   
6       89   f1 m2 f 2      36  4      1       0    1           6   
6       171  f1 m2 f 2      36  4      0       0    2           6   
7       7    f1 m2 f 2      23  1      0       1    0           7   
7       90   f1 m2 f 2      23  1      0       0    0           7   
7       172  f1 m2 f 2      23  1      0       0    0           7   
8       8    f1 m2 f 2      32  2      0       1    0           8   
8       91   f1 m2 f 2      32  2      2       0    0           8   
8       173  f1 m2 f 2      32  2      0       0    1           8   
9       9    f1 m2 f 2      21  1      0       1    1           9   
9       92   f1 m2 f 2      21  1      0       0    1           9   
9       174  f1 m2 f 2      21  1      1       0    0           9   
10      10   f1 m2 f 2      28  2      0       1    0           10  
10      93   f1 m2 f 2      28  2      0       0    1           10  
10      175  f1 m2 f 2      28  2      0       0    1           10  
11      11   f1 m2 f 2      29  2      1       1    0           11  
11      94   f1 m2 f 2      29  2      0       0    0           11  
11      176  f1 m2 f 2      29  2      0       0    1           11  
12      12   f1 m2 f 2      37  4      2       1    1           12  
12      95   f1 m2 f 2      37  4      1       0    1           12  
12      177  f1 m2 f 2      37  4      0       0    1           12  
13      13   f1 m2 f 2      31  1      1       1    0           13  
13      96   f1 m2 f 2      31  1      0       0    0           13  
13      178  f1 m2 f 2      31  1      0       0    0           13  
14      14   f1 m2 f 2      29  3      2       1    0           14  
14      97   f1 m2 f 2      29  3      0       0    1           14  
14      179  f1 m2 f 2      29  3      0       0    2           14  
15      15   f1 m2 f 2      31  2      1       1    1           15  
15      98   f1 m2 f 2      31  2      1       0    0           15  
15      180  f1 m2 f 2      31  2      1       0    0           15  
16      16   f1 m2 f 2      27  2      2       1    0           16  
16      99   f1 m2 f 2      27  2      1       0    0           16  
16      181  f1 m2 f 2      27  2      0       0    0           16  
17      17   f1 m2 f 2      30  5      2       1    1           17  
17      100  f1 m2 f 2      30  5      0       0    2           17  
17      182  f1 m2 f 2      30  5      1       0    2           17  
18      18   f1 m2 f 2      26  1      0       1    1           18  
18      101  f1 m2 f 2      26  1      0       0    0           18  
18      183  f1 m2 f 2      26  1      1       0    0           18  
19      19   f1 m2 f 2      25  3      2       1    1           19  
19      102  f1 m2 f 2      25  3      0       0    1           19  
19      184  f1 m2 f 2      25  3      1       0    1           19  
20      20   f1 m2 f 2      44  1      0       1    1           20  
20      103  f1 m2 f 2      44  1      0       0    0           20  
20      185  f1 m2 f 2      44  1      1       0    0           20  
21      21   f1 m2 f 2      40  1      0       1    1           21  
21      104  f1 m2 f 2      40  1      0       0    0           21  
21      186  f1 m2 f 2      40  1      0       0    0           21  
22      22   f1 m2 f 2      35  2      2       1    0           22  
22      105  f1 m2 f 2      35  2      0       0    0           22  
22      187  f1 m2 f 2      35  2      0       0    0           22  
23      23   f1 m2 f 2      28  2      0       1    2           23  
23      106  f1 m2 f 2      28  2      0       0    0           23  
23      188  f1 m2 f 2      28  2      2       0    0           23  
24      24   f1 m2 f 2      36  1      0       1    1           24  
24      107  f1 m2 f 2      36  1      0       0    0           24  
24      189  f1 m2 f 2      36  1      0       0    1           24  
25      25   f1 m2 f 2      27  2      1       1    1           25  
25      108  f1 m2 f 2      27  2      0       0    1           25  
25      190  f1 m2 f 2      27  2      0       0    2           25  
26      26   f1 m2 f 2      40  2      0       1    2           26  
26      109  f1 m2 f 2      40  2      0       0    0           26  
26      191  f1 m2 f 2      40  2      0       0    0           26  
27      27   f1 m2 f 2      38  2      0       1    2           27  
27      110  f1 m2 f 2      38  2      0       0    0           27  
27      192  f1 m2 f 2      38  2      0       0    0           27  
28      28   f1 m2 f 2      34  3      0       1    2           28  
28      111  f1 m2 f 2      34  3      0       0    0           28  
28      193  f1 m2 f 2      34  3      0       0    0           28  
29      29   f1 m2 f 2      28  4      1       1    2           29  
29      112  f1 m2 f 2      28  4      0       0    2           29  
29      194  f1 m2 f 2      28  4      2       0    1           29  
30      30   f1 m2 f 2      30  4      2       1    0           30  
30      113  f1 m2 f 2      30  4      1       0    1           30  
30      195  f1 m2 f 2      30  4      1       0    1           30  
31      31   f1 m2 f 2      32  1      0       1    1           31  
31      114  f1 m2 f 2      32  1      0       0    0           31  
31      196  f1 m2 f 2      32  1      0       0    0           31  
32      32   f1 m2 f 2      34  2      1       1    0           32  
32      115  f1 m2 f 2      34  2      1       0    0           32  
32      197  f1 m2 f 2      34  2      0       0    0           32  
33      33   f1 m2 f 2      42  1      1       1    0           33  
33      116  f1 m2 f 2      42  1      1       0    0           33  
33      198  f1 m2 f 2      42  1      0       0    0           33  
34      34   f1 m2 f 2      32  2      0       1    2           34  
34      117  f1 m2 f 2      32  2      0       0    0           34  
34      199  f1 m2 f 2      32  2      2       0    0           34  
35      35   f1 m2 f 2      39  1      1       1    0           35  
35      118  f1 m2 f 2      39  1      0       0    0           35  
35      200  f1 m2 f 2      39  1      0       0    0           35  
36      36   f1 m2 f 2      35  2      0       1    2           36  
36      119  f1 m2 f 2      35  2      0       0    0           36  
36      201  f1 m2 f 2      35  2      0       0    0           36  
37      37   f1 m2 f 2      36  1      0       1    1           37  
37      120  f1 m2 f 2      36  1      0       0    0           37  
37      202  f1 m2 f 2      36  1      0       0    0           37  
38      38   f1 m2 f 2      34  3      1       1    2           38  
38      121  f1 m2 f 2      34  3      2       0    0           38  
38      203  f1 m2 f 2      34  3      2       0    0           38  
39      39   f1 m2 f 2      30  3      0       1    0           39  
39      122  f1 m2 f 2      30  3      0       0    2           39  
39      204  f1 m2 f 2      30  3      0       0    1           39  
40      40   f1 m2 f 2      28  1      0       1    1           40  
40      123  f1 m2 f 2      28  1      1       0    0           40  
40      205  f1 m2 f 2      28  1      0       0    0           40  
41      41   f1 m2 f 2      39  3      0       1    2           41  
41      124  f1 m2 f 2      39  3      1       0    0           41  
41      206  f1 m2 f 2      39  3      0       0    0           41  
42      42   f1 m2 f 2      35  1      0       1    0           42  
42      125  f1 m2 f 2      35  1      0       0    0           42  
42      207  f1 m2 f 2      35  1      0       0    0           42  
43      43   f1 m2 f 2      41  1      0       1    0           43  
43      126  f1 m2 f 2      41  1      0       0    0           43  
43      208  f1 m2 f 2      41  1      0       0    0           43  
44      44   f1 m2 f 2      37  2      1       1    1           44  
44      127  f1 m2 f 2      37  2      0       0    0           44  
44      209  f1 m2 f 2      37  2      0       0    0           44  
45      45   f1 m2 f 3      30  1      0       1    0           45  
45      128  f1 m2 f 3      30  1      1       0    0           45  
45      210  f1 m2 f 3      30  1      0       0    0           45  
46      46   f1 m2 f 3      37  1      1       1    0           46  
46      129  f1 m2 f 3      37  1      0       0    0           46  
46      211  f1 m2 f 3      37  1      0       0    1           46  
47      47   f1 m2 f 3      28  2      0       1    2           47  
47      130  f1 m2 f 3      28  2      1       0    0           47  
47      212  f1 m2 f 3      28  2      1       0    0           47  
48      48   f1 m2 f 3      27  4      2       1    0           48  
48      131  f1 m2 f 3      27  4      2       0    1           48  
48      213  f1 m2 f 3      27  4      2       0    0           48  
49      49   f1 m2 f 3      26  2      2       1    0           49  
49      132  f1 m2 f 3      26  2      1       0    0           49  
49      214  f1 m2 f 3      26  2      1       0    0           49  
50      50   f1 m2 f 3      38  3      0       1    2           50  
50      133  f1 m2 f 3      38  3      1       0    0           50  
50      215  f1 m2 f 3      38  3      1       0    0           50  
51      51   f1 m2 f 3      24  3      1       1    2           51  
51      134  f1 m2 f 3      24  3      2       0    1           51  
51      216  f1 m2 f 3      24  3      2       0    0           51  
52      52   f1 m2 f 3      36  5      1       1    2           52  
52      135  f1 m2 f 3      36  5      1       0    1           52  
52      217  f1 m2 f 3      36  5      2       0    1           52  
53      53   f1 m2 f 3      27  3      1       1    1           53  
53      136  f1 m2 f 3      27  3      1       0    1           53  
53      218  f1 m2 f 3      27  3      2       0    0           53  
54      54   f1 m2 f 3      28  1      0       1    1           54  
54      137  f1 m2 f 3      28  1      1       0    0           54  
54      219  f1 m2 f 3      28  1      0       0    1           54  
55      55   f1 m2 f 3      29  2      0       1    2           55  
55      138  f1 m2 f 3      29  2      1       0    0           55  
55      220  f1 m2 f 3      29  2      1       0    1           55  
56      56   f1 m2 f 3      36  2      0       1    2           56  
56      139  f1 m2 f 3      36  2      1       0    0           56  
56      221  f1 m2 f 3      36  2      0       0    1           56  
57      57   f1 m2 f 3      28  2      1       1    0           57  
57      140  f1 m2 f 3      28  2      1       0    1           57  
57      222  f1 m2 f 3      28  2      2       0    0           57  
58      58   f1 m2 f 3      28  2      0       1    2           58  
58      141  f1 m2 f 3      28  2      2       0    0           58  
58      223  f1 m2 f 3      28  2      1       0    0           58  
59      59   f1 m2 f 3      28  1      0       1    1           59  
59      142  f1 m2 f 3      28  1      1       0    0           59  
59      224  f1 m2 f 3      28  1      0       0    0           59  
60      60   f1 m2 f 3      27  2      0       1    2           60  
60      143  f1 m2 f 3      27  2      1       0    0           60  
60      225  f1 m2 f 3      27  2      1       0    0           60  
61      61   f1 m2 f 3      35  2      0       1    2           61  
61      144  f1 m2 f 3      35  2      2       0    0           61  
61      226  f1 m2 f 3      35  2      1       0    0           61  
62      62   f1 m2 f 3      25  1      0       1    1           62  
62      145  f1 m2 f 3      25  1      1       0    0           62  
62      227  f1 m2 f 3      25  1      1       0    0           62  
63      63   f1 m2 f 3      34  1      0       1    1           63  
63      146  f1 m2 f 3      34  1      0       0    0           63  
63      228  f1 m2 f 3      34  1      0       0    0           63  
64      64   f1 m2 f 3      31  2      0       1    2           64  
64      147  f1 m2 f 3      31  2      0       0    0           64  
64      229  f1 m2 f 3      31  2      1       0    0           64  
65      65   f1 m2 f 3      26  2      1       1    0           65  
65      148  f1 m2 f 3      26  2      0       0    1           65  
65      230  f1 m2 f 3      26  2      0       0    2           65  
66      66   f1 m2 f 3      32  1      0       1    1           66  
66      149  f1 m2 f 3      32  1      0       0    0           66  
66      231  f1 m2 f 3      32  1      1       0    0           66  
67      67   f1 m2 f 3      21  1      0       1    1           67  
67      150  f1 m2 f 3      21  1      0       0    1           67  
67      232  f1 m2 f 3      21  1      0       0    0           67  
68      68   f1 m2 f 3      28  3      1       1    2           68  
68      151  f1 m2 f 3      28  3      2       0    0           68  
68      233  f1 m2 f 3      28  3      2       0    0           68  
69      69   f1 m2 f 3      37  3      0       1    2           69  
69      152  f1 m2 f 3      37  3      1       0    1           69  
69      234  f1 m2 f 3      37  3      0       0    2           69  
70      70   f1 m2 f 3      25  1      1       1    0           70  
70      153  f1 m2 f 3      25  1      0       0    0           70  
70      235  f1 m2 f 3      25  1      1       0    0           70  
71      71   f1 m2 f 3      32  1      1       1    0           71  
71      154  f1 m2 f 3      32  1      1       0    0           71  
71      236  f1 m2 f 3      32  1      0       0    0           71  
72      72   f1 m2 f 3      25  1      0       1    1           72  
72      155  f1 m2 f 3      25  1      0       0    0           72  
72      237  f1 m2 f 3      25  1      1       0    0           72  
73      73   f1 m2 f 3      31  1      0       1    1           73  
73      156  f1 m2 f 3      31  1      0       0    1           73  
73      238  f1 m2 f 3      31  1      0       0    0           73  
74      74   f1 m2 f 3      38  6      0       1    2           74  
74      239  f1 m2 f 3      38  6      0       0    2           74  
75      75   f1 m2 f 3      26  2      0       1    2           75  
75      157  f1 m2 f 3      26  2      0       0    2           75  
75      240  f1 m2 f 3      26  2      1       0    1           75  
76      76   f1 m2 f 3      31  1      0       1    1           76  
76      158  f1 m2 f 3      31  1      0       0    0           76  
76      241  f1 m2 f 3      31  1      1       0    0           76  
77      77   f1 m2 f 3      31  2      0       1    1           77  
77      159  f1 m2 f 3      31  2      2       0    0           77  
77      242  f1 m2 f 3      31  2      0       0    1           77  
78      78   f1 m2 f 3      25  1      1       1    0           78  
78      160  f1 m2 f 3      25  1      0       0    0           78  
78      243  f1 m2 f 3      25  1      0       0    1           78  
79      79   f1 m2 f 3      31  1      0       1    1           79  
79      161  f1 m2 f 3      31  1      0       0    0           79  
79      244  f1 m2 f 3      31  1      0       0    1           79  
80      80   f1 m2 f 3      34  1      0       1    1           80  
80      162  f1 m2 f 3      34  1      0       0    0           80  
80      245  f1 m2 f 3      34  1      0       0    0           80  
81      81   f1 m2 f 3      35  2      2       1    0           81  
81      163  f1 m2 f 3      35  2      0       0    0           81  
81      246  f1 m2 f 3      35  2      2       0    0           81  
82      82   f1 m2 f 3      29  1      0       1    1           82  
82      164  f1 m2 f 3      29  1      0       0    1           82  
82      247  f1 m2 f 3      29  1      0       0    1           82  
83      83   f1 m2 f 3      23  1      0       1    1           83  
83      165  f1 m2 f 3      23  1      0       0    1           83  
83      248  f1 m2 f 3      23  1      0       0    1           83  
;;;;
run
recode case
head
clr case = spontaneous induced
