/*variable name in program Variable name in Gator diffcohh MAR_IMPUTE_FIRST_COHAB=1 1st cohabitation date does not match up the MAR_COHABITATION array diffmar MAR_IMPUTE_FIRST_MARRY=1 1st marriage data does not match up the MAR_COHABITATION array The date of first divorce: divorceM CVC_FIRST_DIVORCE_DATE~M Month of first divorce divorceY CVC_FIRST_DIVORCE_DATE~Y year of first divorce divor_CM CVC_FIRST_DIVORCE_MONTH continuous month of first divorce end_mar CVC_FIRST_MARRY_END FIRST MARRIAGE END REASON The CVC_FIRST_MARRY_END classifies the reason to end the first marriage. The variable is classified in three basic marital categories: 1) Divorced 2) Widowed 3) Unknown. marstat CV_MARSTAT cmarstat CV_MARSTAT_COllAPSED cohm CVC_FIRST_COHAB_DATE_M cohy CVC_FIRST_COHAB_DATE_Y marm CVC_FIRST_MARRY_DATE_M mary CVC_FIRST_MARRY_DATE_Y cohcm CVC_FIRST_COHAB_MONTH marcm CVC_FIRST_MARRY_MONTH cohh_ttl CVC_COHAB_TTl mar_ttl CVC_MARRIAGES_TTl Codes for Created Variables Marital/Cohabitation Status 1 = never married, cohabiting 2 = never married, not cohabiting 3 = married, spouse present 4 = married, spouse absent 5 = separated, cohabiting 6 = separated, not cohabiting 7 = divorced, cohabiting 8 = divorced, not cohabiting 9 = widowed, cohabiting 10 = widowed, not cohabiting Collapsed Marital Status 0 = never married 1 = married 2 = separated 3 = divorced 4 = widowed Marital Status in event history 0 = never married, not cohabiting 1 = never married, cohabiting 2 = married 3 = separated 4 = divorced 5 = widowed Variables Used /************************************************************************** YMAR-10100.01.01.01~M M10100111M YMAR-10100.01.01.01~Y M10100111Y YMAR-10800.01.01 M1080011 YMAR-10800.02.01 M1080021 YMAR-10900.01.01 M1090011 YMAR-11400.01 M114001 YMAR-11400.02 M114002 YMAR-11400.03 M114003 YMAR-11700.01~M M117001M YMAR-11700.01~Y M117001Y YMAR-11700.02~M M117002M YMAR-11700.02~Y M117002Y YMAR-13150~M M13150M YMAR-13150~Y M13150Y YMAR-13250~M M13250M YMAR-13250~Y M13250Y YMAR-3100.01~M M31001M YMAR-3100.01~Y M31001Y YMAR-3100.02~M M31002M YMAR-3100.02~Y M31002Y YMAR-3100.03~M M31003M YMAR-3100.03~Y M31003Y YMAR-3100A.01~M M3100A1M YMAR-3100A.01~Y M3100A1Y YMAR-3100A.02~M M3100A2M YMAR-3100A.02~Y M3100A2Y YMAR-3100C.01~M M3100C1M YMAR-3100C.01~Y M3100C1Y YMAR-3100C.02~M M3100C2M YMAR-3100C.02~Y M3100C2Y YMAR-4350A.01 M4350A1 YMAR-4350A.02 M4350A2 YMAR-4360A.01~M M4360A1M YMAR-4360A.01~Y M4360A1Y YMAR-4500.01 M45001 YMAR-4500.02 M45002 YMAR-4500.03 M45003 YMAR-4650.01 M46501 YMAR-4650.02 M46502 YMAR-4650.03 M46503 YMAR-4670.01 M46701 YMAR-4670.02 M46702 YMAR-4700.01 M47001 YMAR-4700.02 M47002 YMAR-4700.03 M47003 YMAR-4800.01~M M48001M YMAR-4800.01~Y M48001Y YMAR-4800.02~M M48002M YMAR-4800.02~Y M48002Y YMAR-4800.03~M M48003M YMAR-4800.03~Y M48003Y YMAR-4900.01 M49001 YMAR-4900.02 M49002 YMAR-4900.03 M49003 YMAR-5400.01.01 M540011 YMAR-5400.01.02 M540012 YMAR-5600.01.01~M M560011M YMAR-5600.01.01~Y M560011Y YMAR-5600.01.02~M M560012M YMAR-5600.01.02~Y M560012Y YMAR-5700.01.01~M M570011M YMAR-5700.01.01~Y M570011Y YMAR-5700.01.02~M M570012M YMAR-5700.01.02~Y M570012Y YMAR-620 M620 YMAR-630 M630 YMAR-6300.01 M63001 YMAR-6300.02 M63002 YMAR-6300.03 M63003 YMAR-640 M640 YMAR-6400.01 M64001 YMAR-6400.02 M64002 YMAR-6400.03 M64003 YMAR-6400A.01 M6400A1 YMAR-6400A.02 M6400A2 YMAR-6400A.03 M6400A3 YMAR-650 M650 YMAR-7000.01.01~M M700011M YMAR-7000.01.01~Y M700011Y YMAR-7000.01.02~M M700012M YMAR-7000.01.02~Y M700012Y YMAR-7000.02.01~M M700021M YMAR-7000.02.01~Y M700021Y YMAR-714J.01 M714J1 YMAR-714JB.01~M M714JB1M YMAR-714JB.01~Y M714JB1Y YMAR-716J.01 M716J1 YMAR-716J.02 M716J2 YMAR-716JB.02~M M716JB2M YMAR-716JB.02~Y M716JB2Y YMAR-717J.01 M717J1 YMAR-717JB.01~M M717JB1M YMAR-717JB.01~Y M717JB1Y YMAR-7300.01.01 M730011 YMAR-7300.02.01 M730021 YMAR-7300.03.01 M730031 YMAR-735 M735 YMAR-736 M736 YMAR-740R M740R YMAR-740RA M740RA YMAR-751 M751 YMAR-7900.01.01.01 M7900111 YMAR-7900.01.01.02 M7900112 YMAR-7900.02.01.01 M7900211 YMAR-7900.03.01.01 M7900311 YMAR-8100.01.01.01~M M8100111M YMAR-8100.01.01.01~Y M8100111Y YMAR-8100.01.01.02~M M8100112M YMAR-8100.01.01.02~Y M8100112Y YMAR-8200.01.01.01~M M8200111M YMAR-8200.01.01.01~Y M8200111Y YMAR-8200.01.01.02~M M8200112M YMAR-8200.01.01.02~Y M8200112Y YMAR-8200.02.01.01~M M8200211M YMAR-8200.02.01.01~Y M8200211Y YMAR-8200.03.01.01~M M8200311M YMAR-8200.03.01.01~Y M8200311Y YMAR-8400.01.01.01 M8400111 YMAR-8400.01.01.02 M8400112 YMAR-8400.02.01.01 M8400211 YMAR-9100.01.01 M910011 YMAR-9100.01.02 M910012 YMAR-9100.02.01 M910021 YMAR-9200.01.01~M M920011M YMAR-9200.01.01~Y M920011Y YMAR-9300.01.01 M930011 YMAR-9300.01.02 M930012 YMAR-9300.02.01 M930021 PARTNERS_UID.01 UID1 PARTNERS_UID.02 UID2 PARTNERS_UID.03 UID3 PARTNERS_ID.01 P1 PARTNERS_ID.02 P2 PARTNERS_ID.03 P3 /* BDATE_M KEY!BDATE_M_1997 */ /* BDATE_Y KEY!BDATE_Y_1997 */ /* P1R1 PARTNERS_ID.01_1997 */ /* P2R1 PARTNERS_ID.02_1997 */ /* P1R3 PARTNERS_ID.01_1999 */ /* P2R3 PARTNERS_ID.02_1999 */ /* P1R2 PARTNERS_ID.01_1998 */ /* P2R2 PARTNERS_ID.02_1998 */ /* P3R2 PARTNERS_ID.03_1998 */ /* P1R4 PARTNERS_ID.01_2000 */ /* P2R4 PARTNERS_ID.02_2000 */ /* P3R4 PARTNERS_ID.03_2000 */ /* P4R4 PARTNERS_ID.04_2000 */ /* P1R5 PARTNERS_ID.01_2001 */ /* P2R5 PARTNERS_ID.02_2001 */ /* P3R5 PARTNERS_ID.03_2001 */ /* P4R5 PARTNERS_ID.04_2001 */ /* P1R6 PARTNERS_ID.01_2002 */ /* P2R6 PARTNERS_ID.02_2002 */ /* P3R6 PARTNERS_ID.03_2002 */ /* P1R7 PARTNERS_ID.01_2003 */ /* P2R7 PARTNERS_ID.02_2003 */ /* P3R7 PARTNERS_ID.03_2003 */ /* P1R8 PARTNERS_ID.01_2004 */ /* P2R8 PARTNERS_ID.02_2004 */ /* P3R8 PARTNERS_ID.03_2004 */ /* P1R9 PARTNERS_ID.01_2005 */ /* P2R9 PARTNERS_ID.02_2005 */ /* P3R9 PARTNERS_ID.03_2005 */ /* P4R9 PARTNERS_ID.04_2005 */ /* P1R10 PARTNERS_ID.01_2006 */ /* P2R10 PARTNERS_ID.02_2006 */ /* P3R10 PARTNERS_ID.03_2006 */ /* P1R11 PARTNERS_ID.01_2007 */ /* P2R11 PARTNERS_ID.02_2007 */ /* P3R11 PARTNERS_ID.03_2007 */ /* P1R12 PARTNERS_ID.01_2008 */ /* P2R12 PARTNERS_ID.02_2008 */ /* P3R12 PARTNERS_ID.03_2008 */ /* P1R13 PARTNERS_ID.01_2009 */ /* P2R13 PARTNERS_ID.02_2009 */ /* P3R13 PARTNERS_ID.03_2009 */ /* P1R14 PARTNERS_ID.01_2010 */ /* P2R14 PARTNERS_ID.02_2010 */ /* P3R14 PARTNERS_ID.03_2010 */ /* P1R15 PARTNERS_ID.01_2011 */ /* P2R15 PARTNERS_ID.02_2011 */ /* P3R15 PARTNERS_ID.03_2011 */ /* P1R16 PARTNERS_ID.01_2013 */ /* P2R16 PARTNERS_ID.02_2013 */ /* P3R16 PARTNERS_ID.03_2013 */ /* P1R17 PARTNERS_ID.01_2015 */ /* P2R17 PARTNERS_ID.02_2015 */ /* P3R17 PARTNERS_ID.03_2015 */ /* P4R17 PARTNERS_ID.04_2015 */ /* P1R18 PARTNERS_ID.01_2017 */ /* P2R18 PARTNERS_ID.02_2017 */ /* P3R18 PARTNERS_ID.03_2017 */ /* P4R18 PARTNERS_ID.04_2017 */ /* int_cmr1 CV_INTERVIEW_CMONTH 1997 */ /* INT_DR1 CV_INTERVIEW_DATE 1997 */ /* INT_MR1 CV_INTERVIEW_DATE 1997 */ /* INT_YR1 CV_INTERVIEW_DATE 1997 */ /* MARSTATR1 CV_MARSTAT 1997 */ /* CMARSTAR1 CV_MARSTAT_COLLAPSED 1997 */ /* int_cmr2 CV_INTERVIEW_CMONTH 1998 */ /* INT_DR2 CV_INTERVIEW_DATE 1998 */ /* INT_MR2 CV_INTERVIEW_DATE 1998 */ /* INT_YR2 CV_INTERVIEW_DATE 1998 */ /* MARSTATR2 CV_MARSTAT 1998 */ /* CMARSTAR2 CV_MARSTAT_COLLAPSED 1998 */ /* int_cmr3 CV_INTERVIEW_CMONTH 1999 */ /* INT_DR3 CV_INTERVIEW_DATE 1999 */ /* INT_MR3 CV_INTERVIEW_DATE 1999 */ /* INT_YR3 CV_INTERVIEW_DATE 1999 */ /* MARSTATR3 CV_MARSTAT 1999 */ /* CMARSTAR3 CV_MARSTAT_COLLAPSED 1999 */ /* int_cmr4 CV_INTERVIEW_CMONTH 2000 */ /* INT_DR4 CV_INTERVIEW_DATE 2000 */ /* INT_MR4 CV_INTERVIEW_DATE 2000 */ /* INT_YR4 CV_INTERVIEW_DATE 2000 */ /* MARSTATR4 CV_MARSTAT 2000 */ /* CMARSTAR4 CV_MARSTAT_COLLAPSED 2000 */ /* int_cmr5 CV_INTERVIEW_CMONTH 2001 */ /* INT_DR5 CV_INTERVIEW_DATE 2001 */ /* INT_MR5 CV_INTERVIEW_DATE 2001 */ /* INT_YR5 CV_INTERVIEW_DATE 2001 */ /* MARSTATR5 CV_MARSTAT 2001 */ /* CMARSTAR5 CV_MARSTAT_COLLAPSED 2001 */ /* int_cmr6 CV_INTERVIEW_CMONTH 2002 */ /* INT_DR6 CV_INTERVIEW_DATE 2002 */ /* INT_MR6 CV_INTERVIEW_DATE 2002 */ /* INT_YR6 CV_INTERVIEW_DATE 2002 */ /* MARSTATR6 CV_MARSTAT 2002 */ /* CMARSTAR6 CV_MARSTAT_COLLAPSED 2002 */ /* int_cmr7 CV_INTERVIEW_CMONTH 2003 */ /* INT_DR7 CV_INTERVIEW_DATE 2003 */ /* INT_MR7 CV_INTERVIEW_DATE 2003 */ /* INT_YR7 CV_INTERVIEW_DATE 2003 */ /* MARSTATR7 CV_MARSTAT 2003 */ /* CMARSTAR7 CV_MARSTAT_COLLAPSED 2003 */ /* int_cmr8 CV_INTERVIEW_CMONTH 2004 */ /* INT_DR8 CV_INTERVIEW_DATE 2004 */ /* INT_MR8 CV_INTERVIEW_DATE 2004 */ /* INT_YR8 CV_INTERVIEW_DATE 2004 */ /* MARSTATR8 CV_MARSTAT 2004 */ /* CMARSTAR8 CV_MARSTAT_COLLAPSED 2004 */ /* int_cmr9 CV_INTERVIEW_CMONTH 2005 */ /* INT_DR9 CV_INTERVIEW_DATE 2005 */ /* INT_MR9 CV_INTERVIEW_DATE 2005 */ /* INT_YR9 CV_INTERVIEW_DATE 2005 */ /* MARSTATR9 CV_MARSTAT 2005; */ /* CMARSTAR9 CV_MARSTAT_COLLAPSED 2005 */ /* int_cmr10 CV_INTERVIEW_CMONTH 2006 */ /* INT_DR10 CV_INTERVIEW_DATE 2006 */ /* INT_MR10 CV_INTERVIEW_DATE 2006 */ /* INT_YR10 CV_INTERVIEW_DATE 2006 */ /* MARSTATR10 CV_MARSTAT 2006 */ /* CMARSTAR10 CV_MARSTAT_COLLAPSED 2006 */ /* int_cmR11 CV_INTERVIEW_CMONTH 2007 */ /* INT_DR11 CV_INTERVIEW_DATE 2007 */ /* INT_MR11 CV_INTERVIEW_DATE 2007 */ /* INT_YR11 CV_INTERVIEW_DATE 2007 */ /* MARSTATR11 CV_MARSTAT 2007 */ /* CMARSTAR11 CV_MARSTAT_COLLAPSED 2007 */ /* int_cmR12 CVC_INTERVIEW_CMONTH 2008 */ /* INT_DR12 CV_INTERVIEW_DATE 2008 */ /* INT_MR12 CV_INTERVIEW_DATE 2008 */ /* INT_YR12 CV_INTERVIEW_DATE 2008 */ /** MARSTATR12 CV_MARSTAT 2008 */ /* CMARSTAR12 CV_MARSTAT_COLLAPSED 2008 */ /* int_cmR13 CVC_INTERVIEW_CMONTH 2009 */ /* INT_DR13 CV_INTERVIEW_DATE 2009 */ /* INT_MR13 CV_INTERVIEW_DATE 2009 */ /* INT_YR13 CV_INTERVIEW_DATE 2009 */ /* MARSTATR13 CV_MARSTAT 2009 */ /* CMARSTAR13 CV_MARSTAT_COLLAPSED 2009 */ /* int_cmR14 CVC_INTERVIEW_CMONTH 2010 */ /* INT_DR14 CV_INTERVIEW_DATE 2010 */ /* INT_MR14 CV_INTERVIEW_DATE 2010 */ /* INT_YR14 CV_INTERVIEW_DATE 2010 */ /* MARSTATR14 CV_MARSTAT 2010 */ /* CMARSTAR14 CV_MARSTAT_COLLAPSED 2010 */ /* int_cmR15 CVC_INTERVIEW_CMONTH 2011 */ /* INT_DR15 CV_INTERVIEW_DATE 2011 */ /* INT_MR15 CV_INTERVIEW_DATE 2011 */ /* INT_YR15 CV_INTERVIEW_DATE 2011 */ /* MARSTATR15 CV_MARSTAT 2011 */ /* CMARSTAR15 CV_MARSTAT_COLLAPSED 2011 */ /* int_cmR16 CVC_INTERVIEW_CMONTH 2013 */ /* INT_DR16 CV_INTERVIEW_DATE 2013 */ /* INT_MR16 CV_INTERVIEW_DATE 2013 */ /* INT_YR16 CV_INTERVIEW_DATE 2013 */ /* MARSTATR16 CV_MARSTAT 2013 */ /* CMARSTAR16 CV_MARSTAT_COLLAPSED 2013 */ /* int_cmR17 CVC_INTERVIEW_CMONTH 2015 */ /* INT_DR17 CV_INTERVIEW_DATE 2015 */ /* INT_MR17 CV_INTERVIEW_DATE 2015 */ /* INT_YR17 CV_INTERVIEW_DATE 2015 */ /* MARSTATR17 CV_MARSTAT 2015 */ /* CMARSTAR17 CV_MARSTAT_COLLAPSED 2015 */ /* int_cmR18 CVC_INTERVIEW_CMONTH 2017 */ /* INT_DR18 CV_INTERVIEW_DATE 2017 */ /* INT_MR18 CV_INTERVIEW_DATE 2017 */ /* INT_YR18 CV_INTERVIEW_DATE 2017 */ /* MARSTATR18 CV_MARSTAT 2017 */ /* CMARSTAR18 CV_MARSTAT_COLLAPSED 2017 */ /* M620R1 YMAR-620_1998 */ /* M620R3 YMAR-620_1999 */ /* M620R4 YMAR-620_2000 */ /* M620R5 YMAR-620_2001 */ /* M620R6 YMAR-620_2002 */ /* M620R7 YMAR-620_2003 */ /* M620R8 YMAR-620_2004 */ /* M620R9 YMAR-620_2005 */ /* M620R10 YMAR-620_2006 */ /* M620R11 YMAR-620_2007 */ /* M620R12 YMAR-620_2008 */ /* M620R13 YMAR-620_2009 */ /* M620R14 YMAR-620_2010 */ /* M620R15 YMAR-620_2011 */ /* M620R16 YMAR-620_2013 */ /* M620R17 YMAR-620_2015 */ /* M620R18 YMAR-620_2017 */ /* cohm CVC_FIRST_COHAB_DATE_M */ /* cohy CVC_FIRST_COHAB_DATE_Y */ /* marm CVC_FIRST_MARRY_DATE_M */ /* mary CVC_FIRST_MARRY_DATE_Y */ /* cohcm CVC_FIRST_COHAB_MONTH */ /* marcm CVC_FIRST_MARRY_MONTH */ /* coh_ttl CVC_COHAB_TTl */ /* mar_ttl CVC_MARRIAGES_TTl */ /* divorceM CVC_first_divorce_date_M */ /* divorceY CVC_first_divorce_date_Y */ /* divor_CM CVC_first_divorce_month */ /* end_mar CVC_FIRST_MARRY_END */ options pagesize=54 linesize=93 nocenter nodate; libname marR18; libname marR19; data one; set marR19.inC_marR19; data two; set marR19.inP_marR19; data five; merge one two; by pubid; if int_y>0 then do; %let m=18; array MARSTATr(&m); array CMARSTAr(&m); array int_mr(&m) ; array int_yr(&m) ; array int_cmr(&m); array m620r(&m); array partner {4,18} p1r1-p1R18 p2r1-p2R18 p3r1-p3R18 p4r1-p4R18; array Hspouse14 {18} spouse1-spouse18; array HpartneR14 {18} partner1-partneR18; array dlipar {4,18} p1r1-p1R18 p2r1-p2R18 p3r1-p3R18 p4r1-p4R18 ; array R18mars (l1) R18mars001-R18mars487; array R18co (l1) R18co001-R18co487; array R18link (l1) R18link001-R18link487; array R19pt (l) R19pt001-R19pt487; array R19sp (l) R19sp001-R19sp487; do i=1 to (&m); if int_cmr(i)>0 then do; dlicm=int_cmr(i); dlimarstat=MARSTATr(i); dlicmarstat=CMARSTAr(i); dlim620=m620r(i); mardli=MARSTATr(i); end; end; drop i; if precoh_ttl<-3 then precoh_ttl=0; if premar_ttl<-3 then premar_ttl=0; if int_y >0 & int_m >0 then do; doicm=(int_y-1980)*12+int_m; end; else if int_y = -5 then do; doicm=-5; end; start=dlicm+1; if m13150y>0 & m13150m>0 then m13150cm=(m13150y -1980)*12 +m13150m; if m13250y>0 & m13250m>0 then m13250cm=(m13250y -1980)*12 +m13250m; if m13150cm>dlicm then do; m7900111=1; m8200111m= m13150m; m8200111y= m13150y; end; if m31001M>0 & m8200111m>0 & m31001M> m8200111m & M45001=1 & m31001y=m8200111y then do; m8200111m=m31001m; flagsame_m=1; end; if m31001y>0 & m8200111y>0 & m31001y > M8200111y & m45001=1 then do; m8200111m=m31001m; m8200111y=m31001y; flagsame_m=2; end; if m31001y>0 & m8200111y>0 & m31001y < M8200111y & m45001=1 then do; m45001=0; flagsame_m=3; end; array m4800m (i) m48001m m48002m m48003m m48004m ; array m4800y (i) m48001y m48002y m48003y m48004y; array m4800cm (i) m48001cm m48002cm m48003cm m48004cm; array m3100m (i) m31001m m31002m m31003m m31004m; array m3100y (i) m31001y m31002y m31003y m31004y; array m3100cm (i) m31001cm m31002cm m31003cm m31004cm; array m3100_CM (i) m3100C1m m3100C2m m3100C3m m3100C4m; array m3100_CY (i) m3100C1y m3100C2y m3100C3y m3100C4y; do i=1 to 4; if m3100_CM>-4 & m3100_CY>-4 then do; m3100m= m3100_CM; m3100y= m3100_CY; end; end; flag4800=0; flag3100=0; do i=1 to 4; if m4800m >0 & m4800y>0 then do; M4800cm= (M4800y-1980)*12+m4800m; if 00 & m3100y>0 then do; M3100cm= (M3100y-1980)*12+m3100m; if m3100cm<=dlicm then flag3100=1; end; end; array par{4} p1-p4; flagp=0; do t=1 to 4; do i=1 to 4; do j=1 to 18; if par(t)=partner{i,j}& partner{i,j}>0 then flagp=1; end; end; end; array maf (l) maf001-maf487; array mas (l) mas001-mas487; array mat (l) mat001-mat487; array mafour (l) mafour001-mafour487; array mbf (l) mbf001-mbf487; array mbs (l) mbs001-mbs487; array mbt (l) mbt001-mbt487; array mbfour (l) mbfour001-mbfour487; array mcf (l) mcf001-mcf487; array mcs (l) mcs001-mcs487; array mct (l) mct001-mct487; array mcfour (l) mcfour001-mcfour487; array mdf (l) mdf001-mdf487; array mds (l) mds001-mds487; array mdt (l) mdt001-mdt487; array mdfour (l) mdfour001-mdfour487; array maa (l) maa001-maa487; array mbb (l) mbb001-mbb487; array mcc (l) mcc001-mcc487; array mdd (l) mdd001-mdd487; array mal (l) mal001-mal487; array cohaf (l) cohaf001-cohaf487; array cohas (l) cohas001-cohas487; array cohat (l) cohat001-cohat487; array cohafour (l) cohafour001-cohafour487; array cohbf (l) cohbf001-cohbf487; array cohbs (l) cohbs001-cohbs487; array cohbt (l) cohbt001-cohbt487; array cohbfour (l) cohbfour001-cohbfour487; array cohcf (l) cohcf001-cohcf487; array cohcs (l) cohcs001-cohcs487; array cohct (l) cohct001-cohct487; array cohcfour (l) cohcfour001-cohcfour487; array cohdf (l) cohdf001-cohdf487; array cohds (l) cohds001-cohds487; array cohdt (l) cohdt001-cohdt487; array cohdfour (l) cohdfour001-cohdfour487; array cohaa (l) cohaa001-cohaa487; array cohbb (l) cohbb001-cohbb487; array cohcc (l) cohcc001-cohcc487; array cohdd (l) cohdd001-cohdd487; array cohal (l) cohal001-cohal487; array mars (l) mars001-mars487; array cohh (l) cohh001-cohh487; array link (l) link001-link487; array dual (l) dual001-dual487; do l=start to 487; cohal=0; cohaf=0; cohas=0; cohat=0; cohafour=0; cohbf=0; cohbs=0; cohbt=0; cohbfour=0; cohcf=0; cohcs=0; cohct=0; cohcfour=0; cohdf=0; cohds=0; cohdt=0; cohdfour=0; mal=0 ; maf=0 ; mas=0 ; mat=0 ; mafour=0 ; mbf=0 ; mbs=0 ; mbt=0 ; mbfour=0 ; mcf=0 ; mcs=0 ; mct=0 ; mcfour=0 ; mdf=0 ; mds=0 ; mdt=0 ; mdfour=0 ; end; array ysca (t) M31001Y M700011Y M700012Y ysca4; array yscb (t) M31002Y M700021Y yscb3 yscb4; array yscc (t) M31003Y yscc2 yscc3 yscc4; array yscd (t) M31004Y yscd2 yscd3 yscd4; array msca (t) M31001M M700011M M700012M msca4; array mscb (t) M31002M M700021M mscb3 mscb4; array mscc (t) M31003M mscc2 mscc3 mscc4; array mscd (t) M31004M mscd2 mscd3 mscd4; array ysc (p) ysca yscb yscc yscd; array msc (p) msca mscb mscc mscd; array ycoa (t) M31001Y M700011Y M700012Y ycoa4; array ycob (t) M31002Y M700021Y ycob3 ycob4; array ycoc (t) M31003Y ycoc2 ycoc3 ycoc4; array ycod (t) M31004Y ycod2 ycod3 ycod4; array mcoa (t) M31001M M700011M M700012M mcoa4; array mcob (t) M31002M M700021M mcob3 mcob4; array mcoc (t) M31003M mcoc2 mcoc3 mcoc4; array mcod (t) M31004M mcod2 mcod3 mcod4; array yco (p) ycoa ycob ycoc ycod; array mco (p) mcoa mcob mcoc mcod; array ymcoa (t) ymcoa1 ymcoa2 ymcoa3 ymcoa4; array ymcob (t) ymcob1 ymcob2 ymcob3 ymcob4; array ymcoc (t) ymcoc1 ymcoc2 ymcoc3 ymcoc4; array ymcod (t) ymcod1 ymcod2 ymcod3 ymcod4; array ymco (p) ymcoa ymcob ymcoc ymcod; array yspa (t) M3100A1Y yspa2 yspa3 yspa4; array yspb (t) M3100A2Y yspb2 yspb3 yspb4; array yspc (t) M3100A3Y yspc2 yspc3 yspc4; array yspd (t) yspd1 yspd2 yspd3 yspd4; array mspa (t) M3100A1M mspa2 mspa3 mspa4; array mspb (t) M3100A2M mspb2 mspb3 mspb4; array mspc (t) M3100A3M mspc2 mspc3 mspc4; array mspd (t) mspd1 mspd2 mspd3 mspd4; array ysp (p) yspa yspb yspc yspd; array msp (p) mspa mspb mspc mspd; array ymspa (t) ymspa1 ymspa2 ymspa3 ymspa4; array ymspb (t) ymspb1 ymspb2 ymspb3 ymspb4; array ymspc (t) ymspc1 ymspc2 ymspc3 ymspc4; array ymspd (t) ymspd1 ymspd2 ymspd3 ymspd4; array ymsp (p) ymspa ymspb ymspc ymspd; array yBsca (t) M3100C1Y yBscb2 yBsca3 yBsca4; array yBscb (t) yBscb1 yBscb2 yBscb3 yBscb4; array yBscc (t) yBscc1 yBscc2 yBscc3 yBscc4; array yBscd (t) yBscd1 yBscd2 yBscd3 yBscd4; array mBsca (t) M3100C1M mBsca2 mBsca3 mBsca4; array mBscb (t) mBscb1 mBscb2 mBscb3 mBscb4; array mBscc (t) mBscc1 mBscc2 mBscc3 mBscc4; array mBscd (t) mBsca1 mBscd2 mBscd3 mBscd4; array yBsc (p) yBsca yBscb yBscc yBscd; array mBsc (p) mBsca mBscb mBscc mBscd; array yeca (t) M48001Y M920011Y M920012Y yeca4; array yecb (t) M48002Y M920021Y yecb3 yecb4; array yecc (t) M48003Y yecc2 yecc3 yecc4; array yecd (t) M48004Y yecd2 yecd3 yecd4; array yec (p) yeca yecb yecc yecd; array meca (t) M48001M M920011M M920012M meca4; array mecb (t) M48002M M920021M mecb3 mecb4; array mecc (t) M48003M mecc2 mecc3 mecc4; array mecd (t) M48004M mecd2 mecd3 mecd4; array mec (p) meca mecb mecc mecd; array csmca (t) csmca1 csmca2 csmca3 csmca4; array csmcb (t) csmcb1 csmcb2 csmcb3 csmcb4; array csmcc (t) csmcc1 csmcc2 csmcc3 csmcc4; array csmcd (t) csmcd1 csmcd2 csmcd3 csmcd4; array csmc (p) csmca csmcb csmcc csmcd; array cemca (t) cemca1 cemca2 cemca3 cemca4; array cemcb (t) cemcb1 cemcb2 cemcb3 cemcb4; array cemcc (t) cemcc1 cemcc2 cemcc3 cemcc4; array cemcd (t) cemcd1 cemcd2 cemcd3 cemcd4; array cemc (p) cemca cemcb cemcc cemcd; array csmma (t) csmma1 csmma2 csmma3 csmma4; array csmmb (t) csmmb1 csmmb2 csmmb3 csmmb4; array csmmc (t) csmmc1 csmmc2 csmmc3 csmmc4; array csmmd (t) csmmd1 csmmd2 csmmd3 csmmd4; array csmm (p) csmma csmmb csmmc csmmd; array csmsa (t) csmsa1 csmsa2 csmsa3 csmsa4; array csmsb (t) csmsb1 csmsb2 csmsb3 csmsb4; array csmsc (t) csmsc1 csmsc2 csmsc3 csmsc4; array csmsd (t) csmsd1 csmsd2 csmsd3 csmsd4; array csms (p) csmsa csmsb csmsc csmsd; array cemma (t) cemma1 cemma2 cemma3 cemma4; array cemmb (t) cemmb1 cemmb2 cemmb3 cemmb4; array cemmc (t) cemmc1 cemmc2 cemmc3 cemmc4; array cemmd (t) cemmd1 cemmd2 cemmd3 cemmd4; array cemm (p) cemma cemmb cemmc cemmd; array cemsa (t) cemsa1 cemsa2 cemsa3 cemsa4; array cemsb (t) cemsb1 cemsb2 cemsb3 cemsb4; array cemsc (t) cemsc1 cemsc2 cemsc3 cemsc4; array cemsd (t) cemsd1 cemsd2 cemsd3 cemsd4; array cems (p) cemsa cemsb cemsc cemsd; array cdlia (t) M46701 cdlia2 cdlia3 cdlia4; array cdlib (t) M46702 cdlib2 cdlib3 cdlib4; array cdlic (t) cdlic1 cdlic2 cdlic3 cdlic4; array cdlid (t) cdlid1 cdlid2 cdlid3 cdlid4; array cdli (p) cdlia cdlib cdlic cdlid; array cbega (t) M47001 cbega2 cbega3 cbega4; array cbegb (t) M47002 cbegb2 cbegb3 cbegb4; array cbegc (t) M47003 cbegc2 cbegc3 cbegc4; array cbegd (t) M47004 cbegd2 cbegd3 cbegd4; array cbeg (p) cbega cbegb cbegc cbegd; array cbeg2a (t) cbeg2a1 M910011 M910012 cbeg2a4; array cbeg2b (t) cbeg2b1 M910021 cbeg2b3 cbeg2b4; array cbeg2c (t) cbeg2c1 cbeg2c2 cbeg2c3 cbeg2c4; array cbeg2d (t) cbeg2d1 cbeg2d2 cbeg2d3 cbeg2d4; array cbeg2 (p) cbeg2a cbeg2b cbeg2c cbeg2d; array mbega (t) M45001 mbega2 mbega3 mbega4; array mbegb (t) M45002 mbegb2 mbegb3 mbegb4; array mbegc (t) M45003 mbegc2 mbegc3 mbegc4; array mbegd (t) M45004 mbegd2 mbegd3 mbegd4; array mbeg (p) mbega mbegb mbegc mbegd; array mdlia (t) mardli mdlia2 mdlia3 mdlia4; array mdlib (t) mdlib1 mdlib2 mdlib3 mdlib4; array mdlic (t) mdlic1 mdlic2 mdlic3 mdlic4; array mdlid (t) mdlid1 mdlid2 mdlid3 mdlid4; array mdli (p) mdlia mdlib mdlic mdlid; array cstata1 (t) m540011 m540012 cstata13 cstata14; array cstata2 (t) m7900111 m7900112 cstata23 cstata24; array cstata3 (t) m9800111 m9800112 cstata33 cstata34; array cstatb1 (t) m7900211 cstatb12 cstatb13 cstatb13; array cstatb2 (t) m540021 cstatb22 cstatb23 cstatb24; array cstatc2 (t) m540031 m540032 m540033 cstatbc4; array cstat1 (p) cstata1 cstatb1 cstatc1 cstatd1; array cstat2 (p) cstata2 cstatb2 cstatc2 cstatd2; array cstat3 (p) cstata3 cstatb3 cstatc3 cstatd3; array cstat4 (p) cstata4 cstatb4 cstatc4 cstatd4; array cstat (i) cstat1 cstat2 cstat3 cstat4 ; array cyma1 (t) M570011y cyma12 cyma13 cyma14; array cmma1 (t) M570011m cmma12 cmma13 cmma14; array cyma2 (t) m8200111y cyma22 cyma23 cyma24; array cmma2 (t) m8200111m cmma22 cmma23 cmma24; array cyma3 (t) m10100111y cyma32 cyma33 cyma34; array cmma3 (t) m10100111m cmma32 cmma33 cmma34; array cymb1 (t) m8200211y cymb12 cymb13 cymb14; array cmmb1 (t) m8200211m cmmb12 cmmb13 cmmb14; array cymb2 (t) m570021y cymb22 cymb23 cymb24; array cmmb2 (t) m570021m cmmb22 cmmb23 cmmb24; array cymc2 (t) cymc21 m570032y cymc23 cymc24; array cmmc2 (t) cmmc21 m570032m cmmc23 cmmc24; array cym1 (p) cyma1 cymb1 cymc1 cymd1; array cmm1 (p) cmma1 cmmb1 cmmc1 cmmd1; array cym2 (p) cyma2 cymb2 cymc2 cymd2; array cmm2 (p) cmma2 cmmb2 cmmc2 cmmd2; array cym3 (p) cyma3 cymb3 cymc3 cymd3; array cmm3 (p) cmma3 cmmb3 cmmc3 cmmd3; array cym4 (p) cyma4 cymb4 cymc4 cymd4; array cmm4 (p) cmma4 cmmb4 cmmc4 cmmd4; array cym (i) cym1 cym2 cym3 cym4; array cmm (i) cmm1 cmm2 cmm3 cmm4; array cyssa1 (t) M560011Y M560012Y cyssa13 cyssa14; array cmssa1 (t) M560011M M560012M cmssa13 cmssa14; array cyssa2 (t) M8100111Y M8100112Y cyssa23 cyssa24; array cmssa2 (t) M8100111M M8100112M cmssa23 cmssa24; array cyssa3 (t) cyssa31 m10000112y cyssa33 cyssa34; array cmssa3 (t) cmssa31 m10000112m cmssa33 cmssa34; array cyssc2 (t) cyssc21 cyssc22 M560033Y cyssc24; array cmssc2 (t) cmssc21 cmssc22 M560033M cmssc24; array cyss1 (p) cyssa1 cyssb1 cyssc1 cyssd1; array cmss1 (p) cmssa1 cmssb1 cmssc1 cmssd1; array cyss2 (p) cyssa2 cyssb2 cyssc2 cyssd2; array cmss2 (p) cmssa2 cmssb2 cmssc2 cmssd2; array cyss3 (p) cyssa3 cyssb3 cyssc3 cyssd3; array cmss3 (p) hmsas3 cmssb3 cmssc3 cmssd3; array cyss4 (p) cyssa4 cyssb4 cyssc4 cyssd4; array cmss4 (p) cmssa4 cmssb4 cmssc4 cmssd4; array cyss (i) cyss1 cyss2 cyss3 cyss4; array cmss (i) cmss1 cmss2 cmss3 cmss4; array fstata (t) M114001 fstata2 fstata3 fstata4; array fstatb (t) M114002 fstatb2 fstatb3 fstatb4; array fstatc (t) M114003 fstatc2 fstatc3 fstatc4; array fstatd (t) M114004 fstatd2 fstatd3 fstatd4; array fstat (p) fstata fstatb fstatc fstatd; array eyra (t) M117001Y eyra2 eyra3 eyra4; array eyrb (t) M117002Y eyrb2 eyrb3 eyrb4; array eyrc (t) M117003Y eyrc2 eyrc3 eyrc4; array emoa (t) M117001M emoa2 emoa3 emoa4; array emob (t) M117002M emob2 emob3 emob4; array emoc (t) M117003M emoc2 emoc3 emoc4; array eyr (p) eyra eyrb eyrc eyrd; array emo (p) emoa emob emoc emod; array divoa (t) M4350A1 divoa2 divoa3 divoa4; array divob (t) M4350A2 divob2 divob3 divob4; array divoc (t) M4350A3 divoc2 divoc3 divoc4; array divo (p) divoa divob divoc divod; array divoya (t) M4360A1y divoya2 divoya3 divoya4; array divoyb (t) M4360A2y divoyb2 divoyb3 divoyb4; array divoyc (t) divoyc1 divoyc2 divoyc3 divoyc4; array divoyd (t) divoyd1 divoyd2 divoyd3 divoyd4; array divoma (t) M4360A1M divoma2 divoma3 divoma4; array divomb (t) M4360A2M divomb2 divomb3 divomb4; array divomc (t) divomc1 divomc2 divomc3 divomc4; array divomd (t) divomd1 divomd2 divomd3 divomd4; array divoy (p) divoya divoyb divoyc divoyd; array divom (p) divoma divomb divomc divomd; array annua (t) M4350A1 annua2 annua3 annua4; array annub (t) M4350A2 annub2 annub3 annub4; array annuc (t) annuc1 annuc2 annuc3 annuc4; array annud (t) annud1 annud2 annud3 annud4; array annu (p) annua annub annuc annud; array annuya (t) m4370A1y annuya2 annuya3 annuya4; array annuyb (t) annuyb1 annuyb2 annuyb3 annuyb4; array annuyc (t) annuyc1 annuyc2 annuyc3 annuyc4; array annuyd (t) annuyd1 annuyd2 annuyd3 annuyd4; array annuma (t) m4370A1M annuma2 annuma3 annuma4; array annumb (t) annumb1 annumb2 annumb3 annumb4; array annumc (t) annumc1 annumc2 annumc3 annumc4; array annumd (t) annumd1 annumd2 annumd3 annumd4; array annuy (p) annuya annuyb annuyc annuyd; array annum (p) annuma annumb annumc annumd; array cmarsa (t) cmarsa1 cmarsa2 cmarsa3 cmarsa4; array cmarsb (t) cmarsb1 cmarsb2 cmarsb3 cmarsb4; array cmarsc (t) cmarsc1 cmarsc2 cmarsc3 cmarsc4; array cmarsd (t) cmarsd1 cmarsd2 cmarsd3 cmarsd4; array cmars (p) cmarsa cmarsb cmarsc cmarsd; array imarsa (t) imarsa1 imarsa2 imarsa3 imarsa4; array imarsb (t) imarsb1 imarsb2 imarsb3 imarsb4; array imarsc (t) imarsc1 imarsc2 imarsc3 imarsc4; array imarsd (t) imarsd1 imarsd2 imarsd3 imarsd4; array imars (p) imarsa imarsb imarsc imarsd; array coha (t) cohaf cohas cohat cohafour; array cohb (t) cohbf cohbs cohbt cohbfour; array cohc (t) cohcf cohcs cohct cohcfour; array cohd (t) cohdf cohds cohdt cohdfour; array coh (p) coha cohb cohc cohd ; array ma (t) maf mas mat mafour; array mb (t) mbf mbs mbt mbfour; array mc (t) mcf mcs mct mcfour; array md (t) mcf mcs mdt mdfour; array m (p) ma mb mc md; do l1=1 to dlicm; if l1=dlicm then do; cohhR18=R18co; marsR18=R18mars; linkR18=R18link; end; end; %let n=12; array slink(&n); array sco(&n); do i=1 to (&n); slink{i}= -4; sco{i}=-4; end; drop i; do i=1 to (&n); do l1=169 to dlicm until (R18mars=2 & R18co=200+i); if R18mars=2 & R18co=200+i then do ; slink{i}=R18link; sco{i}=R18co; end; end; end; do i=1 to (&n); flaglink=0; if slink{i}>0 then flaglink=1; if flaglink=0 then do l1=169 to dlicm; if R18mars=2 & R18co in (. -4) & R18link>-4 then do ; if R18link ne slink1 & R18link ne slink2 & R18link ne slink3 & R18link ne slink4 & R18link ne slink5 & R18link ne slink6 & R18link ne slink7 & R18link ne slink8 & R18link ne slink9 & R18link ne slink10 & R18link ne slink11 & R18link ne slink12 then do; slink{i}=R18link; sco{i}=R18co; end; if cohhR18 in (. -4) & flagp=1 & 21 then cohhR18=premar_ttl+200; end; end; end; end; do i=1 to 11; if slink(i)>0 then do j=i+1 to 12; if slink(i)=slink(j) then slink(j)=-4; end; end; flags=0; do t=1 to 4; do i=1 to (&n); if slink{i} =par{t} & slink{i}>0 & flagp=1 then flags=1; end; end; if (cohhR18<=-4 or cohhR18=. or 1000 then do ; if -33 & flags=0 then cohhR18=100+precoh_ttl; end; if cohhR18<=-4 or cohhR18=.& flags=1 & flagp=1 then do; if (dlimarstat>3 | dlimarstat=-3 ) then cohhR18=200+premar_ttl; end; do p=1 to 4; do t=1 to 4; if 100-4 & cemm in (-4, .) & m620 ~=3 then do; cohhR18=200+premar_ttl; end; end; end; do p=1 to 4; do t=1 to 4; if cdli=1 then do; csmc=dlicm; end; if -40 then do; cmars=1; csmm=dlicm; end; if ysc>0 & msc>0 then csmc=(ysc-1980)*12+msc; if (-3 le ysc le 0 | -3 le msc le 0) then csmc=dlicm; if ysc > 0 & -3 <= msc <= 0 then csmc=(ysc-1980)*12+1; if yec>0 & mec>0 then cemc=(yec-1980)*12+mec-1; if csmc>=cemc+1 & cemc>0 then do; x0f1=1; cemc=csmc; end; if (-3 <=yec <= 0 | -3 <= mec <= 0) & (cbeg ^= 1 & cdli^= 1 & cbeg2 ^= 1) then do; x0f2=1; cemc=doicm-1; end; if (-3 <=yec <= 0 | -3 <= mec <= 0) & (fstat=3 | fstat=4) & eyr >0 & emo >0 then cemc=(eyr-1980)*12+emo-1; if yec>0 & -4=doicm then cemc=doicm-1; end; if (cdli=1 | cbeg=1) | (cbeg2=1) then cemc=doicm; do j=1 to (&n); if cdli=1 & m620=4 & slink(j)=par(p) then cemc=doicm-1; if dlimarstat in(3 4) & csmc>-4 & slink(j)=p1 & p2=-4 then do; if mbeg=0 then mbeg=1; end; end; if mbeg=1 then do; cmars=1; imars=1; csmm=csmc; ymsc=.; do j=1 to (&n); if (dlimarstat=3 | dlimarstat=4) & (cbeg=1 | yec>-4) & slink(j)=par(p) then do; cmars=1; imars=1; csmm=start; cemm=doicm; ymsc=.; end; end; if(dlimarstat=7 | dlimarstat=8) & premar_ttl>=1 then do; cmars=3; imars=3; csms=start; cems=csmma1; ymsc=.; end; end; if -40 & cemc>0 & csmc LE C LE cemc then cohal=1; if csmc>0 & cemc>0 & csmc LE C LE cemc then coh=1; if cmars=1 & csmm>0 & cemm>0 & csmm le c le cemm then mal=1; if imars=1 & csmm>0 & cemm>0 & csmm le c le cemm then m=1; if 2 le cmars le 3 & csms>0 & cems>0 & csms le c le cems then mal=cmars; if 2 le imars le 3 & csms>0 & cems>0 & csms le c le cems then m=imars; if cmars =4 & csms>0 & cems>0 & csms le c le cems then mal=cmars; if imars =4 & csms>0 & cems>0 & csms le c le cems then m=imars; if cmars=0 then mal=cmars; if -2 le cmars le -1 then mal=cmars; if -2 le imars le -1 then m=imars; end; if ysp>0 & msp>0 then do; csmm = (ysp-1980)*12+msp; cemm=doicm; cmars=1; imars=1; marcm=csmm; flagsp=1; end; if -40 & -40 & divom>0 then do; edit4360a(i)=(divoy -1980)*12 +divom; end; end; if divo=1 then do; if divoy>0 & divom>0 then do; cmars=3; imars=3; csms=(divoy-1980)*12 +divom; cems=doicm; end; if (-40 & -4cems then csms=cems; end; end; if (divo=2 ) & annuy>0 & annum>0 then do; cemc=(annuy-1980)*12 +annum; cmars=0; imars=0; cemm=cemc; csmm=start; end; do i=1 to 4; if cstat=1 & cym>0 & cmm>0 then do; cmars=1; imars=1; csmm=(cym-1980)*12+cmm; end; if -2 le cstat le -1 & (cym>-4 & cmm>-4 )then do; cmars=cstat; imars=cstat; end; if cstat=1 & (-40 & cmss>0 then do; cmars=2; imars=2; csms=(cyss-1980)*12+cmss; end; if cstat=3 & (-40 & cmss>0 then do; cmars=3; imars=3; csms=(cyss-1980)*12+cmss; end; if cstat=4 & (-40 & cmss>0 then do; cmars=0; imars=0; cemc=(cyss-1980)*12+cmss; end; if csmc>0 & cemc in(-4 .) then cemc=doicm; if csmm>0 & cemm in(-4 .) then cemm=doicm; end; if fstat=2 & 00 & emo>0 then do; cmars=2; imars=2; csms=(eyr-1980)*12 + emo; cems=doicm; end; if fstat=3 & (-40 & -4doicm then csms=doicm; cems=doicm; end; if fstat=4 & eyr>0 & emo>0 then do; cmars=3; imars=3; csms=(eyr-1980)*12 + emo; cems=doicm; end; if fstat=4 & (-40 & -4doicm then csms=doicm ; cems=doicm; end; if fstat=5 & eyr>0 & emo>0 & csmm>0 then do; cmars=4; imars=4; *10-31-03 widowed; csms=(eyr-1980)*12 + emo; widow=1; end; if fstat=5 & (-40 then do; cmars=4; imars=4; csms=doicm; widow=1; end; if csms>0 then cems=doicm; if -4=1 then cems=csmm-1; if mbeg=1 & (dlimarstat=7 | dlimarstat=8) & premar_ttl>=1 & p1>1900 then cemm=doicm; C=0; do l=1 to doicm; C=C+1; if csmc>0 & cemc>0 & csmc LE C LE cemc then cohal=1; if csmc>0 & cemc>0 & csmc LE C LE cemc then coh=1; if cmars>=1 & csmm>0 & cemm>0 & csmm le c le cemm then mal=1; if imars>=1 & csmm>0 & cemm>0 & csmm le c le cemm then m=1; if cmars = 2 & csms>0 & cems>0 & csms le c le cems then mal=2; if imars = 2 & csms>0 & cems>0 & csms le c le cems then m=2; if cmars = 3 & csms>0 & cems>0 & csms le c le cems then mal=3; if imars = 3 & csms>0 & cems>0 & csms le c le cems then m=3; if cmars = 4 & csms>0 & cems>0 & csms le c le cems then mal=4; if imars = 4 & csms>0 & cems>0 & csms le c le cems then m=4; if -2 le cmars le -1 then mal=cmars; if -2 le imars le -1 then m=imars; if cmars=0 then mal=cmars; if c=doicm then do; if mal in (-1,-2) then mars=-3; if mal=0 & cohal=1 then mars=1; if mal=0 & cohal=0 then mars=2; if mal=1 & cohal=1 then mars=3; if mal=1 & cohal=0 then mars=4; if mal=2 & cohal=1 then mars=5; if mal=2 & cohal=0 then mars=6; if mal=3 & cohal=1 then mars=7; if mal=3 & cohal=0 then mars=8; if mal=4 & cohal=1 then mars=9; if mal=4 & cohal=0 then mars=10; marstat=mars; end; end; end; end; if 1 le marstat le 2 then cmarstat=0; if 3 le marstat le 4 then cmarstat=1; if 5 le marstat le 6 then cmarstat=2; if 7 le marstat le 8 then cmarstat=3; if 9 le marstat le 10 then cmarstat=4; if -2 le marstat le -1 then cmarstat=marstat; if marstat=-3 then cmarstat=-3; do p=1 to 4; do t=1 to 4; do i=1 to 4; if yco>0 & mco>0 then do; ymco=((yco-1980)*12)+mco; cohcm= min (ymcoa1, ymcoa2, ymcoa3, ymcoa4, ymcob1, ymcob2, ymcob3, ymcob4, ymcoc1, ymcoc2, ymcoc3, ymcoc4, ymcod1, ymcod2, ymcod3, ymcod4); end; if cym >0 & cmm>0 & csmm>0 then do; marcm= min(csmma1, csmma2, csmma3, csmma4, csmmb1, csmmb2, csmmb3, csmmb4, csmmc1, csmmc2, csmmc3, csmmc4, csmmd1, csmmd2, csmmd3, csmmd4); end; if mbeg=1 then do; if dlicohcm<=-4 & dlimarcm<=-4 then do; marcm=cohcm; cohcm=-4; end; else if dlicohcm>=-3 and dlimarcm<=-4 then do; marcm=cohcm; end; end; end; end; end; array cvcm cohcm marcm; array cvy cohy mary; array cvm cohm marm; array fir firstc492-firstc157; do over cvcm; do i=492 to 157 by -12; fir=i; if fir-11<=cvcm<=fir then do; cvm=cvcm-(fir-12); cvy=(cvcm-cvm)/12+1980; end; end; if cvcm in (. -4) then do; cvcm=-4; cvy=-4; cvm=-4; end; end; if cohcm<0 then do; do p=1 to 4; do t=1 to 4; if -40 then do; cohm=msc; cohy=-3; cohcm=-3; end; if ysc>0 & -40 & marcm<0 then do; marm=cmm; mary=-3; marcm=-3; end; if cym>0 & -40 & -30 & dlicohcm=-4 & dlimarm=-4 then do; cohm=-4; cohcm=-4; cohy=-4; marcm=-3; mary=-3; marm=msc; end; end; end; end; do p=1 to 4; do t=1 to 4; if mbeg=1 then do; if -3-4 then do; marcm=-3; mary=-3; marm=-3; end; if ysc>0 & -3-4 then do; marcm=-3; mary=ysc; marm=-3; end; if -30 & dlimarm=-4 & dlicohcm>-4 then do; marcm=-3; mary=-3; marm=msc; end; end; end; end; do p=1 to 4; do t=1 to 4; if M45001=0 & M45002=1 & dlicohcm = -4 then do; if m31001m>0 & m31001y>0 then do; cohm= m31001m; cohy= m31001y; cohcm= m31001cm; end; end; end; end; if dlicohm >-4 then do; cohcm= dlicohcm; cohy= dlicohy; cohm= dlicohm; end; if dlimarcm>-4 then do; marcm=dlimarcm; mary=dlimary; marm=dlimarm; end; do l=start to doicm; if -2 le mal le -1 then mars=-3; if mal=0 & cohal=1 then mars=1; if mal=0 & cohal=0 then mars=0; if mal=1 & cohal=1 then mars=2; if mal=1 & cohal=0 then mars=2; if mal=2 & cohal=1 then mars=3; if mal=2 & cohal=0 then mars=3; if mal=3 & cohal=1 then mars=4; if mal=3 & cohal=0 then mars=4; if mal=4 & cohal=1 then mars=5; if mal=4 & cohal=0 then mars=5; end; do l=start to doicm; if cohaf=1 | cohas=1 | cohat=1 | cohafour=1 then cohaa=1; if cohbf=1 | cohbs=1 | cohbt=1 | cohbfour=1 then cohbb=1; if cohcf=1 | cohcs=1 | cohct=1 | cohcfour=1 then cohcc=1; if cohdf=1 | cohds=1 | cohdt=1 | cohdfour=1 then cohdd=1; if maf=1 | mas=1 | mat=1 | mafour=1 then maa=1; if mbf=1 | mbs=1 | mbt=1 | mbfour=1 then mbb=1; if mcf=1 | mcs=1 | mct=1 | mcfour=1 then mcc=1; if mdf=1 | mds=1 | mdt=1 | mdfour=1 then mdd=1; end; do l=1 to doicm; if dual <0 then do; if cohdd=1 | mdd=1 then link=p4; if cohcc=1 | mcc=1 then link=p3; if cohbb=1 | mbb=1 then link=p2; if cohaa=1 | maa=1 then link=p1; end; else if dual>0 then do; if cohaa=1 | maa=1 then link=p1; if cohbb=1 | mbb=1 then link=p2; if cohcc=1 | mcc=1 then link=p3; if cohdd=1 | mdd=1 then link=p4; end; end; do l=start to doicm; if cohdd=1 then do; cohh=104; link=p4; end; if cohcc=1 then do; cohh=103; link=p3; end; if cohbb=1 then do; cohh=102; link=p2; end; if cohaa=1 then do; cohh=101; link=p1; end; if mdd=1 & cohdd=1 then do; cohh=204; link=p4; end; if mcc=1 & cohcc=1 then do; cohh=203; link=p3; end; if mbb=1 & cohbb=1 then do; cohh=202; link=p2; end; if maa=1 & cohaa=1 then do; cohh=201; link=p1; end; end; do l=start to doicm; if cohh>0 & cohh^=202 & cohaa=1 & cohbb=1 then dual=102; if dual=102 & mbb=1 then dual=202; if cohh>0 & cohh ^= 203 & cohh ^= 103 & cohcc=1 then dual=103; if dual=103 & mcc=1 then dual=203; if cohh>0 & cohh ^= 204 & cohh ^= 104 & cohdd=1 then dual=104; if dual=104 & mdd=1 then dual=204; if dual>0 & cohh>0 & cohh ^=dual then do; cohh=dual; if cohh in (102 202) then link=p2; else if cohh in (103 203) then link=p3; else if cohh in (104 204) then link=p4; end; if dual>0 then chkdual=1; end; countP=0; if p1>1900 & p2=-4 then countP=1; else if 0 1900 & 01900 & p2>1900 & p3=-4 then countP=4; else if (01900) & p3=-4 then countP=5; else if 01900 & p2>1900 & 01900 & p3>1900 | p1>1900 & 01900) then countP=7; else if ( p1>1900 & 01900 & 01900 then countP=9; else if 01900 & p2>1900 & p3>1900 then countP=11; nump18=0; nums18=0; nump18=precoh_ttl; nums18=premar_ttl; if (nump18=. | nump18<0) then nump18=0; if (nums18=. | nums18<0) then nums18=0; numR18=max(nump18,nums18); flagpcount=0; if countP=5 & cemca1200 then cohh=cohh-1; end; end; if min(of r19sp209-r19sp487)=3 then do; x48f1=22; do l=start to doicm; if cohh>200 then cohh=cohh-2; end; end; if min(of r19sp209-r19sp487)=4 then do; x48f1=23; do l=start to doicm; if cohh>200 then cohh=cohh-3; end; end; do l=start to doicm; if 100 lt cohh lt 200 then r19pt=cohh-100; if cohh gt 200 then r19sp=cohh-200; end; if max(of r19pt209-r19pt487)=4 then do; x48f2_14=1; do l=start to doicm; if cohh=103 then x48f2_143=1; if cohh=102 then x48f2_142=1; end; if x48f2_143 ne 1 & x48f2_142=1 then do l=start to doicm; if cohh=104 then cohh=103; end; else if x48f2_143=1 & x48f2_142 ne 1 then do l=start to doicm; if cohh in (103 104) then cohh=cohh-1; end; else if x48f2_143 ne 1 & x48f2_142 ne 1 then do l=start to doicm; if cohh=104 then cohh=102; end; end; if max(of r19pt209-r19pt487)=3 then do; x48f2_13=1; do l=start to doicm; if cohh=102 then x48f2_132=1; end; if x48f2_132 ne 1 then do l=start to doicm; if cohh=103 then cohh=102; end; end; if max(of r19sp209-r19sp487)=4 then do; x48f2_24=1; do l=start to doicm; if cohh=203 then x48f2_243=1; if cohh=202 then x48f2_242=1; end; if x48f2_243 ne 1 & x48f2_242=1 then do l=start to doicm; if cohh=204 then cohh=203; end; else if x48f2_243=1 & x48f2_242 ne 1 then do l=start to doicm; if cohh in (203 204) then cohh=cohh-1; end; else if x48f2_243 ne 1 & x48f2_242 ne 1 then do l=start to doicm; if cohh=204 then cohh=202; end; end; if max(of r19sp209-r19sp487)=3 then do; x48f2_23=1; do l=start to doicm; if cohh=202 then x48f2_232=1; end; if x48f2_232 ne 1 then do l=start to doicm; if cohh=203 then cohh=202; end; end; c1minus=0; c2minus=0; if flagp=1 then do l=start to doicm; if cohh=101 & 0200 then do; c2minus=c2minus+1; c201=R18co; l1=487; end; end; if cohh=202 & 0200 then do; c2minus=c2minus+1; c202=R18co; l1=487; end; end; if cohh=203 & 0200 then do; c2minus=c2minus+1; c203=R18co; l1=487; end; end; if cohh=204 & 0200 then do; c2minus=c2minus+1; c204=R18co; l1=487; end; end; if cohh=201 & 0200 then do; cohh=c201; cohh_p=cohh_p+1; end; else if cohh=202 & c202>200 then do; cohh=c202; cohh_p=cohh_p+1; end; else if cohh=203 & c203>200 then do; cohh=c203; cohh_p=cohh_p+1; end; else if cohh=204 & c204>200 then do; cohh=c204; cohh_p=cohh_p+1; end; end; do l=start to doicm; if cohh>0 & cohh_p=0 then do; if 100200 then cohh=cohh+nums18-c2minus; end; end; do l =start to doicm; if (dual=101 | dual=102 | dual=103 | dual=104 ) then r19dual=1; else if cohh>0 & R19dual~=1 then R19dual=0; else if R19dual ^= 1 & R19dual ^= 0 then R19dual=-4; end; if R19dual=. then R19dual=-4; do l=start to 487; if 100 lt cohh lt 200 then r19pt=cohh-100; if cohh gt 200 then r19sp=cohh-200; end; nump19=max(of R19pt209-R19pt487); nums19=max(of R19sp209-R19sp487); numR19=max(nump19,nums19); coh_ttl=nump19; mar_ttl=nums19; if nump19=. then coh_ttl=precoh_ttl; if nums19=. then mar_ttl=premar_ttl; if p1=-4 & p2=-4 & p3=-4 then do; coh_ttl=precoh_ttl; mar_ttl=premar_ttl; end; do l=start to doicm; do i=2 to 3; if cohh>=200+i & premar_ttl=0 & mbb=1 then do; cohh=cohh-1; mar_ttl=i; flagmbb=l; end; if cohh>=200+i & premar_ttl=0 & mcc=1 then do; cohh=cohh-1; mar_ttl=i; flagmcc=l; end; end; end; do l=start to doicm; do i=1 to 3; if cohh>=200+i & premar_ttl= i-1 & flagmcc=. then do; cohh=200+i; mar_ttl=i; flagmbb1=l; end; end; end; flagcohh101=0; do l=start to doicm; if cohh=101 then flagcohh101=1; end; if flagcohh101=0 & precoh_ttl<1 & countP in (4 11) then x2f1=1; do p=1 to 4; do t=1 to 4; / if mbeg=1 & Mar_ttl>0 & precoh_ttl<=0 & dlicohcm<-3 & p2=-4 then do; coh_ttl=0; cohcm=-4; cohy=-4; cohm=-4; end; if ysp>-4 & msp>-4 & premar_ttl=0 then mar_ttl=1; if ysp>-4 & msp>-4 & premar_ttl>0 & marstat=4 & dlimarstat>5 & (p1>1900 | p2>1900) then do; x2f2=1; mar_ttl=premar_ttl+1; end; end; end; if mar_ttl>0 & coh_ttl=0 & cohcm>-4 then do; cohy=-4; cohm=-4; cohcm=-4; end; if precoh_ttl>coh_ttl & coh_ttl>0 then coh_ttl=precoh_ttl; do p=1 to 4; do t=1 to 4; if M4350a1=3 & mar_ttl>0 & m620=1 & (dlimarstat=5 | dlimarstat=6) & m114001 not in (3,4) then do; x2f3=1; do l=start to doicm; mars=2; if mar_ttl=1 & cohh>0 & p2=-4 then cohh=200+1; if link=-4 | link=. then link=slink1; marstat=4; cmarstat=1; if M700011Y>-4 & marstat=4 then marstat=3; end; end; end; end; if x2f3=1 then do; do l=start to 487; if 100 lt cohh lt 200 then r19pt=cohh-100; if cohh gt 200 then r19sp=cohh-200; end; nump19=max(of R19pt209-R19pt487); nums19=max(of R19sp209-R19sp487); if nump19>0 and nump19 ne coh_ttl then do; coh_ttl=nump19; x2f3p=1;end; if nums19>0 and nums19 ne mar_ttl then do; mar_ttl=nums19; x2f3s=1;end; if nump19=. then coh_ttl=precoh_ttl; if nums19=. then mar_ttl=premar_ttl; end; do p=1 to 4; do t=1 to 4; if dlimarstat=-3 and premar_ttl>0 and m620 in (-1, -2) and fstat=-4 then do; x2f4=1; marstat=dlimarstat; cmarstat=dlimarstat; do l=start to doicm; mars=dlimarstat; end; end; end; end; matchlink=0; flagslink=0; do i=1 to (&n); if slink(i)<=-4 then do l=start to doicm; if mars=2 & (cohh>200 or cohh in (. -4)) & link>-4 & (link ne slink1 & link ne slink2 & link ne slink3 & link ne slink4 & link ne slink5 & link ne slink6 & link ne slink7 & link ne slink8 & link ne slink9 & link ne slink10 & link ne slink11 & link ne slink12) then do; slink{i}=link; sco{i}=cohh; end; end; end; do i=1 to 11; if slink(i)>0 then do j=i+1 to 12; if slink(i)=slink(j) then slink(j)=-4; end; end; do l=doicm; if marstat=3 & cohh>200 then r19link=link; end; do i=1 to (&n); if slink{i}=r19link & r19link>0 then matchlink=1; end; if matchlink ne 1 & marstat=3 then flagslink=1; if flagslink=1 then marstat=4 ; if flagslink=1 & (dlim620=0 & m620=0 ) then do; marstat=-3; cmarstat=-3; do l=start to doicm; mars=-3; end; end; if flagslink=1 & dlim620=2 & m620=2 & m4350a1 ne 3 & m114001 ne 1 then do; marstat=5; cmarstat=2; do l=start to doicm; mars=3; if link=slink1 | link=slink2 | link=slink3 | link=slink4 | link=slink5 | link=slink6 | link=slink7 | link=slink8 | link=slink9 | link=slink10 | link=slink11 | link=slink12 then link=-4; end; end; spouse=0; cspouse=0; do t=1 to 4; do i=1 to (&n); if slink{i}>0 & par{t}>0 then do; if par{t}=slink{i} & marstat>2 then spouse=1 ; end; end; end; if mar_ttl>0 & 20 then cspouse=1; if cspouse=1 and (dlim620=3 & m620=3 ) & p1=-4 then do; do l=start to doicm; mars=4; marstat=8; cmarstat=3; end; end; if cspouse=1 and ((dlim620=3 )& m620=3 & p1>0 & marstat<=6) then do; do l=start to doicm; mars=4; marstat=7; cmarstat=3; end; end; if cspouse=1 and (dlim620=2 & m620=3 ) & p1=-4 then do; do l=start to doicm; mars=3; marstat=6; cmarstat=2; end; end; if cspouse=1 and (dlim620=2 & m620=2) & (p1=-4 | (m46701 ne 1 & m47001 ne 1 & m114001 ne 1)) then do; do l=start to doicm; mars=3; marstat=6; cmarstat=2; end; end; if cspouse=1 and (dlim620=2 & m620=3 & p1>0 & marstat<=6) then do; x3f1=1; do l=start to doicm; mars=3; marstat=5; cmarstat=2; end; end; if( cspouse=1 and ((dlim620<2 & m620=3) | m620=0 ) & dlimarstat ne 5 and dlimarstat ne 6 | cspouse=1 & m620 in (-1,-2) | cspouse=1 & dlimarstat=-3 & m620 =2 & p1=-4) then do; x3f2=1; do l=start to doicm; mars=-3; marstat=-3; cmarstat=-3; if link=<-3 then link=-3; end; end; slink_end=-4; do i=1 to (&n); if slink(i)>0 then slink_end=slink(i); end; if (3<=marstat<=4) & (3<=dlimarstat<=4) then do; do l=start to doicm; if link=. then do; x3f3=l; end; end; end; if (((M540011=3 | m7900111=3) & M114001=4)|((M540021=3 | m7900211=3) & M114002=4)) & csmsa1-1>=start then do; x3f4=1; end; if (m620=0 ) & mar_ttl=0 & premar_ttl=0 then do; x3f5=1; if marstat=6 then do; marstat=2; cmarstat=0; do l=start to doicm; if mars=3 then mars=0; end; end; end; count620=0; do i=11 to 18; if m620r(i)=1 then do; count620=count620+1; end; end; if count620>=2 & m620>=1 & mar_ttl<=0 then flag620=1; if m620=1 & mar_ttl=0 & count620>=1 then case620=1; if dlimarstat=4 & flagp=1 & marstat<5 then do; do l=start to doicm; if mars=2 & 100=-4 | cohh>-4) then chk2=1; if mars in(1,3,5,7) & link in(. 4) then chk3=1; end; end; if mar_ttl>premar_ttl+1 then chk4=1; if dlimarstat in (7,8) & spouse=1 & mar_ttl=premar_ttl & marstat<7 &( m7900111=1 | m45001=1 | m7900112=1) then do; x3f6=1; do l=start to doicm; if cohh>200 then cohh=cohh+1; end; mar_ttl=mar_ttl+1; end; max_slink=0; do i=1 to (&n); if slink(i) ne -4 then max_slink=i; end; if max_slink>mar_ttl>=0 then do; mtl_slink=1; mar_ttl=max_slink; end; checkfir=0; checkfir_M=0; checkmar=0; checkcohh=0; type1=0; type2=0; if dlicohcm<-3 & cohcm>-4 then do; do l=start to doicm until (1000 & checkcohh>0 & cohcm ne checkcohh then checkfir=1; if checkfir=1 then do; if checkcohh = cohcm+1 then type1=1; if checkcohh > cohcm+1 then type1=2; if checkcohh < cohcm then type1=3; end; end; n=15; if dlimarcm<-3 & marcm>-4 then do; do l=start to doicm until (mars=2); if mars=2 then checkmar=l; end; if marcm>0 & checkmar>0 & marcm ne checkmar then checkfir_M=1; if checkfir_M=1 then do; if checkmar = marcm+1 then type2=1; if checkmar > marcm+1 then type2=2; if checkmar < marcm then type2=3; end; end; diffcohh=-4; diffmar=-4; if diffcohhR18 =1 | checkfir=1 then diffcohh=1; if diffMARR18 =1 | checkfir_M=1 then diffmar=1; if m716jb1M >-4 then m716jb1=4; if m716jb2M >-4 then m716jb2=4; if m716jb3M >-4 then m716jb3=4; if m716wb1m>-4 then m716w1 =4; if m716wb2m>-4 then m716w2 =4; if m716wb3m>-4 then m716w3 =4; if m717j1=3 then m717j1 =4; if m4350a1=1 then m4350a1=4; if m4350a2=1 then m4350a2=4; if m13250m>-4 then m13225 =4; divo_11=m716jb1; divo_12=m716jb2; divo_13=m716jb3; divo_14=m716w1; divo_15=m716w3; divo_16=m717j1; divoM_11=m716jb1m; divoM_12=m716jb2m; divoM_13=m716jb3m; divoM_14=m716wb1m; divoM_15=m716wb3m; divoM_16=m717jb1m; divoY_11=m716jb1y; divoY_12=m716jb2y; divoY_13=m716jb3y; divoY_14=m716wb1y; divoY_15=m716wb3y; divoY_16=m717jb1y; divo_21=m540011; divo_22=m540012; divo_23=m540021; divo_24=m540031; divo_25=m540032; divo_26=m540033; divoM_21=m560011m; divoM_22=m560012m; divoM_23=m560021m; divoM_24=m560031m; divoM_25=m560032m; divoM_26=m560033m; divoY_21=m560011y; divoY_22=m560012y; divoY_23=m560021y; divoY_24=m560031y; divoY_25=m560032y; divoY_26=m560033y; divo_31=m7900111; divo_32=m7900112; divo_33=m7900113; divo_35 =m7900211; divo_36 =m7900311; divoM_31=m8100111m; divoM_32=m8100112m; divoM_33=m8100113m; divoM_35=m8100211m; divoM_36=m8100311m; divoY_31=m8100111y; divoY_32=m8100112y; divoY_33=m8100113y; divoY_35=m8100211y; divoY_36=m8100311y; divo_41 =m9800111; divo_42=M9800112; divo_43=m114001; divo_44=m114002; divo_45=m114003; divo_46=m13225; divoM_41=m10000111m; divoM_42=M10000112M; divoM_43=m117001m; divoM_44=m117002m; divoM_45=m117003m; divoM_46=m13250m; divoY_41=M10000111Y; divoY_42=M10000112Y; divoY_43=m117001y; divoY_44=m117002y; divoY_45=m117003y; divoY_46=m13250y; divo_51=m4350a1; divo_52=m4350a2; divoM_51=m4360a1m; divoM_52=m4360a2m; divoY_51=m4360a1y; divoY_52=m4360a2y; array divo_ (5,6) divo_11 - divo_16 divo_21- divo_26 divo_31- divo_36 divo_41- divo_46 divo_51- divo_56; array divo_M (5,6) divoM_11 - divoM_16 divoM_21 - divoM_26 divoM_31 - divoM_36 divoM_41 - divoM_46 divoM_51 - divoM_56; array divo_Y (5,6) divoY_11- divoY_16 divoY_21 - divoY_26 divoY_31 - divoY_36 divoY_41 - divoY_46 divoY_51 - divoY_56; array divo_cm (5,6) divocm_11 - divoCM_16 divocm_21 - divocm_26 divocm_31 - divocm_36 divocm_41 - divocm_46 divocm_51 - divocm_56; divo_cmR19=-4; divo_MR19 =-4; divo_YR19 =-4; end_marR19 =-4; divorceM=-4; divorceY=-4; divor_CM=-4; do k=1 to 5; do i=1 to 6; if divo_(k,i)=4 & divo_M(k,i)>0 & divo_Y(k,i)>0 then do; divo_cm(k,i) = (divo_Y(k,i)-1980)*12+ divo_M(k,i); divo_cmR19=min(of divocm_11 - divoCM_16 divocm_21 - divocm_26 divocm_31 - divocm_36 divocm_41 - divocm_46 divocm_51 - divocm_56); end; end; end; if divo_cmR19>0 then do; do i=492 to 204 by -12; fir=i; if fir-11<=divo_CMR19 <=fir then do; divo_MR19=divo_cmR19-(fir-12); divo_YR19=(divo_cmR19-divo_MR19)/12+1980; end; end; end; if divo_cmR19=-4 then do; do k=1 to 5; do i=1 to 6; if divo_(k,i)=4 & (divo_M(k,i) in ( -1, -2,-3 ) | divo_y(k,i) in (-1,-2,-3)) then do; divo_cm(k,i)=-3; divo_mR19=divo_m(k,i); divo_YR19=divo_Y(k,i); if divo_mR19 in (-1 -2,-3) then divo_mR19=-3; if divo_YR19 in (-1 -2,-3) then divo_YR19=-3; divo_cmR19=-3; end; end; end; end; count_w=0; do i=1 to 18; if m620r(i)=4 then count_w=count_w+1; end; if cspouse=1 and (dlim620=3 & m620=3 ) & marstat>6 & divo_CMR18 in (-4 ,.)then do; x3f7=1; divo_cmR19=-3; divo_mR19=-3; divo_YR19=-3; end_marR19=1; end_mar=1; end; if divo_cmR19>-4 & divo_CMR18<=-4 then do; divor_CM= divo_cmR19; divorceM= divo_mR19; divorceY= divo_YR19; end_marR19=1; end_mar=1; end; if (premar_ttl=2 | mar_ttl=2) & divo_CMR18=-4 & divo_cmR19=-4 & count_w=0 then do; divor_CM=-3; divorceM=-3; divorceY=-3; end_mar=3; x4f1=1; end; if divo_cmR18 >-4 & ( m4350a1 =2 | m4350a2 =2) then end_mar=3; IF (premar_ttl=2 | mar_ttl=2) & divo_cmR18>-4 & divo_CMR19>-4 & divo_CMR19 ~= divo_cmR18 then check_divo=1; if check_divo=1 then do; divor_CM= divo_CMR18; divorceM=divo_MR18; divorceY=divo_YR18; end_mar=end_marR18; end; if divo_cmR19=-4 & divo_CMR18>-4 then do; divor_CM= divo_CMR18; divorceM=divo_MR18; divorceY=divo_YR18; end_mar=end_marR18; end; if divo_CMR18>-4 then do; divor_CM= divo_CMR18; divorceM=divo_MR18; divorceY=divo_YR18; end_mar=end_marR18; end; if end_marR19=3 & end_mar=. then end_mar=end_marR19; if divo_CMR19>-4 & divo_cmR18>-4 & divo_CMR19 ~= divo_cmR18 then check_divo=2; if check_divo=2 then do; divor_CM= divo_CMR18; divorceM=divo_MR18; divorceY=divo_YR18; end_mar=end_marR18; if( 0-4 then do; divor_CM= divo_cmR19; divorceM= divo_mR19; divorceY= divo_YR19; end; end; if divor_CM<-3 then do; divor_CM=-4; divorceM=-4; divorceY=-4; end_mar=-4; end; if divor_CM>-4 & (mar_ttl=0 | marcm=-4) then flagM=1; do l=start to doicm; if mars=0 & (dlimarstat >2 | dlimarstat=-3) then flagmars1=l; if flagmars1=l & m4350a1 ^=2 then mars=2; if (mars=1 or mars=2) & link=-4 then flagmars2=l; if cohh=100 | cohh=200 then flagcohh=l; end; end; if int_y in(-5 -4 .) then do; cohm=dlicohm; cohy=dlicohy; cohcm=dlicohcm; marm=dlimarm; mary=dlimary; marcm=dlimarcm; mar_ttl=premar_ttl; coh_ttl=precoh_ttl; diffcohh=diffcohhR18; diffmar=diffmarR18; divorceM = divo_mR18 ; divorceY = divo_YR18 ; divor_CM = divo_cmR18; end_mar = end_marR18; end; if R19dual = . then R19dual = -4; if marstat =. then marstat = -5; if cmarstat=. then cmarstat = -5; do l= 208 to 487; if mars=. then mars=-4; if cohh=. then cohh=-4; if link=. then link=-4; end; if r19dual=-5 then r19dual=-4; if marstat=4 then do; do l=doicm; r19mars_a=mars; r19cohh_a=cohh; r19link_a=link; end; do l1=dlicm; if R18mars=2 then do; dli_cohh_a=R18cohh; dli_link_a=r18link; end; end; end; array cmarscvc (*) CMARSTAr1-cmarstar18 cmarstat; cvc_cmars=-4; do i=1 to dim(cmarscvc); if cmarscvc(i)>-4 then cvc_cmars=cmarscvc(i); end;