label var so2 "Sulphur dioxide, mg/m3"

do m_mycuecrit.mata

do m_mycuecrit.mata

type m_mygmm2s.mata

do m_mygmm2s.mata

do m_mygmm2s.mata

do m_mygmmcue.mata

do m_mygmmcue.mata

type m_myomega.mata

do m_myomega.mata

do m_myomega.mata

savedresults compare ivhomo 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

generate byte newengland = inlist(state, "Connecticut", "Maine", ///
"Massachusetts", "New Hampshire", "Rhode Island", "Vermont")

keep make model price

l make model

list make price weight length in 1/10, sep(0) subvarname

	local temp = make[`i']

  cap noi spinst_itsp_wrk makematrix http://fmwww.bc.edu/repec/bocode/m

ssc install makematrix

do makevec.mata

do makevec.mata

file write `hh' "\hline & Mean & Mean & \$N\$ & \$t\$ test & Mann--Whitney \\" _newline

label var marr "Marriages, '000"

outtable using ch9.03t, mat(allbeta) replace ///
format(%9.4f) center caption("Grunfeld company estimates")

pctrange regday specneed bilingua occupday tot_day tchratio, mat(MCAS) noprint

mkmat hw mt final total, mat(grades)

  cap noi spinst_itsp_wrk mat2txt http://fmwww.bc.edu/repec/bocode/m

type calcextrema.mata

type avgper.mata

type gbpval.mata

type gbpval3.mata

type makevec.mata

mata: pointer2(2, 3)

type shufstr.mata

type aggreg2.mata

type loadcty.mata

type compcty2.mata

type mm_suregub.mata

type m_myomega.mata

type m_mycuecrit.mata

do centerv.mata

/*
program aggreg, rclass
	version 10.1
	syntax varlist(max=1 numeric) [if] [in], per(integer) [func(string)]
	marksample touse
	quietly count if `touse'
	if `r(N)' == 0 {
	        error 2000
	}
* validate per versus selected sample
	if `per' <= 0 | `per' >= `r(N)' {
		display as error "per must be > 0 and < N of observations."
		error 198
	}
	if mod(`r(N)',`per' != 0) {
		display as error "N of observations must be a multiple of per."
		error 198
	}
* validate func option; default is average (code A)
    local ops A S F L
    local opnames average sum first last
    if "`func'" == "" {
    	local op "A"
    }
    else {
    	local nop : list posof "`func'" in opnames
    	if !`nop' {
    		display as err "Error: func must be chosen from `opnames'"
    		error 198
    	}
    	local op : word `nop' of `ops'
    }
* validate the new varname
	local newvar = "`varlist'`op'`per'"
     quietly generate `newvar' = .
* pass the varname and newvarname to Mata
	mata: aggreg("`varlist'", "`newvar'",  `per', "`op'", "`touse'")
	end
*/
do aggreg.mata

do gbpval2.mata

do myvecstr.mata

do myvecsub.mata

do makevec.mata

do mf_flipmat.mata

do indcorr.mata

do country.mata

do compcty.mata

do loadcty.mata

do compcty2.mata

do mf_nneighbor.mata

* webuse grunfeld, clear
do mm_suregub.mata

do m_myomega.mata

do m_mycuecrit.mata

do m_myomega.mata

do m_mycuecrit.mata

mata: mata clear

mata: mata clear

mata: mata clear

mata: mata clear

mata: mata clear

set matastrict on

nneighbor pop, y(so2) matchobs(mo1) matchval(mv1)

nneighbor pop temp wind precip days, y(so2) matchobs(mo5) matchval(mv5)

nneighbor pop temp wind, y(so2) matchobs(mo3) matchval(mv3)

use mathpnl1992, clear

append using mathpnl1994

append using mathpnl1998

use mathpnl_long, clear

use mathpnl_long, clear

use mathpnl_long, clear

save mathpnl_long2, replace

merge distid using mathpnl_state, sort uniqusing 

use mathpnl_wide, clear

use mathpnl_wide, clear