larsvilhuber’s gists

larsvilhuber / TEST.md

Last active October 8, 2026 19:38

Miniconda for AEAREP-9960

  source "$HOME/miniconda3/etc/profile.d/conda.sh"
  conda create -n rmc-py314-mkl -c defaults python=3.14 numpy scipy "blas=*=mkl" -y

larsvilhuber / part1_of_main.do

Last active December 19, 2025 15:47

Generic minimal AEA-specific config

	// insert this near the top of your main/master/run.do

	======== AEA specific parts =============*
	local scenario "B" // Scenario B: main is in project top-level directory
	local pwd : pwd // This always captures the current directory

	if "`scenario'" == "A" { // If in Scenario A, we need to change directory first
	cd ..
	}
	// Now capture the directory to use as rootdir. Do not change this in the submitted version!

larsvilhuber / create_synthetic_sdv.py

Created October 17, 2025 15:35

Create synthetic data with SDV

	import pandas as pd
	from sdv.single_table import GaussianCopulaSynthesizer
	from sdv.metadata import Metadata
	import os

	# created based on simple example at https://docs.sdv.dev/sdv/explore/sdv-community

	# Read my Stata format data from the same directory as this code
	current_dir = os.path.dirname(os.path.abspath(__file__))
	data_path = os.path.join(current_dir, 'Individual.dta')

larsvilhuber / analysis.do

Created May 26, 2025 02:57

Day 1 Figure

	use "data/derived/analysis.dta", clear

	generate same_country = (iso_o == iso_d )
	generate ln_value = ln(VALUE_EURO)
	generate ln_N = ln(N)
	generate ln_dist = ln(dist)

	tabulate same_country [fw = N]

	regress ln_value ln_dist, robust

larsvilhuber / create_analysis_sample.do

Created May 26, 2025 02:50

Day 1 Cleaning code

	import delimited "data/country-codes.csv", varnames(1) case(preserve) encoding("utf-8") clear
	keep ISO31661Alpha2 ISO31661Alpha3
	save "data/country-codes.dta", replace

	import delimited "data/ted-sample.csv", varnames(1) case(preserve) encoding("utf-8") clear
	keep ID_NOTICE_CAN YEAR CAE_NAME ISO_COUNTRY_CODE VALUE_EURO WIN_NAME WIN_COUNTRY_CODE

	rename ISO_COUNTRY_CODE ISO31661Alpha2
	merge m:1 ISO31661Alpha2 using "data/derived/country-codes.dta", nogenerate keep(master match)
	rename ISO31661Alpha3 iso_d

larsvilhuber / download_ipeds.do

Created June 7, 2024 19:53

Downloading files from IPEDS (untested)

	// Create a program to download IPDS Completion data files
	// The program will check if the file exists and if it does, it will download the file
	// If the file does not exist, it will print a message to the user



	capture program drop download_ipeds_completion
	program define download_ipeds_completion
	version 17.0
	args year

larsvilhuber / checks_n.do

Created March 29, 2024 13:22

Checks and controls in Stata data flow

	. *Check sample sizes
	. qui desc using ${scratch}/simple2010
	. assert r(N)==1941473
	. qui desc using ${scratch}/simple2018
	. assert r(N)==2011482

larsvilhuber / download.R

Last active April 29, 2025 16:18

Downloading from Zenodo (2025)

	# this has worked as of 2024-03-01
	library(jsonlite)
	library(dplyr)
	library(tidyr)
	library(stringr)

	# Specific DOI - resolves to a fixed version
	zenodo.id <- "6568295"
	# We will recover the rest from Zenodo API
	zenodo.api = "https://zenodo.org/api/records/"

larsvilhuber / main.m

Last active June 7, 2023 19:36

Sample main.m file for repeated strategies

	# Instructions to author:
	# - Remove ALL paths in each of the files called. They are not needed.
	# - Change all lines related to the bootstrap to read:
	# [x_M, vc_bs] = bs_abf(match, id_all, strg, K, start_bs, session, bootstrapsize);
	#
	# This replaces the entire section from the old README:
	# Once the MATLAB data is generated, then need to run the MATLAB script est_abf_.m (depending on whether interested in PD, BoSLI, BOSHI and on the shorter or expanded versions) and to get the results after this script is finished run the corresponding est_test_.m script. In all these scripts the appropriate path needs to be updated as indicated. Additionally, in all est_abf_*.m scripts the desired bootstrapping needs to be adjusted for each separate estimation at the parts of code starting with [x_M, vc_bs] (for example at line 87 of est_abf_PD_twohalves.m). The MATLAB scripts call on functions bs_abf.m, chis_prb_gf.m, lik_i_abf.m, ll_abf.m, opt_abf.m, and waldtest_gf.m which are all included in the strategies fol

larsvilhuber / 01_download_sipp1992.do

Created March 3, 2023 03:18

A bunch of Stata do files to download SIPP files from NBER

	copy https://data.nber.org/sipp/1992/sipp92w1.zip $extractdir/92w1.zip
	copy https://data.nber.org/sipp/1992/sip92w1.dct $extractdir/92w1.dct
	copy https://data.nber.org/sipp/1992/sipp92w2.zip $extractdir/92w2.zip
	copy https://data.nber.org/sipp/1992/sip92w2.dct $extractdir/92w2.dct
	copy https://data.nber.org/sipp/1992/sipp92w3.zip $extractdir/92w3.zip
	copy https://data.nber.org/sipp/1992/sip92w3.dct $extractdir/92w3.dct
	copy https://data.nber.org/sipp/1992/sipp92w4.zip $extractdir/92w4.zip
	copy https://data.nber.org/sipp/1992/sip92w4.dct $extractdir/92w4.dct
	copy https://data.nber.org/sipp/1992/sipp92w5.zip $extractdir/92w5.zip
	copy https://data.nber.org/sipp/1992/sip92w5.dct $extractdir/92w5.dct

Lars Vilhuber larsvilhuber