ltc-code 0.2.0__tar.gz → 0.2.1__tar.gz

{ltc_code-0.2.0 → ltc_code-0.2.1}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.3
 Name: ltc-code
-Version: 0.2.0
+Version: 0.2.1
 Summary: Add your description here
 Requires-Python: >=3.9
 Description-Content-Type: text/markdown

{ltc_code-0.2.0 → ltc_code-0.2.1}/pyproject.toml

@@ -1,6 +1,6 @@
 [project]
 name = "ltc-code"
-version = "0.2.0"
+version = "0.2.1"
 description = "Add your description here"
 readme = "README.md"
 requires-python = ">=3.9"

ltc_code-0.2.1/src/ltc_code/20260706_ceprscripts/BALANCE.do

@@ -0,0 +1,123 @@
+use "P:\Long_Term_Charter\data\clean\new_build\cmo_analysis\cmo_analysis.dta", clear
+cap postclose handle
+clear matrix
+scalar drop _all
+
+* Generate numeric cohort
+destring cohort_lottery, gen(cohort_numeric)
+
+* Create tempfile to store results
+tempfile results
+postfile handle str50 sample str30 depvar str30 instrument double coef double se double N_obs double N_ppl double control_mean double N_ppl_control using `results', replace
+
+
+* Store local for depvars and indep vars
+local depvars p_his
+local instruments offer initial_offer
+
+* Overall and CMO specific samples
+levelsof cmo_name, local(cmos)
+
+local cmo_census
+foreach c of local cmos {
+    local cmo_census `"`cmo_census' "`c' (1990-2005)""'
+}
+
+local samples `" "No Harmony" "Overall" "Census Sample" "KIPP Sample" `cmos' `cmo_census' "'
+* local samples `" "No Harmony" "Overall" "Census Sample"  "KIPP Sample" "'
+*local samples `" `cmo_census' "'
+
+
+foreach s of local samples {
+    foreach dep of local depvars {
+        foreach z of local instruments {
+            preserve
+
+            di "SAMPLE: `s'"
+            di "INSTRUMENT: `z'"
+            di "DEPVAR: `dep'"
+
+            * Filter to relevant sample
+            if "`s'" == "No Harmony" {
+                drop if cmo_name == "Harmony"
+            }
+			else if "`s'" == "Census Sample" {
+				keep if inrange(cohort_numeric, 1990, 2005)
+				drop if cmo_name == "Harmony"
+            }
+			else if "`s'" == "KIPP Sample" {
+				keep if inrange(cohort_numeric, 1990, 2003)
+				keep if inlist(cmo_name, "KIPP_NJ", "KIPP TX", "KIPP NC", "KIPP MA")
+            }
+            else if "`s'" == "Overall" {
+
+            }
+			else if strpos("`s'", " (1990-2005)") {
+				keep if cmo_name == subinstr("`s'", " (1990-2005)", "", .)
+				keep if inrange(cohort_numeric, 1990, 2005)
+            }
+            else {
+                keep if cmo_name == "`s'"
+            }
+          
+
+            * Require all non missing regression vars
+            drop if mi(`dep') | mi(`z') | mi(cmo_name) | mi(risk_id_name) | mi(sid_cepr)
+
+            * Drop risk sets with no offer variation
+            bys risk_id_name: egen min_z = min(`z')
+            bys risk_id_name: egen max_z = max(`z')
+            keep if min_z != max_z
+            drop min_z max_z
+            count
+            di "N Remaining: `r(N)'"
+
+            if `r(N)' == 0 {
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+
+            * Generate weights
+            bys sid_cepr: gen ww = _N
+            gen w = 1/ww
+
+            cap areg `dep' `z' [aw=w], a(risk_id_name) cluster(sid_cepr)
+
+            * If insufficient obs post this
+            if _rc{
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+
+            di _b[`z']
+
+            * Store the coefficient, se, and N
+            scalar coef = _b[`z']
+            scalar se_store = _se[`z']
+            scalar N_obs = e(N)
+
+            egen tag = tag(sid_cepr) if e(sample)
+            qui count if tag
+            scalar N_ppl = r(N)
+            qui su `dep' [aw=w] if `z' == 0 & e(sample)
+            scalar control_mean = r(mean)
+            egen tag2 = tag(sid_cepr) if offer == 0 & e(sample)
+            qui count if tag2
+            scalar N_ppl_control = r(N)
+
+            post handle ("`s'") ("`dep'") ("`z'") (coef) (se_store) (N_obs) (N_ppl) (control_mean) (N_ppl_control)
+
+            restore
+        }
+    }
+}
+
+postclose handle
+
+use `results', clear
+sort instrument sample depvar
+order sample instrument depvar
+
+save "P:\Long_Term_Charter\data\clean\new_build\cmo_analysis\BALANCE\balance.dta", replace

ltc_code-0.2.1/src/ltc_code/20260706_ceprscripts/FS.do

@@ -0,0 +1,129 @@
+use "P:\Long_Term_Charter\data\clean\new_build\cmo_apps\stata_collapsed\cmo_apps_sid.dta", clear
+cap postclose handle
+clear matrix
+scalar drop _all
+
+* Generate numeric cohort
+destring cohort_lottery, gen(cohort_numeric)
+
+* Create tempfile to store results
+tempfile results
+postfile handle str50 sample str30 depvar str30 instrument double coef double se double N_obs double N_ppl double control_mean double N_ppl_control using `results', replace
+
+
+* Store local for depvars and indep vars
+local depvars enroll enrollment_years
+local instruments offer initial_offer
+
+* Overall and CMO specific samples
+levelsof cmo_name, local(cmos)
+
+local cmo_census
+foreach c of local cmos {
+    local cmo_census `"`cmo_census' "`c' (1990-2005)""'
+}
+
+local samples `" "No Harmony" "Overall" "Census Sample" "KIPP Sample" `cmos' `cmo_census' "'
+* local samples `" "No Harmony" "Overall" "Census Sample"  "KIPP Sample" "'
+*local samples `" `cmo_census' "'
+
+
+
+foreach s of local samples {
+    foreach dep of local depvars {
+        foreach z of local instruments {
+            preserve
+
+            di "SAMPLE: `s'"
+            di "INSTRUMENT: `z'"
+            di "DEPVAR: `dep'"
+
+            * Filter to relevant sample
+            if "`s'" == "No Harmony" {
+                drop if cmo_name == "Harmony"
+            }
+			else if "`s'" == "Census Sample" {
+				keep if inrange(cohort_numeric, 1990, 2005)
+				drop if cmo_name == "Harmony"
+            }
+			else if "`s'" == "KIPP Sample" {
+				keep if inrange(cohort_numeric, 1990, 2003)
+				keep if inlist(cmo_name, "KIPP_NJ", "KIPP TX", "KIPP NC", "KIPP MA")
+            }
+            else if "`s'" == "Overall" {
+
+            }
+			else if strpos("`s'", " (1990-2005)") {
+				keep if cmo_name == subinstr("`s'", " (1990-2005)", "", .)
+				keep if inrange(cohort_numeric, 1990, 2005)
+            }
+            else {
+                keep if cmo_name == "`s'"
+            }
+			
+            * Require all non missing regression vars
+            drop if mi(`dep') | mi(`z') | mi(cmo_name) | mi(risk_set) | mi(sid_cepr)
+			count
+			if `r(N)' == 0 {
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+			
+            * Drop risk sets with no offer variation
+            bys risk_set: egen min_z = min(`z')
+            bys risk_set: egen max_z = max(`z')
+            keep if min_z != max_z
+            drop min_z max_z
+            count
+            di "N Remaining: `r(N)'"
+
+            if `r(N)' == 0 {
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+
+            * Generate weights
+            bys sid_cepr: gen ww = _N
+            gen w = 1/ww
+
+            cap areg `dep' `z' [aw=w], a(risk_set) vce(robust)
+
+            * If insufficient obs post this
+            if _rc{
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+			
+            di _b[`z']
+
+            * Store the coefficient, se, and N
+            scalar coef = _b[`z']
+            scalar se_store = _se[`z']
+            scalar N_obs = e(N)
+
+            egen tag = tag(sid_cepr) if e(sample)
+            qui count if tag
+            scalar N_ppl = r(N)
+            qui su `dep' [aw=w] if `z' == 0 & e(sample)
+            scalar control_mean = r(mean)
+            egen tag2 = tag(sid_cepr) if offer == 0 & e(sample)
+            qui count if tag2
+            scalar N_ppl_control = r(N)
+
+            post handle ("`s'") ("`dep'") ("`z'") (coef) (se_store) (N_obs) (N_ppl) (control_mean) (N_ppl_control)
+
+            restore
+        }
+    }
+}
+
+postclose handle
+
+use `results', clear
+sort instrument sample depvar
+order sample instrument depvar
+
+save "P:\Long_Term_Charter\data\clean\new_build\cmo_analysis\FS\FS.dta", replace

ltc_code-0.2.1/src/ltc_code/20260706_ceprscripts/ITT.do

@@ -0,0 +1,123 @@
+use "P:\Long_Term_Charter\data\clean\new_build\cmo_analysis\cmo_analysis.dta", clear
+cap postclose handle
+clear matrix
+scalar drop _all
+
+* Generate numeric cohort
+destring cohort_lottery, gen(cohort_numeric)
+
+* Create tempfile to store results
+tempfile results
+postfile handle str50 sample str30 depvar str30 instrument double coef double se double N_obs double N_ppl double control_mean double N_ppl_control using `results', replace
+
+
+* Store local for depvars and indep vars
+local depvars att_2yr_byY2 att_4yr_byY2 att_any_byY2 
+*cmp_AA_byY6 cmp_BA_byY6 cmp_any_byY6
+local instruments offer initial_offer
+
+* Overall and CMO specific samples
+levelsof cmo_name, local(cmos)
+
+local cmo_census
+foreach c of local cmos {
+    local cmo_census `"`cmo_census' "`c' (1990-2005)""'
+}
+
+local samples `" "No Harmony" "Overall" "Census Sample" "KIPP Sample" `cmos' `cmo_census' "'
+* local samples `" "No Harmony" "Overall" "Census Sample"  "KIPP Sample" "'
+*local samples `" `cmo_census' "'
+
+
+foreach s of local samples {
+    foreach dep of local depvars {
+        foreach z of local instruments {
+            preserve
+
+            di "SAMPLE: `s'"
+            di "INSTRUMENT: `z'"
+            di "DEPVAR: `dep'"
+
+            * Filter to relevant sample
+            if "`s'" == "No Harmony" {
+                drop if cmo_name == "Harmony"
+            }
+			else if "`s'" == "Census Sample" {
+				keep if inrange(cohort_numeric, 1990, 2005)
+				drop if cmo_name == "Harmony"
+            }
+            else if "`s'" == "KIPP Sample" {
+				keep if inrange(cohort_numeric, 1990, 2003)
+				keep if inlist(cmo_name, "KIPP_NJ", "KIPP TX", "KIPP NC", "KIPP MA")
+            }
+            else if "`s'" == "Overall" {
+
+            }
+			else if strpos("`s'", " (1990-2005)") {
+				keep if cmo_name == subinstr("`s'", " (1990-2005)", "", .)
+				keep if inrange(cohort_numeric, 1990, 2005)
+            }
+            else {
+                keep if cmo_name == "`s'"
+            }
+
+            * Require all non missing regression vars
+            drop if mi(`dep') | mi(`z') | mi(cmo_name) | mi(risk_id_name) | mi(sid_cepr)
+
+            * Drop risk sets with no offer variation
+            bys risk_id_name: egen min_z = min(`z')
+            bys risk_id_name: egen max_z = max(`z')
+            keep if min_z != max_z
+            drop min_z max_z
+            count
+            di "N Remaining: `r(N)'"
+
+            if `r(N)' == 0 {
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+
+            * Generate weights
+            bys sid_cepr: gen ww = _N
+            gen w = 1/ww
+
+            cap areg `dep' `z' [aw=w], a(risk_id_name) cluster(sid_cepr)
+
+            * If insufficient obs post this
+            if _rc{
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+
+            di _b[`z']
+
+            * Store the coefficient, se, and N
+            scalar coef = _b[`z']
+            scalar se_store = _se[`z']
+            scalar N_obs = e(N)
+
+            egen tag = tag(sid_cepr) if e(sample)
+            qui count if tag
+            scalar N_ppl = r(N)
+            qui su `dep' [aw=w] if `z' == 0 & e(sample)
+            scalar control_mean = r(mean)
+            egen tag2 = tag(sid_cepr) if offer == 0 & e(sample)
+            qui count if tag2
+            scalar N_ppl_control = r(N)
+
+            post handle ("`s'") ("`dep'") ("`z'") (coef) (se_store) (N_obs) (N_ppl) (control_mean) (N_ppl_control)
+
+            restore
+        }
+    }
+}
+
+postclose handle
+
+use `results', clear
+sort instrument sample depvar
+order sample instrument depvar
+
+save "P:\Long_Term_Charter\data\clean\new_build\cmo_analysis\ITT\ITT.dta", replace

ltc_code-0.2.1/src/ltc_code/20260706_ceprscripts/TOT.do

@@ -0,0 +1,123 @@
+use "P:\Long_Term_Charter\data\clean\new_build\cmo_analysis\cmo_analysis_wide.dta", clear
+cap postclose handle
+clear matrix
+scalar drop _all
+
+* Generate numeric cohort
+destring cohort_lottery, gen(cohort_numeric)
+
+* Create tempfile to store results
+tempfile results
+postfile handle str50 sample str30 depvar str30 instrument double coef double se double N_obs double N_ppl double control_mean double N_ppl_control using `results', replace
+
+
+* Store local for depvars and indep vars
+local depvars att_2yr_byY2 att_4yr_byY2 att_any_byY2 
+*cmp_AA_byY6 cmp_BA_byY6 cmp_any_byY6
+local instruments offer initial_offer
+
+* Overall and CMO specific samples
+levelsof cmo_name, local(cmos)
+
+local cmo_census
+foreach c of local cmos {
+    local cmo_census `"`cmo_census' "`c' (1990-2005)""'
+}
+
+local samples `" "No Harmony" "Overall" "Census Sample" "KIPP Sample" `cmos' `cmo_census' "'
+* local samples `" "No Harmony" "Overall" "Census Sample"  "KIPP Sample" "'
+*local samples `" `cmo_census' "'
+
+
+foreach s of local samples {
+    foreach dep of local depvars {
+        foreach z of local instruments {
+            preserve
+
+            di "SAMPLE: `s'"
+            di "INSTRUMENT: `z'"
+            di "DEPVAR: `dep'"
+
+            * Filter to relevant sample
+            if "`s'" == "No Harmony" {
+                drop if cmo_name == "Harmony"
+            }
+			else if "`s'" == "Census Sample" {
+				keep if inrange(cohort_numeric, 1990, 2005)
+				drop if cmo_name == "Harmony"
+            }
+            else if "`s'" == "KIPP Sample" {
+				keep if inrange(cohort_numeric, 1990, 2003)
+				keep if inlist(cmo_name, "KIPP_NJ", "KIPP TX", "KIPP NC", "KIPP MA")
+            }
+            else if "`s'" == "Overall" {
+
+            }
+			else if strpos("`s'", " (1990-2005)") {
+				keep if cmo_name == subinstr("`s'", " (1990-2005)", "", .)
+				keep if inrange(cohort_numeric, 1990, 2005)
+            }
+            else {
+                keep if cmo_name == "`s'"
+            }
+
+            * Require all non missing regression vars
+            drop if mi(`dep') | mi(`z') | mi(cmo_name) | mi(risk_set) | mi(sid_cepr)
+
+            * Drop risk sets with no offer variation
+            bys risk_set: egen min_z = min(`z')
+            bys risk_set: egen max_z = max(`z')
+            keep if min_z != max_z
+            drop min_z max_z
+            count
+            di "N Remaining: `r(N)'"
+
+            if `r(N)' == 0 {
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+
+            * Generate weights
+            bys sid_cepr: gen ww = _N
+            gen w = 1/ww
+
+            cap ivreghdfe `dep' (enroll = `z') [aw=w], a(risk_set) robust
+
+            * If insufficient obs post this
+            if _rc{
+                post handle ("`s'") ("`dep'") ("`z'") (.) (.) (0) (.) (.) (.)
+                restore
+                continue
+            }
+
+            di _b[enroll]
+
+            * Store the coefficient, se, and N
+            scalar coef = _b[enroll]
+            scalar se_store = _se[enroll]
+            scalar N_obs = e(N)
+
+            egen tag = tag(sid_cepr) if e(sample)
+            qui count if tag
+            scalar N_ppl = r(N)
+            qui su enroll [aw=w] if enroll == 0 & e(sample)
+            scalar control_mean = r(mean)
+            egen tag2 = tag(sid_cepr) if offer == 0 & e(sample)
+            qui count if tag2
+            scalar N_ppl_control = r(N)
+
+            post handle ("`s'") ("`dep'") ("`z'") (coef) (se_store) (N_obs) (N_ppl) (control_mean) (N_ppl_control)
+
+            restore
+        }
+    }
+}
+
+postclose handle
+
+use `results', clear
+sort instrument sample depvar
+order sample instrument depvar
+
+save "P:\Long_Term_Charter\data\clean\new_build\cmo_analysis\TOT\TOT.dta", replace