ml model d0 mynormal_lf_d0_c1 (a:) (mu: price = weight turn) /lnsigma, ///
maximize nolog diparm(lnsigma, exp label("sigma")) 

* label var reg "Census region"
label data "1980 Census data for NE and NC states"

note: Subset of Census data, prepared on TS for Chapter 2 

label data ""

display _newline "Datafile: `c(filename)' (N=`c(N)', k=`c(k)')  as of `c(filedate)'"

rename daten caldate

quietly by month: generate adjchange = sum(dDTB3/sqrt(dcal))

label var dDTB3 "Daily change in 3-mo Treasury rate"

tsline dDTB3, yline(0) scheme(s2mono)

local country US UK DE FR

local country US UK DE FR

local country US UK DE FR

label def satis 1 "Poor" 2 "Good" 3 "Better" 4 "Best"

label def tlab 0 "below mean temp" 1 "above mean temp"

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

local inv2 "Defendant (0/1)"

label define netmarr2xc 1 "marr > 2 divr" 2 "marr <= 2 divr"

label define popsize 0 "<1m" 1 "1-4m" 2 "4-8m" 3 ">8m"

note popsize: value label popsize defined for this variable

qui generate deld = date - date[_n-1]

nl ces @ lnoutput capital labor, parameters(b0 rho delta) ///
 initial(b0 0 rho 1 delta 0.5)

local inv Audit Defendant Ins\_Trading Departed Age Tenure Voting\_Share

	display _newline "Dependent variable `c':"

savedresults compare ivrobust e(), include(macros: depvar scalar: N j matrix: b V) /// 
             tol(1e-7) verbose

savedresults compare ivreg2cue e(), include(macros: depvar scalar: N j matrix: b V) ///
             tol(1e-4) verbose

// summarize, detail computes percentiles of stratio
quietly summarize stratio if e(sample), detail 

qui summarize so2, detail

		summarize invest if company == `c', detail

    di in smcl "{p}{err}At least one of the packages was not able to be " ///
               "installed.  Try typing {cmd:spinst_itsp} again.  If that " ///
               "fails, look at the output of the error message above to " ///
               "diagnose the problem.  Perhaps you have an out-of-date " ///
               "version of the package already installed.  In that case, " ///
               "type {cmd:adoupdate} to update it."

bootstrap diff=(r(mu_1)-r(mu_2)), nodots: ttest bp, by(sex)

label var so2 "Sulphur dioxide, mg/m3"

serset dir

serset dir

display _newline "Current working directory: `cwd'"

ml display

scatter expp1992 expp1998 sal, m(oh x) scheme(s2mono) ///
  title("School District Expenditures vs. Average Salaries")  ///
  xtitle("Average salaries") ytitle("Expenditures") ///
  legend(position(5) ring(0) label(1 "1992") label(2 "1998"))

type testall.do

do appB.9.do

do ch10.1.do

do ch10.4.do

do ch11.01.do

do ch12.2.do

do ch12.4.do

do ch12.6.do

do ch13.09a.do

do ch13.9.do

do ch14.1.do

do ch14.3.do

do ch14.8.do

do ch2.0.do

do ch2.4.do

do ch2.6.do

do ch3.01.do

do ch4.1.do

do ch4.4.do

do ch6.5.do

do ch7.00.do

do ch8.2.do

do ch9.06.do

do ch9.09.do

do make_bpress.do

do teststruct.do

* get the month number
generate mnr = month(dofm(ym))

generate qtr = qofd(dofm(ym))

textbarplot compname beta_k, recast(dot) xline(0) ///
scheme(s2mono) name(fig9_2d,replace)

	generate double lngdp`i' = log(gdp`i')

	generate double lngdp`c' = log(`c'gdp)

	generate double `c'gdp = gdp`cc'

generate double rootage = age * root2

bysort qm: generate double futret = log(settle) - log(sprev) if qmdy == rollover

svmat double retcorr

qui generate ds2 = dscorrd^2