scalar ub2 = 97

	tabulate anxiety if inrange(iq, lb`i', ub`i') 

foreach c in US UK DE FR  {

local country US UK DE FR

local country US UK DE FR

use un, clear

lrtest unconstr constr

	qui g child`i' = cond(nchild >= `i',int(30 * uniform()) + 1, .)

g p = 10*(uniform()+1) + 40

g final = int(4*hw + 7*invnorm(uniform()))

	merge distid using mathpnl`y'm, sort unique

merge tabrv using f10_4_2 f10_4_3, unique 

merge tabrv using f10_4_3a, update unique

merge rep78 using repstats, uniqusing

merge tabrv using f10_4_2a, update unique

histogram medagebrack, discrete frequency ///
lcolor(black) fcolor(gs15) addlabels ///
addlabopts(mlabposition(6)) xtitle(Upper limits of median age) ///
title(Northeast and North Central States: Median Age) 

label def ind 1 Kappa 2 Lambda 3 Nu 4 Xi 5 Tau 6 Upsilon 7 Phi 8 Chi 9 Psi

label def ind 1 Kappa 2 Lambda 3 Nu 4 Xi 5 Tau 6 Upsilon 7 Phi 8 Chi 9 Psi

use urates, clear

use urates, clear

label variable urbanized "Population in urban areas, %"

replace urbanized = 100 * urbanized

summarize urbanized

describe pop smallpop largepop popsize urbanized

label data "1980 US Census data with population size indicators"

local country US UK DE FR

local country US UK DE FR

use usairquality, clear

use usairquality, clear

serset use

serset use

serset use