data list free file="ont_males_by_csd2.prn" /
 geog ttl ag65 to ag75.

/* The following commands assume a file consisting of CD-level  */
/* records with the exception of 3528, 3548, 3557, 3558, and 3560 */
/* and CSD-level records for all CSDs in                        */
/* 3528, 3548, 3557, 3558, and 3560                 */

compute phu06=0.
if (any(geog,3501,3502)) phu06 eq 3558.
if (geog eq 3506) phu06 eq 3551.
if (any(geog,3507,3509)) phu06 eq 3543.
if (any(geog,3510,3511)) phu06 eq 3541.
if (any(geog,3512,3513)) phu06 eq 3538.
if (any(geog,3514,3516,3546)) phu06 eq 3535.
if (geog eq 3515) phu06 eq 3555.
if (geog eq 3518) phu06 eq 3530.
if (geog eq 3519) phu06 eq 3570.
if (geog eq 3520) phu06 eq 3595.
if (geog eq 3521) phu06 eq 3553.
if (any(geog,3522,3523)) phu06 eq 3566.
if (geog eq 3524) phu06 eq 3536.
if (geog eq 3525) phu06 eq 3537.
if (geog eq 3526) phu06 eq 3546.
if (geog eq 3529) phu06 eq 3527.
if (geog eq 3530) phu06 eq 3565.
if (geog eq 3531) phu06 eq 3554.
if (geog eq 3532) phu06 eq 3552.
if (geog eq 3534) phu06 eq 3531.
if (geog eq 3536) phu06 eq 3540.
if (geog eq 3537) phu06 eq 3568.
if (geog eq 3538) phu06 eq 3542.
if (geog eq 3539) phu06 eq 3544.
if (geog eq 3540) phu06 eq 3539.
if (any(geog,3541,3542)) phu06 eq 3533.
if (any(geog,3543,3544)) phu06 eq 3560.
if (geog eq 3547) phu06 eq 3557.
if (geog eq 3549) phu06 eq 3547.
if (any(geog,3551,3552,3553)) phu06 eq 3561.
if (geog eq 3554) phu06 eq 3563.
if (geog eq 3556) phu06 eq 3556.
if (geog eq 3559) phu06 eq 3549.

if (any(geog,3528018,3528052)) phu06 eq 3534.
if (any(geog,3528035,23528037)) phu06 eq 3527.

if (any(geog,3548001,3548091)) phu06 eq 3557.
if (any(geog,3548069,3548072)) phu06 eq 3563.
if (any(geog,3548013,3548019,3548021,3548022)) phu06 eq 3547.
if (any(geog,3548027,3548031,3548034,3548044)) phu06 eq 3547.
if (any(geog,3548055,3548073,3548094)) phu06 eq 3547.

if ((geog eq 3557096)) phu06 eq 3556.
if (any(geog,3557001,3557004,3557006,3557008)) phu06 eq 3526.
if (any(geog,3557011,3557014,3557016,3557019)) phu06 eq 3526.
if (any(geog,3557021,3557026,3557028,3557035)) phu06 eq 3526.
if (any(geog,3557038,3557039,3557040,3557041)) phu06 eq 3526.
if (any(geog,3557051,3557061,3557066,3557071)) phu06 eq 3526.
if (any(geog,3557072,3557073,3557074,3557075)) phu06 eq 3526.
if (any(geog,3557076,3557077,3557078,3557079)) phu06 eq 3526.
if (any(geog,3557082,3557091,3557094,3557095)) phu06 eq 3526. 

if (any(geog,3558001,3558003,3558004,3558011)) phu06 eq 3562.
if (any(geog,3558012,3558016,3558019,3558028)) phu06 eq 3562.
if (any(geog,3558034,3558041,3558044,3558051)) phu06 eq 3562.
if (any(geog,3558054,3558059,3558060,3558061)) phu06 eq 3562.
if (any(geog,3558062,3558063,3558064,3558065)) phu06 eq 3562.
if (any(geog,3558066,3558067,3558068,3558069)) phu06 eq 3562.
if (any(geog,3558075,3558076,3558090,3558095)) phu06 eq 3562.
if (any(geog,3558097,3558100)) phu06 eq 3562.
if (any(geog,3558080,3558085)) phu06 eq 3549.

if (any(geog,3560053,3560072,3560074,3560075)) phu06 eq 3562.
if (any(geog,3560078,3560079,3560081,3560085)) phu06 eq 3562.
if (any(geog,3560086,3560088,3560090,3560093)) phu06 eq 3562.
if (any(geog,3560096,3560098,3560100)) phu06 eq 3562.
if (any(geog,3560050,3560051,3560052,3560091)) phu06 eq 3556.
if (any(geog,3560001,3560004,3560005,3560007)) phu06 eq 3549.
if (any(geog,3560008,3560010,3560021,3560024)) phu06 eq 3549.
if (any(geog,3560027,3560032,3560034,3560042)) phu06 eq 3549.
if (any(geog,3560046,3560049)) phu06 eq 3549.
if (any(geog,3560054,3560055,3560056,3560057)) phu06 eq 3549.
if (any(geog,3560058,3560059,3560060,3560061)) phu06 eq 3549.
if (any(geog,3560062,3560063,3560064,3560065)) phu06 eq 3549.
if (any(geog,3560066,3560067,3560068,3560069)) phu06 eq 3549.
if (any(geog,3560070,3560071,3560076,3560077)) phu06 eq 3549.
if (any(geog,3560080,3560082,3560083,3560084)) phu06 eq 3549.
if (any(geog,3560087,3560089,3560090,3560095)) phu06 eq 3549.
if (any(geog,3560097,3560102,3560104)) phu06 eq 3549.
frequencies variables=phu06.
descriptives variables=geog to ag75.
execute.
temporary
select if (phu06 eq 0).
list variables=geog
execute.
save outfile="ont_by_phu06.sav" / keep=all.
execute
finish.