
/* Program combines main data with parents' and spouses' data into one dataset to be used in the analyses in STATA */

libname lib100 "W:\Tutkimusprojektit\Kirjanpito\Dokumentaatio\JAR\Datat\";
run;


data Kirjanpito_analyysi;
set lib100.kirjanpitodata_tekovuosi1_220925;   			

Ens_rikos_tvuosi = min(Ens_muu_rikos_tvuosi, Ens_talousrikos_tvuosi);

if Ens_rikos_tvuosi = . then Ens_rikos_tvuosi = 9999;

Alle_18v_rikos = 0;
if Ens_rikos_tvuosi - BIRTHYEAR < 18 then Alle_18v_rikos = 1;

if BUSINESS_SOCIALDEG = . then BUSINESS_SOCIALDEG = 0;
if EDULEVEL = . then EDULEVEL = 0;
if TOTALINCOME = . then TOTALINCOME = 0;
if LN_TOTALINCOME = . then LN_TOTALINCOME = 0;

OMAVARAISUUS3 = OMAVARAISUUS/100;
ROA3 = ROA/100;

SUM_PERCENTAGEOWN2 = SUM_PERCENTAGEOWN - PERCENTAGEOWN;

if MEAN_D_CONVICTIONS = . then MEAN_D_CONVICTIONS = 0;
if N_Board = . then N_Board = 1;

if D_Mean_D_muu_rikos_4 = . then D_Mean_D_muu_rikos_4 = 0;
if Mean_D_Muu_rikos_4 = . then Mean_D_Muu_rikos_4 = 0;
if D_Mean_D_Talous_kaikki = . then D_Mean_D_Talous_kaikki = 0;
if Mean_D_Talous_kaikki = . then Mean_D_Talous_kaikki = 0;

D_FMOB = 0;
D_FMCEO = 0;

if Mean_fmOB > 0 then D_FMOB = 1;
if Mean_fmCEO > 0 then D_FMCEO = 1;

D_FMOB_firm = 0;
D_FMCEO_firm = 0;

if Mean_fmOB_firm > 0 then D_FMOB_firm = 1;
if Mean_fmCEO_firm > 0 then D_FMCEO_firm = 1;

LN_AGE = log(AGE);
LN_EDULEVEL = log(EDULEVEL+1);
LN_Lifetime_SUM_DIVORCE = log(Lifetime_SUM_DIVORCE+1);
LN_Lifetime_SUM_Married = log(Lifetime_SUM_Married+1);

LN_CUM_SUM_DIVORCE = log(CUM_SUM_DIVORCE+1);
LN_CUM_SUM_Married = log(CUM_SUM_Married+1);

LN_NUM_CHILDREN_UNDER_18 = log(NUM_CHILDREN_UNDER_18+1);

LN_Boardsize = log(N_Board);

Lifetime_tuomio = 0;
if D_lifetime_muu_rikos_4_tuomio = 1 or D_lifetime_talous_kaikki_tuomio = 1 then Lifetime_tuomio = 1;

if talous_tuomio ^= . then Vuosi_dummy_talous_tuomio = 1;
if muu_rikos_4_tuomio ^= . then Vuosi_dummy_muu_rikos_4_tuomio = 1;
if talous_kaikki_tuomio ^= . then Vuosi_dummy_talous_kaikki_tuomio = 1;

LN_CUM_Muu_rikos_4_tuomio_2 = log(CUM_Muu_rikos_4_tuomio_2 + 1);
LN_CUM_Talous_kaikki_tuomio_2 = log(CUM_Talous_kaikki_tuomio_2 + 1);

LN_CUM_Muu_rikos_4_2 = log(CUM_Muu_rikos_4_2 + 1);
LN_CUM_Talous_kaikki_2 = log(CUM_Talous_kaikki_2 + 1);

if 1970 <= Birthyear <= 1979 then Birthyear_class = 1;
if 1980 <= Birthyear <= 1989 then Birthyear_class = 2;
if Birthyear >= 1990 then Birthyear_class = 3;

AGE_FIRST_MAAKUNTA = YEAR_FIRST - BIRTHYEAR;

if D_MARITAL_STATUS = . then delete;

if Current_ratio = . then delete;

if MAAKUNTA ^= 21 and MAAKUNTA_FIRST ^= 21;

Industry_2002_2 = Industry_2002 * 1;

I001 = 0; I002 = 0; I005 = 0; I010 = 0; I013 = 0; I014 = 0; I015 = 0; I016 = 0; I017 = 0; I018 = 0;
I019 = 0; I020 = 0; I021 = 0; I022 = 0; I023 = 0; I024 = 0; I025 = 0; I026 = 0; I027 = 0; I028 = 0;
I029 = 0; I030 = 0; I031 = 0; I032 = 0; I033 = 0; I034 = 0; I035 = 0; I036 = 0; I037 = 0; I040 = 0;
I041 = 0; I045 = 0; I050 = 0; I051 = 0; I052 = 0; I055 = 0; I060 = 0; I061 = 0; I062 = 0; I063 = 0;
I064 = 0; I065 = 0; I066 = 0; I067 = 0; I070 = 0; I071 = 0; I072 = 0; I073 = 0; I074 = 0; I075 = 0;
I080 = 0; I085 = 0; I090 = 0; I091 = 0; I092 = 0; I093 = 0; 

if Industry_2002_2 = 1 then I001 = 1; if Industry_2002_2 = 2 then I002 = 1; if Industry_2002_2 = 5 then I005 = 1; if Industry_2002_2 = 10 then I010 = 1; if Industry_2002_2 = 13 then I013 = 1;
if Industry_2002_2 = 14 then I014 = 1; if Industry_2002_2 = 15 then I015 = 1; if Industry_2002_2 = 16 then I016 = 1; if Industry_2002_2 = 17 then I017 = 1; if Industry_2002_2 = 18 then I018 = 1;
if Industry_2002_2 = 19 then I019 = 1; if Industry_2002_2 = 20 then I020 = 1; if Industry_2002_2 = 21 then I021 = 1; if Industry_2002_2 = 22 then I022 = 1; if Industry_2002_2 = 23 then I023 = 1;
if Industry_2002_2 = 24 then I024 = 1; if Industry_2002_2 = 25 then I025 = 1; if Industry_2002_2 = 26 then I026 = 1; if Industry_2002_2 = 27 then I027 = 1; if Industry_2002_2 = 28 then I028 = 1;
if Industry_2002_2 = 29 then I029 = 1; if Industry_2002_2 = 30 then I030 = 1; if Industry_2002_2 = 31 then I031 = 1; if Industry_2002_2 = 32 then I032 = 1; if Industry_2002_2 = 33 then I033 = 1;
if Industry_2002_2 = 34 then I034 = 1; if Industry_2002_2 = 35 then I035 = 1; if Industry_2002_2 = 36 then I036 = 1; if Industry_2002_2 = 37 then I037 = 1; if Industry_2002_2 = 40 then I040 = 1;
if Industry_2002_2 = 41 then I041 = 1; if Industry_2002_2 = 45 then I045 = 1; if Industry_2002_2 = 50 then I050 = 1; if Industry_2002_2 = 51 then I051 = 1; if Industry_2002_2 = 52 then I052 = 1;
if Industry_2002_2 = 55 then I055 = 1; if Industry_2002_2 = 60 then I060 = 1; if Industry_2002_2 = 61 then I061 = 1; if Industry_2002_2 = 62 then I062 = 1; if Industry_2002_2 = 63 then I063 = 1;
if Industry_2002_2 = 64 then I064 = 1; if Industry_2002_2 = 65 then I065 = 1; if Industry_2002_2 = 66 then I066 = 1; if Industry_2002_2 = 67 then I067 = 1; if Industry_2002_2 = 70 then I070 = 1;
if Industry_2002_2 = 71 then I071 = 1; if Industry_2002_2 = 72 then I072 = 1; if Industry_2002_2 = 73 then I073 = 1; if Industry_2002_2 = 74 then I074 = 1; if Industry_2002_2 = 75 then I075 = 1;
if Industry_2002_2 = 80 then I080 = 1; if Industry_2002_2 = 85 then I085 = 1; if Industry_2002_2 = 90 then I090 = 1; if Industry_2002_2 = 91 then I091 = 1; if Industry_2002_2 = 92 then I092 = 1;
if Industry_2002_2 = 93 then I093 = 1; 

/*--- Selection ---*/

*if D_FMOB_firm = 0;
*if year >= 2006 and omistus_firma_mukaan = 1;

drop ROA OMAVARAISUUS;
rename ROA3 = ROA OMAVARAISUUS3 = OMAVARAISUUS N_Board = Boardsize;

run;

proc sort data = kirjanpito_analyysi;
by ID YEAR firm descending asema;
run;

* Removes multiple positions in one firm;

data kirjanpito_analyysi_0;
set kirjanpito_analyysi;

if id = lag(id) and year = lag(year) and firm = lag(firm) then delete;

run;

* Calculates number of positions in different firms per person per year;

proc means noprint data = Kirjanpito_analyysi_0;

by ID YEAR;
var totalincome;
output out = Sum_years sum = sum_totalincome;

run;

* Parents' data;

data Kirjanpito_vanhemmat;
set lib100.kaikki_vanhemmat_tvuosi_250925; 

if m_sos_vanh_dummy = . then m_sos_vanh_dummy = 0;
if f_sos_vanh_dummy = . then f_sos_vanh_dummy = 0;
if p_sos_vanh_dummy = . then p_sos_vanh_dummy = 0;

m_biol_vanh_dummy = 1;
f_biol_vanh_dummy = 1;
p_biol_vanh_dummy = 1;

if m_sos_vanh_dummy = 1 then m_biol_vanh_dummy = 0;
if f_sos_vanh_dummy = 1 then f_biol_vanh_dummy = 0;
if m_sos_vanh_dummy = 1 or f_sos_vanh_dummy = 1 then p_biol_vanh_dummy = 0;

M_D_Talous_kaikki_3 = M_D_Talous_kaikki_2;
F_D_Talous_kaikki_3 = F_D_Talous_kaikki_2;

P_D_Talous_kaikki_3 = 0;
if M_D_Talous_kaikki_2 = 1 and F_D_Talous_kaikki_2 = 1 then P_D_Talous_kaikki_3 = 1;
if P_D_Talous_kaikki_3 = 1 then M_D_Talous_kaikki_3 = 0;
if P_D_Talous_kaikki_3 = 1 then F_D_Talous_kaikki_3 = 0;

P_D_Talous_kaikki_4 = 0;
if M_D_Talous_kaikki_2 = 1 or F_D_Talous_kaikki_2 = 1 then P_D_Talous_kaikki_4 = 1;
if P_D_Talous_kaikki_3 = 1 then P_D_Talous_kaikki_4 = 0;

M_D_muu_rikos_4_3 = M_D_muu_rikos_4_2;
F_D_muu_rikos_4_3 = F_D_muu_rikos_4_2;

P_D_muu_rikos_4_3 = 0;
if M_D_muu_rikos_4_2 = 1 and F_D_muu_rikos_4_2 = 1 then P_D_muu_rikos_4_3 = 1;
if P_D_muu_rikos_4_3 = 1 then M_D_muu_rikos_4_3 = 0;
if P_D_muu_rikos_4_3 = 1 then F_D_muu_rikos_4_3 = 0;

P_D_muu_rikos_4_4 = 0;
if M_D_muu_rikos_4_2 = 1 or F_D_muu_rikos_4_2 = 1 then P_D_muu_rikos_4_4 = 1;
if P_D_muu_rikos_4_3 = 1 then P_D_muu_rikos_4_4 = 0;

LN_M_CUM_Muu3_tuomio = log(M_CUM_Muu3_tuomio + 1);
LN_M_CUM_Muu4_tuomio = log(M_CUM_Muu4_tuomio + 1);
LN_M_CUM_Muu5_tuomio = log(M_CUM_Muu5_tuomio + 1);
LN_M_CUM_Talous_tuomio = log(M_CUM_Talous_tuomio + 1);
LN_M_CUM_Talous_kaikki_tuomio = log(M_CUM_Talous_kaikki_tuomio + 1);
LN_M_CUM_Kirjanpitorikos_tuomio = log(M_CUM_Kirjanpitorikos_tuomio + 1);

LN_F_CUM_Muu3_tuomio = log(F_CUM_Muu3_tuomio + 1);
LN_F_CUM_Muu4_tuomio = log(F_CUM_Muu4_tuomio + 1);
LN_F_CUM_Muu5_tuomio = log(F_CUM_Muu5_tuomio + 1);
LN_F_CUM_Talous_tuomio = log(F_CUM_Talous_tuomio + 1);
LN_F_CUM_Talous_kaikki_tuomio = log(F_CUM_Talous_kaikki_tuomio + 1);
LN_F_CUM_Kirjanpitorikos_tuomio = log(F_CUM_Kirjanpitorikos_tuomio + 1);

LN_P_CUM_Muu3_tuomio = log(P_CUM_Muu3_tuomio + 1);
LN_P_CUM_Muu4_tuomio = log(P_CUM_Muu4_tuomio + 1);
LN_P_CUM_Muu5_tuomio = log(P_CUM_Muu5_tuomio + 1);
LN_P_CUM_Talous_tuomio = log(P_CUM_Talous_tuomio + 1);
LN_P_CUM_Talous_kaikki_tuomio = log(P_CUM_Talous_kaikki_tuomio + 1);
LN_P_CUM_Kirjanpitorikos_tuomio = log(P_CUM_Kirjanpitorikos_tuomio + 1);

LN_M_CUM_Muu3 = log(M_CUM_Muu3 + 1);
LN_M_CUM_Muu4 = log(M_CUM_Muu4 + 1);
LN_M_CUM_Muu5 = log(M_CUM_Muu5 + 1);
LN_M_CUM_Talous = log(M_CUM_Talous + 1);
LN_M_CUM_Talous_kaikki = log(M_CUM_Talous_kaikki + 1);
LN_M_CUM_Kirjanpitorikos_kaikki = log(M_CUM_Kirjanpitorikos_kaikki + 1);

LN_F_CUM_Muu3 = log(F_CUM_Muu3 + 1);
LN_F_CUM_Muu4 = log(F_CUM_Muu4 + 1);
LN_F_CUM_Muu5 = log(F_CUM_Muu5 + 1);
LN_F_CUM_Talous = log(F_CUM_Talous + 1);
LN_F_CUM_Talous_kaikki = log(F_CUM_Talous_kaikki + 1);
LN_F_CUM_Kirjanpitorikos_kaikki = log(F_CUM_Kirjanpitorikos_kaikki + 1);

LN_P_CUM_Muu3 = log(P_CUM_Muu3 + 1);
LN_P_CUM_Muu4 = log(P_CUM_Muu4 + 1);
LN_P_CUM_Muu5 = log(P_CUM_Muu5 + 1);
LN_P_CUM_Talous = log(P_CUM_Talous + 1);
LN_P_CUM_Talous_kaikki = log(P_CUM_Talous_kaikki + 1);
LN_P_CUM_Kirjanpitorikos_kaikki = log(P_CUM_Kirjanpitorikos_kaikki + 1);

M_LN_EDULEVEL = log(M_EDULEVEL+1);
F_LN_EDULEVEL = log(F_EDULEVEL+1);
P_LN_EDULEVEL = log(P_EDULEVEL+1);

if P_TOTALINCOME = . then P_TOTALINCOME = 0;
if P_SALARYTOTAL = . then P_SALARYTOTAL = 0;
if P_CAPITALTOTAL = . then P_CAPITALTOTAL = 0;

P_LN_TOTALINCOME = log(P_TOTALINCOME+1);
P_LN_SALARYTOTAL = log(P_SALARYTOTAL+1);
P_LN_CAPITALTOTAL = log(P_CAPITALTOTAL+1);

P_AGE = ((F_AGE + M_AGE) / 2);

if M_AGE = . then P_AGE = F_AGE;
if F_AGE = . then P_AGE = M_AGE;

P_LN_AGE = log(P_AGE);

M_LN_Lifetime_SUM_DIVORCE = log(M_Lifetime_SUM_DIVORCE+1);
F_LN_Lifetime_SUM_DIVORCE = log(F_Lifetime_SUM_DIVORCE+1);
P_LN_Lifetime_SUM_DIVORCE = log(P_Lifetime_SUM_DIVORCE+1);

M_LN_CUM_SUM_DIVORCE = log(M_CUM_SUM_DIVORCE+1);
F_LN_CUM_SUM_DIVORCE = log(F_CUM_SUM_DIVORCE+1);
P_LN_CUM_SUM_DIVORCE = log(P_CUM_SUM_DIVORCE+1);

if M_D_lifetime_muu_rikos_4_tuomio = 1 or M_D_lifetime_KPrikos_tuomio = 1 or M_D_lifetime_talous_tuomio = 1 then M_Lifetime_tuomio = 1;
if F_D_lifetime_muu_rikos_4_tuomio = 1 or F_D_lifetime_KPrikos_tuomio = 1 or F_D_lifetime_talous_tuomio = 1 then F_Lifetime_tuomio = 1;
if P_D_lifetime_muu_rikos_4_tuomio = 1 or P_D_lifetime_KPrikos_tuomio = 1 or P_D_lifetime_talous_tuomio = 1 then P_Lifetime_tuomio = 1;

if P_D_talous_kaikki_tuomio = 0 and P_D_muu_rikos_4_tuomio = 0 then P_D_kaikki_tuomiot = 0;
if P_D_talous_kaikki_tuomio = 1 or P_D_muu_rikos_4_tuomio = 1 then P_D_kaikki_tuomiot = 1;

if P_D_talous_kaikki_2 = 0 and P_D_muu_rikos_4_2 = 0 then P_D_kaikki_rikos = 0;
if P_D_talous_kaikki_2 = 1 or P_D_muu_rikos_4_2 = 1 then P_D_kaikki_rikos = 1;

P_D_Finan_rep_fraud = 0;
P_D_Tax_relat_fraud = 0;
P_D_Other_finan_crime = 0;

if P_CUM_Finan_rep_fraud > 0 then P_D_Finan_rep_fraud = 1;
if P_CUM_Tax_relat_fraud > 0 then P_D_Tax_relat_fraud = 1;
if P_CUM_other_finan_crime > 0 then P_D_Other_finan_crime = 1;

P_D_Talous_kaikki_sakko_2 = 0;
P_D_muu_rikos_4_sakko_2 = 0;

P_D_Talous_kaikki_sakko_tuomio_2 = 0;
P_D_muu_rikos_4_sakko_tuomio_2 = 0;

P_D_Talous_kaik_ehdol_vankeus_3 = P_D_Talous_kaik_ehdol_vankeus_2;
if P_D_Talous_kaik_ehdot_vankeus_2 = 1 then P_D_Talous_kaik_ehdol_vankeus_3 = 0;

P_D_muu_rikos_4_ehdol_vankeus_3 = P_D_muu_rikos_4_ehdol_vankeus_2;
if P_D_muu_rikos_4_ehdot_vankeus_2 = 1 then P_D_muu_rikos_4_ehdol_vankeus_3 = 0;

if P_D_Talous_kaikki_2 = 1 and P_D_Talous_kaikki_vankeus_2 = 0 then P_D_Talous_kaikki_sakko_2 = 1;
if P_D_muu_rikos_4_2 = 1 and P_D_muu_rikos_4_vankeus_2 = 0 then P_D_muu_rikos_4_sakko_2 = 1;

if P_D_Talous_kaikki_tuomio_2 = 1 and P_D_Talous_kaikki_vankeus_tuom_2 = 0 then P_D_Talous_kaikki_sakko_tuomio_2 = 1;
if P_D_muu_rikos_4_tuomio_2 = 1 and P_D_muu_rikos_4_vankeus_tuomio_2 = 0 then P_D_muu_rikos_4_sakko_tuomio_2 = 1;

P_D_talous_kaikki_3kat = 0;
if P_D_Talous_kaikki_sakko_2 = 1 then P_D_Talous_kaikki_3kat = 1;
if P_D_Talous_kaikki_vankeus_2 = 1 then P_D_Talous_kaikki_3kat = 2;

P_D_muu_rikos_4_3kat = 0;
if P_D_muu_rikos_4_sakko_2 = 1 then P_D_muu_rikos_4_3kat = 1;
if P_D_muu_rikos_4_vankeus_2 = 1 then P_D_muu_rikos_4_3kat = 2;

if P_D_MARITAL_STATUS = . then delete;

vah_toinen_vanh_mm = 0; if m_skunta = '200' or f_skunta = '200' then vah_toinen_vanh_mm = 1;
molemmat_vanh_mm = 0; if m_skunta = '200' and f_skunta = '200' then molemmat_vanh_mm = 1;

run;

proc sort data = kirjanpito_vanhemmat;
by id year;
run;

* Spouses' data;

data puolisot;
set lib100.puolisoiden_datat_011025;
run;

proc sort data = puolisot;
by id year;
run;

* Combines main data with parents' and spouses' data;

data Kirjanpito_analyysi_2;
merge kirjanpito_analyysi_0 kirjanpito_vanhemmat puolisot sum_years;
by id year;

LN_Pa_Length = log(Pa_Length + 1);

Pa_Length_5v = 0;
if Pa_Length >= 5 then Pa_Length_5v = 1;

if firm = '' then delete;
if m_id = '' and f_id = '' then delete;

if D_lifetime_Talous_tuomio = 0 then vaste_talous = 0;
if D_lifetime_Muu_rikos_4_tuomio = 0 then vaste_muu_rikos_4 = 0;
if D_lifetime_Talous_kaikki_tuomio = 0 then vaste_talous_kaikki = 0;

if Vuosi_dummy_talous_tuomio = 1 then vaste_talous = 1;
if Vuosi_dummy_muu_rikos_4_tuomio = 1 then vaste_muu_rikos_4 = 1;
if Vuosi_dummy_talous_kaikki_tuomio = 1 then vaste_talous_kaikki = 1;

if kirjanpitorikos_tuomio = . then kirjanpitorikos_tuomio = 0;
if talous_tuomio = . then talous_tuomio = 0;
if muu_rikos_4_tuomio = . then muu_rikos_4_tuomio = 0;
if muu_rikos_5_tuomio = . then muu_rikos_5_tuomio = 0;
if talous_kaikki_tuomio = . then talous_kaikki_tuomio = 0;
if liikenne_tuomio = . then liikenne_tuomio = 0;

if kirjanpitorikos_kaikki = . then kirjanpitorikos_kaikki = 0;
if talous = . then talous = 0;
if muu_rikos_4 = . then muu_rikos_4 = 0;
if muu_rikos_5 = . then muu_rikos_5 = 0;
if talous_kaikki = . then talous_kaikki = 0;

D_Vaste_kirjanpitorikos_tuomio = 0;
D_Vaste_talous_tuomio = 0;
D_Vaste_muu_rikos_4_tuomio = 0;
D_Vaste_muu_rikos_5_tuomio = 0;
D_Vaste_talous_kaikki_tuomio = 0;
D_Vaste_liikenne_tuomio = 0;

D_Vaste_Finan_rep_fraud_tuomio = 0;
D_Vaste_Tax_relat_fraud_tuomio = 0;
D_Vaste_Other_finan_crime_tuomio = 0;

if kirjanpitorikos_tuomio > 0 then D_Vaste_kirjanpitorikos_tuomio = 1;
if talous_tuomio > 0 then D_Vaste_talous_tuomio = 1;
if muu_rikos_4_tuomio > 0 then D_Vaste_muu_rikos_4_tuomio = 1;
if muu_rikos_5_tuomio > 0 then D_Vaste_muu_rikos_5_tuomio = 1;
if talous_kaikki_tuomio > 0 then D_Vaste_talous_kaikki_tuomio = 1;
if liikenne_tuomio > 0 then D_Vaste_liikenne_tuomio = 1;

if Finan_rep_fraud_tuomio > 0 then D_Vaste_Finan_rep_fraud_tuomio = 1;
if Tax_relat_fraud_tuomio > 0 then D_Vaste_Tax_relat_fraud_tuomio = 1;
if Other_finan_crime_tuomio > 0 then D_Vaste_Other_finan_crime_tuomio = 1;

D_Vaste_kirjanpitorikos_kaikki = 0;
D_Vaste_kirjpitorikos_rikoslaki = 0;
D_Vaste_talous = 0;
D_Vaste_muu_rikos_4 = 0;
D_Vaste_muu_rikos_5 = 0;
D_Vaste_talous_kaikki = 0;

Vaste_talous_kaikki_num = 0;
if talous_kaikki ^= . then Vaste_talous_kaikki_num = talous_kaikki;
LN_Vaste_talous_kaikki_num = log(Vaste_talous_kaikki_num + 1);

Vaste_Finan_rep_fraud_num = 0;
if Finan_rep_fraud ^= . then Vaste_Finan_rep_fraud_num = Finan_rep_fraud;
LN_Vaste_Finan_rep_fraud_num = log(Vaste_Finan_rep_fraud_num + 1);

Vaste_Tax_relat_fraud_num = 0;
if Tax_relat_fraud ^= . then Vaste_Tax_relat_fraud_num = Tax_relat_fraud;
LN_Vaste_Tax_relat_fraud_num = log(Vaste_Tax_relat_fraud_num + 1);

Vaste_Other_finan_crime_num = 0;
if Other_finan_crime ^= . then Vaste_Other_finan_crime_num = Other_finan_crime;
LN_Vaste_Other_finan_crime_num = log(Vaste_Other_finan_crime_num + 1);

D_Vaste_talous_kaikki_vankeus = 0;
D_Vaste_talous_kaikki_3kat = 0;

D_Vaste_Finan_rep_fraud = 0;
D_Vaste_Tax_relat_fraud = 0;
D_Vaste_Other_finan_crime = 0;
D_Vaste_Velallisen_rikos = 0;

if kirjanpitorikos_kaikki > 0 then D_Vaste_kirjanpitorikos_kaikki = 1;
if talous > 0 then D_Vaste_talous = 1;
if muu_rikos_4 > 0 then D_Vaste_muu_rikos_4 = 1;
if muu_rikos_5 > 0 then D_Vaste_muu_rikos_5 = 1;
if talous_kaikki > 0 then D_Vaste_talous_kaikki = 1;
if talous_kaikki_vankeus > 0 then D_Vaste_talous_kaikki_vankeus = 1;

if talous_kaikki > 0 and talous_kaikki_vankeus = . then D_Vaste_talous_kaikki_3kat = 1;
if talous_kaikki > 0 and talous_kaikki_vankeus > 0 then D_Vaste_talous_kaikki_3kat = 2;

if Finan_rep_fraud > 0 then D_Vaste_Finan_rep_fraud = 1;
if Tax_relat_fraud > 0 then D_Vaste_Tax_relat_fraud = 1;
if Other_finan_crime > 0 then D_Vaste_Other_finan_crime = 1;
if Velallisen_rikos > 0 then D_Vaste_Velallisen_rikos = 1;

if D_Vaste_talous_kaikki_tuomio = 0 and D_Vaste_talous_kaikki = 1 then D_Vaste_talous_kaikki_tuomio = .;
if D_Vaste_talous_kaikki_tuomio = . then D_Vaste_Finan_rep_fraud_tuomio = .;
if D_Vaste_talous_kaikki_tuomio = . then D_Vaste_Tax_relat_fraud_tuomio = .;
if D_Vaste_talous_kaikki_tuomio = . then D_Vaste_Other_finan_crime_tuomio = .;
if D_Vaste_muu_rikos_4_tuomio = 0 and D_Vaste_muu_rikos_4 = 1 then D_Vaste_muu_rikos_4_tuomio = .;

if D_Vaste_talous_kaikki = 0 then Vaste_sakot_talous_kaikki = 0;
if D_Vaste_talous_kaikki = 0 then Vaste_vankeus_talous_kaikki = 0;

if D_Vaste_talous_kaikki = 0 then Vaste_sakot_finan_rep_fraud = 0;
if D_Vaste_talous_kaikki = 0 then Vaste_vankeus_finan_rep_fraud = 0;

if D_Vaste_talous_kaikki = 0 then Vaste_sakot_tax_relat_fraud = 0;
if D_Vaste_talous_kaikki = 0 then Vaste_vankeus_tax_relat_fraud = 0;

if D_Vaste_talous_kaikki = 0 then Vaste_sakot_other_finan = 0;
if D_Vaste_talous_kaikki = 0 then Vaste_vankeus_other_finan = 0;

if sakot_talous_kaikki > 0 then Vaste_sakot_talous_kaikki = sakot_talous_kaikki;
if vankeus_talous_kaikki > 0 then Vaste_vankeus_talous_kaikki = vankeus_talous_kaikki;

if sakot_finan_rep_fraud > 0 then Vaste_sakot_finan_rep_fraud = sakot_finan_rep_fraud;
if vankeus_finan_rep_fraud > 0 then Vaste_vankeus_finan_rep_fraud = vankeus_finan_rep_fraud;

if sakot_tax_relat_fraud > 0 then Vaste_sakot_tax_relat_fraud = sakot_tax_relat_fraud;
if vankeus_tax_relat_fraud > 0 then Vaste_vankeus_tax_relat_fraud = vankeus_tax_relat_fraud;

if sakot_other_finan > 0 then Vaste_sakot_other_finan = sakot_other_finan;
if vankeus_other_finan > 0 then Vaste_vankeus_other_finan = vankeus_other_finan;

Vaste_LN_sakot_talous_kaikki = log(Vaste_sakot_talous_kaikki + 1);
Vaste_LN_vankeus_talous_kaikki = log(Vaste_vankeus_talous_kaikki + 1);

Vaste_LN_sakot_finan_rep_fraud = log(Vaste_sakot_finan_rep_fraud + 1);
Vaste_LN_vankeus_finan_rep_fraud = log(Vaste_vankeus_finan_rep_fraud + 1); 

Vaste_LN_sakot_tax_relat_fraud = log(Vaste_sakot_tax_relat_fraud + 1); 
Vaste_LN_vankeus_tax_relat_fraud = log(Vaste_vankeus_tax_relat_fraud + 1); 

Vaste_LN_sakot_other_finan = log(Vaste_sakot_other_finan + 1); 
Vaste_LN_vankeus_other_finan = log(Vaste_vankeus_other_finan + 1); 

if D_Vaste_talous_kaikki = 0 or D_Vaste_Finan_rep_fraud = 1 then Finan_mukaan = 1;
if D_Vaste_talous_kaikki = 0 or D_Vaste_Tax_relat_fraud = 1 then Tax_mukaan = 1;
if D_Vaste_talous_kaikki = 0 or D_Vaste_Other_finan_crime = 1 then Other_mukaan = 1;

Tuomio = 0;

if muu_rikos_4_tuomio = 1 or talous_kaikki_tuomio = 1 then tuomio = 1;

YEAR = YEAR + 1;

if YEAR < 2020;

Vanh_omistus = 0;

if firm = m_firm1 or firm = m_firm2 or firm = m_firm3 or firm = m_firm4 or firm = m_firm5 or firm = m_firm6 or firm = m_firm7 or firm = m_firm8 or firm = m_firm9 or 
firm = m_firm10 or firm = m_firm11 or firm = m_firm12 or firm = m_firm13 or firm = m_firm14 or firm = m_firm15 or firm = m_firm16 or firm = m_firm17 or firm = m_firm18 or 
firm = m_firm19 or firm = m_firm20 or firm = m_firm21 or firm = m_firm22 or firm = m_firm23 or 
firm = f_firm1 or firm = f_firm2 or firm = f_firm3 or firm = f_firm4 or firm = f_firm5 or firm = f_firm6 or firm = f_firm7 or firm = f_firm8 or firm = f_firm9 or 
firm = f_firm10 or firm = f_firm11 or firm = f_firm12 or firm = f_firm13 or firm = f_firm14 or firm = f_firm15 or firm = f_firm16 or firm = f_firm17 or firm = f_firm18 or
firm = f_firm19 or firm = f_firm20 or firm = f_firm21 or firm = f_firm22 or firm = f_firm23 or firm = f_firm24 or firm = f_firm25 or firm = f_firm26 or firm = f_firm27 or
firm = f_firm28 or firm = f_firm29 or firm = f_firm30 or firm = f_firm31 or firm = f_firm32 or firm = f_firm33 then vanh_omistus = 1;

drop sum_totalincome _type_ _freq_ m_firm1 m_firm2 m_firm3 m_firm4 m_firm5 m_firm6 m_firm7 m_firm8 m_firm9 
m_firm10  m_firm11  m_firm12  m_firm13  m_firm14 m_firm15 m_firm16 m_firm17 m_firm18
m_firm19  m_firm20  m_firm21  m_firm22  m_firm23
f_firm1  f_firm2  f_firm3  f_firm4  f_firm5  f_firm6  f_firm7  f_firm8  f_firm9 
f_firm10  f_firm11  f_firm12  f_firm13  f_firm14  f_firm15  f_firm16  f_firm17  f_firm18 
f_firm19  f_firm20  f_firm21  f_firm22  f_firm23  f_firm24  f_firm25  f_firm26  f_firm27 
f_firm28  f_firm29  f_firm30  f_firm31  f_firm32  f_firm33;
 
run; 

* Winsorize firm variables;

proc means P1 P5 P95 P99 data = Kirjanpito_analyysi_2 noprint;
var TASE LNTASE ROA Current_Ratio;
output out = kvanKirjanpito P1 = P1_TASE P1_LNTASE P1_ROA P1_Current_Ratio P5 = P5_TASE P5_LNTASE P5_ROA P5_Current_Ratio
P95 = P95_TASE P95_LNTASE P95_ROA P95_Current_Ratio P99 = P99_TASE P99_LNTASE P99_ROA P99_Current_Ratio;
run;

data Kirjanpito_analyysi_4;

if _n_= 1 then set kvanKirjanpito;

set Kirjanpito_analyysi_2;

TASE2 = TASE;
LNTASE2 = LNTASE;
ROA2 = ROA;
Current_Ratio2 = Current_Ratio;

if TASE > P99_TASE then TASE2=P99_TASE;
if TASE < P1_TASE then TASE2=P1_TASE;

if LNTASE > P99_LNTASE then LNTASE2=P99_LNTASE;
if LNTASE < P1_LNTASE then LNTASE2=P1_LNTASE;

if ROA > P99_ROA then ROA2=P99_ROA;
if ROA < P1_ROA then ROA2=P1_ROA;

if Current_Ratio > P95_Current_Ratio then Current_Ratio2=P95_Current_Ratio;
if Current_Ratio < P5_Current_Ratio then Current_Ratio2=P5_Current_Ratio;

if Current_Ratio = . then Current_Ratio2 = .;

OMAVARAISUUS2 = OMAVARAISUUS;
if OMAVARAISUUS < 0 then OMAVARAISUUS2 = 0;

drop _TYPE_ _FREQ_ P1_TASE P1_LNTASE P99_TASE P99_LNTASE P1_ROA P99_ROA P1_Current_Ratio P99_Current_Ratio P5_TASE
P5_LNTASE P5_ROA P5_Current_Ratio P95_TASE P95_LNTASE P95_ROA P95_Current_Ratio;

run;

* Winsorize personal variables;

proc means P1 P99 data = Kirjanpito_analyysi_4 noprint;
var SALARYTOTAL LN_TOTALINCOME TOTALINCOME M_LN_TOTALINCOME F_LN_TOTALINCOME P_LN_TOTALINCOME P_TOTALINCOME Pa_LN_TOTALINCOME Pa_TOTALINCOME;
output out = kvanKirjanpito P1 = P1_SALARYTOTAL P1_LN_TOTALINCOME P1_TOTALINCOME P1_M_LN_TOTALINCOME P1_F_LN_TOTALINCOME 
P1_P_LN_TOTALINCOME P1_P_TOTALINCOME P1_Pa_LN_TOTALINCOME P1_Pa_TOTALINCOME P99 = P99_SALARYTOTAL P99_LN_TOTALINCOME P99_TOTALINCOME P99_M_LN_TOTALINCOME
P99_F_LN_TOTALINCOME P99_P_LN_TOTALINCOME P99_P_TOTALINCOME P99_Pa_LN_TOTALINCOME P99_Pa_TOTALINCOME;
run;

data Kirjanpito_analyysi_5;

if _n_= 1 then set kvanKirjanpito;

set Kirjanpito_analyysi_4;

SALARYTOTAL2 = SALARYTOTAL;

LN_TOTALINCOME2 = LN_TOTALINCOME;
TOTALINCOME2 = TOTALINCOME;

M_LN_TOTALINCOME2 = M_LN_TOTALINCOME;
F_LN_TOTALINCOME2 = F_LN_TOTALINCOME;

P_LN_TOTALINCOME2 = P_LN_TOTALINCOME;
P_TOTALINCOME2 = P_TOTALINCOME;

Pa_LN_TOTALINCOME2 = Pa_LN_TOTALINCOME;
Pa_TOTALINCOME2 = Pa_TOTALINCOME;


if SALARYTOTAL > P99_SALARYTOTAL and SALARYTOTAL ne . then SALARYTOTAL2=P99_SALARYTOTAL;
if SALARYTOTAL < P1_SALARYTOTAL and SALARYTOTAL ne . then SALARYTOTAL2=P1_SALARYTOTAL;

if LN_TOTALINCOME > P99_LN_TOTALINCOME and LN_TOTALINCOME ne . then LN_TOTALINCOME2=P99_LN_TOTALINCOME;
if LN_TOTALINCOME < P1_LN_TOTALINCOME and LN_TOTALINCOME ne . then LN_TOTALINCOME2=P1_LN_TOTALINCOME;

if TOTALINCOME > P99_TOTALINCOME and TOTALINCOME ne . then TOTALINCOME2=P99_TOTALINCOME;
if TOTALINCOME < P1_TOTALINCOME and TOTALINCOME ne . then TOTALINCOME2=P1_TOTALINCOME;

if M_LN_TOTALINCOME > P99_M_LN_TOTALINCOME and M_LN_TOTALINCOME ne . then M_LN_TOTALINCOME2=P99_M_LN_TOTALINCOME;
if M_LN_TOTALINCOME < P1_M_LN_TOTALINCOME and M_LN_TOTALINCOME ne . then M_LN_TOTALINCOME2=P1_M_LN_TOTALINCOME;

if F_LN_TOTALINCOME > P99_F_LN_TOTALINCOME and F_LN_TOTALINCOME ne . then F_LN_TOTALINCOME2=P99_F_LN_TOTALINCOME;
if F_LN_TOTALINCOME < P1_F_LN_TOTALINCOME and F_LN_TOTALINCOME ne . then F_LN_TOTALINCOME2=P1_F_LN_TOTALINCOME;

if P_LN_TOTALINCOME > P99_P_LN_TOTALINCOME and P_LN_TOTALINCOME ne . then P_LN_TOTALINCOME2=P99_P_LN_TOTALINCOME;
if P_LN_TOTALINCOME < P1_P_LN_TOTALINCOME and P_LN_TOTALINCOME ne . then P_LN_TOTALINCOME2=P1_P_LN_TOTALINCOME;

if P_TOTALINCOME > P99_P_TOTALINCOME and P_TOTALINCOME ne . then P_TOTALINCOME2=P99_P_TOTALINCOME;
if P_TOTALINCOME < P1_P_TOTALINCOME and P_TOTALINCOME ne . then P_TOTALINCOME2=P1_P_TOTALINCOME;

if Pa_LN_TOTALINCOME > P99_Pa_LN_TOTALINCOME and Pa_LN_TOTALINCOME ne . then Pa_LN_TOTALINCOME2=P99_Pa_LN_TOTALINCOME;
if Pa_LN_TOTALINCOME < P1_Pa_LN_TOTALINCOME and Pa_LN_TOTALINCOME ne . then Pa_LN_TOTALINCOME2=P1_Pa_LN_TOTALINCOME;

if Pa_TOTALINCOME > P99_Pa_TOTALINCOME and Pa_TOTALINCOME ne . then Pa_TOTALINCOME2=P99_Pa_TOTALINCOME;
if Pa_TOTALINCOME < P1_Pa_TOTALINCOME and Pa_TOTALINCOME ne . then Pa_TOTALINCOME2=P1_Pa_TOTALINCOME;

drop _TYPE_ _FREQ_ P1_SALARYTOTAL P1_LN_TOTALINCOME P1_TOTALINCOME  P1_P_LN_TOTALINCOME P1_P_TOTALINCOME P1_Pa_LN_TOTALINCOME
P1_Pa_TOTALINCOME P99_SALARYTOTAL P99_LN_TOTALINCOME P99_TOTALINCOME
P99_P_LN_TOTALINCOME P99_P_TOTALINCOME P99_Pa_LN_TOTALINCOME P99_Pa_TOTALINCOME ;

run;

*  Data on all individuals;

data kirjanpito_kaikki;
set Kirjanpito_analyysi_5;
run;

proc sort data = kirjanpito_kaikki;
by ID YEAR firm descending asema;
run;

data kirjanpito_kaikki_5;
set kirjanpito_kaikki;

if id = lag(id) and year = lag(year) and firm = lag(firm) then delete;

run;

* Identifies persons with crimes during the sample period;

proc means noprint data = kirjanpito_kaikki_5;
by id;
var D_Vaste_talous_kaikki D_Vaste_muu_rikos_4;
output out = rikokset mean = Talous_rikos_henk Muu_rikos_henk;
run;

* Identifies firms where executives' parents are owners;

proc means noprint data = kirjanpito_kaikki_5;
by id;
var vanh_omistus;
output out = vanhempien_omistus mean = Mean_vanh_omistus;
run;

proc sort data = kirjanpito_kaikki_5;
by id year;
run;

* Adds crime rate data from childhood municipalities;

data Rikokset_lapsuus;
set lib100.rikokset_lapsuus_030624;
run;

data Rikokset_lapsuus_maakunta;
set lib100.rikokset_lapsuus_maakunta_030624;
run;

data Syntymakunta_rikos;
set lib100.Syntymakunta_rikos_110624;
run;

proc sort data = Syntymakunta_rikos;
by id;
run;

data Muuttokunta_rikos;
set lib100.Muuttokunta_rikos_011025;
run;

proc sort data = Muuttokunta_rikos;
by id;
run;

* Bias score variables;

data demog_lisamuuttujat;
set lib100.demog_lisamuuttujat_011025;
run;

proc sort data = demog_lisamuuttujat;
by id;
run;

proc means noprint data = kirjanpito_kaikki_5;
by id;
var vah_toinen_vanh_mm;
output out = mm_vanh mean = mean_mm_vanh;
run;

data kirjanpito_kaikki_6;
merge kirjanpito_kaikki_5 rikokset Rikokset_lapsuus Rikokset_lapsuus_maakunta Syntymakunta_rikos Muuttokunta_rikos vanhempien_omistus demog_lisamuuttujat mm_vanh;
by id;

if firm = '' then delete;

if Talous_rikos_henk > 0 and D_Vaste_talous_kaikki = 0 then delete;
*if Muu_rikos_henk > 0 and D_Vaste_Muu_rikos_4 = 0 then delete;

* Identifies executives from immigrant families;

D_mean_mm_vanh = 0;
if mean_mm_vanh > 0 then D_mean_mm_vanh = 1;

drop _type_ _freq_;

run;

* Executives may have positions in several firms in a given year (calculate average firm variables per year);

proc means noprint data = kirjanpito_kaikki_6;

by ID YEAR;
var OMAVARAISUUS2 TASE2 LNTASE2 ROA2 Current_Ratio2 Boardsize LN_Boardsize Mean_D_talous_kaikki Mean_D_muu_rikos_4
D_Mean_D_talous_kaikki D_Mean_D_muu_rikos_4 CEO_Duality;

output out = Means_kaikki mean = Mean_OMAVARAISUUS2 Mean_TASE2 Mean_LNTASE2 Mean_ROA2 Mean_Current_Ratio2 Mean_Boardsize Mean_LN_Boardsize 
Mean_Mean_D_talous_kaikki Mean_Mean_D_muu_rikos_4 Mean_D_Mean_D_talous_kaikki Mean_D_Mean_D_muu_rikos_4 Mean_CEO_Duality;

run;

proc means noprint data = kirjanpito_kaikki_6;
by ID YEAR;
where omistus_firma_mukaan = 1;
var PERCENTAGEOWN SUM_PERCENTAGEOWN2;

output out = Means_kaikki_2 mean = MEAN_PERCENTAGEOWN MEAN_SUM_PERCENTAGEOWN2;

run;

proc means noprint data = kirjanpito_kaikki_6;

by ID YEAR;
var I001 I002 I005 I010 I013 I014 I015 I016 I017 I018 I019 I020 I021 I022 I023 I024 I025 I026 I027 I028 
I029 I030 I031 I032 I033 I034 I035 I036 I037 I040 I041 I045 I050 I051 I052 I055 I060 I061 I062 I063 
I064 I065 I066 I067 I070 I071 I072 I073 I074 I075 I080 I085 I090 I091 I092 I093;

output out = Sum_kaikki sum = sum_I001 sum_I002 sum_I005 sum_I010 sum_I013 sum_I014 sum_I015 sum_I016 sum_I017
sum_I018 sum_I019 sum_I020 sum_I021 sum_I022 sum_I023 sum_I024 sum_I025 sum_I026 sum_I027 sum_I028 
sum_I029 sum_I030 sum_I031 sum_I032 sum_I033 sum_I034 sum_I035 sum_I036 sum_I037 sum_I040 sum_I041
sum_I045 sum_I050 sum_I051 sum_I052 sum_I055 sum_I060 sum_I061 sum_I062 sum_I063 
sum_I064 sum_I065 sum_I066 sum_I067 sum_I070 sum_I071 sum_I072 sum_I073 sum_I074
sum_I075 sum_I080 sum_I085 sum_I090 sum_I091 sum_I092 sum_I093;

run;

proc sort nodupkey data = kirjanpito_kaikki_6;
by ID YEAR;
run;

data kirjanpito_kaikki_22;
merge kirjanpito_kaikki_6 Means_kaikki Means_kaikki_2 Sum_kaikki;
by ID YEAR;

if LV = . then delete;
if ID = '' then delete;

Score_muuttuja_2 = 0;

if Vuokra_17v = 1 then Score_muuttuja_2 + 1;
if Isan_sose_asema_17v = 1 then Score_muuttuja_2 + 1;
if Vanhempien_tulot_17v = 1 then Score_muuttuja_2 + 1;
if D_mean_mm_vanh = 1 then Score_muuttuja_2 + 1;
if M_EDULEVEL = 0 or F_EDULEVEL = 0 then Score_muuttuja_2 + 1;

if Score_muuttuja_2 = 0 then Score_muuttuja_2_3lk = 0;
if Score_muuttuja_2 = 1 then Score_muuttuja_2_3lk = 1;
if Score_muuttuja_2 >= 2 then Score_muuttuja_2_3lk = 2;

Score_vuokra_17v = 0;
Score_isan_sose_asema_17v = 0;
Score_vanh_tulot_17v = 0;
Score_vanh_maahanmuuttaja = 0;
Score_vanh_koulutus = 0;

if Vuokra_17v = 1 then Score_vuokra_17v = 1;
if Isan_sose_asema_17v = 1 then Score_isan_sose_asema_17v = 1;
if Vanhempien_tulot_17v = 1 then Score_vanh_tulot_17v = 1;
if D_mean_mm_vanh = 1 then Score_vanh_maahanmuuttaja = 1;
if M_EDULEVEL = 0 or F_EDULEVEL = 0 then Score_vanh_koulutus = 1;

Max_Industry = max(sum_I001, sum_I002, sum_I005, sum_I010, sum_I013, sum_I014, sum_I015, sum_I016, sum_I017,
sum_I018, sum_I019, sum_I020, sum_I021, sum_I022, sum_I023, sum_I024, sum_I025, sum_I026, sum_I027, sum_I028,
sum_I029, sum_I030, sum_I031, sum_I032, sum_I033, sum_I034, sum_I035, sum_I036, sum_I037, sum_I040, sum_I041,
sum_I045, sum_I050, sum_I051, sum_I052, sum_I055, sum_I060, sum_I061, sum_I062, sum_I063, 
sum_I064, sum_I065, sum_I066, sum_I067, sum_I070, sum_I071, sum_I072, sum_I073, sum_I074,
sum_I075, sum_I080, sum_I085, sum_I090, sum_I091, sum_I092, sum_I093);

if Max_Industry = sum_I001 then Industry_2002_3 = 1;
if Max_Industry = sum_I002 then Industry_2002_3 = 2;
if Max_Industry = sum_I005 then Industry_2002_3 = 5;
if Max_Industry = sum_I010 then Industry_2002_3 = 10;
if Max_Industry = sum_I013 then Industry_2002_3 = 13;
if Max_Industry = sum_I014 then Industry_2002_3 = 14;
if Max_Industry = sum_I015 then Industry_2002_3 = 15;
if Max_Industry = sum_I016 then Industry_2002_3 = 16;
if Max_Industry = sum_I017 then Industry_2002_3 = 17;
if Max_Industry = sum_I018 then Industry_2002_3 = 18;
if Max_Industry = sum_I019 then Industry_2002_3 = 19;
if Max_Industry = sum_I020 then Industry_2002_3 = 20;
if Max_Industry = sum_I021 then Industry_2002_3 = 21;
if Max_Industry = sum_I022 then Industry_2002_3 = 22;
if Max_Industry = sum_I023 then Industry_2002_3 = 23;
if Max_Industry = sum_I024 then Industry_2002_3 = 24;
if Max_Industry = sum_I025 then Industry_2002_3 = 25;
if Max_Industry = sum_I026 then Industry_2002_3 = 26;
if Max_Industry = sum_I027 then Industry_2002_3 = 27;
if Max_Industry = sum_I028 then Industry_2002_3 = 28;
if Max_Industry = sum_I029 then Industry_2002_3 = 29;
if Max_Industry = sum_I030 then Industry_2002_3 = 30;
if Max_Industry = sum_I031 then Industry_2002_3 = 31;
if Max_Industry = sum_I032 then Industry_2002_3 = 32;
if Max_Industry = sum_I033 then Industry_2002_3 = 33;
if Max_Industry = sum_I034 then Industry_2002_3 = 34;
if Max_Industry = sum_I035 then Industry_2002_3 = 35;
if Max_Industry = sum_I036 then Industry_2002_3 = 36;
if Max_Industry = sum_I037 then Industry_2002_3 = 37;
if Max_Industry = sum_I040 then Industry_2002_3 = 40;
if Max_Industry = sum_I041 then Industry_2002_3 = 41;
if Max_Industry = sum_I045 then Industry_2002_3 = 45;
if Max_Industry = sum_I050 then Industry_2002_3 = 50;
if Max_Industry = sum_I051 then Industry_2002_3 = 51; 
if Max_Industry = sum_I052 then Industry_2002_3 = 52;
if Max_Industry = sum_I055 then Industry_2002_3 = 55;
if Max_Industry = sum_I060 then Industry_2002_3 = 60;
if Max_Industry = sum_I061 then Industry_2002_3 = 61;
if Max_Industry = sum_I062 then Industry_2002_3 = 62;
if Max_Industry = sum_I063 then Industry_2002_3 = 63;
if Max_Industry = sum_I064 then Industry_2002_3 = 64;
if Max_Industry = sum_I065 then Industry_2002_3 = 65;
if Max_Industry = sum_I066 then Industry_2002_3 = 66;
if Max_Industry = sum_I067 then Industry_2002_3 = 67;
if Max_Industry = sum_I070 then Industry_2002_3 = 70;
if Max_Industry = sum_I071 then Industry_2002_3 = 71;
if Max_Industry = sum_I072 then Industry_2002_3 = 72;
if Max_Industry = sum_I073 then Industry_2002_3 = 73;
if Max_Industry = sum_I074 then Industry_2002_3 = 74;
if Max_Industry = sum_I075 then Industry_2002_3 = 75;
if Max_Industry = sum_I080 then Industry_2002_3 = 80;
if Max_Industry = sum_I085 then Industry_2002_3 = 85;
if Max_Industry = sum_I090 then Industry_2002_3 = 90;
if Max_Industry = sum_I091 then Industry_2002_3 = 91;
if Max_Industry = sum_I092 then Industry_2002_3 = 92;
if Max_Industry = sum_I093 then Industry_2002_3 = 93;

D_Mean_D_Mean_D_talous_kaikki = 0;
if Mean_D_Mean_D_talous_kaikki > 0 then D_Mean_D_Mean_D_talous_kaikki = 1;

D_Mean_D_Mean_D_muu_rikos_4 = 0;
if Mean_D_Mean_D_muu_rikos_4 > 0 then D_Mean_D_Mean_D_muu_rikos_4 = 1;

NUM_POSITIONS = _freq_;
LN_NUM_POSITIONS = log(_freq_);

drop  _type_ _freq_ omavaraisuus2 lntase2 roa2 Current_ratio2 ROE Current_ratio LV TASE LNTASE OMAVARAISUUS ROA LN_Boardsize
I001 I002 I005 I010 I013 I014 I015 I016 I017 I018 I019 I020 I021 I022 I023 I024 I025 I026 I027 I028 
I029 I030 I031 I032 I033 I034 I035 I036 I037 I040 I041 I045 I050 I051 I052 I055 I060 I061 I062 I063 
I064 I065 I066 I067 I070 I071 I072 I073 I074 I075 I080 I085 I090 I091 I092 I093;

run;

proc means noprint data = kirjanpito_kaikki_22;
var Mean_LNTASE2 MEAN_PERCENTAGEOWN MEAN_SUM_PERCENTAGEOWN2 Mean_OMAVARAISUUS2 Mean_LN_Boardsize;
output out = Median_tase median = Median_Mean_LNTASE2 Median_MEAN_PERCENTAGEOWN Median_MEAN_SUM_PERCENTAGEOWN2
Median_Mean_OMAVARAISUUS2 Median_Mean_LN_Boardsize;
run;

proc univariate noprint data = kirjanpito_kaikki_22;
var Mean_LNTASE2;
output out = Assets pctlpts = 33.3 66.7 pctlpre = P_Assets_;
run;

data kirjanpito_kaikki_2222;

if _n_= 1 then set Median_tase;
if _n_= 1 then set Assets;

set kirjanpito_kaikki_22;

FS_A_Median = 0;
if Mean_LNTASE2 < Median_Mean_LNTASE2 then FS_A_Median = 1;

Yli_Med_MEAN_PERCENTAGEOWN = 0;
if MEAN_PERCENTAGEOWN >= Median_MEAN_PERCENTAGEOWN then Yli_Med_MEAN_PERCENTAGEOWN = 1;

Yli_Med_MEAN_SUM_PERCENTAGEOWN2 = 0;
if MEAN_SUM_PERCENTAGEOWN2 >= Median_MEAN_SUM_PERCENTAGEOWN2 then Yli_Med_MEAN_SUM_PERCENTAGEOWN2 = 1;

EQ_Yli_Median = 0;
if Mean_OMAVARAISUUS2 >= Median_Mean_OMAVARAISUUS2 then EQ_Yli_Median = 1;

BoardSize_Yli_Median = 0;
if Mean_LN_Boardsize >= Median_Mean_LN_Boardsize then BoardSize_Yli_Median = 1;

MEAN_SUM_PERCENTAGEOWN_yli_0 = 0;
if MEAN_SUM_PERCENTAGEOWN2 > 0 then MEAN_SUM_PERCENTAGEOWN_yli_0 = 1;

Assets_3lk = 0;

if P_Assets_33_3 <= Mean_LNTASE2 <= P_Assets_66_7 then Assets_3lk = 1;
if Mean_LNTASE2 > P_Assets_66_7 then Assets_3lk = 2;


drop _freq_ _type_;

run;

proc sort data = kirjanpito_kaikki_2222;
by year Industry_2002_3 kunta;
run;

data kirjanpito_kaikki_22222;
set kirjanpito_kaikki_2222;

if year ne lag(year) or Industry_2002_3 ne lag(Industry_2002_3) or kunta ne lag(kunta) then YearIndustryMunicipality+1;

run;

proc sort data = kirjanpito_kaikki_22222;
by year Industry_2002_3 maakunta;
run;

data kirjanpito_kaikki_222222;
set kirjanpito_kaikki_22222;

if year ne lag(year) or Industry_2002_3 ne lag(Industry_2002_3) or maakunta ne lag(maakunta) then YearIndustryProvince+1;

run;

proc sort data = kirjanpito_kaikki_222222;
by year Industry_2002_3;
run;

data kirjanpito_kaikki_2222222;
set kirjanpito_kaikki_222222;

if year ne lag(year) or Industry_2002_3 ne lag(Industry_2002_3) then YearIndustry+1;

run;

proc sort data = kirjanpito_kaikki_2222222;
by year Kunta;
run;

data kirjanpito_kaikki_22222222;
set kirjanpito_kaikki_2222222;

if year ne lag(year) or Kunta ne lag(Kunta) then YearMunicipality+1;

run;

proc sort data = kirjanpito_kaikki_22222222;
by year Maakunta;
run;

data kirjanpito_kaikki_222222222;
set kirjanpito_kaikki_22222222;

if year ne lag(year) or Maakunta ne lag(Maakunta) then YearProvince+1;

run;

proc sort data = kirjanpito_kaikki_222222222;
by id year;
run;

proc means P1 P99 data = kirjanpito_kaikki_222222222 noprint;
var Talous_rikokset_lapsuus;
output out = Talous_rikokset_lapsuus_win P1 = P1_Talous_rikokset_lapsuus P99 = P99_Talous_rikokset_lapsuus;
run;

data kirjanpito_kaikki_22222222222;

if _n_= 1 then set Talous_rikokset_lapsuus_win;

set Kirjanpito_kaikki_222222222;

Talous_rikokset_lapsuus2 = Talous_rikokset_lapsuus;

if Talous_rikokset_lapsuus > P99_Talous_rikokset_lapsuus then Talous_rikokset_lapsuus2=P99_Talous_rikokset_lapsuus;

if Kaikki_rikokset_17 = . and YEAR_FIRST < BIRTHYEAR + 17 then Kaikki_rikokset_17 = Kaikki_rikokset_ensim;
if Talous_rikokset_17 = . and YEAR_FIRST < BIRTHYEAR + 17 then Talous_rikokset_17 = Talous_rikokset_ensim;
if Muut_rikokset_17 = . and YEAR_FIRST < BIRTHYEAR + 17 then Muut_rikokset_17 = Muut_rikokset_ensim;

LN_Kaikki_rikokset_vuosi = log(Kaikki_rikokset_vuosi + 1);
LN_Talous_rikokset_vuosi = log(Talous_rikokset_vuosi + 1);
LN_Muut_rikokset_vuosi = log(Muut_rikokset_vuosi + 1);
LN_Kaikki_rikokset_17 = log(Kaikki_rikokset_17 + 1);

LN_Kaikki_rikokset_m_vuosi = log(Kaikki_rikokset_maakunta_vuosi + 1);
LN_Talous_rikokset_m_vuosi = log(Talous_rikokset_maakunta_vuosi + 1);
LN_Muut_rikokset_m_vuosi = log(Muut_rikokset_maakunta_vuosi + 1);

LN_Kaikki_rikokset_lapsuus = log(Kaikki_rikokset_lapsuus + 1);
LN_Talous_rikokset_lapsuus = log(Talous_rikokset_lapsuus + 1);
LN_Talous_rikokset_lapsuus2 = log(Talous_rikokset_lapsuus2 + 1);
LN_Muut_rikokset_lapsuus = log(Muut_rikokset_lapsuus + 1);

Pa_LN_Kaikki_rikokset_lapsuus = log(Pa_Kaikki_rikokset_lapsuus + 1);
Pa_LN_Talous_rikokset_lapsuus = log(Pa_Talous_rikokset_lapsuus + 1);
Pa_LN_Muut_rikokset_lapsuus = log(Pa_Muut_rikokset_lapsuus + 1);

LN_Kaikki_rikokset_m_lapsuus = log(Kaikki_rikokset_lapsuus_maakunta + 1);
LN_Talous_rikokset_m_lapsuus = log(Talous_rikokset_lapsuus_maakunta + 1);
LN_Muut_rikokset_m_lapsuus = log(Muut_rikokset_lapsuus_maakunta + 1);

Finan_filter = 0; if D_Vaste_talous_kaikki = 1 and D_Vaste_Finan_rep_fraud = 0 then Finan_filter = 1;
Tax_filter = 0; if D_Vaste_talous_kaikki = 1 and D_Vaste_Tax_relat_fraud = 0 then Tax_filter = 1;
Other_filter = 0; if D_Vaste_talous_kaikki = 1 and D_Vaste_Other_finan_crime = 0 then Other_filter = 1;

run;

proc univariate data = kirjanpito_kaikki_22222222222;
var LN_Talous_rikokset_lapsuus;
output out = rikokset_kaikki_l pctlpts = 33.3 66.7 pctlpre = P_rikokset_;
run;

proc univariate data = kirjanpito_kaikki_22222222222;
var Kaikki_rikokset_17;
output out = rikokset_kaikki_p pctlpts = 33.3 66.7 pctlpre = P_rikokset_;
run;

data kirjanpito_kaikki_2;
set kirjanpito_kaikki_22222222222;

Kaikki_rikokset_17_3lk = 0;
if 128.67 < Kaikki_rikokset_17 <= 164.56 then Kaikki_rikokset_17_3lk = 1;
if Kaikki_rikokset_17 > 164.56 then Kaikki_rikokset_17_3lk = 2;

if Kaikki_rikokset_17 = . then Kaikki_rikokset_17_3lk = .;

Talous_rikokset_lapsuus_Q = 1;
if 1.7173950539 < LN_Talous_rikokset_lapsuus <= 2.2619714101 then Talous_rikokset_lapsuus_Q = 2;
if 2.2619714101 < LN_Talous_rikokset_lapsuus then Talous_rikokset_lapsuus_Q = 3;

Talous_rikokset_lapsuus_High = 0;
if 2.2619714101 < LN_Talous_rikokset_lapsuus then Talous_rikokset_lapsuus_High = 1;

run;

proc sort nodupkey data = kirjanpito_kaikki_2;
by id year;
run;

*  Data for analyses in panel B of Table 6;

/*proc means noprint data = kirjanpito_kaikki_2;
by id;
var D_Vaste_talous_kaikki P_D_talous_kaikki_2 P_D_muu_rikos_4_2 P_BUSINESS_SOCIALDEG P_EDULEVEL P_CUM_SUM_DIVORCE P_TOTALINCOME2 P_D_CUM_SUBSTANCE P_D_CUM_ALL_BUT_SUBS
D_Muu_rikos_4 MALE BUSINESS_SOCIALDEG EDULEVEL CUM_SUM_DIVORCE TOTALINCOME2 D_CUM_SUBSTANCE D_CUM_ALL_BUT_SUBS D_FMOB LN_NUM_POSITIONS Mean_OMAVARAISUUS2 Mean_LNTASE2 Mean_ROA2 Mean_LN_Boardsize Mean_Current_Ratio2;
output out = henkilo_data mean = M_D_Vaste_talous_kaikki M_P_D_talous_kaikki_2 M_P_D_muu_rikos_4_2 M_P_BUSINESS_SOCIALDEG M_P_EDULEVEL M_P_CUM_SUM_DIVORCE M_P_TOTALINCOME2 M_P_D_CUM_SUBSTANCE M_P_D_CUM_ALL_BUT_SUBS
M_D_Muu_rikos_4 M_MALE M_BUSINESS_SOCIALDEG M_EDULEVEL M_CUM_SUM_DIVORCE M_TOTALINCOME2 M_D_CUM_SUBSTANCE M_D_CUM_ALL_BUT_SUBS M_D_FMOB M_LN_NUM_POSITIONS M_Mean_OMAVARAISUUS2 M_Mean_LNTASE2 M_Mean_ROA2
M_Mean_LN_Boardsize M_Mean_Current_Ratio2;
run;

data kirjanpito_kaikki_3;
merge kirjanpito_kaikki_2 henkilo_data;

by id;

if muuttokunta = . then delete;
if skunta = 200 then delete;

LN_Talous_rikokset_mkunta_80 = log(Talous_rikokset_mkunta_80 + 1);
LN_Talous_rikokset_mkunta_90 = log(Talous_rikokset_mkunta_90 + 1);
LN_Talous_rikokset_mkunta_00 = log(Talous_rikokset_mkunta_00 + 1);
LN_Talous_rikokset_mkunta_10 = log(Talous_rikokset_mkunta_10 + 1);

LN_Talous_rikokset_skunta_80 = log(Talous_rikokset_skunta_80 + 1);
LN_Talous_rikokset_skunta_90 = log(Talous_rikokset_skunta_90 + 1);
LN_Talous_rikokset_skunta_00 = log(Talous_rikokset_skunta_00 + 1);
LN_Talous_rikokset_skunta_10 = log(Talous_rikokset_skunta_10 + 1);

if 1970 <= muuttovuosi <= 1979 then Talous_rikokset_muutos = Talous_rikokset_mkunta_80 - Talous_rikokset_skunta_80;
if 1980 <= muuttovuosi <= 1989 then Talous_rikokset_muutos = Talous_rikokset_mkunta_80 - Talous_rikokset_skunta_80;
if 1990 <= muuttovuosi <= 1999 then Talous_rikokset_muutos = Talous_rikokset_mkunta_90 - Talous_rikokset_skunta_90;
if 2000 <= muuttovuosi <= 2009 then Talous_rikokset_muutos = Talous_rikokset_mkunta_00 - Talous_rikokset_skunta_00;
if 2010 <= muuttovuosi <= 2018 then Talous_rikokset_muutos = Talous_rikokset_mkunta_10 - Talous_rikokset_skunta_10;

if 1970 <= muuttovuosi <= 1979 then Talous_rikokset_muutos_LN = LN_Talous_rikokset_mkunta_80 - LN_Talous_rikokset_skunta_80;
if 1980 <= muuttovuosi <= 1989 then Talous_rikokset_muutos_LN = LN_Talous_rikokset_mkunta_80 - LN_Talous_rikokset_skunta_80;
if 1990 <= muuttovuosi <= 1999 then Talous_rikokset_muutos_LN = LN_Talous_rikokset_mkunta_90 - LN_Talous_rikokset_skunta_90;
if 2000 <= muuttovuosi <= 2009 then Talous_rikokset_muutos_LN = LN_Talous_rikokset_mkunta_00 - LN_Talous_rikokset_skunta_00;
if 2010 <= muuttovuosi <= 2018 then Talous_rikokset_muutos_LN = LN_Talous_rikokset_mkunta_10 - LN_Talous_rikokset_skunta_10;

LN_muuttoika = log(muuttoika + 1);

D_D_Vaste_talous_kaikki = 0;
if M_D_Vaste_talous_kaikki > 0 then D_D_Vaste_talous_kaikki = 1;

D_P_D_talous_kaikki_2 = 0;
if M_P_D_talous_kaikki_2 > 0 then D_P_D_talous_kaikki_2 = 1;

D_P_D_muu_rikos_4_2 = 0;
if M_P_D_muu_rikos_4_2 > 0 then D_P_D_muu_rikos_4_2 = 1;

D_P_BUSINESS_SOCIALDEG = 0;
if M_P_BUSINESS_SOCIALDEG > 0 then D_P_BUSINESS_SOCIALDEG = 1;

M_P_LN_EDULEVEL = log(M_P_EDULEVEL + 1);

M_P_LN_CUM_SUM_DIVORCE = log(M_P_CUM_SUM_DIVORCE + 1);

M_P_LN_TOTALINCOME2 = log(M_P_TOTALINCOME2 + 1);

D_P_D_CUM_SUBSTANCE = 0;
if M_P_D_CUM_SUBSTANCE > 0 then D_P_D_CUM_SUBSTANCE = 1;

D_P_D_CUM_ALL_BUT_SUBS = 0;
if M_P_D_CUM_ALL_BUT_SUBS then D_P_D_CUM_ALL_BUT_SUBS = 1;

D_D_Muu_rikos_4 = 0;
if M_D_Muu_rikos_4 > 0 then D_D_Muu_rikos_4 = 1;

D_BUSINESS_SOCIALDEG = 0;
if M_BUSINESS_SOCIALDEG > 0 then D_BUSINESS_SOCIALDEG = 1;

M_LN_EDULEVEL = log(M_EDULEVEL + 1);

M_LN_CUM_SUM_DIVORCE = log(M_CUM_SUM_DIVORCE + 1);

M_LN_TOTALINCOME2 = log(M_TOTALINCOME2 + 1); 

D_D_CUM_SUBSTANCE = 0;
if M_D_CUM_SUBSTANCE > 0 then D_D_CUM_SUBSTANCE = 1;

D_D_CUM_ALL_BUT_SUBS = 0;
if M_D_CUM_ALL_BUT_SUBS > 0 then D_D_CUM_ALL_BUT_SUBS = 1;

D_Talous_rikokset_muutos_Q = 0;

if Talous_rikokset_muutos > 0.922 then D_Talous_rikokset_muutos_Q = 1;

T_rikokset_D_muutos_ika_inter = D_Talous_rikokset_muutos_Q * muuttoika;

run;

proc sort nodupkey data = kirjanpito_kaikki_3;
by id;
run;

proc univariate noprint data = kirjanpito_kaikki_3;
var Talous_rikokset_muutos;
output out = rikokset_muutos_muutto pctlpts = 33.3 66.7 pctlpre = P_rikokset_;
run;

proc means P1 P99 data = Kirjanpito_kaikki_3 noprint;
var Talous_rikokset_muutos Talous_rikokset_mkunta_80 Talous_rikokset_mkunta_90 Talous_rikokset_mkunta_00 Talous_rikokset_mkunta_10
Talous_rikokset_skunta_80 Talous_rikokset_skunta_90 Talous_rikokset_skunta_00 Talous_rikokset_skunta_10;
output out = Talous_rikokset_muutos_win P1 = P1_Talous_rikokset_muutos P1_Talous_rikokset_mkunta_80 P1_Talous_rikokset_mkunta_90 P1_Talous_rikokset_mkunta_00
P1_Talous_rikokset_mkunta_10 P1_Talous_rikokset_skunta_80 P1_Talous_rikokset_skunta_90 P1_Talous_rikokset_skunta_00 P1_Talous_rikokset_skunta_10
P99 = P99_Talous_rikokset_muutos P99_Talous_rikokset_mkunta_80 P99_Talous_rikokset_mkunta_90 P99_Talous_rikokset_mkunta_00
P99_Talous_rikokset_mkunta_10 P99_Talous_rikokset_skunta_80 P99_Talous_rikokset_skunta_90 P99_Talous_rikokset_skunta_00 P99_Talous_rikokset_skunta_10;
run;

data Kirjanpito_kaikki_33;

if _n_= 1 then set Talous_rikokset_muutos_win;

set Kirjanpito_kaikki_3;

Talous_rikokset_muutos2 = Talous_rikokset_muutos;

if Talous_rikokset_muutos > P99_Talous_rikokset_muutos then Talous_rikokset_muutos2=P99_Talous_rikokset_muutos;
if Talous_rikokset_muutos < P1_Talous_rikokset_muutos then Talous_rikokset_muutos2=P1_Talous_rikokset_muutos;

Talous_rikokset_muutos_ika_inter = Talous_rikokset_muutos2 * muuttoika;

Talous_rikokset_mkunta_80_2 = Talous_rikokset_mkunta_80;
Talous_rikokset_mkunta_90_2 = Talous_rikokset_mkunta_90;
Talous_rikokset_mkunta_00_2 = Talous_rikokset_mkunta_00;
Talous_rikokset_mkunta_10_2 = Talous_rikokset_mkunta_10;
Talous_rikokset_skunta_80_2 = Talous_rikokset_skunta_80;
Talous_rikokset_skunta_90_2 = Talous_rikokset_skunta_90;
Talous_rikokset_skunta_00_2 = Talous_rikokset_skunta_00;
Talous_rikokset_skunta_10_2 = Talous_rikokset_skunta_10;

if Talous_rikokset_mkunta_80 > P99_Talous_rikokset_mkunta_80 then Talous_rikokset_mkunta_80_2=P99_Talous_rikokset_mkunta_80;
if Talous_rikokset_mkunta_90 > P99_Talous_rikokset_mkunta_90 then Talous_rikokset_mkunta_90_2=P99_Talous_rikokset_mkunta_90;
if Talous_rikokset_mkunta_00 > P99_Talous_rikokset_mkunta_00 then Talous_rikokset_mkunta_00_2=P99_Talous_rikokset_mkunta_00;
if Talous_rikokset_mkunta_10 > P99_Talous_rikokset_mkunta_10 then Talous_rikokset_mkunta_10_2=P99_Talous_rikokset_mkunta_10;
if Talous_rikokset_skunta_80 > P99_Talous_rikokset_skunta_80 then Talous_rikokset_skunta_80_2=P99_Talous_rikokset_skunta_80;
if Talous_rikokset_skunta_90 > P99_Talous_rikokset_skunta_90 then Talous_rikokset_skunta_90_2=P99_Talous_rikokset_skunta_90;
if Talous_rikokset_skunta_00 > P99_Talous_rikokset_skunta_00 then Talous_rikokset_skunta_00_2=P99_Talous_rikokset_skunta_00;
if Talous_rikokset_skunta_10 > P99_Talous_rikokset_skunta_10 then Talous_rikokset_skunta_10_2=P99_Talous_rikokset_skunta_10;

T_rikokset_muutos_LN_ika_inter = Talous_rikokset_muutos_LN * muuttoika;

LN_Talous_rikokset_mkunta_80_2 = log(Talous_rikokset_mkunta_80_2 + 1);
LN_Talous_rikokset_mkunta_90_2 = log(Talous_rikokset_mkunta_90_2 + 1);
LN_Talous_rikokset_mkunta_00_2 = log(Talous_rikokset_mkunta_00_2 + 1);
LN_Talous_rikokset_mkunta_10_2 = log(Talous_rikokset_mkunta_10_2 + 1);

LN_Talous_rikokset_skunta_80_2 = log(Talous_rikokset_skunta_80_2 + 1);
LN_Talous_rikokset_skunta_90_2 = log(Talous_rikokset_skunta_90_2 + 1);
LN_Talous_rikokset_skunta_00_2 = log(Talous_rikokset_skunta_00_2 + 1);
LN_Talous_rikokset_skunta_10_2 = log(Talous_rikokset_skunta_10_2 + 1);

if 1970 <= muuttovuosi <= 1979 then Talous_rikokset_muutos_win = Talous_rikokset_mkunta_80_2 - Talous_rikokset_skunta_80_2;
if 1980 <= muuttovuosi <= 1989 then Talous_rikokset_muutos_win = Talous_rikokset_mkunta_80_2 - Talous_rikokset_skunta_80_2;
if 1990 <= muuttovuosi <= 1999 then Talous_rikokset_muutos_win = Talous_rikokset_mkunta_90_2 - Talous_rikokset_skunta_90_2;
if 2000 <= muuttovuosi <= 2009 then Talous_rikokset_muutos_win = Talous_rikokset_mkunta_00_2 - Talous_rikokset_skunta_00_2;
if 2010 <= muuttovuosi <= 2018 then Talous_rikokset_muutos_win = Talous_rikokset_mkunta_10_2 - Talous_rikokset_skunta_10_2;

if 1970 <= muuttovuosi <= 1979 then Talous_rikokset_muutos_LN_win = LN_Talous_rikokset_mkunta_80_2 - LN_Talous_rikokset_skunta_80_2;
if 1980 <= muuttovuosi <= 1989 then Talous_rikokset_muutos_LN_win = LN_Talous_rikokset_mkunta_80_2 - LN_Talous_rikokset_skunta_80_2;
if 1990 <= muuttovuosi <= 1999 then Talous_rikokset_muutos_LN_win = LN_Talous_rikokset_mkunta_90_2 - LN_Talous_rikokset_skunta_90_2;
if 2000 <= muuttovuosi <= 2009 then Talous_rikokset_muutos_LN_win = LN_Talous_rikokset_mkunta_00_2 - LN_Talous_rikokset_skunta_00_2;
if 2010 <= muuttovuosi <= 2018 then Talous_rikokset_muutos_LN_win = LN_Talous_rikokset_mkunta_10_2 - LN_Talous_rikokset_skunta_10_2;

T_rikokset_muutos_LN_ika_int_win = Talous_rikokset_muutos_LN_win * muuttoika;
T_rikokset_muutos_LN_ika_inter = Talous_rikokset_muutos_LN * muuttoika;
T_rikokset_muutos_ika_int_win = Talous_rikokset_muutos_win * muuttoika;

run;

proc sort nodupkey data = kirjanpito_kaikki_33;
by id;
run;*/

* Save final datasets to be used in analyses in STATA;

data lib100.stata_talous_011025 /*lib100.Statadata_fmob_firm_0_011025 lib100.Statadata_owner_011025 lib100.stata_t_muuttaneet_1k_011025*/;
 
set Kirjanpito_kaikki_2 /*Kirjanpito_kaikki_33*/;

if Pa_MALE ^= . then Puoliso = 1;

keep ID YEAR FIRM D_Vaste_kirjanpitorikos_tuomio D_Vaste_talous_tuomio D_Vaste_talous_kaikki_tuomio D_Vaste_muu_rikos_4_tuomio D_Vaste_muu_rikos_5_tuomio 
D_Vaste_kirjanpitorikos_kaikki D_Vaste_talous D_Vaste_talous_kaikki D_Vaste_talous_kaikki_vankeus D_Vaste_talous_kaikki_3kat D_Vaste_muu_rikos_4 D_Vaste_muu_rikos_5
D_Vaste_Liikenne_tuomio D_Vaste_Finan_rep_fraud D_Vaste_Tax_relat_fraud D_Vaste_Other_finan_crime D_Vaste_Finan_rep_fraud_tuomio D_Vaste_Tax_relat_fraud_tuomio D_Vaste_Other_finan_crime_tuomio
MALE LN_AGE BUSINESS_SOCIALDEG LN_Edulevel D_Marital_status LN_Cum_Sum_Divorce D_FMOB D_Vaste_talous_kaikki_yks D_Vaste_Finan_rep_fraud_yks D_Vaste_Tax_relat_fraud_yks D_Vaste_Other_finan_crime_yks
LN_TOTALINCOME2 TOTALINCOME2 LN_TOTALINCOME4 LN_SALARYTOTAL2 SALARYTOTAL2 D_talous_kaikki D_talous_kaikki_tuomio D_muu_rikos_4 D_muu_rikos_4_tuomio D_muu_rikos_5_tuomio D_Liikenne_tuomio LN_NUM_POSITIONS NUM_POSITIONS MEAN_PERCENTAGEOWN CEO_PERCENTAGEOWN
LN_CUM_Muu_rikos_4_tuomio_2 LN_CUM_Talous_kaikki_tuomio_2 LN_CUM_Muu_rikos_4_2 LN_CUM_Talous_kaikki_2 SALARYTOTAL LN_SALARYTOTAL TOTALINCOME LN_TOTALINCOME
LN_P_CUM_Muu4_tuomio LN_P_CUM_talous_kaikki_tuomio LN_P_CUM_Muu4 LN_P_CUM_talous_kaikki
P_BUSINESS_SOCIALDEG P_LN_EDULEVEL P_D_MARITAL_STATUS P_LN_CUM_SUM_DIVORCE P_D_Kaikki_rikos
P_biol_D_talous_kaikki_tuomio_2 P_sos_D_talous_kaikki_tuomio_2 P_biol_D_muu_rikos_4_tuomio_2 P_sos_D_muu_rikos_4_tuomio_2
P_D_talous_kaikki_tuomio_2 P_D_muu_rikos_4_tuomio_2 P_D_talous_kaikki_2 P_D_muu_rikos_4_2 P_D_muu_rikos_5_tuomio_2 P_D_Liikenne_tuomio_2 P_LN_TOTALINCOME2 P_LN_TOTALINCOME4
M_D_Talous_kaikki_2 M_D_muu_rikos_4_2 F_D_Talous_kaikki_2 F_D_muu_rikos_4_2 P_biol_D_talous_kaikki_2 P_sos_D_talous_kaikki_2 P_biol_D_muu_rikos_4_2 P_sos_D_muu_rikos_4_2
Mean_Mean_D_talous_kaikki Mean_Mean_D_muu_rikos_4 D_Mean_D_Mean_D_talous_kaikki D_Mean_D_Mean_D_muu_rikos_4 Mean_D_Talous_kaikki Mean_D_Muu_rikos_4
Mean_OMAVARAISUUS2 Mean_LNTASE2 Mean_ROA2 Mean_LN_Boardsize Mean_Current_Ratio2 MEAN_SUM_PERCENTAGEOWN2 MAAKUNTA MAAKUNTA_FIRST KUNTA KUNTA_FIRST OMAVARAISUUS2 LNTASE2 LNTASE4 ROA2 LN_Boardsize Current_Ratio2
Industry_2002_3 Industry_2002_2 YearIndustryProvince YearIndustryMunicipality YearIndustry YearMunicipality YearProvince
FMOB_0_P_D_talous_kaikki_2 FMOB_0_P_D_muu_rikos_4_2 FMOB_1_P_D_talous_kaikki_2 FMOB_1_P_D_muu_rikos_4_2 FMOB_0_P_D_talous_kaikki_2_2 FMOB_0_P_D_muu_rikos_4_2_2 FMOB_1_P_D_talous_kaikki_2_2 FMOB_1_P_D_muu_rikos_4_2_2 PEP_0_P_D_talous_kaikki_2 PEP_0_P_D_muu_rikos_4_2 PEP_1_P_D_talous_kaikki_2 PEP_1_P_D_muu_rikos_4_2
VR_a16v_P_D_talous_kaikki_2 VR_a16v_P_D_muu_rikos_4_2 VR_y16v_P_D_talous_kaikki_2 VR_y16v_P_D_muu_rikos_4_2 Vanh_rikos_j_16v_ei_sita_ennen Vanhempien_rikos_ennen_16v VR_5_15v_P_D_talous_kaikki_2 VR_5_15v_P_D_muu_rikos_4_2
VR_a5v_P_D_talous_kaikki_2 VR_a5v_P_D_muu_rikos_4_2 Vuosia_vanh_viim_talousrikos_2 Vuosia_vanh_viim_muu_rikos_2 FS_A_Median
FS_A_median_P_D_talous_kaikki_2 FS_B_median_P_D_talous_kaikki_2 FS_A_median_P_D_muu_rikos_4_2 FS_B_median_P_D_muu_rikos_4_2 Big_Four_0_P_D_talous_kaikki_2 Big_Four_0_P_D_muu_rikos_4_2 Big_Four_1_P_D_talous_kaikki_2 Big_Four_1_P_D_muu_rikos_4_2 D_Big_four
P_D_talous_kaikki_vankeus_2 P_D_muu_rikos_4_vankeus_2 P_D_talous_kaikki_sakko_2 P_D_muu_rikos_4_sakko_2 F_D_talous_kaikki_3 F_D_muu_rikos_4_3 M_D_talous_kaikki_3 M_D_muu_rikos_4_3 P_D_talous_kaikki_3 P_D_muu_rikos_4_3
P_D_Talous_kaikki_2_sama_mk P_D_Talous_kaikki_2_eri_mk P_D_muu_rikos_4_2_sama_mk P_D_muu_rikos_4_2_eri_mk P_D_Talous_kaikki_2_sama_kunta P_D_Talous_kaikki_2_eri_kunta P_D_muu_rikos_4_2_sama_kunta P_D_muu_rikos_4_2_eri_kunta
P_D_talous_kaikki_4 P_D_muu_rikos_4_4 P_D_talous_kaik_ehdol_vankeus_2 P_D_talous_kaik_ehdot_vankeus_2 P_D_Talous_kaik_ehdol_vankeus_3
P_D_muu_rikos_4_ehdol_vankeus_3 P_D_muu_rikos_4_ehdol_vankeus_2 P_D_muu_rikos_4_ehdot_vankeus_2 P_D_talous_kaikki_3kat P_D_muu_rikos_4_3kat
Y_MED_OWN_0_P_D_talous_kaik_2 Y_MED_OWN_1_P_D_talous_kaik_2 Y_MED_OWN_0_P_D_muu_rikos_4_2 Y_MED_OWN_1_P_D_muu_rikos_4_2  
Y_M_SUMOWN_0_P_D_talous_kaik_2 Y_M_SUMOWN_1_P_D_talous_kaik_2 Y_M_SUMOWN_0_P_D_muu_rikos_4_2 Y_M_SUMOWN_1_P_D_muu_rikos_4_2
EQ_Y_Median_0_P_D_talous_kaik_2 EQ_Y_Median_1_P_D_talous_kaik_2 EQ_Y_Median_0_P_D_muu_rikos_4_2 EQ_Y_Median_1_P_D_muu_rikos_4_2
BS_Y_Median_0_P_D_talous_kaik_2 BS_Y_Median_1_P_D_talous_kaik_2 BS_Y_Median_0_P_D_muu_rikos_4_2 BS_Y_Median_1_P_D_muu_rikos_4_2
Yli_Med_MEAN_PERCENTAGEOWN Yli_Med_MEAN_SUM_PERCENTAGEOWN2 SUMOWN_yli_0_0_P_D_talous_kaik_2 SUMOWN_yli_0_0_P_D_muu_rikos_4_2 SUMOWN_yli_0_1_P_D_talous_kaik_2 SUMOWN_yli_0_1_P_D_muu_rikos_4_2 MEAN_SUM_PERCENTAGEOWN_yli_0 familymember_on_board
P_D_CUM_SUBSTANCE P_D_CUM_ALL_BUT_SUBS D_CUM_SUBSTANCE D_CUM_ALL_BUT_SUBS Assets_3lk Assets_3lk_2 Vaste_sakot_talous_kaikki Vaste_vankeus_talous_kaikki
Vaste_sakot_finan_rep_fraud Vaste_vankeus_finan_rep_fraud Vaste_sakot_tax_relat_fraud Vaste_vankeus_tax_relat_fraud Vaste_sakot_other_finan Vaste_vankeus_other_finan
Vaste_LN_sakot_talous_kaikki Vaste_LN_vankeus_talous_kaikki Vaste_LN_sakot_finan_rep_fraud Vaste_LN_vankeus_finan_rep_fraud
Vaste_LN_sakot_tax_relat_fraud Vaste_LN_vankeus_tax_relat_fraud Vaste_LN_sakot_other_finan Vaste_LN_vankeus_other_finan
Alle_18v_rikos Alle_18v_P_D_talous_kaikki_2 Alle_18v_P_D_muu_rikos_4_2 Yli_18v_P_D_talous_kaikki_2 Yli_18v_P_D_muu_rikos_4_2 D_muu_rikos_4_yli_18 m_id f_id
P_D_Talous_kaikki_vankeus_tuom_2 P_D_Talous_kaikki_sakko_tuomio_2 P_D_muu_rikos_4_vankeus_tuomio_2 P_D_muu_rikos_4_sakko_tuomio_2 LN_Vaste_talous_kaikki_num LN_Vaste_Finan_rep_fraud_num LN_Vaste_Tax_relat_fraud_num LN_Vaste_Other_finan_crime_num
Mean_CEO_Duality pa_id Pa_D_Muu_rikos_4_2 Pa_D_Talous_kaikki_2 Pa_D_muu_rikos_4_vankeus_2 Pa_D_Talous_kaikki_vankeus_2
Pa_D_CUM_SUBSTANCE Pa_D_CUM_ALL_BUT_SUBS Pa_CUM_NALL_BUT_SUBS Pa_CUM_NSUBSTANCE Pa_LN_TOTALINCOME2 Pa_TOTALINCOME2 Pa_EDULEVEL Pa_LN_EDULEVEL Pa_BUSINESS_SOCIALDEG Pa_LN_CUM_Sum_divorce LN_Pa_Length Pa_MALE Puoliso Pa_Length_yli_median
Kaikki_rikokset_17_3lk LN_Kaikki_rikokset_vuosi LN_Talous_rikokset_vuosi LN_Muut_rikokset_vuosi LN_Kaikki_rikokset_17 Kaikki_rikokset_vuosi Kaikki_rikokset_17 Kaikki_rikokset_lapsuus Talous_rikokset_lapsuus Muut_rikokset_lapsuus 
LN_Kaikki_rikokset_lapsuus LN_Talous_rikokset_lapsuus LN_Muut_rikokset_lapsuus Lapsuusvuosia LN_Kaikki_rikokset_m_lapsuus LN_Talous_rikokset_m_lapsuus LN_Muut_rikokset_m_lapsuus LN_Kaikki_rikokset_m_vuosi LN_Talous_rikokset_m_vuosi LN_Muut_rikokset_m_vuosi
D_D_Vaste_talous_kaikki D_P_D_talous_kaikki_2 D_P_D_muu_rikos_4_2 D_P_BUSINESS_SOCIALDEG M_P_LN_EDULEVEL M_P_LN_CUM_SUM_DIVORCE M_P_LN_TOTALINCOME2 D_P_D_CUM_SUBSTANCE D_P_D_CUM_ALL_BUT_SUBS
D_D_Muu_rikos_4 M_MALE D_BUSINESS_SOCIALDEG M_LN_EDULEVEL M_LN_CUM_SUM_DIVORCE M_LN_TOTALINCOME2 D_D_CUM_SUBSTANCE D_D_CUM_ALL_BUT_SUBS Talous_rikokset_muutos Talous_rikokset_muutos_ika_inter Talous_rikokset_muutos_ika_i_24 muuttoika muuttoik_24 skunta Talous_rikokset_lapsuus_Q
P_Pa_LN_EDULEVEL P_Pa_LN_CUM_Sum_divorce P_Pa_LN_TOTALINCOME2 P_Pa_D_muu_rikos_4_2 P_Pa_D_talous_kaikki_2 P_Pa_D_CUM_SUBSTANCE P_Pa_D_CUM_ALL_BUT_SUBS P_Pa_BUSINESS_SOCIALDEG P_Pa_D_MARITAL_STATUS
Pa_LN_Kaikki_rikokset_lapsuus Pa_LN_Talous_rikokset_lapsuus Pa_LN_Muut_rikokset_lapsuus Pa_LN_AGE Puolison_vanhemmat Pa_Length_5 M_D_FMOB M_LN_NUM_POSITIONS M_Mean_OMAVARAISUUS2 M_Mean_LNTASE2 M_Mean_ROA2
M_Mean_LN_Boardsize M_Mean_Current_Ratio2 Finan_filter Tax_filter Other_filter Talous_rikokset_lapsuus_24_Q LN_Kaikki_rikokset_lapsuus_24 LN_Talous_rikokset_lapsuus_24 LN_Muut_rikokset_lapsuus_24 Talous_rikokset_lapsuus_High Talous_rikokset_muutos2
Testimuuttuja D_Talous_rikokset_muutos_Q T_rikokset_D_muutos_ika_inter LN_muuttoika LN_Talous_rikokset_muutos LN_T_rikokset_muutos_ika_int LN_T_rikokset_muutos_LN_ika_int T_rikokset_muutos_LN_ika_inter Talous_rikokset_muutos_LN
T_rikokset_muutos_LN_ika_int_win Talous_rikokset_muutos_LN_win T_rikokset_muutos_ika_int_win Talous_rikokset_muutos_win LN_Talous_rikokset_lapsuus2 Talous_rikokset_lapsuus2 Finan_mukaan Tax_mukaan Other_mukaan Tax_mukaan_2 Other_mukaan_2
Finan_mukaan_yks Tax_mukaan_yks Other_mukaan_yks vah_toinen_vanh_mm molemmat_vanh_mm sama_tuomioistuin D_sama_tuomioistuin BIRTHYEAR P_D_Lifetime_talous_kaikki P_D_Lifetime_talous_kaikki P_D_Lifetime_muu_rikos_4 Eri_kunta Score_muuttuja
LN_Score_muuttuja Score_muuttuja_3lk Score_vuokra_17v Score_isan_sose_asema_17v Score_vanh_tulot_17v Score_vanh_maahanmuuttaja Score_vanh_kieli Score_vanh_koulutus Score_muuttuja_2_3lk Score_muuttuja_2 D_mean_mm_vanh
m_sos_vanh_dummy f_sos_vanh_dummy p_sos_vanh_dummy;

run;

*  Statistics in Tables 1 and 2;

proc means noprint data = kirjanpito_kaikki_2;
by id;
var AGE TotalIncome2 NUM_Positions P_TOTALINCOME2 D_Vaste_talous_kaikki D_Vaste_Finan_rep_fraud D_Vaste_tax_relat_fraud D_Vaste_Other_finan_crime
D_Vaste_muu_rikos_4 Sakot_talous_kaikki Vankeus_talous_kaikki;
output out = Means_Table1 mean = Mean_AGE Mean_TotalIncome2 Mean_NUM_Positions Mean_P_TOTALINCOME2 Mean_D_Vaste_talous_kaikki Mean_D_Vaste_Finan_rep_fraud
Mean_D_Vaste_tax_relat_fraud Mean_D_Vaste_Other_finan_crime Mean_D_Vaste_muu_rikos_4 Mean_Sakot_talous_kaikki Mean_Vankeus_talous_kaikki;
run;

proc means noprint data = kirjanpito_kaikki_2;
by id;
var Talous_kaikki Finan_rep_fraud Tax_relat_fraud Other_finan_crime muu_rikos_4 Talous_kaikki_tuomio Finan_rep_fraud_tuomio Tax_relat_fraud_tuomio Other_finan_crime_tuomio muu_rikos_4_tuomio;
output out = Sums_Table1 sum = Sum_Talous_kaikki Sum_Finan_rep_fraud Sum_Tax_relat_fraud Sum_Other_finan_crime Sum_muu_rikos_4 Sum_Talous_kaikki_tuomio Sum_Finan_rep_fraud_tuomio Sum_Tax_relat_fraud_tuomio Sum_Other_finan_crime_tuomio Sum_muu_rikos_4_tuomio;
run;

data Table1;
merge kirjanpito_kaikki_2 Means_Table1 Sums_Table1;
by id;

D_D_Vaste_talous_kaikki = 0;
D_D_Vaste_Finan_rep_fraud = 0;
D_D_Vaste_tax_relat_fraud = 0;
D_D_Vaste_Other_finan_crime = 0;
D_D_Vaste_muu_rikos_4 = 0;

D_Sakot_talous_kaikki = 0;
D_Vankeus_talous_kaikki = 0;

if Mean_D_Vaste_talous_kaikki > 0 then D_D_Vaste_talous_kaikki = 1;
if Mean_D_Vaste_Finan_rep_fraud > 0 then D_D_Vaste_Finan_rep_fraud = 1;
if Mean_D_Vaste_tax_relat_fraud > 0 then D_D_Vaste_tax_relat_fraud = 1;
if Mean_D_Vaste_Other_finan_crime > 0 then D_D_Vaste_Other_finan_crime = 1;
if Mean_D_Vaste_muu_rikos_4 > 0 then D_D_Vaste_muu_rikos_4 = 1;

if Mean_Sakot_talous_kaikki > 0 then D_Sakot_talous_kaikki = 1;
if Mean_Vankeus_talous_kaikki > 0 then D_Vankeus_talous_kaikki = 1; 

if Sum_Talous_kaikki = . then Sum_Talous_kaikki = 0;
if Sum_Finan_rep_fraud = . then Sum_Finan_rep_fraud = 0;
if Sum_Tax_relat_fraud = . then Sum_Tax_relat_fraud = 0;
if Sum_Other_finan_crime = . then Sum_Other_finan_crime = 0;
if Sum_muu_rikos_4 = . then Sum_muu_rikos_4 = 0;

D_Sum_talous_kaikki = 0;
if Sum_Talous_kaikki > 0 then D_Sum_Talous_kaikki = 1;

D_Sum_talous_kaikki_tuomio = 0;
if Sum_Talous_kaikki_tuomio > 0 then D_Sum_Talous_kaikki_tuomio = 1;

D_Sum_muu_rikos_4 = 0;
if Sum_muu_rikos_4 > 0 then D_Sum_muu_rikos_4 = 1;

if Sum_talous_kaikki = 1 then D3_sum_talous_kaikki = 1; if Sum_talous_kaikki = 2 then D3_sum_talous_kaikki = 2; if Sum_talous_kaikki >= 3 then D3_sum_talous_kaikki = 3;
if Lifetime_CUM_Talous_kaikki = 1 then D3_Lifetime_CUM_Talous_kaikki = 1; if Lifetime_CUM_Talous_kaikki = 2 then D3_Lifetime_CUM_Talous_kaikki = 2; if Lifetime_CUM_Talous_kaikki >= 3 then D3_Lifetime_CUM_Talous_kaikki = 3;
if F_Lifetime_CUM_Talous_kaikki = 1 then D3_F_Lifetime_CUM_Talous_kaikki = 1; if F_Lifetime_CUM_Talous_kaikki = 2 then D3_F_Lifetime_CUM_Talous_kaikki = 2; if F_Lifetime_CUM_Talous_kaikki >= 3 then D3_F_Lifetime_CUM_Talous_kaikki = 3;
if M_Lifetime_CUM_Talous_kaikki = 1 then D3_M_Lifetime_CUM_Talous_kaikki = 1; if M_Lifetime_CUM_Talous_kaikki = 2 then D3_M_Lifetime_CUM_Talous_kaikki = 2; if M_Lifetime_CUM_Talous_kaikki >= 3 then D3_M_Lifetime_CUM_Talous_kaikki = 3;

if Sum_muu_rikos_4 = 1 then D3_sum_muu_rikos_4 = 1; if Sum_muu_rikos_4 = 2 then D3_sum_muu_rikos_4 = 2; if Sum_muu_rikos_4 >= 3 then D3_sum_muu_rikos_4 = 3;
if Lifetime_CUM_muu4 = 1 then D3_Lifetime_CUM_muu4 = 1; if Lifetime_CUM_muu4 = 2 then D3_Lifetime_CUM_muu4 = 2; if Lifetime_CUM_muu4 >= 3 then D3_Lifetime_CUM_muu4 = 3;
if F_Lifetime_CUM_muu4 = 1 then D3_F_Lifetime_CUM_muu4 = 1; if F_Lifetime_CUM_muu4 = 2 then D3_F_Lifetime_CUM_muu4 = 2; if F_Lifetime_CUM_muu4 >= 3 then D3_F_Lifetime_CUM_muu4 = 3;
if M_Lifetime_CUM_muu4 = 1 then D3_M_Lifetime_CUM_muu4 = 1; if M_Lifetime_CUM_muu4 = 2 then D3_M_Lifetime_CUM_muu4 = 2; if M_Lifetime_CUM_muu4 >= 3 then D3_M_Lifetime_CUM_muu4 = 3;

run;

proc sort data = Table1;
by id descending year;
run;

* Data for mothers and fathers;

data mothers;
set Table1;
if M_MALE = . then delete;
run;

proc sort nodupkey data = mothers;
by m_id;
run;

data fathers;
set Table1;
if F_MALE = . then delete;
run;

proc sort nodupkey data = fathers;
by f_id;
run;

proc sort nodupkey data = Table1;
by id;
run;

* Table 1, panels A and B;

proc freq data = Table1;
tables D_Sum_Talous_kaikki*D_Sum_Talous_kaikki_tuomio D_Lifetime_Talous_kaikki*D_Lifetime_Talous_kaikki_tuomio D3_Sum_Talous_kaikki D3_Lifetime_CUM_Talous_kaikki;
run;

proc freq data = mothers;
tables M_D_lifetime_Talous_kaikki*M_D_lifetime_Talous_kaik_tuomio D3_M_Lifetime_CUM_Talous_kaikki;
run;

proc freq data = fathers;
tables F_D_lifetime_Talous_kaikki*F_D_lifetime_Talous_kaik_tuomio D3_F_Lifetime_CUM_Talous_kaikki;
run;

* Table 1, panel C;

proc freq data = Table1;
tables D_D_Vaste_Finan_rep_fraud D_D_Vaste_tax_relat_fraud D_D_Vaste_Other_finan_crime;
run;

proc freq data = Table1;
tables D_Lifetime_Finan_fraud D_Lifetime_Tax_fraud D_Lifetime_Other_fraud;
run;

proc freq data = mothers;
tables M_D_Lifetime_Finan_fraud M_D_Lifetime_Tax_fraud M_D_Lifetime_Other_fraud;
run;

proc freq data = fathers;
tables F_D_Lifetime_Finan_fraud F_D_Lifetime_Tax_fraud F_D_Lifetime_Other_fraud;
run;

* Table 1, panel D; 

proc freq data = Table1; 
tables D_Sakot_talous_kaikki D_Lifetime_Talous_kaikki_s D_Vankeus_talous_kaikki D_Lifetime_Talous_kaikki_v;
run;

proc freq data = mothers;
tables M_D_Lifetime_Talous_kaikki_s M_D_Lifetime_Talous_kaikki_v;
run;

proc freq data = fathers;
tables F_D_Lifetime_Talous_kaikki_s F_D_Lifetime_Talous_kaikki_v;
run;

* Statistics for Table 2;

proc means n mean median std noprint data = kirjanpito_kaikki_2;
by id;
var D_Vaste_talous_kaikki EDULEVEL BUSINESS_SOCIALDEG D_MARITAL_STATUS CUM_SUM_Divorce TOTALINCOME2 D_CUM_SUBSTANCE D_CUM_ALL_BUT_SUBS D_muu_rikos_4 MALE AGE NUM_POSITIONS D_FMOB
P_D_talous_kaikki_2 P_D_muu_rikos_4_2 P_EDULEVEL P_BUSINESS_SOCIALDEG P_D_MARITAL_STATUS P_CUM_Sum_Divorce P_TOTALINCOME2 P_D_CUM_SUBSTANCE P_D_CUM_ALL_BUT_SUBS;
output out = Table2 mean = Mean_D_Vaste_talous_kaikki Mean_EDULEVEL Mean_BUSINESS_SOCIALDEG Mean_D_MARITAL_STATUS Mean_CUM_SUM_Divorce Mean_TOTALINCOME2 Mean_D_CUM_SUBSTANCE
Mean_D_CUM_ALL_BUT_SUBS Mean_D_muu_rikos_4 Mean_MALE Mean_AGE Mean_NUM_POSITIONS Mean_D_FMOB Mean_P_D_talous_kaikki_2 Mean_P_D_muu_rikos_4_2 Mean_P_EDULEVEL Mean_P_BUSINESS_SOCIALDEG
Mean_P_D_MARITAL_STATUS Mean_P_CUM_Sum_Divorce Mean_P_TOTALINCOME2 Mean_P_D_CUM_SUBSTANCE Mean_P_D_CUM_ALL_BUT_SUBS;
run;

data Table2_2;
merge kirjanpito_kaikki_6 Table2;
by id;

D_D_Vaste_talous_kaikki = 0;

if Mean_D_Vaste_talous_kaikki > 0 then D_D_Vaste_talous_kaikki = 1;

run;

proc sort nodupkey data = Table2_2;
by id;
run;

* Table 2, panel A;

proc means mean std p25 median p75 n data = Table2_2;
var Mean_D_Vaste_talous_kaikki Mean_EDULEVEL Mean_BUSINESS_SOCIALDEG Mean_D_MARITAL_STATUS Mean_CUM_SUM_Divorce Mean_TOTALINCOME2 Mean_D_CUM_SUBSTANCE
Mean_D_CUM_ALL_BUT_SUBS Mean_D_muu_rikos_4 Mean_MALE Mean_AGE Mean_NUM_POSITIONS Mean_D_FMOB Mean_P_D_talous_kaikki_2 Mean_P_D_muu_rikos_4_2 Mean_P_EDULEVEL Mean_P_BUSINESS_SOCIALDEG
Mean_P_D_MARITAL_STATUS Mean_P_CUM_Sum_Divorce Mean_P_TOTALINCOME2 Mean_P_D_CUM_SUBSTANCE Mean_P_D_CUM_ALL_BUT_SUBS;
run;

proc sort data = table2_2;
by D_D_Vaste_talous_kaikki;
run; 

* Table 2, panel B;

proc means mean std p25 median p75 n data = Table2_2;
by D_D_Vaste_talous_kaikki;
var Mean_EDULEVEL Mean_BUSINESS_SOCIALDEG Mean_D_MARITAL_STATUS Mean_CUM_SUM_Divorce Mean_TOTALINCOME2 Mean_D_CUM_SUBSTANCE
Mean_D_CUM_ALL_BUT_SUBS Mean_D_muu_rikos_4 Mean_MALE Mean_AGE Mean_NUM_POSITIONS Mean_D_FMOB Mean_P_D_talous_kaikki_2 Mean_P_D_muu_rikos_4_2 Mean_P_EDULEVEL Mean_P_BUSINESS_SOCIALDEG
Mean_P_D_MARITAL_STATUS Mean_P_CUM_Sum_Divorce Mean_P_TOTALINCOME2 Mean_P_D_CUM_SUBSTANCE Mean_P_D_CUM_ALL_BUT_SUBS;
run;

proc ttest plots=none data = Table2_2;

class D_D_Vaste_talous_kaikki;

var Mean_EDULEVEL Mean_BUSINESS_SOCIALDEG Mean_D_MARITAL_STATUS Mean_CUM_SUM_Divorce Mean_TOTALINCOME2 Mean_D_CUM_SUBSTANCE
Mean_D_CUM_ALL_BUT_SUBS Mean_D_muu_rikos_4 Mean_MALE Mean_AGE Mean_NUM_POSITIONS Mean_D_FMOB Mean_P_D_talous_kaikki_2 Mean_P_D_muu_rikos_4_2 Mean_P_EDULEVEL Mean_P_BUSINESS_SOCIALDEG
Mean_P_D_MARITAL_STATUS Mean_P_CUM_Sum_Divorce Mean_P_TOTALINCOME2 Mean_P_D_CUM_SUBSTANCE Mean_P_D_CUM_ALL_BUT_SUBS;

run;

* Table 2, panel C;

proc means mean std p25 median p75 n data = kirjanpito_kaikki_2;

var D_Vaste_talous_kaikki D_Vaste_Finan_rep_fraud D_Vaste_Tax_relat_fraud D_Vaste_Other_finan_crime

P_D_talous_kaikki_2 P_D_muu_rikos_4_2 P_EDULEVEL P_BUSINESS_SOCIALDEG P_D_MARITAL_STATUS P_CUM_Sum_Divorce P_TOTALINCOME2 P_D_CUM_SUBSTANCE P_D_CUM_ALL_BUT_SUBS

EDULEVEL BUSINESS_SOCIALDEG D_MARITAL_STATUS CUM_SUM_Divorce TOTALINCOME2 D_CUM_SUBSTANCE D_CUM_ALL_BUT_SUBS D_muu_rikos_4 MALE AGE NUM_POSITIONS D_FMOB

Mean_OMAVARAISUUS2 Mean_TASE2 Mean_ROA2 Mean_Boardsize Mean_Current_ratio2;

run;

* Table 2, panel D tests;

proc ttest plots=none data = kirjanpito_kaikki_2;

class P_D_talous_kaikki_2;

var P_D_muu_rikos_4_2 P_BUSINESS_SOCIALDEG P_LN_EDULEVEL P_D_MARITAL_STATUS P_LN_CUM_SUM_DIVORCE P_LN_TOTALINCOME2 P_D_CUM_SUBSTANCE P_D_CUM_ALL_BUT_SUBS D_Muu_rikos_4
MALE LN_AGE BUSINESS_SOCIALDEG LN_EDULEVEL D_MARITAL_STATUS LN_CUM_SUM_DIVORCE LN_TOTALINCOME2 D_CUM_SUBSTANCE D_CUM_ALL_BUT_SUBS D_FMOB LN_NUM_POSITIONS Mean_OMAVARAISUUS2
Mean_LNTASE2 Mean_ROA2 Mean_LN_Boardsize Mean_Current_Ratio2;

run;

* Statistics for Table 7, panel A;

proc means n mean median std noprint data = kirjanpito_kaikki_2;
by id;
where Pa_MALE ^= .;
var Pa_D_talous_kaikki_2;
output out = table7a mean = mean_Pa_D_talous_kaikki_2;
run;

data table7a_2;
set table7a;

puoliso_talousrikos = 0;

if mean_Pa_D_talous_kaikki_2 > 0 then puoliso_talousrikos = 1;

run;

data table7a_3;
set kirjanpito_kaikki_2;

if Pa_MALE ^= .;

D_Talous_rikos_henk = 0;
if Talous_rikos_henk > 0 then D_Talous_rikos_henk = 1;

keep id D_Talous_rikos_henk;

run;

proc sort nodupkey data = table7a_3;
by id;
run;

data table7a_4;
merge table7a_3 table7a_2;
by id;
run;

proc freq data = table7a_4;
tables puoliso_talousrikos * D_Talous_rikos_henk / chisq;
run;
