type gbpval.mata

type gbpval2.mata

do gbpval2.mata

do gbpval3.mata

do gbpval.mata

use gdp4cty, clear

use gdp4cty, clear

use gdp4cty, clear

use gdp4cty, clear

describe gdp*_*  

	generate double `c'gdp = gdp`cc'

		rename gdp`i'_`y' gdp`cname'_`y'

		rename gdp`i'_`y' gdp`cname'_`y'

		summarize gdp`i'_`y' 

levelsof cty, local(gdplist)

	qui replace gdptot = gdptot + `c'

summarize gdpUS*

recode cc `rc', gen(newcc)

tscollap (last) aireoq=air (sum) airsum=air (mean) airavg=air, to(q) gen(yq)

centervars price mpg if foreign, gen(cf_) double

encode state, generate(stateid)

mvsumm LD.linvest, generate(invrate) stat(mean) window(3) end

quietly ivreg2 lw s expr tenure rns smsa  (iq=med kww age mrt), gmm2s robust

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

flipmat grades

mat rownames grades = `idnr'

mat list grades, noheader

*mata: mata clear
*program drop _all
// use http://fmwww.bc.edu/ec-p/data/hayashi/griliches76.dta, clear
use griliches76, clear

use grun1, clear

mat list table1, format(%9.3f) ti("Grunfeld company statistics")

/*
*! onespell 1.1.1  CFBaum  13jan2005
* locate units with internal gaps in varlist and zap all but longest spell
program onespell, rclass
	version 8.2
	syntax varlist(numeric) [if] [in], Saving(string) [ REPLACE NOIsily]
	preserve
	quietly tsset
	local pv  "`r(panelvar)'"
	local tv "`r(timevar)'"
	summarize `pv', meanonly
	local n1 = r(N)
	tsfill
	marksample touse
	tempvar testgap spell seq end maxspell keepspell wantspell

	local sss "quietly"
	if "`noisily'" == "noisily" {
        local sss "noisily"
     }

	`sss' {
* testgap is panelvar if obs is usable, 0 otherwise	
	generate `testgap' = cond(`touse',`pv',.)
	tsspell `testgap' if `testgap'<., spell(`spell') seq(`seq') end(`end')
	drop if `spell'==0 | `touse'==0
* if `spell' > 1 for a unit, there are gaps in usable data
* calculate max length spell for each unit and identify
* that spell as the one to be retained
	egen `maxspell'= max(`seq'), by(`pv')
	bysort `pv': generate `keepspell' = cond(`seq'==`maxspell',`spell',0)
	egen `wantspell' = max(`keepspell'), by(`pv')
* in case of ties, latest spell of max length is selected
	list `pv' `tv' `spell' `seq' `maxspell' `keepspell' `wantspell', sepby(`pv') 
	summarize `spell' `wantspell'
	keep if `wantspell' == `spell'
	summarize `pv', meanonly
	local n2 = r(N)
	drop __*
	}
	display _n "Observations removed: " `n1'-`n2'
	save `saving',`replace'
	restore
end
*/
webuse grunfeld, clear

use grunfeldavg

use grunfeldavg, clear

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)