/*
File: 	17th_pums_household.do
Date:	Feb 12, 2004
Desc:	Creates a consistent set of selected household variables from 
		the 1980, 1990, and 2000 Decennial Census PUMS 1%
*/

/* Determine data year */

local year=`1' /* get data year from 17th_pums_master.do */

/* geography */

	/* region */

if `year'==1980 | `year'==1990 {
gen byte region=.
replace region=1 if division==1 | division==2
replace region=2 if division==3 | division==4
replace region=3 if 5<=division & division<=7
replace region=4 if division==8 | division==9
}
if `year'==2000{
capture lab drop region
}
lab var region "Region"
lab def region 1 Northeast 2 Midwest 3 South 4 West
lab val region region
notes region: Census region
notes region: Missing if region is not identifiable

	/* division */

if `year'==1980 | `year'==1990 | `year'==2000 {
replace division=. if division==0
}
lab var division "Census division"
#delimit ;
lab def div
1 "New England"
2 "Middle Atlantic"
3 "East North Central"
4 "West North Central"
5 "South Atlantic"
6 "East South Central"
7 "West South Central"
8 "Mountain"
9 "Pacific";
#delimit cr
lab val division div
notes division: Census division
notes division: Missing if division not identifiable

	/* state */
	
if `year'==1980 | `year'==1990 | `year'==2000 {
capture lab drop state
#delimit ;
lab def state
 1 AL
 2 AK
 4 AZ
 5 AR
 6 CA
 8 CO
 9 CT
10 DE
11 DC
12 FL
13 GA
15 HI
16 ID
17 IL
18 IN
19 IA
20 KS
21 KY
22 LA
23 ME
24 MD
25 MA
26 MI
27 MN
28 MS
29 MO
30 MT
31 NE
32 NV
33 NH
34 NJ
35 NM
36 NY
37 NC
38 ND
39 OH
40 OK
41 OR
42 PA
44 RI
45 SC
46 SD
47 TN
48 TX
49 UT
50 VT
51 VA
53 WA
54 WV
55 WI
56 WY
72 PR
99 99
;
#delimit cr
lab val state state
}
notes state: PUMS 1% sample does not allow for unique identification of /*
*/ state for households residing in PUMAs that overlap state boundaries
notes state: Households in PUMAs that overlap state boundaries have state /*
*/ set equal to 99
	
	/* puma */
	
gen puma1980=.
gen puma1990=.
gen puma2000=.

if `year'==1980{

	/* 
	PUMAs in 1980 identified uniquely using state, country group, and
	SMSA; current version of this program has not yet done this matching
	*/

}
if `year'==1990 {
replace puma1990=puma
}
if `year'==2000 {
replace puma2000=puma1
}
lab var puma1980 "Public Use Microdata Area"
notes puma1980: Coding is state dependent
notes puma1980: Not yet coded for 17TH Street PUMS extract

lab var puma1990 "Public Use Microdata Area"
notes puma1990: Coding is state dependent
notes puma1990: Not consistent with puma2000

lab var puma2000 "Public Use Microdata Area"
notes puma2000: Super-Public Use Microdata Area code
notes puma2000: Designates area of 400,000 or more population
notes puma2000: Coding is state dependent
notes puma2000: Not consistent with puma1990

	/* msa-pmsa */

gen msa=.

if `year'==1980 {
replace msa=smsa
}
if `year'==1990 {
replace msa=msapmsa
}
if `year'==2000 {
replace msa=msapmsa1
}

	/* modifications for consistency across 1980-2000 */

/* The MSA variable used here is based on the IPUMS <www.ipums.org>
metarea variable and is essentially consistent over the period
1980-2000. For a full discussion of metropolitan areas see the IPUMS
web page <http://www.ipums.org/usa/hgeographic/metareaa.html>.

Some key issues raised in the IPUMS documentation:

(1) In the 1970 and 1980 censuses the variable used here was described as 
"Standard Metropolitan Area" (SMA). In 1990 and 2000, the variable was 
described as "Metropolitan Area" (MA), which included: free-standing 
"Metropolitan Statistical Areas" (MSA); and, "Consolidated Metropolitan 
Statistical Areas" (CMSA), which consisted of two or more economically 
and socially linked metropolitan areas, called "Primary Metropolitan 
Statistical Areas" (PMSA). Many PMSAs were separate SMSAs or SMAs before 
1990.

(2)"Despite the terminological shifts, the general concept is the same for 
all years: A metropolitan area (whether an SMA, SMSA, MSA, or a CMSA 
containing several PMSAs) is an area consisting of a large population center
and adjacent communities (usually counties) that have a high degree of 
economic and social interaction with that center. Metropolitan areas often 
cross state lines. Some metropolitan areas contain more than one central city.
In previous years, those cities may either have been the central cities of 
separate metropolitan areas or not part of any metropolitan areas."

(3) In 1980, SMSA definitions were essentially identical to those used in
1970. The 1980 1% sample identifies 282 SMSAs (the 5% sample only identifies
180 SMSAs that do not cross state lines). For confidentiality reasons, the 
1980 1% sample pools 36 additional SMSAs that have fewer than 100,000 people
into 18 pairs of SMSAs with combined populations above 100,000. The 18 paired
SMSAs are:

    * Bangor, ME and Lewiston-Auburn, ME 
    * Tioga County, NY (part of the Binghampton, NY/PA SMSA) and Elmira, NY 
    * Bismark, ND and Grand Forks, ND/MN 
    * Bloomington, IN and Owensboro, KY 
    * Bristol, CT and Meriden, CT 
    * Bryan-College Station, TX and Sherman-Dennison, TX 
    * Guilford County, NC (excluding Greensboro city and High Point city; 
      the entire county is part of the Greensboro-Winston-Salem-High Point, 
      NC SMSA) and Burlington, NC 
    * Casper, WY and Great Falls, MT 
    * Dubuque, IA and Iowa City, IA 
    * El Paso, TX (excluding El Paso city) and Las Cruces, NM 
    * Lawton, OK and Enid, OK 
    * Fitchburg-Leominster, MA and Pittsfield, MA 
    * Fort Walton Beach, FL and Panama City, FL 
    * La Crosse, WI and Rochester, MN 
    * Laredo, TX and Victoria, TX 
    * Topeka, KS (excluding Topeka City) and Lawrence, KS 
    * Saline County, AR (part of Little Rock-N. Little Rock, AR SMSA) and 
      Pine Bluff, AR 
    * Midland, TX and San Angelo, TX 

(4) Metropolitan areas are only partially identified in the 1980-2000 
decennial censuses. "The lowest level of geography available in [the PUMS] 
files is the PUMA (in 1990 files and the 2000 5% file), the Super-PUMA 
(in the 2000 1% file), and the County Group (in the 1980 5% file). In most 
cases, the boundaries for PUMAs (as well as SuperPUMAs and County Groups) 
were drawn so that metropolitan areas would completely overlap with one or 
more PUMAs. In some cases, however, a metropolitan area only partly overlapped 
with a given PUMA or set of PUMAs." Following the IPUMS procedure, the msa 
variable "...assigned the appropriate metropolitan area code to persons 
residing within PUMAs that the metropolitan area completely overlapped, but 
NOT to individuals residing within PUMAs that the metropolitan area only 
partly overlapped. ...[Don't] assume that the identifiable portion of 
a partly-identified metropolitan area is a representative sample of the 
entire metropolitan area... weighted population counts for these 
metropolitan areas will be low by amounts ranging from 5%-60% (since 
the unidentified individuals will not be counted as living in the area)."

(5) For consistency across decennial censuses and other reasons, nearby cities
are sometimes covered by the same code:

Anaheim-Santa Ana (See Los Angeles)
Anderson, SC (See Greenville)
Aurora-Elgin, IL (See Chicago)
Bay City, Michigan (See Saginaw)
Beaver County, PA (See Pittsburgh)
Bergen-Passaic (See New York)
Boulder, CO (See Denver)
Brazoria, TX (See Houston)
Bristol, CT (See Hartford)
Fall River, MA (See Providence)
Fort Worth, TX (See Dallas)
Gary, IN (See Chicago)
Jersey City, NJ (See New York)
Joliet, IL (See Chicago)
Lake County, IL (See Chicago)
Lowell, MA (See Boston)
Lawrence-Haverhill (See Boston)
Middlesex-Somerset-Hunterdon (See New York)
Middleton, CT (See Hartford)
Nassau-Suffolk (See New York)
New Britian, CT (See Hartford)
Niagara Falls, NY (See Buffalo)
Oakland, CA (See San Francisco)
Ogden, Utah (See Salt Lake)
Paterson-Clifton (See New York)
Pawtuckett, RI (See Providence)
Petersburg-Col.Heights, VA (See Richmond)
Rock Hill, SC (See Charlotte)
Springfield, OH (See Dayton)
Salem-Gloucester, MA (See Boston)
San Bernadino, CA (see Riverside)
Springfield, OH (See Dayton)
Vallejo-Fairfield-Napa, CA (See San Francico)
Vancouver, WA (see Portland)
Wilkes-Barre, PN (See Scranton)
Winston-Salem, NC (see Greensboro)

The variable msac below attempts to collapse the codes for the preceding cities
on a consistent basis.

*/

/*
Some additional city codes for consolidated metropolitan areas from the 
Missouri Census Data Center
http://mcdc2.missouri.edu/sastools/sas_formats/Smetro.sas

1122 "Boston-Worcester-Lawrence, MA-NH-ME-C(C)"
1602	"Chicago-Gary-Kenosha, IL-IN-WI (C)"
1642 "Cincinnati-Hamilton, OH-KY-IN (C)"
1692 "Cleveland-Akron, OH (C)"
1922 "Dallas-Fort Worth, TX (C)"
2082 "Denver-Boulder-Greeley, CO (C)"
2162 "Detroit-Ann Arbor-Flint, MI (C)"
3362 "Houston-Galveston-Brazoria, TX (C)"
4472 "Los Angeles-Riverside-Orange County, (C)"
4992 "Miami-Fort Lauderdale, FL (C)"
5082 "Milwaukee-Racine, WI (C)"
5602 "New York-Northern New Jersey-Long Isl(C)"
6162 "Philadelphia-Wilmington-Atlantic City, PA-NJ-DE-MD (C)"
6442 "Portland-Salem, OR-WA (C)"
6922 "Sacramento-Yolo, CA (C)"
7362 "San Francisco-Oakland-San Jose, CA (C)"

*/

/* some tidying */

if `year'==1980 | `year'==1990 | `year'==2000 {
replace msa=6960 /* Saginaw */ if msa==800 /* Bay City, Michigan */
replace msa=1310 /* Burlington, VT */ if msa==1305 /* Burlington, VT */
replace msa=2660 /* Florence, SC */ if msa==2655 /* Florence, SC */
replace msa=4640 /* Lynchburg, VA */ if msa==4645 /* unidentified */
replace msa=5720 /* "Norfolk-VA Beach-Newport News, VA" */ if msa==5680 /*
				Newport News-Hampton, VA */
replace msa=5640 /* Newark, NJ */ if msa==5645 /* mislabeled "Newark, OH" */
replace msa=7560 /* "Scranton-Wilkes-Barre, PA" */ if msa==5745 /* 
				"Northeastern Pennsylvania" */
replace msa=5600 /* "New York-Northeastern, NJ" */ if msa==6040 /*
				"Paterson-Clifton-Passaic, NJ" */
replace msa=6760 /*	"Richmond-Petersburg, VA" */ if msa==6140 /*
				"Petersburg-Colonial Heights-Hopewell, VA" */
replace msa=1520 /* "Charlotte-Gastonia-Rock Hill, SC" */ if msa==6885 /*
				Rock Hill, SC */
replace msa=2000 /* "Dayton-Springfield, OH" */ if msa==7960 /* 
				Springfield, OH */
}
lab var msa "MSA"
#delimit ;
lab def msa
40	"Abilene, TX"
80	"Akron, OH"
120	"Albany, GA"
160	"Albany-Schenectady-Troy, NY"
200	"Albuquerque, NM"
220	"Alexandria, LA"
240	"Allentown-Bethlehem-Easton, PA/NJ"
360  "Anaheim-Santa Ana (See Los Angeles)"
280	"Altoona, PA"
320	"Amarillo, TX"
380	"Anchorage, AK"
400	"Anderson, IN"
405  "Anderson, SC (See Greenville)"
440	"Ann Arbor, MI"
450	"Anniston, AL"
460	"Appleton-Oskosh-Neenah, WI"
480	"Asheville, NC"
500	"Athens, GA"
520	"Atlanta, GA"
560	"Atlantic City, NJ"
580	"Auburn-Opelika, AL"
600	"Augusta-Aiken, GA-SC"
620  "Aurora-Elgin, IL (See Chicago)"
640	"Austin, TX"
680	"Bakersfield, CA"
720	"Baltimore, MD"
730	"Bangor, ME"
740	"Barnstable-Yarmouth, MA"
760	"Baton Rouge, LA"
780	"Battle Creek, MI"
800	"Bay City, MI"
840	"Beaumont-Port Arthur-Orange, TX"
845	"Beaver County, PA (See Pittsburgh)"
860	"Bellingham, WA"
870	"Benton Harbor, MI"
875	"Bergen-Passaic (See New York)"
880	"Billings, MT"
920	"Biloxi-Gulfport, MS"
960	"Binghamton, NY"
1000	"Birmingham, AL"
1010	"Bismarck, ND"
1020	"Bloomington, IN"
1040	"Bloomington-Normal, IL"
1080 "Boise City, ID"
1120	"Boston, MA"
1125	"Boulder, CO (See Denver)"
1140	"Bradenton, FL"
1145	"Brazoria, TX (See Houston)"
1150	"Bremerton, WA"
1160	"Bridgeport, CT"
1200	"Brockton"
1240	"Brownsville - Harlingen-San Benito, TX"
1260	"Bryan-College Station, TX"
1280	"Buffalo-Niagara Falls, NY"
1300	"Burlington, NC"
1310	"Burlington, VT"
1320	"Canton, OH"
1350	"Casper, WY"
1360	"Cedar Rapids, IA"
1400	"Champaign-Urbana-Rantoul, IL"
1440	"Charleston-N.Charleston,SC"
1480	"Charleston, WV"
1520	"Charlotte-Gastonia-Rock Hill, SC"
1540	"Charlottesville, VA"
1560	"Chattanooga, TN/GA"
1580	"Cheyenne, WY"
1600	"Chicago-Gary-Lake IL"
1620	"Chico, CA"
1640	"Cincinnati, OH/KY/IN"
1660	"Clarksville- Hopkinsville, TN/KY"
1680	"Cleveland, OH"
1720	"Colorado Springs, CO"
1740	"Columbia, MO"
1760	"Columbia, SC"
1800	"Columbus, GA/AL"
1840	"Columbus, OH"
1880	"Corpus Christi, TX"
1900	"Cumberland, MD/WV"
1920	"Dallas-Fort Worth, TX"
1930	"Danbury, CT"
1950	"Danville, VA"
1960	"Davenport, IA-Rock Island-Moline, IL"
2000	"Dayton-Springfield, OH"
2020	"Daytona Beach, FL"
2030	"Decatur, AL"
2040	"Decatur, IL"
2080	"Denver-Boulder-Longmont, CO"
2120	"Des Moines, IA"
2160	"Detroit, MI"
2180	"Dothan, AL"
2190	"Dover, DE"
2200	"Dubuque, IA"
2240	"Duluth-Superior, MN/WI"
2280	"Dutchess County, NY"
2290	"Eau Claire, WI"
2310	"El Paso, TX"
2320	"Elkhart-Goshen, IN"
2330	"Elmira, NY"
2340	"Enid, OK"
2360	"Erie, PA"
2400	"Eugene-Springfield, OR"
2440	"Evansville, IN/KY"
2480	"Fall River, MA (See Providence)"
2520	"Fargo-Moorhead, ND/MN"
2560	"Fayetteville, NC"
2580	"Fayetteville-Springdale, AR"
2600	"Fitchburg-Leominster, MA"
2620	"Flagstaff, AZ"
2640	"Flint, MI"
2650	"Florence, AL"
2655 "Florence, SC"
2660	"Florence, SC"
2670	"Fort Collins-Loveland, CO"
2680	"Fort Lauderdale-Hollywood-Pompano Beach, FL"
2700	"Fort Myers-Cape Coral, FL"
2710	"Fort Pierce, FL"
2720	"Fort Smith, AR/OK"
2750	"Fort Walton Beach, FL"
2760	"Fort Wayne, IN"
2800 "Fort Worth, TX (See Dallas)"
2840	"Fresno, CA"
2880	"Gadsden, AL"
2900	"Gainesville, FL"
2920	"Galveston-Texas City, TX"
2960 "Gary, IN (See Chicago)"
2970	"Glens Falls, NY"
2975	"Warren/Washington Counties, NY"
2980	"Goldsboro, NC"
2990	"Grand Forks, ND/MN"
3000	"Grand Rapids, MI"
3040	"Great Falls, MT"
3060	"Greeley, CO"
3080	"Green Bay, WI"
3120	"Greensboro-Winston Salem"
3160	"Greenville-Spartenburg-Anderson, SC"
3180	"Hagerstown, MD"
3200	"Hamilton-Middleton, OH"
3240	"Harrisburg-Lebanon-Carlisle, PA"
3280	"Hartford-Bristol-Middleton-New Britian, CT"
3290	"Hickory-Morgantown, NC"
3320	"Honolulu, HI"
3350	"Houma-Thibodoux, LA"
3360	"Houston-Brazoria, TX"
3400	"Huntington-Ashland, WV/KY/OH"
3440	"Huntsville, AL"
3480	"Indianapolis, IN"
3500	"Iowa City, IA"
3520	"Jackson, MI"
3560	"Jackson, MS"
3580	"Jackson, TN"
3590	"Jacksonville, FL"
3600	"Jacksonville, NC"
3605	"Onslow County, NC"
3610	"Jamestown-Dunkirk, NY"
3620	"Janesville-Beloit, WI"
3640	"Jersey City, NJ (See New York)"
3660	"Johnson City-Kingsport-Bristol, TN/VA"
3680	"Johnstown, PA"
3690 "Joliet, IL (See Chicago)"
3710	"Joplin, MO"
3720	"Kalamazoo-Portage, MI"
3740	"Kankakee, IL"
3760	"Kansas City, MO-KS"
3800	"Kenosha, WI"
3810	"Killeen-Temple, TX"
3840	"Knoxville, TN"
3850	"Kokomo, IN"
3870	"LaCrosse, WI"
3880	"Lafayette, LA"
3920	"Lafayette-W. Lafayette,IN"
3960	"Lake Charles, LA"
3965 "Lake County, IL (See Chicago)"
3980	"Lakeland-Winterhaven, FL"
4000	"Lancaster, PA"
4040	"Lansing-E. Lansing, MI"
4080	"Laredo, TX"
4100	"Las Cruces, NM"
4120	"Las Vegas, NV"
4150	"Lawrence, KS"
4160	"Lawrence-Haverhill (See Boston)"
4200	"Lawton, OK"
4240	"Lewiston-Auburn, ME"
4280	"Lexington-Fayette, KY"
4320	"Lima, OH"
4360	"Lincoln, NE"
4400	"Little Rock-North Little Rock, AR"
4410	"Long Branch-Asbury Park, NJ"
4420	"Longview-Marshall, TX"
4440	"Lorain-Elyria, OH"
4480	"Los Angeles-Long Beach, CA"
4520	"Louisville, KY/IN"
4560	"Lowell, MA (See Boston)"
4600	"Lubbock, TX"
4640	"Lynchburg, VA"
4680	"Macon-Warner Robins, GA"
4720	"Madison, WI"
4760	"Manchester, NH"
4800	"Mansfield, OH"
4880	"McAllen-Edinburg-Pharr-Mission, TX"
4890	"Medford, OR"
4900	"Melbourne-Titusville-Cocoa-Palm Bay, FL"
4920	"Memphis, TN/AR/MS"
4940	"Merced, CA"
5000	"Miami-Hialeah, FL"
5015 "Middlesex-Somerset-Hunterdon, NJ (P)" /* To New York, NY */
5020	"Middleton, CT (See Hartford)"
5040	"Midland, TX"
5080	"Milwaukee Co., WI"
5120	"Minneapolis-St. Paul, MN"
5160	"Mobile, AL"
5170	"Modesto, CA"
5190	"Monmouth-Ocean, NJ"
5200	"Monroe, LA"
5240	"Montgomery, AL"
5280	"Muncie, IN"
5320	"Muskegon-Norton Shores-Muskegon Heights, MI"
5330	"Myrtle Beach, SC"
5340	"Naples, FL"
5345	"Collier County, FL"
5350	"Nashua, NH"
5360	"Nashville, TN"
5380 "Nassau-Suffolk, NY (P)" /* To New York, NY */
5400	"New Bedford, MA"
5440 "New Britian, CT (See Hartford)"
5460	"New Brunswick-Perth Amboy-Sayreville, NJ"
5480	"New Haven-Meriden, CT"
5520	"New London-Norwich, CT/RI"
5560	"New Orleans, LA"
5600	"New York-Northeastern, NJ"
5640	"Newark, NJ"
5660	"Newburgh-Middletown, NY"
5700	"Niagara Falls, NY (See Buffalo)"
5720	"Norfolk-VA Beach-Newport News, VA"
5760	"Norwalk, CT"
5775	"Oakland, CA (See San Francisco)"
5790	"Ocala, FL"
5800	"Odessa, TX"
5880	"Oklahoma City, OK"
5910	"Olympia, WA"
5920	"Omaha, NE/IA"
5945 "Orange County, CA (P)" /* To Los Angeles [JTS] */
5950	"Orange, NY"
5960	"Orlando, FL"
5990	"Owensboro, KY"
6000	"Oxnard-Ventura, CA"
6010	"Panama City, FL"
6020	"Parkersburg-Marietta,WV/OH"
6025	"Jackson County, MS"
6030	"Pascagoula-Moss Point, MS"
6060	"Pawtuckett, RI (See Providence)"
6080	"Pensacola, FL"
6120	"Peoria, IL"
6160	"Philadelphia, PA/NJ"
6200	"Phoenix, AZ"
6240	"Pine Bluff, AR"
6280	"Pittsburgh-Beaver Valley, PA"
6320	"Pittsfield, MA"
6400	"Portland, ME"
6440	"Portland-Vancouver, OR"
6450	"Portsmouth-Dover-Rochester, NH/ME"
6460	"Poughkeepsie, NY"
6480	"Providence-Fall River-Pawtuckett, MA"
6520	"Provo-Orem, UT"
6560	"Pueblo, CO"
6580	"Punta Gorda, FL"
6600	"Racine, WI"
6640	"Raleigh-Durham, NC"
6660	"Rapid City, SD"
6680	"Reading, PA"
6690	"Redding, CA"
6720	"Reno, NV"
6740	"Richland-Kennewick-Pasco, WA"
6760	"Richmond-Petersburg, VA"
6780	"Riverside-San Bernadino,CA"
6800	"Roanoke, VA"
6820	"Rochester, MN"
6840	"Rochester, NY"
6880	"Rockford, IL"
6890	"Rocky Mount, NC"
6920	"Sacramento, CA"
6960	"Saginaw-Bay City-Midland, MI"
6980	"St. Cloud, MN"
7000	"St. Joseph, MO"
7040	"St. Louis, MO-IL"
7080	"Salem, OR"
7090	"Salem-Gloucester, MA"
7120	"Salinas-Sea Side-Monterey, CA"
7140	"Salisbury-Concord, NC"
7160	"Salt Lake City-Ogden, UT"
7200	"San Angelo, TX"
7240	"San Antonio, TX"
7320	"San Diego, CA"
7360	"San Francisco-Oakland-Vallejo, CA"
7400	"San Jose, CA"
7460	"San Luis Obispo-Atascad-P Robles, CA"
7470	"Santa Barbara-Santa Maria-Lompoc, CA"
7480	"Santa Cruz, CA"
7485	"Santa Cruz, CA" /* coding change? */
7490	"Santa Fe, NM"
7500	"Santa Rosa-Petaluma, CA"
7510	"Sarasota, FL"
7520	"Savannah, GA"
7560	"Scranton-Wilkes-Barre, PA"
7600	"Seattle-Everett, WA"
7610	"Sharon, PA"
7620	"Sheboygan, WI"
7640	"Sherman-Denison, TX"
7680	"Shreveport, LA"
7720	"Sioux City, IA/NE"
7760	"Sioux Falls, SD"
7800	"South Bend-Mishawaka, IN"
7840	"Spokane, WA"
7880	"Springfield, IL"
7920	"Springfield, MO"
8000	"Springfield-Holyoke-Chicopee, MA"
8040	"Stamford, CT"
8050	"State College, PA"
8080	"Steubenville-Weirton,OH/WV"
8120	"Stockton, CA"
8140	"Sumter, SC"
8160	"Syracuse, NY"
8200	"Tacoma, WA"
8240	"Tallahassee, FL"
8280	"Tampa-St. Petersburg-Clearwater, FL"
8320	"Terre Haute, IN"
8360	"Texarkana, TX/AR"
8400	"Toledo, OH/MI"
8440	"Topeka, KS"
8480	"Trenton, NJ"
8520	"Tucson, AZ"
8560	"Tulsa, OK"
8600	"Tuscaloosa, AL"
8640	"Tyler, TX"
8680	"Utica-Rome, NY"
8720	"Vallejo-Fairfield-Napa, CA (See San Francisco)"
8725	"Vancouver, WA (see Portland)"
8730	"Ventura-Oxnard-Simi Valley, CA"
8735 "Ventura, CA (P)" /* To Ventura-Oxnard-Simi Valley [JTS] */
8750	"Victoria, TX"
8760	"Vineland-Milville-Bridgetown, NJ"
8780	"Visalia-Tulare -Porterville, CA"
8800	"Waco, TX"
8840	"Washington, DC/MD/VA"
8880	"Waterbury, CT"
8920	"Waterloo-Cedar Falls, IA"
8940	"Wausau, WI"
8960	"West Palm Beach-Boca Raton -Delray Beach, FL"
9000	"Wheeling, WV/OH"
9040	"Wichita, KS"
9080	"Wichita Falls, TX"
9140	"Williamsport, PA"
9160	"Wilmington, DE/NJ/MD"
9200	"Wilmington, NC"
9240	"Worcester, MA"
9260	"Yakima, WA"
9270	"Yolo, CA"
9280	"York, PA"
9320	"Youngstown-Warren, OH-PA"
9340	"Yuba City, CA"
9360	"Yuma, AZ"
9997	"Mixed MSA/PMSA NONMSA/PMSA area"
9998	"2 or more MSAs"
9999	"Not in MA"
; 
#delimit cr
lab val msa msa
notes msa: Primary metropolitan statistical area
notes msa: See 17th_pums_household.do for further details

	/* collapsed version of msa variable */

gen msac=msa
replace msac=4480 if msa==360 /* Anaheim-Santa Ana, CA to Los Angeles, CA */
replace msac=3160 if msa==405 /* Anderson, SC to Greenville, SC */
replace msac=1600 if msa==620 /* Aurora-Elgin, IL to Chicago, IL */
replace msac=6280 if msa==845 /* Beaver County, PA to Pittsburgh, PA */
replace msac=5600 if msa==875 /* Bergen-Passaic, NJ to New York, NY */
replace msac=2080 if msa==1125 /* Boulder, CO to Denver, CO */
replace msac=3360 if msa==1145 /* Brazoria, TX to Houston, TX */
replace msac=6480 if msa==2480 /* Fall River, MA to Providence, RI */
replace msac=1920 if msa==2800 /* Fort Worth, TX to Dallas, TX */
replace msac=1600 if msa==2960 /* Gary, IN to Chicago, IL */
replace msac=5600 if msa==3640 /* Jersey City, NJ to New York, NY */
replace msac=1600 if msa==3690 /* Joliet, IL to Chicago, IL */
replace msac=1600 if msa==3965 /* Lake County, IL to Chicago, IL */
replace msac=1120 if msa==4560 /* Lowell, MA to Boston, MA */
replace msac=1120 if msa==4160 /* Lawrence-Haverhill, MA to Boston, MA */
replace msac=5600 if msa==5015 /* Middlesex-Somerset-Hunterdon,	NJ 
				  		    to New York, NY */
replace msac=3280 if msa==5020 /* Middleton, CT to Hartford, CT */
replace msac=5600 if msa==5380 /* Nassau-Suffolk, NY to New York, NY */
replace msac=3280 if msa==5440 /* New Britian, CT to Hartford, CT */
replace msac=1280 if msa==5700 /* Niagara Falls, NY to Buffalo, NY */
replace msac=7360 if msa==5775 /* Oakland, CA to San Francisco, CA */
replace msac=4480 if msa==5945 /* Orange County, CA to Los Angeles, CA [JTS] */
replace msac=6480 if msa==6060 /* Pawtuckett, RI to Providence, RI */
replace msac=1120 if msa==7090 /* Salem-Gloucester, MA to Boston, MA */
replace msac=7480 if msa==7485 /* Santa Cruz, CA to Santa Cruz, CA [JTS] */
replace msac=7360 if msa==8720 /* Vallejo-Fairfield-Napa, CA 
						    to San Francisco, CA */
replace msac=6440 if msa==8725 /* Vancouver, WA to Portland, OR */
replace msac=8730 if msa==8735 /* Ventura, CA to Ventura-Oxnard-Simi, CA [JTS]*/
lab var msac "MSA"
lab val msac msa
notes msac: See 17th_pums_household.do for complete description

/* no relevant codes for 1990, 2000; merge may have taken place earlier
Bay City, Michigan (See Saginaw)
Bristol, CT (See Hartford)
Ogden, Utah (See Salt Lake)
Paterson-Clifton (See New York)
Petersburg-Col.Heights, VA (See Richmond)
Rock Hill, SC (See Charlotte)
Springfield, OH (See Dayton)
San Bernadino, CA (see Riverside)
Wilkes-Barre, PA (See Scranton)
Wilkes-Barre, PA (See Scranton)
*/

/* housing */

	/* number of persons in household */

notes persons: Households with zero persons are vacant and have been /*
*/ dropped from sample

	/* institutionalized group housing */

gen byte gq=.

if `year'==1980 {
replace gq=1 if 1<=gqtype & gqtype<=8
replace gq=0 if gqtype==0
}
if `year'==1990 {
replace gq=1 if gqinst==1 | gqinst==2
replace gq=0 if gqinst==0
}
if `year'==2000 {
replace gq=1 if unittype==1 | unittype==2
replace gq=0 if unittype==0
}
lab var gq "Group quarters"

gen byte instl=.

if `year'==1980 {
replace instl=0 if gqtype==0 | (5<=gqtype & gqtype<=8)
replace instl=1 if 1<=gqtype & gqtype<=4
}
if `year'==1990 {
replace instl=0 if gqinst==0 | gqinst==2
replace instl=1 if gqinst==1 /* institutional group quarters */
}
if `year'==2000 {
replace instl=0 if unittype==0 | unittype==2
replace instl=1 if unittype==1 /* institutional group quarters */
}
lab var instl "Institutional group quarters"
notes instl: Includes residents of mental hospitals, nursing homes, /*
*/ correctional institutions, and other institutions
notes instl: Excludes residents of military quarters, college dormitories, /*
*/ boarding houses and non-inmates living in institutions

	/* housing tenure */

gen byte ownhouse=.

if `year'==1980 {
replace ownhouse=0 if tenure==2 | tenure==3
replace ownhouse=1 if tenure==1
}
if `year'==1990 | `year'==2000 {
replace tenure=. if tenure==0
replace ownhouse=0 if tenure==3 | tenure==4
replace ownhouse=1 if tenure==1 | tenure==2
}

gen byte ownclear=.

if `year'==1980 {
replace ownclear=.
}
if `year'==1990 | `year'==2000 {
replace ownclear=0 if tenure~=.
replace ownclear=1 if tenure==2
}
lab var ownclear "Own home, no mortgage"

if `year'==1980 {
replace tenure=. /*	1980 tenure variable not consistent with 1990, 2000
				use ownhouse and ownclear */
}
capture lab drop tenure
#delimit ;
lab def tenure
1  "Own, mortgage"
2  "Own, no mortgage"
3  "Rent"
4  "No cash rent";
#delimit cr
lab val tenure tenure
notes tenure: Set to missing for 1980 because 1980 code differs from 1990, 2000
notes tenure: See also ownhouse, ownclear
notes tenure: Missing if group quarters or housing vacant

	/* property value */

gen byte pv1980=.
gen byte pv1990=.
gen byte pv2000=.
gen byte pv1=.

if `year'==1980 {

replace pv1980=value if 1<=value & value<=24

replace pv1=1 if value==1
replace pv1=2 if value==2
replace pv1=3 if value==3 | value==4
replace pv1=4 if value==5 | value==6
replace pv1=5 if value==7 | value==8
replace pv1=6 if value==9
replace pv1=7 if value==10
replace pv1=8 if value==11 | value==12
replace pv1=9 if value==13 | value==14
replace pv1=10 if value==15 | value==16
replace pv1=11 if value==17 | value==18
replace pv1=12 if value==19
replace pv1=13 if value==20
replace pv1=14 if value==21
replace pv1=15 if value==22
replace pv1=16 if value==23
replace pv1=17 if value==24

}
if `year'==1990 {

replace pv1990=value if 1<=value & value<=25

replace pv1=1 if value==1
replace pv1=2 if value==2
replace pv1=3 if value==3
replace pv1=4 if value==4
replace pv1=5 if value==5
replace pv1=6 if value==6
replace pv1=7 if value==7
replace pv1=8 if value==8 | value==9
replace pv1=9 if value==10 | value==11
replace pv1=10 if value==12 | value==13
replace pv1=11 if value==14 | value==15
replace pv1=12 if value==16
replace pv1=13 if value==17
replace pv1=14 if value==18
replace pv1=15 if value==19
replace pv1=16 if value==20 | value==21
replace pv1=17 if 22<=value & value<=25

}
if `year'==2000 {

replace pv2000=value if 1<=value & value<=24

replace pv1=1 if value==1
replace pv1=2 if value==2
replace pv1=3 if value==3
replace pv1=4 if value==4
replace pv1=5 if value==5
replace pv1=6 if value==6
replace pv1=7 if value==7
replace pv1=8 if value==8
replace pv1=9 if value==9
replace pv1=10 if value==10
replace pv1=11 if value==11
replace pv1=12 if value==12
replace pv1=13 if value==13
replace pv1=14 if value==14
replace pv1=15 if value==15
replace pv1=16 if value==16 | value==17
replace pv1=17 if 18<=value & value<=24

}

lab var pv1980 "Property value, 1980"
#delimit ;
lab def pv1980
1 " - $10,000"			
2 "$10,000 to $14,999"	
3 "$15,000 to $17,499"	
4 "$17,500 to $19,999"	
5 "$20,000 to $22,499"	
6 "$22,500 to $24,999"	
7 "$25,000 to $27,499"	
8 "$27,500 to $29,999"	
9 "$30,000 to $34,999"	
10 "$35,000 to $39,999"	
11 "$40,000 to $44,999"	
12 "$45,000 to $49,999"	
13 "$50,000 to $54,999"	 
14 "$55,000 to $59,999"
15 "$60,000 to $64,999"
16 "$65,000 to $69,999"	 
17 "$70,000 to $74,999"	 
18 "$75,000 to $79,999"	 
19 "$80,000 to $89,999"	 
20 "$90,000 to $99,999"	
21 "$100,000 to $124,999"	
22 "$125,000 to $149,999"	
23 "$150,000 to $199,999"	 
24 "$200,000 - ";
#delimit cr
lab val pv1980 pv1980
notes pv1980: Missing if group quarters, rental unit, vacant, or /*
*/ not for sale only
notes pv1980: Not available in 1990, 2000

lab var pv1990 "Property value, 1990"
#delimit ;
lab def pv1990
1 " - $10,000"			
2 "$10,000 to $14,999"	
3 "$15,000 to $19,999"	
4 "$20,000 to $24,999"	
5 "$25,000 to $29,999"	
6 "$30,000 to $34,999"	
7 "$35,000 to $39,999"	
8 "$40,000 to $44,999"	
9 "$45,000 to $49,999"	
10 "$50,000 to $54,999"	
11 "$55,000 to $59,999"	
12 "$60,000 to $64,999"	
13 "$65,000 to $69,999"	 
14 "$70,000 to $74,999"
15 "$75,000 to $79,999"
16 "$80,000 to $89,999"	 
17 "$90,000 to $99,999"	 
18 "$100,000 to $124,999"	 
19 "$125,000 to $149,999"	 
20 "$150,000 to $174,999"	
21 "$175,000 to $199,999"	
22 "$200,000 to $249,999"	
23 "$250,000 to $299,999"	 
24 "$300,000 to $399,999"
25 "$400,000 -";
#delimit cr
lab val pv1990 pv1990
notes pv1990: Missing if group quarters, rental unit, vacant, or /*
*/ not for sale only
notes pv1990: Not available in 1980, 2000

lab var pv2000 "Property value, 2000"
#delimit ;
lab def pv2000
1 " - $10,000"			
2 "$10,000 to $14,999"	
3 "$15,000 to $19,999"	
4 "$20,000 to $24,999"	
5 "$25,000 to $29,999"	
6 "$30,000 to $34,999"	
7 "$35,000 to $39,999"	
8 "$40,000 to $49,999"	
9 "$50,000 to $59,999"	
10 "$60,000 to $69,999"	
11 "$70,000 to $79,999"	
12 "$80,000 to $89,999"	
13 "$90,000 to $99,999"	 
14 "$100,000 to $124,999"
15 "$125,000 to $149,999"
16 "$150,000 to $174,999"	 
17 "$175,000 to $199,999"	 
18 "$200,000 to $249,999"	 
19 "$250,000 to $299,999"	 
20 "$300,000 to $399,999"	
21 "$400,000 to $499,999"	
22 "$500,000 to $749,999"	
23 "$750,000 to $999,999"	 
24 "$1,000,000 -";
#delimit cr
lab val pv2000 pv2000
notes pv2000: Missing if group quarters, rental unit, vacant, or /*
*/ not for sale only
notes pv2000: Not available in 1980, 1990

lab var pv1 "Nom. prop. value, 1980-2000"
#delimit ;
lab def pv1
1 " - $10,000"			
2 "$10,000 to $14,999"	
3 "$15,000 to $19,999"	
4 "$20,000 to $24,999"	
5 "$25,000 to $29,999"	
6 "$30,000 to $34,999"	
7 "$35,000 to $39,999"	
8 "$40,000 to $49,999"	
9 "$50,000 to $59,999"	
10 "$60,000 to $69,999"	
11 "$70,000 to $79,999"	
12 "$80,000 to $89,999"	
13 "$90,000 to $99,999"	 
14 "$100,000 to $124,999"
15 "$125,000 to $149,999"
16 "$150,000 to $199,999"	 
17 "$200,000 -";
#delimit cr
lab val pv1 pv1
notes pv1: Consistent nominal bands for 1980-2000
notes pv1: Top-coding: 1980, $200,000; 1990, $400,000; 2000, $1,000,000
notes pv1: Variables pv1980, pv1990, and pv2000 have finer valuations for /*
*/ each year, but are not consistent over time
notes pv1: Missing if group quarters, rental unit, vacant, or /*
*/ not for sale only

	/* monthly mortgage payment */

gen piti=. /* combined principal, interest, taxes, insurance */

if `year'==1980 {

	/* if not owned or being bought; not single-family house, mobile home or 
		condo; group quarters; vacant */
replace piti=. if mortgag3==0
	/* 	if no regular payment required */
replace piti=0 if mortgag3==0 & ownhouse==1 /* for consistency w/ 1990, 2000 */
	/* 	if tax and insurance are included in mortgage payment */
replace piti=mortgag3 if (2<=mortgag3 & mortgag3~=.) & taxincl==1 & /*
*/ insincl==1
	/* if taxes or insurance are paid separately */
replace piti=mortgag3+(taxinsur/12) /* reported annually 
*/ if (2<=mortgag3 & mortgag3~=.) & (taxincl==2 | insincl==2)
	/* if owned free and clear, set piti to taxes and insurance */
replace piti=(taxinsur/12) if ownhouse==1 & mortgag1==3 /* no mortgage */

	/*	if household has second mortgage, should be included in mortgag3 */

}
if `year'==1990 {

	/* 	for observations where taxes are not included in mortgage payments,
		recode tax variable; use midpoint of nominal bands; convert	annual 
		taxes to monthly rates
	*/
gen taxmnth=0 if rtaxamt==1
replace taxmnth=25/12 if rtaxamt==2
replace taxmnth=75/12 if rtaxamt==3
replace taxmnth=125/12 if rtaxamt==4
replace taxmnth=175/12 if rtaxamt==5
replace taxmnth=225/12 if rtaxamt==6
replace taxmnth=275/12 if rtaxamt==7
replace taxmnth=325/12 if rtaxamt==8
replace taxmnth=375/12 if rtaxamt==9
replace taxmnth=425/12 if rtaxamt==10
replace taxmnth=475/12 if rtaxamt==11
replace taxmnth=525/12 if rtaxamt==12
replace taxmnth=575/12 if rtaxamt==13
replace taxmnth=625/12 if rtaxamt==14
replace taxmnth=675/12 if rtaxamt==15
replace taxmnth=725/12 if rtaxamt==16
replace taxmnth=775/12 if rtaxamt==17
replace taxmnth=825/12 if rtaxamt==18
replace taxmnth=875/12 if rtaxamt==19
replace taxmnth=925/12 if rtaxamt==20
replace taxmnth=975/12 if rtaxamt==21
replace taxmnth=1050/12 if rtaxamt==22
replace taxmnth=1150/12 if rtaxamt==23
replace taxmnth=1250/12 if rtaxamt==24
replace taxmnth=1350/12 if rtaxamt==25
replace taxmnth=1450/12 if rtaxamt==26
replace taxmnth=1550/12 if rtaxamt==27
replace taxmnth=1650/12 if rtaxamt==28
replace taxmnth=1750/12 if rtaxamt==29
replace taxmnth=1850/12 if rtaxamt==30
replace taxmnth=1950/12 if rtaxamt==31
replace taxmnth=2050/12 if rtaxamt==32
replace taxmnth=2150/12 if rtaxamt==33
replace taxmnth=2250/12 if rtaxamt==34
replace taxmnth=2350/12 if rtaxamt==35
replace taxmnth=2450/12 if rtaxamt==36
replace taxmnth=2550/12 if rtaxamt==37
replace taxmnth=2650/12 if rtaxamt==38
replace taxmnth=2750/12 if rtaxamt==39
replace taxmnth=2850/12 if rtaxamt==40
replace taxmnth=2950/12 if rtaxamt==41
replace taxmnth=3050/12 if rtaxamt==42
replace taxmnth=3150/12 if rtaxamt==43
replace taxmnth=3250/12 if rtaxamt==44
replace taxmnth=3350/12 if rtaxamt==45
replace taxmnth=3450/12 if rtaxamt==46
replace taxmnth=3550/12 if rtaxamt==47
replace taxmnth=3650/12 if rtaxamt==48
replace taxmnth=3750/12 if rtaxamt==49
replace taxmnth=3850/12 if rtaxamt==50
replace taxmnth=3950/12 if rtaxamt==51
replace taxmnth=4050/12 if rtaxamt==52
replace taxmnth=4150/12 if rtaxamt==53
replace taxmnth=4250/12 if rtaxamt==54
replace taxmnth=4350/12 if rtaxamt==55
replace taxmnth=4450/12 if rtaxamt==56
replace taxmnth=4500/12 if rtaxamt==57
replace taxmnth=5000/12 if rtaxamt==58
replace taxmnth=6500/12 if rtaxamt==59
replace taxmnth=7500/12 if rtaxamt==60

	/* 	if not owned or being bought; not single-family house, mobile home or 
		condo; group quarters; vacant */
replace piti=. if mortgag3==0
	/* 	if no regular payment required */
replace piti=0 if mortgag3==1 
	/* 	if tax and insurance included in mortgage payment */
replace piti=mortgag3 if (2<=mortgag3 & mortgag3~=.) & taxincl==1 & /*
*/ insincl==1
	/* if tax included, but insurance paid separately */
replace piti=mortgag3+(insamt/12) /* reported annually 
*/ if (2<=mortgag3 & mortgag3~=.) & taxincl==1 & insincl==2
	/* if insurance included, but tax paid separately */
replace piti=mortgag3+(taxmnth) /*
*/ if (2<=mortgag3 & mortgag3~=.) & taxincl==2 & insincl==1
	/* if insurance and tax both paid separately */
replace piti=mortgag3+(taxmnth)+(insamt/12) /* reported annually 
*/ if (2<=mortgag3 & mortgag3~=.) & taxincl==2 & insincl==2
	/* if owned free and clear, set piti to taxes and insurance  */
replace piti=(insamt/12)+taxmnth if ownhouse==1 & mortgag==3 /* no mortgage */

	/* if household has second mortgage */
* replace piti=piti+mortamt2 if mortgag2==1 & (2<=mortamt2 & mortamt2~=.)
replace piti=piti+mortamt if mortgag2==1 & (2<=mortamt & mortamt~=.)
}
if `year'==2000 {

	/* 	for observations where taxes are not included in mortgage payments,
		recode tax variable; use midpoint of nominal bands; convert	annual 
		taxes to monthly rates
	*/
gen taxmnth=0 if taxamt==1
replace taxmnth=25/12 if taxamt==2
replace taxmnth=75/12 if taxamt==3
replace taxmnth=125/12 if taxamt==4
replace taxmnth=175/12 if taxamt==5
replace taxmnth=225/12 if taxamt==6
replace taxmnth=275/12 if taxamt==7
replace taxmnth=325/12 if taxamt==8
replace taxmnth=375/12 if taxamt==9
replace taxmnth=425/12 if taxamt==10
replace taxmnth=475/12 if taxamt==11
replace taxmnth=525/12 if taxamt==12
replace taxmnth=575/12 if taxamt==13
replace taxmnth=625/12 if taxamt==14
replace taxmnth=675/12 if taxamt==15
replace taxmnth=725/12 if taxamt==16
replace taxmnth=775/12 if taxamt==17
replace taxmnth=825/12 if taxamt==18
replace taxmnth=875/12 if taxamt==19
replace taxmnth=925/12 if taxamt==20
replace taxmnth=975/12 if taxamt==21
replace taxmnth=1050/12 if taxamt==22
replace taxmnth=1150/12 if taxamt==23
replace taxmnth=1250/12 if taxamt==24
replace taxmnth=1350/12 if taxamt==25
replace taxmnth=1450/12 if taxamt==26
replace taxmnth=1550/12 if taxamt==27
replace taxmnth=1650/12 if taxamt==28
replace taxmnth=1750/12 if taxamt==29
replace taxmnth=1850/12 if taxamt==30
replace taxmnth=1950/12 if taxamt==31
replace taxmnth=2050/12 if taxamt==32
replace taxmnth=2150/12 if taxamt==33
replace taxmnth=2250/12 if taxamt==34
replace taxmnth=2350/12 if taxamt==35
replace taxmnth=2450/12 if taxamt==36
replace taxmnth=2550/12 if taxamt==37
replace taxmnth=2650/12 if taxamt==38
replace taxmnth=2750/12 if taxamt==39
replace taxmnth=2850/12 if taxamt==40
replace taxmnth=2950/12 if taxamt==41
replace taxmnth=3050/12 if taxamt==42
replace taxmnth=3150/12 if taxamt==43
replace taxmnth=3250/12 if taxamt==44
replace taxmnth=3350/12 if taxamt==45
replace taxmnth=3450/12 if taxamt==46
replace taxmnth=3550/12 if taxamt==47
replace taxmnth=3650/12 if taxamt==48
replace taxmnth=3750/12 if taxamt==49
replace taxmnth=3850/12 if taxamt==50
replace taxmnth=3950/12 if taxamt==51
replace taxmnth=4050/12 if taxamt==52
replace taxmnth=4150/12 if taxamt==53
replace taxmnth=4250/12 if taxamt==54
replace taxmnth=4350/12 if taxamt==55
replace taxmnth=4450/12 if taxamt==56
replace taxmnth=4550/12 if taxamt==57
replace taxmnth=4650/12 if taxamt==58
replace taxmnth=4750/12 if taxamt==59
replace taxmnth=4850/12 if taxamt==60
replace taxmnth=4950/12 if taxamt==61
replace taxmnth=5250/12 if taxamt==62
replace taxmnth=5750/12 if taxamt==63
replace taxmnth=6500/12 if taxamt==64
replace taxmnth=7500/12 if taxamt==65
replace taxmnth=8550/12 if taxamt==66
replace taxmnth=9100/12 if taxamt==64

	/* 	if vacant; group quarters; renter-occupied; owner-occupied and no 
   		mortgage */
replace piti=. if mrt1amt==. 
	/* 	if no regular payment required */
replace piti=0 if mrt1amt==0 
	/* 	if tax and insurance included in mortgage payment */
replace piti=mrt1amt if (1<=mrt1amt & mrt1amt~=.) & (taxincl==1 & insincl==1)
	/* if tax included, but insurance paid separately */
replace piti=mrt1amt+(insamt/12) /* reported annually 
*/ if (1<=mrt1amt & mrt1amt~=.) & taxincl==1 & insincl==2
	/* if insurance included, but tax paid separately */
replace piti=mrt1amt+(taxmnth) /*
*/ if (1<=mrt1amt & mrt1amt~=.) & taxincl==2 & insincl==1
	/* if insurance and tax both paid separately */
replace piti=mrt1amt+(taxmnth)+(insamt/12) /* reported annually 
*/ if (2<=mrt1amt & mrt1amt~=.) & taxincl==2 & insincl==2
	/* if owned free and clear, set piti to taxes and insurance  */
replace piti=(insamt/12)+taxmnth if ownhouse==1 & mortg1==3 /* no mortgage */

	/* if household has second mortgage */
replace piti=piti+mrt2amt if (mortg2==1 | mortg2==2 | mortg2==4) /*
*/  & (1<=mrt2amt & mrt2amt~=.)
}
lab var piti "Principal, interest, taxes, insurance"
notes piti: Includes second mortgage payments
notes piti: Top-code for mortgage: 1980, $1,500; 1990, $2,000; 2000, $3,000
notes piti: Top-code for annual taxes: 1980, $3,000 (combined taxes and /*
*/ insurance); 1990, $7,500; 2000, $9,100
notes piti: Top-code for annual insurance: 1980, $3,000 (combined taxes and /*
*/ insurance); 1990, $1,301; 2000, $2,500
notes piti: Top-code for second mortgages: 1980 (included in total mortgage /*
*/ variable); 1990, $1,001; 2000, $1,100
notes piti: Households above top-codes in 1990 and 2000 assigned state median /*
*/ above the topcode by the Census Bureau

	/* monthly cost of house ownership as a percent of household income */

if `year'==1980 {

	/* calc. annualized ownership cost as percent of annual household income */
gen owncost=100*(ownercst*12)/hhincome if (1<=ownercst & ownercst<=2000) & /*
*/ (1<=hhincome & hhincome<=75000) & ownhouse==1

gen byte mrtg_29=0 if ownhouse==1 & owncost~=. /* owners only */
replace mrtg_29=1 if (0<owncost & owncost<30) & ownhouse==1 /* 1-29% */
gen byte mrtg3049=0 if ownhouse==1 & owncost~=.
replace mrtg3049=1 if (30<=owncost & owncost<50) & ownhouse==1 /* 30-49% */
gen byte mrtg50_=0 if ownhouse==1 & owncost~=.
replace mrtg50_=1 if (50<=owncost & owncost~=.) & ownhouse==1 /* 50-100+% */
}
if `year'==1990 {
gen byte mrtg_29=0 if ownhouse==1 & rnsmocpi~=0 /* owners only */
replace mrtg_29=1 if ownhouse==1 & (0<rnsmocpi & rnsmocpi<30) /* 1-29% */
gen byte mrtg3049=0 if ownhouse==1 & rnsmocpi~=0
replace mrtg3049=1 if ownhouse==1 & (30<=rnsmocpi & rnsmocpi<50) /* 30-49% */
gen byte mrtg50_=0 if ownhouse==1 & rnsmocpi~=0
replace mrtg50_=1 if ownhouse==1 & (50<=rnsmocpi & rnsmocpi~=.) /* 50-100+% */
}
if `year'==2000 {
gen byte mrtg_29=0 if ownhouse==1 & smocapi~=0 /* owners only */
replace mrtg_29=1 if ownhouse==1 & (0<smocapi & smocapi<30) /* 1-29% */
gen byte mrtg3049=0 if ownhouse==1 & smocapi~=0
replace mrtg3049=1 if ownhouse==1 & (30<=smocapi & smocapi<50) /* 30-49% */
gen byte mrtg50_=0 if ownhouse==1 & smocapi~=0
replace mrtg50_=1 if ownhouse==1 & (50<=smocapi & smocapi~=.) /* 50-100+% */
}
lab var mrtg_29 "Mortgage/income <30%"
lab var mrtg3049 "Mortgage/income 30-49%"
lab var mrtg50_ "Mortgage/income 50%+"
notes mrtg_29: Based on Census Bureau's "selected monthly owner costs," /*
*/ which is the sum of real estate taxes, fire and hazard insurance, /*
*/ utilities, fuels, and mortgage (including second mortgages)
notes mrtg_29: For one-family houses only on less than 10 acres and without /*
*/ a commercial establishment or medical office
notes mrtg_29: Excludes owner-occupied condominiums, mobile homes, and trailers
notes mrtg_29: Excludes maintenance and repairs
notes mrtg_29: Comparable to monthly gross rent data
notes mrtg3049: Based on Census Bureau's "selected monthly owner costs," /*
*/ which is the sum of real estate taxes, fire and hazard insurance, /*
*/ utilities, fuels, and mortgage (including second mortgages)
notes mrtg3049: For one-family houses only on less than 10 acres and without /*
*/ a commercial establishment or medical office
notes mrtg3049: Excludes owner-occupied condominiums, mobile homes, and trailers
notes mrtg3049: Excludes maintenance and repairs
notes mrtg3049: Comparable to monthly gross rent data
notes mrtg50_: Based on Census Bureau's "selected monthly owner costs," /*
*/ which is the sum of real estate taxes, fire and hazard insurance, /*
*/ utilities, fuels, and mortgage (including second mortgages)
notes mrtg50_: For one-family houses only on less than 10 acres and without /*
*/ a commercial establishment or medical office
notes mrtg50_: Excludes owner-occupied condominiums, mobile homes, and trailers
notes mrtg50_: Excludes maintenance and repairs
notes mrtg50_: Comparable to monthly gross rent data

	/* piti plus utilities and fuel */

gen pitiplus=.

if `year'==1980 {
gen water12=wtrcost/12 if 1<=wtrcost & wtrcost<=500
gen fuel12=fuelcost/12 if 1<=fuelcost & fuelcost<=2000
egen utilfuel=rsum(eleccost gascost water12 fuel12)
replace pitiplus=piti+utilfuel if ownhouse==1
replace pitiplus=round(pitiplus,1) if ownhouse==1
drop water12 fuel12
}
if `year'==1990 {
gen elec12=eleccost/12 if 3<=eleccost & eleccost~=.
gen gas12=gascost/12 if 3<=gascost & gascost~=.
gen water12=watrcost/12 if 3<=watrcost & watrcost~=.
gen fuel12=fuelcost/12 if 3<=fuelcost & fuelcost~=.
egen utilfuel=rsum(elec12 gas12 water12 fuel12)
replace pitiplus=piti+utilfuel if ownhouse==1
replace pitiplus=round(pitiplus,1) if ownhouse==1
drop elec12 gas12 water12 fuel12
}
if `year'==2000 {
gen elec12=elec/12 if 2<=elec & elec~=.
gen gas12=gas/12 if 2<=gas & gas~=.
gen water12=water/12 if 2<=water & water~=.
gen fuel12=oil/12 if 2<=oil & oil~=.
egen utilfuel=rsum(elec12 gas12 water12 fuel12)
replace pitiplus=piti+utilfuel if ownhouse==1
replace pitiplus=round(pitiplus,1) if ownhouse==1
drop elec12 gas12 water12 fuel12
}
lab var pitiplus "PITI + utilities, fuel"
notes pitiplus: Top-coding in 1980: electricity, $200; gas, $150; water, /*
*/ $500 per year; fuel, $2,000 per year
notes pitiplus: Top-coding in 1990: electricity, $3,100 per year; gas, /*
*/ $2,100 per year; water, $1,000 per year; fuel, $1,900 per year
notes pitiplus: Top-coding in 2000:  electricity, $4,800 per year; gas, /*
*/ $3,000 per year; water, $2,000 per year; fuel, $2,100 per year
notes pitiplus: Top-coded values in 1990 and 2000 replaced with state means /*
*/ above the top-code by the Census Bureau 

	/* monthly rent */

gen rentgr=. /* 	gross rent includes monthly cost of utilities (water,
				electricity, gas) and fuels (oil, coal, kerosene, wood,
				etc.
			*/

if `year'==1980 {
replace rentgr=rent2 if (1<=rent2 & rent2<=999) & ownhouse==0
}
if `year'==1990 {
replace rentgr=rgrent if (1<=rgrent & rgrent~=.) & ownhouse==0
}
if `year'==2000 { 
replace rentgr=grent if (1<=grent & grent<=9999) & ownhouse==0
}
capture lab drop rentgr
lab var rentgr "Gross rent (month)"
notes rentgr: Gross rent includes average monthly cost for utilities /*
*/ (excluding telephone) and fuel
notes rentgr: Not directly comparable with piti, which excludes utility /*
*/ and fuel costs
notes rentgr: Comparable with pitiplus, which includes average monthly /*
*/ utility and fuel costs
notes rentgr: Missing if group quarters or not a rental unit

	/* monthly rent as percent of household income */

if `year'==1980 {

	/* calc. annualized rental cost as percent of annual household income */
gen rentcost=100*(rent2*12)/hhincome if (1<=rent2 & rent2<=999) & /*
*/ (1<=hhincome & hhincome<=75000) & ownhouse==0

gen byte rent_29=0 if ownhouse==0 & rentcost~=. /* renters only */
replace rent_29=1 if (0<rentcost & rentcost<30) & ownhouse==0 /* 1-29% */
gen byte rent3049=0 if ownhouse==0 & rentcost~=.
replace rent3049=1 if ownhouse==0 & (30<=rentcost & rentcost<50) /* 30-49% */
gen byte rent50_=0 if ownhouse==0 & rentcost~=.
replace rent50_=1 if ownhouse==0 & (50<=rentcost & rentcost~=.) /* 50-100+% */
}
if `year'==1990 {
gen byte rent_29=0 if ownhouse==0 & rgrapi~=0 /* renters only */
replace rent_29=1 if (1<=rgrapi & rgrapi<=5) & ownhouse==0 /* 1-29% */
gen byte rent3049=0 if ownhouse==0 & rgrapi~=0
replace rent3049=1 if (6<=rgrapi & rgrapi<=8) & ownhouse==0 /* 30-49% */
gen byte rent50_=0 if ownhouse==0 & rgrapi~=0
replace rent50_=1 if (9<=rgrapi & rgrapi<=14) & ownhouse==0 /* 50-100+% */
}
if `year'==2000 {
gen byte rent_29=0 if ownhouse==0 & grapi~=0 /* renters only */
replace rent_29=1 if (1<=grapi & grapi<30) & ownhouse==0 /* 1-29% */
gen byte rent3049=0 if ownhouse==0 & grapi~=0
replace rent3049=1 if (30<=grapi & grapi<50) & ownhouse==0 /* 30-49% */
gen byte rent50_=0 if ownhouse==0 & grapi~=0
replace rent50_=1 if (50<=grapi & grapi<=101) & ownhouse==0 /* 50-100+% */
}
lab var rent_29 "Rent/income <30%"
lab var rent3049 "Rent/income 30-49%"
lab var rent50_ "Rent/income 50%+"
notes rent_29: Rental cost includes contract rent plus average monthly cost /*
*/ of utilities (excluding telephone) and fuel
notes rent3049: Rental cost includes contract rent plus average monthly cost /*
*/ of utilities (excluding telephone) and fuel
notes rent50_: Rental cost includes contract rent plus average monthly cost /*
*/ of utilities (excluding telephone) and fuel

	/* house, apartment, mobile home or other */

gen byte hamo=.

if `year'==1980 | `year'==1990 {
replace hamo=1 if units1==2 | units1==3 /* house */
replace hamo=2 if 4<=units1 & units1<=9 /* apartment */
replace hamo=3 if units1==1 /* mobile home or trailer */
replace hamo=4 if units1==10 /* other */
replace hamo=. if units1==0 /* group quarters */
drop units1
}
if `year'==2000 {
replace hamo=1 if bldgsz==2 | bldgsz==3 /* house */
replace hamo=2 if 4<=bldgsz & bldgsz<=9 /* apartment */
replace hamo=3 if bldgsz==1 /* mobile home */
replace hamo=4 if bldgsz==10 /* other */
replace hamo=. if bldgsz==0 /* group quarters */
drop bldgsz
}
lab var hamo "Type of housing"
lab def hamo 1 House 2 Apartment 3 "Mobile Home" 4 Other
lab val hamo hamo
notes hamo: Based on units1 in 1990; bldgsz in 2000
notes hamo: Missing if live in group quarters

	/* number of cars */

gen byte numveh=.

if `year'==1980 {
	/* convert autos and trucks variables to numbers */
tempvar c80 t80
gen byte `c80'=0 if autos==1
replace `c80'=1 if autos==2
replace `c80'=2 if autos==3
replace `c80'=3 if autos==4 /* 3 or more cars */
gen byte `t80'=0 if trucks==1
replace `t80'=1 if trucks==2
replace `t80'=2 if trucks==3
replace `t80'=3 if trucks==4 /* 3 or more trucks */
	/* add variables for autos and trucks together */
replace numveh=`c80'+`t80'
}
if `year'==1990 {
replace numveh=0 if autos==1
replace numveh=1 if autos==2
replace numveh=2 if autos==3
replace numveh=3 if autos==4
replace numveh=4 if autos==5
replace numveh=5 if autos==6
replace numveh=6 if autos==7 /* 6 vehicles */
replace numveh=6 if autos==8 /* 7 or more vehicles */
}
if `year'==2000 {
replace numveh=0 if vehicl==0
replace numveh=1 if vehicl==1
replace numveh=2 if vehicl==2
replace numveh=3 if vehicl==3 
replace numveh=4 if vehicl==4 
replace numveh=5 if vehicl==5 
replace numveh=6 if vehicl==6 /* 6 or more vehicles */
}
lab var numveh "Number of vehicles"
#delimit ;
lab def numveh
0 "None"
1 "One"
2 "Two"
3 "Three"
4 "Four"
5 "Five"
6 "Six or more";
#delimit cr
lab val numveh numveh
notes numveh: Top-coding in underlying variables: 1980, 3 cars, 3 trucks; /*
*/ 1990, 7; 2000, 6
notes numveh: Data on vehicles in 1990 and 2000 comparable with combined /*
*/ data for cars and trucks in 1980 (Census Bureau)
notes numveh: Missing if live in group quarters or vacant

/* family structure and employment */

	/* married couple family */

gen byte mcf=.

if `year'==1980 {
rename hhtype hhtype1 /* to avoid variable name conflict below */
replace mcf=0 if (1<=hhtype1 & hhtype1<=4) & (0<=children & children<=4)
replace mcf=1 if hhtype1==1 & (1<=children & children<=4)
}
if `year'==1990 {
replace mcf=0 if 0<=rhhfamtp & rhhfamtp<=22 /* includes group quarters and
									  nonfamily households */
replace mcf=1 if rhhfamtp==1
}
if `year'==2000 {
replace mcf=0 if 0<=hht & hht<=7 /* includes group quarters and
							 nonfamily households */
replace mcf=1 if hht==1
}
lab var mcf "Married couple family"
notes mcf: Based on rhhfamtp in 1990; on hht in 2000
notes mcf: Population includes those living in group quarters and /*
*/ nonfamily households

	/* household type */

gen byte hhtype=.

if `year'==1980 {
	/* not available */
}
if `year'==1990 {
	* other
replace hhtype=7
	* mcf, kids
replace hhtype=1 if rhhfamptp==1 /* mcf */ & /*
*/ & ((1<=rnatadpt & rnatadpt<=28) | (1<=rstpchld & rstpchld<=28))
	* mcf, no kids
replace hhtype=2 if rhhfamptp==1 /* mcf */ & /*
*/ (rnatadpt==0 & rstpchld==0)
	* single female, children (natural, adopted, step)
replace hhtype=3 if (rhhfamtp==3 /* single female w/ family 
*/ | rhhfamtp==21 /* single female, no family, alone
*/ | rhhfamtp==22) /* single female, no family, not alone
*/  & ((1<=rnatadpt & rnatadpt<=28) | (1<=rstpchld & rstpchld<=28))
	* single female, no children
replace hhtype=4 if (rhhfamtp==3 /* single female w/ family 
*/ | rhhfamtp==21 /* single female, no family, alone
*/ | rhhfamtp==22) /* single female, no family, not alone
*/  & (rnatadpt==0 & rstpchld==0)
	* single male, children
replace hhtype=3 if (rhhfamtp==2 /* single male w/ family 
*/ | rhhfamtp==11 /* single male, no family, alone
*/ | rhhfamtp==12) /* single male, no family, not alone
*/  & ((1<=rnatadpt & rnatadpt<=28) | (1<=rstpchld & rstpchld<=28))
	* single male, no children
replace hhtype=4 if (rhhfamtp==2 /* single male w/ family 
*/ | rhhfamtp==11 /* single male, no family, alone
*/ | rhhfamtp==12) /* single male, no family, not alone
*/  & (rnatadpt==0 & rstpchld==0)
}
if `year'==2000 {
	* other
replace hhtype=7
	* mcf, kids (own children, related children
replace hhtype=1 if hht==1 /* mcf */ & /*
*/ & ((1<=noc & noc<=96) | (1<=nrc & nrc<=96))
	* mcf, no kids
replace hhtype=2 if hht==1 /* mcf */ & /*
*/ (noc==0 & nrc==0)
	* single female, children (natural, adopted, step)
replace hhtype=3 if (hht==3 /* single female w/ family 
*/ | hht==6 /* single female, no family, alone
*/ | hht==7) /* single female, no family, not alone
*/  & ((1<=noc & noc<=96) | (1<=nrc & nrc<=96))
	* single female, no children
replace hhtype=3 if (hht==3 /* single female w/ family 
*/ | hht==6 /* single female, no family, alone
*/ | hht==7) /* single female, no family, not alone
*/  & (noc==0 & nrc==0)
	* single male, children
replace hhtype=3 if (hht==2 /* single male w/ family 
*/ | hht==4 /* single male, no family, alone
*/ | hht==5) /* single male, no family, not alone
*/ & ((1<=noc & noc<=96) | (1<=nrc & nrc<=96))
	* single male, no children
replace hhtype=4 if (hht==2 /* single male w/ family 
*/ | hht==4 /* single male, no family, alone
*/ | hht==5) /* single male, no family, not alone
*/ & (noc==0 & nrc==0)
}
lab var hhtype "Household type"
#delimit ;
lab def hhtype
1 "MCF w/ children"
2 "MCF w/o children"
3 "SF w/ children"
4 "SF w/o children"
5 "SM w/ children"
6 "SM w/o children"
7 "Other";
#delimit cr
lab val hhtype hhtype

	/* child or children in household */

gen byte perlt18=.

if `year'==1980 {
replace perlt18=0 /* includes group quarters */
replace perlt18=1 if 1<=children & children<=3
replace perlt18=1 if gq==1 & (0<=age & age<=17) /* children in gq */
}
if `year'==1990 {
replace perlt18=0 if r18undr==0 /* includes group quarters */
replace perlt18=1 if r18undr==1
replace perlt18=1 if gq==1 & (0<=age & age<=17) /* children in gq */
}
if `year'==2000 {
replace perlt18=0 if p18==0 /* includes group quarters */
replace perlt18=1 if 1<=p18 & p18<=97
replace perlt18=1 if gq==1 & (0<=age & age<=17) /* children in gq */
}
lab var perlt18 "Child present"
notes perlt18: Based on r18undr in 1990; on p18 in 2000

	/* senior citizen in household */

gen byte pergt64=.

if `year'==1980 {
/* see 17th_pums_mixed.do */
}
if `year'==1990 {
replace pergt64=0 if r65over==0 /* includes group quarters */
replace pergt64=1 if r65over==1
}
if `year'==2000 {
replace pergt64=0 if p65==0 /* includes group quarters */
replace pergt64=1 if 1<=p65 & p65<=97
}
lab var pergt64 "Senior citizen present"
notes pergt64: Based on r65over in 1990; on p65 in 2000

	/* non-nuclear family members in household */

gen byte nonnuc=.

if `year'==1980 {
/* under construction : in 1980, may need person-level data to generate */
}
if `year'==1990 {
gen nn=1 if relat1==6 | relat1==7 /* grandchild, other relative */
sort state serialno
by state serialno: egen knn=max(nn) /* takes value missing or 1 */
replace knn=0 if knn==.
replace nonnuc=knn
drop nn knn
}
if `year'==2000 {
gen nn=1 if 8<=relate & relate<=16 
/* grandchild; parent-in-law; son-in-law, daughter-in-law;
other relative; brother-in-law,sister-in-law; nephew,niece; 
grandparent; uncle, aunt; cousin */
sort state serialno
by state serialno: egen knn=max(nn) /* takes value missing or 1 */
replace knn=0 if knn==.
replace nonnuc=knn
drop nn knn
}
lab var nonnuc "Household includes non-nuclear members"
notes nonnuc: Based on relat1 in 1990; on TK in 2000

	/* non-relatives in household */

gen byte nonrel=.

if `year'==1980 {
/* see 17th_pums_mixed.do */
}
if `year'==1990 {
gen nr=1 if (8<=relat1 & relat1<=11) | relat1==13
/* roomer, boarder, foster child; housemate, roommate; unmarried 
partner; other nonrelative; other persons in group quarters */
sort state serialno
by state serialno: egen knr=max(nr) /* takes value missing or 1 */
replace knr=0 if knr==.
replace nonrel=knr
drop nr knr
}
if `year'==2000 {
gen nr=1 if (17<=relate & relate<=21) | relate==23
/* roomer/boarder; housemate,roommate; unmarried partner;
foster child; other nonrelative; noninstitutionalized general 
quarter person */
sort state serialno
by state serialno: egen knr=max(nr) /* takes value missing or 1 */
replace knr=0 if knr==.
replace nonrel=knr
drop nr knr
}
lab var nonrel "Household includes non-relatives"
notes nonrel: Based on relat1 in 1990; on TK in 2000

/* income */

	/* total household income */

gen double hhinc=.

if `year'==1980 {
replace hhinc=hhincome
replace hhinc=. if gq==1 /* Census assigns zero to group quarters */ | /*
*/ (1<=vacancy1 & vacancy1<=2) | (1<=vacancy2 & vacancy2<=5) | /*
*/ (1<=vacancy3 & vacancy3<=2) /* Census assigns zero to vacant units */
}
if `year'==1990 {
replace hhinc=rhhinc
replace hhinc=. if gq==1 /* Census assigns zero to group quarters */ | /*
*/ persons==0 /* Census assigns zero to vacant units */
lab var hhinc "Household total income, 1989"
}
if `year'==2000 {
replace hhinc=hinc
replace hhinc=. if gq==1 /* Census assigns zero to group quarters */ | /*
*/ persons==0 /* Census assigns zero to vacant units */
lab var hhinc "Household total income, 1999"
}
notes hhinc: Top-code in 1980: $75,000; 1990: $9,999,999; in 2000: $99,999,999
notes hhinc: Top-coded values assigned state-specific median above /*
*/ the topcode in 1990, 2000
notes hhinc: Bottom-code in 1980: -$9,995; 1990: -$999,999; in 2000: -$59,999
notes hhinc: Based on hhincome in 1980; rhhinc in 1990; on hinc in 2000

	/* total family income */
	
gen double faminc=.

if `year'==1980 {
replace faminc=famincom
replace faminc=. if gq==1 /* Census assigns zero to group quarters */ | /*
*/ (1<=vacancy1 & vacancy1<=2) | (1<=vacancy2 & vacancy2<=5) | /*
*/ (1<=vacancy3 & vacancy3<=2) /* Census assigns zero to vacant units */
}
if `year'==1990 {
replace faminc=rfaminc
replace rfaminc=. if gq==1 /* Census assigns zero to group quarters */ | /*
*/ persons==0 /* Census assigns zero to vacant units */
lab var faminc "Family total income, 1989"
}
if `year'==2000 {
replace faminc=finc
replace faminc=. if gq==1 /* Census assigns zero to group quarters */ | /*
*/ persons==0 /* Census assigns zero to vacant units */
lab var faminc "Family total income, 1999"
}
notes faminc: Top-code in 1980: $75,000; 1990: $9,999,999; in 2000: $99,999,999
notes faminc: Top-coded values assigned state-specific median above /*
*/ the topcode in 1990 and 2000
notes faminc: Bottom-code in 1980: -$9,995; 1990: -$999,999; in 2000: -$59,999
notes faminc: Based on faminc in 1980; rfaminc in 1990; on finc in 2000

	/* house weight */

if `year'==1980 {
gen hhwgt=1 /* no weights in 1980 */
}
if `year'==1990 {
rename houswgt hhwgt
}
if `year'==2000 {
rename hweight hhwgt
}
lab var hhwgt "Housing weight"
notes hhwgt: Integer weight of housing unit
notes hhwgt: PUMS data are unweighted in 1980

/* 
Copyright 2004, 2005 Sarah Gammage and John Schmitt

This program and all programs referenced in it are free software. You
can redistribute the program or modify it under the terms of the GNU
General Public License as published by the Free Software Foundation;
either version 2 of the License, or (at your option) any later version.

This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
GNU General Public License for more details.

You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307
USA.
*/