/******************************************************************************
Project: Immigration Data Linkage Project (ircc_link)

File: S:\ircc_link\prog\gillesd\concept_samples\outmigration.sas

Programmer: Gilles R. Detillieux

First Created: Sep. 8, 2022

Revisions:
  - Dec. 20, 2022: - adapted to updated "v2" datasets & split crosswalks
  - Dec. 22, 2022: - updated end dates to cover new data
  - Jan. 12, 2023: - adapted to most recent "v2" variable names
  - Mar. 5, 2025: - adapted to 2024 "V2" dataset & variable names, acqdt values
  - Jul. 22, 2025: - use get_ircc_pr macro to fetch data

Input datasets:
registry.ircc_crosswalk_1985jan_v2: IRCC linkage dataset, 1985-2023
registry.ircc_legacy_landing_1985jan_v2: IRCC landing data from 1985 to 1999
registry.ircc_landing_2000jan_v2: IRCC landing data from 2000 to 2023
registry.mhmrs_1970jun: MH insurance registry

Description: This program finds landing date of immigrants living in Manitoba,
for anyone listed in the IRCC immigration landing datasets and crosswalk,
then gets the first MB health insurance registry coverage episode ending after
landing, to calculate rates of outmigration of permanent residents in Manitoba.


******************************************************************************/

options nocenter;


*** Get a linked immigrant cohort from IRCC crosswalk and landing tables,;
*** keep one record per valid PHIN.;

%get_ircc_pr(cohort=, out=landing, dup=first);


*** coverage episodes pulled from registry using regcov macro;
*** ignore coverage gaps of 360 days or less;
%regcov(dsname=landing,
        startdt=01jan1985,
        enddt=30nov2023,
        id=scrphin,
        outds=regcov,
        gapdays=360);

*** sort by person so we can merge with grouped landing data;
proc sort data=regcov;
   by scrphin covdt enddt;
run;

*** merge in landing data to find coverage episodes during or after landing;
data landing;
   merge landing (in=in_la)
         regcov (in=in_co);
   by scrphin;
   if in_la and in_co and covdt > . and enddt >= landing_dt;
run;

*** keep one record per person, categorizing reason for leaving;
data landing;
   set landing;
   by scrphin covdt enddt;
   retain firstcanccode dropphin;
   if first.scrphin then do;
      firstcanccode = canccode;
      dropphin = 0;
   end;
   *** drop dupl. or registered in error (may show up under another phin);
   if canccode in ('W','5','6') then dropphin = 1;
   if last.scrphin & ^dropphin;
   canccode = firstcanccode;
   pop = 1;
   *** categorize outmigration;
   * Note: canccode A, C-R = Left Manitoba to some other known destination;
   length outmigration $18;
   if canccode = '0' then outmigration = 'Remained';
   else if canccode = 'V' then outmigration = 'Remained';
   else if canccode = '2' then outmigration = 'Remained';   * or 'Deceased';
   else if canccode = '7' then outmigration = 'Cannot Locate';
   else if canccode = 'A' or 'C' <= canccode <= 'R' then outmigration = 'Moved';
   else outmigration = 'Other';
   leftmb = (outmigration in ('Moved','Cannot Locate'));
   lyear = year(landing_dt);
   label lyear = 'Landing Year';
   label outmigration = 'Outmigration by Reason for Coverage Cancellation';
   label leftmb = 'Status of outmigration (Moved or Cannot Locate)';
run;


*** summarize outmigration counts by landing year and reason for cancellation;
proc tabulate data=landing;
   class lyear outmigration / order=formatted;
   var pop;
   table lyear, pop='Immigrants'*sum='N'*f=6.0
                outmigration*pop=''*(sum='N'*f=6.0 rowpctsum<pop>='%')
         / nocellmerge;
   title 'Immigrants leaving Manitoba, by landing year and reason';
run;


*** a little cleanup;
proc delete data=landing regcov;
run;