/*
-------------------------------------------------------------------------------------
Description of variables in nickel dataset:
ID (identification number)
ICDCODE (160=death due to nasal sinus cancer)
EXPOSURE (exposure level)
DOB (Date of birth)
AGE_EMP (Age at first employment)
AGE_STRT (Age at start of study)
AGE_STOP (Age at death or end of study, whichever is earlier)
SUBCO20 (1=included in 20% subcohort, 0=not included)
Data are from Breslow and Day "Statistical Methods in Cancer
Research Volume 2: The Design and Analysis of Cohort Studies",
1987, IARC Publications Number 82, Lyon France, p 375-390.
--------------------------------------------------------------------------------------
*/
data nickel20;
input id icdcode exposure dob age_emp age_strt age_stop subco20;
yr_emp=floor(dob+age_emp);
if icdcode=160 then nasal=1; else nasal=2;
logafe=log(age_emp-10); yfe10=(dob+age_emp-1915)/10;
yfe100=(dob+age_emp-1915)**2/100;
logexp=log(exposure+1);
/* Time = years since first employment */
timestrt=age_strt-age_emp; timestop=min(age_stop,85)-age_emp;
if id=0 & icdcode=490 then id=990;
if id=0 & icdcode=177 then id=991;
if id=0 & icdcode=430 then id=992;
if id=0 & icdcode=0 then id=993;
cards;
3 0 5 1889.0192 17.4808 45.2273 92.9808 0
4 162 5 1885.978 23.1864 48.2684 63.2712 0
6 163 10 1881.2548 25.2452 52.9917 54.1644 0
8 527 9 1886.3397 24.7206 47.9067 69.6794 0
9 150 0 1879.5 29.9575 54.7465 76.8442 0
10 163 2 1889.915 21.2877 44.3314 62.5413 0
15 334 0 1890.5 23.2836 43.7465 62 0
16 160 0.5 1874.3315 50.3569 59.9149 65.5835 0
17 420 0 1909.5 15.4618 34.7465 50.5137 0
18 12 0 1892.5 24.1394 51.7465 57.5931 1
19 160 10 1881.7261 33.452 52.5204 58.4897 0
21 14 0 1877.8 30.2026 56.4464 65.5698 0
22 177 2.5 1879.5 28.6694 54.7465 67.4808 1
23 162 0 1900.5 22.6479 33.7465 52.7028 0
26 999 0 1887.474 28.1973 46.7725 77.056 0
27 420 0 1893.8494 28.9615 40.3971 68.5205 0
28 420 0 1905.8494 14.052 28.3971 61.5287 0
29 434 0 1873.5 36.7439 60.7465 63.0082 0
30 420 1 1899.4 15.4 34.8464 66.5808 0
32 153 0 1905.5 14.4014 28.7465 53.7684 0
33 151 1.5 1883.5 27.7959 50.7465 61.8096 0
34 154 4 1886.5 20.9274 47.7465 53.2561 0
35 420 0.5 1886.3945 23.5233 47.8519 82.8274 0
37 162 0 1889.0931 34.4905 45.1533 63.0435 0
40 162 1 1894.6383 28.7507 39.6082 74.5425 0
41 331 2.5 1886.9672 27.8191 47.2793 65.1094 1
42 420 0 1899.4987 24.1945 34.7478 83.0493 0
44 998 0 1889.5 34.959 44.7465 56.478 0
47 160 1 1893.3124 18.3945 45.9341 47.1139 1
48 592 4 1884.5 26.8342 49.7465 63.6339 0
49 162 2.5 1890.9698 20.2274 43.2766 62.7891 0
50 420 3.5 1896.5 26.3602 37.7465 54.9685 0
51 151 1 1883.5 22.6973 50.7465 70 0
52 420 0.5 1894.7479 20.7521 39.4985 86.6877 0
53 160 1 1900.7103 22.8925 33.5361 72.7609 1
54 160 11 1889.3972 25.2877 44.8492 77.8082 0
55 0 0 1900.9863 22.726 33.2601 81.0137 1
56 502 2 1904.5 18.2917 29.7465 51.5847 0
57 27 0.5 1882.5 34.5603 51.7465 58.4344 1
58 241 0 1872.5 45.6398 61.7465 66.7684 0
59 420 0 1901.2246 22.1644 33.0219 81.0193 0
60 422 5 1889.5 21.626 44.7465 60.5713 0
61 420 0 1879.5 35.8835 54.7465 76.2014 0
62 592 6.5 1879.5 34.1521 54.7465 73.4316 0
63 160 14 1878.0302 30.4098 56.2163 58.0654 0
64 334 6 1888.5 18.8069 45.7465 64.8124 0
65 502 6 1883.5 28.2946 50.7465 64.4808 0
68 162 5 1897.8494 16.6506 36.3971 79.8301 0
70 162 0 1901.2739 21.5699 32.9725 54.7479 1
71 163 0.5 1895.989 17.6466 38.2574 61.4932 1
73 163 8.5 1887.2931 25.1877 46.9534 78.7151 1
75 0 0 1905.4082 16.8795 28.8383 76.5918 0
76 332 7.5 1879.5 33.4208 54.7465 72.0562 0
77 420 0 1890.5 32.3521 43.7465 74.2732 0
78 153 0.5 1894.915 28.0986 39.3314 72.389 0
79 334 1 1902.863 21.5386 31.3834 71.0876 1
80 2 2 1889.5 33.3109 44.7465 64.1383 0
82 527 3 1876.5 32.235 57.7465 60.9795 1
83 332 2 1880.5 30.7767 53.7465 76.0219 0
84 0 0 1895.4795 27.5425 38.767 86.5205 0
85 151 6.5 1881.411 33.3917 52.8354 86.3807 1
86 163 21 1885.5 22.2671 48.7465 59.6123 0
87 163 4 1882.5 20.7849 51.7465 58.489 1
88 502 0 1904.0684 14.222 30.1781 71.4495 0
89 502 3 1875.5 43.4315 58.7465 68.5055 1
90 491 0 1900.9535 19.5957 33.293 76.8821 1
91 154 3.5 1886.5 23.8809 47.7465 54.5055 0
92 163 5.5 1895.1753 24.1563 39.0712 52.4768 0
93 443 0 1888.3032 20.4509 45.9432 81.0885 0
94 420 5.5 1881.5 25.5192 52.7465 75.1093 0
95 410 0 1874.5 35.7576 59.7465 65.8716 0
96 420 0 1890.022 32.7972 44.2245 89.9452 0
97 0 0 1902.5 21.7897 31.7465 79.5 1
98 162 3.5 1894.1589 28.8082 40.0875 61.8219 0
99 420 0 1883.5 31.7684 50.7465 62.3932 1
100 331 5 1874.5 31.6589 59.7465 60.5493 0
101 420 0 1904.5273 14.0316 44.7191 75.5901 1
103 491 0 1901.7343 21.9342 32.5122 56.3315 0
104 420 17.5 1883.5 26.0205 50.7465 73.5548 0
106 160 11 1880.3115 28.4426 53.9349 58.072 0
107 420 0 1897.6548 25.4055 36.5917 79.8466 0
108 153 2.5 1880.5 35.653 53.7465 65.4452 0
109 592 5.5 1899.5 20.6666 34.7465 45.2595 0
110 0 0 1909.2466 14.0302 24.9999 72.7534 0
111 0 1 1905.2219 17.4247 29.0245 76.7781 0
112 491 0 1897.452 25.5726 36.7944 78.578 1
113 998 7.5 1880.5 27.1 53.7465 61.111 0
115 331 5.5 1888.5 25.1521 45.7465 66.6561 0
116 162 1.5 1901.1918 22.3672 33.0547 55.0896 0
117 502 1 1891.6521 24.3643 42.5944 66.3369 0
118 493 7 1887.5 26.1493 46.7465 57.7849 0
119 501 11.5 1881.5 26.3356 52.7465 58.2836 0
120 502 7 1873.5 47.1311 60.7465 69.4342 0
121 163 14 1882.011 27.5562 52.2355 68.578 1
122 162 0 1893.2384 30.3013 41.0081 82.8053 0
124 331 0 1871.5 36 62.7465 84.724 0
126 162 18 1883.5 24.7103 50.7465 57.2787 0
127 792 1 1900.5 22.3192 33.7465 62 1
128 163 0 1886.5 20.9247 47.7465 53.6366 0
129 450 2 1897.126 25.8302 37.1205 74.7014 1
131 163 18 1881.9589 26.675 52.2876 67.1973 0
132 332 0 1900.929 22.0875 33.3175 69.745 0
135 491 0 1892.1093 31.2085 42.1372 84.8934 0
137 465 0 1905.3864 19.0808 28.8601 74.2794 0
138 181 3.5 1882.5 24 51.7465 71.8726 0
140 443 0 1887.5 36.2616 46.7465 71.2369 1
141 163 0.5 1891.5 18.5027 42.7465 62.2534 0
142 160 7.5 1883.915 26.7124 50.3314 57.4493 0
143 420 3 1887.5 32.9508 46.7465 77.3032 1
144 160 1.5 1884.5055 28.5958 49.741 67.2288 0
145 177 10.5 1882.7288 27.0383 51.5177 91.8274 0
147 420 4 1878.5 30 55.7465 87.911 0
148 420 0 1902.5 21.2041 31.7465 65.6584 0
149 163 2 1889.9973 33.3068 44.2491 69.0355 1
150 160 9.5 1871.9315 38.6904 62.3149 64.7324 0
151 160 10.5 1875.3644 32.8296 58.8821 77.7452 0
152 420 4.5 1874.5 40.2616 59.7465 83.0726 0
154 463 2.5 1900.1503 22.6552 44.0962 70.0004 0
155 434 0 1880.5 32.7192 53.7465 58.2369 0
156 196 12.5 1887.5 27.1603 46.7465 58.837 0
157 332 0.5 1905.7507 16.3314 28.4957 74.752 0
158 163 5.5 1885.5 22.5 48.7465 54.5819 0
159 502 3.5 1886.5 23.7466 47.7465 55.6123 0
161 162 2.5 1891.5726 31.2767 42.6738 60.5531 0
162 157 0 1897.5 26.089 36.7465 46.9181 0
164 163 3 1896.5 24.4126 37.7465 52.6096 1
165 502 2 1891.5 31.437 42.7465 50.7822 0
166 451 2.5 1901.063 21.6411 33.1835 67.5107 0
170 420 0 1893.6191 29.4329 40.6273 78.7415 0
172 163 3 1888.6257 26.8743 45.6207 58.5222 1
173 0 3 1895.9287 26.9207 38.3177 86.0713 0
175 332 0 1887.5 26 46.7465 76.2946 0
176 198 4.5 1882.7754 33.3093 51.4711 68.4329 0
177 420 0 1902.4137 10.7836 31.8328 67.8494 0
180 163 3 1892.5 21.3109 41.7465 48.1174 0
181 526 2 1881.5 24.7576 52.7465 62.6011 0
186 163 1.5 1897.9287 25.6521 36.3177 61.0959 0
187 999 0 1906.5 18.0984 27.7465 34.5 0
190 163 5 1888.4727 25.0273 45.7738 57.9877 0
191 181 4 1886.5 35.0918 47.7465 69.8142 0
194 446 1 1894.0411 28.8849 40.2053 71.715 0
195 334 0 1892.5 30.3384 41.7465 70 1
196 491 3.5 1892.5629 27.9371 41.6836 81.8317 0
197 356 2.5 1902.5 20.4727 31.7465 53.9399 0
199 163 0 1880.3606 24.1257 53.8859 68.1257 1
200 0 7.5 1885.8657 20.548 48.3807 96.1343 0
202 541 0.5 1898.5 24.3575 35.7465 42.1011 0
203 163 1 1886.6493 36.4713 47.5972 71.8795 0
204 502 6.5 1876.5 30.9137 57.7465 66.2507 0
206 162 0 1892.5 30.5795 41.7465 64.5192 1
207 162 4 1895.6904 18.074 38.556 60.3123 0
208 163 1.5 1887.0273 26.4727 47.2191 67.2111 0
210 163 0 1890.5068 30.0287 43.7396 55.011 0
211 163 5 1886.5 31.9219 47.7465 50.03 0
212 163 1.5 1895.5 18 38.7465 49.6233 0
213 162 0 1910.1288 14.2018 24.1177 63.0493 0
214 420 1.5 1886.6548 29.1425 47.5917 78.559 1
215 434 0 1885.5 25.111 48.7465 54.3657 0
216 420 0 1881.5 23.4645 52.7465 66.6694 1
218 491 10 1883.5 30.5027 50.7465 78.7794 0
219 151 0.5 1884.3907 38.7299 49.8557 83.1819 1
221 163 11 1888.9181 21.8271 45.3284 68.0272 0
222 160 3.5 1876.0847 37.7372 58.1617 70.9098 0
223 162 0 1889.5 33.8671 44.7465 63.7576 1
224 434 0 1903.5 13.7739 30.7465 37.6453 0
225 177 0 1883.8494 25.6082 50.3971 83.7698 1
226 260 0 1883.5 21 50.7465 67.9987 0
227 0 0 1904.3442 12.5847 29.9022 77.6558 0
228 502 0 1906.9124 15.6713 27.3341 71.6603 0
229 420 0 1901.5 23.3634 32.7465 62.735 1
230 0 0 1903.5 13.0629 30.7465 78.5 0
231 162 3 1896.5 21.3301 37.7465 59.5874 0
234 420 2 1882.5 27.6506 51.7465 71.4562 0
236 161 0 1892.5 21.889 41.7465 68.2814 1
237 163 5 1895.5 27.3219 38.7465 52.4205 0
239 420 0 1878.5 26.0684 55.7465 80.8315 0
240 0 1 1902.5 20.5328 31.7465 79.5 0
243 998 2 1887.5 26.1603 46.7465 51.5055 0
245 422 1.5 1886.5 23.9302 47.7465 73.5874 0
246 163 0 1894.5 20.7877 39.7465 55.6425 0
248 420 1 1895.6219 28.4354 38.6245 69.5397 0
250 163 10 1885.5 24.8424 48.7465 59.8972 0
251 153 3.5 1877.5 34.7732 56.7465 72.9712 0
252 420 1 1888.5 17.5795 45.7465 58.4534 0
255 0 0 1899.1616 23.9425 35.0848 82.8384 0
256 223 2 1880.5 39.9563 53.7465 68.235 0
257 490 0 1899.5 23.3685 34.7465 62.9904 0
259 154 2.5 1885.5 37.426 48.7465 76.489 0
261 0 0 1901.3835 21.9288 32.8629 80.6165 1
265 999 0 1902.5 13.4945 31.7465 55.8917 1
266 331 0 1887.5 34.5137 46.7465 61.0601 0
267 2 3 1881.5 38.9344 52.7465 65.9302 0
268 420 1 1894.5 20.7904 39.7465 88.4124 0
270 163 2 1895.8 27.0876 38.4464 53.4301 1
271 160 9 1884.5192 30.5493 49.7273 64.9905 0
272 420 5.5 1882.5 38.4153 51.7465 68.6343 0
273 154 0 1909.5 14.6913 24.7465 55.9302 0
274 491 2 1883.4384 36.9769 50.8081 75.0658 0
275 331 0 1880.5 44.1121 53.7465 54.1548 0
276 153 1 1900.3524 22.6696 33.894 69.464 1
277 163 0 1903.5 21.4016 30.7465 65 1
279 163 1 1880.5 27.653 53.7465 63.3082 0
281 420 0 1893.2849 26.9391 40.9615 76.5315 0
282 162 6 1891.5 22.5055 42.7465 54.0535 0
283 500 0 1894.7151 29 39.5314 57.5609 1
284 431 6.5 1880.5 26.4507 53.7465 56.3771 0
285 420 0 1899.063 15.7206 35.1835 68.3342 1
286 421 5 1886.5 23.6836 47.7465 63.8315 0
287 331 6 1878.5 28 55.7465 84.0643 0
289 162 0 1884.9508 35.7705 49.2957 76.9424 0
291 502 8 1874.5 28.3411 59.7465 66.0847 0
292 422 3.5 1871.5 37.1366 62.7465 86.2233 0
293 610 8.5 1882.5 28.1219 51.7465 56.9302 0
294 160 5.5 1885.2986 35.2698 48.9479 62.0959 0
295 196 0 1889.0192 35.448 45.2273 74.937 0
297 420 5 1886.5 27.1329 47.7465 70.8232 0
298 160 12 1885.8219 25.2301 48.4246 54.7 0
299 160 0 1898.5 23.4315 35.7465 68.6534 0
300 163 5 1900.1039 20.295 34.1426 46.7509 0
301 525 8 1889.5 22.5819 44.7465 51.3224 1
302 593 0 1900.5082 15.2316 33.7383 75.4205 0
303 422 7.5 1879.5 31.7959 54.7465 67.6534 0
305 334 3 1899.8712 16.6288 34.3752 73.9069 1
306 999 0 1897.5 25.3219 36.7465 37.3 0
307 422 0 1877.5 27.3524 56.7465 85.2151 1
309 331 3.5 1900.8907 14.5311 33.3557 68.5394 0
310 422 12 1878.5 27.6013 55.7465 76.4342 0
313 0 0.5 1896.5629 19.196 37.6836 85.4371 0
315 420 0 1879.5 24.1849 54.7465 79.6206 0
316 491 0 1889.2246 33.7974 45.0219 91.9152 0
318 541 0.5 1899.5 17 34.7465 43.7466 0
319 160 0 1900.7377 22.1226 33.5088 66.3198 1
321 0 0 1897.4329 25.5836 36.8136 84.5671 0
322 592 0 1900.5 22.6753 33.7465 55.059 0
324 332 0 1885.2493 37.5863 48.9972 83.8685 0
325 434 0 1910.5 14.0737 23.7465 43.0343 0
326 502 3 1879.5 43.3109 54.7465 64.7814 1
327 465 1.5 1903.9617 13.6548 30.2848 73.0575 1
328 331 1.5 1884.5 20.7959 49.7465 66.4452 1
329 163 6 1888.5 23.7897 45.7465 51.3439 0
330 160 0 1881.5781 28.8986 52.6683 70.0795 0
332 420 0 1884.3279 38.5598 49.9186 69.3104 0
334 162 6 1892.5 28.3989 46.7465 71.9426 0
335 440 6 1884.5 26.6781 54.7465 78.4672 0
336 420 12 1876.5 28 62.7465 83.8087 1
337 420 0 1885.5 36.8151 48.7465 67.9932 0
338 0 0 1894.9287 29.4019 49.3177 87.0713 1
339 331 5 1882.5 28.4589 51.7465 72.448 0
340 502 0 1894.2849 28.7261 39.9615 66.5675 0
341 154 0 1888.5 31.6858 45.7465 51.4808 0
344 420 1 1908.8907 14.5586 25.3557 68.29 0
345 162 2 1900.8279 22.0269 33.4186 54.3502 1
346 163 4.5 1894.2028 28.7917 40.0437 58.2017 0
348 0 2 1901.3287 21.7919 32.9177 80.6713 0
349 332 0 1899.2493 23.9672 34.9972 75.4137 0
351 331 6 1873.5 49.3137 65.7465 73.2261 0
353 204 0 1901.5 23.1694 32.7465 48.7631 0
354 420 0 1877.1151 31.3439 57.1313 97.7781 0
357 420 0 1896.8005 22.6995 37.4459 75.6011 0
360 450 0 1883.9069 30.052 50.3396 86.2411 0
361 443 0 1887.5836 31.6301 46.6628 81.8383 0
362 331 0.5 1889.1808 21.4329 45.0657 88.2576 0
363 160 10 1883.4219 26.5973 50.8246 57.7891 0
364 345 3 1886.5 21.6968 47.7465 60.6644 0
365 420 0 1900.2021 22.6691 34.0443 60.4153 1
368 420 2 1888.5629 22.9371 50.6836 81.2865 0
369 420 0 1904.3279 19.3872 29.9186 70.3762 0
370 420 1 1899.6191 14.5946 34.6273 76.9628 1
371 600 8 1886.5 28.5876 47.7465 57.765 0
375 592 1 1885.5 26.7897 48.7465 60.2261 0
376 332 0 1907.0466 15.0355 27.1998 68.8383 0
377 153 9.5 1881.5 27.6644 52.7465 80.2123 0
378 160 3.5 1893.5 19.911 40.7465 63.974 1
380 411 0 1907.1534 14.7397 27.093 53.6635 0
381 160 13 1884.5956 29.4592 49.6509 78.0867 0
383 204 0 1903.0959 19.926 31.1505 72.5753 0
384 431 7 1877.5 37.4617 56.7465 64.5328 0
385 0 0 1903.5343 16.5996 30.7122 78.4657 0
386 451 5 1901.674 15.3287 32.5725 78.3752 0
388 434 4 1892.5 30.3164 41.7465 49.5466 1
389 420 0.5 1903 20.3041 31.2465 58.8547 1
390 420 2.5 1886.5042 22.9958 52.7423 85.7909 0
391 162 13 1881.5 28.8917 52.7465 60.6918 0
392 160 15 1891.0192 21.5054 43.2273 57.5847 0
394 237 0 1904.571 14.0043 29.6754 65.829 0
395 160 3 1896.1639 22.8552 38.0825 70.8059 1
397 332 0 1890.5 32.5603 43.7465 46.4098 1
400 162 0 1903.3699 19.4247 30.8766 77.0673 0
403 502 7 1887.5 35.3301 46.7465 48.5956 0
404 155 0 1902.9041 20.5781 31.3424 77.5549 0
406 163 5.5 1892.5 21.1986 41.7465 66.7083 1
407 422 3 1883.5 39.3602 50.7465 66.5658 0
408 443 1 1882.5 25.5055 51.7465 74.4808 1
409 163 2 1874.5 38.4044 59.7465 70.8069 0
410 162 0 1900.5 22.5247 33.7465 51.6613 0
412 420 0 1880.5 35.1823 53.7465 78.2151 1
416 162 0 1893.8109 29.7699 40.4355 61.6849 0
417 162 7 1875.5 38.8069 58.7465 69.459 0
419 155 0 1902.5397 20.4548 31.7068 70.8082 0
420 420 3 1894.5 24.8232 39.7465 71.3384 1
421 160 1 1884.5 30 49.7465 69.7822 0
423 502 0 1874.5 31.7 59.7465 74.522 0
424 332 0 1891.7507 31.041 42.4957 72.7958 1
425 163 3 1897.5 13.3164 36.7465 55.5137 1
427 422 3.5 1888.3361 26.9708 45.9104 67.4174 0
428 162 0 1896.5 26.8781 37.7465 45.2123 0