we use this function to get prices for all available PRAs in the country, i.e. not matching to our 200/100 city samples. we use those prices for the neighborhood IV exercise and for plotting purposes.

all_pra_prices(yr1950 = 1950, readdisk = TRUE)