Chang ,JX

рассказ и язык

Archive for October, 2009

stata time gap

Posted by changjx on October 30, 2009

seb obs 30
gen t=_n+1
tsset t,d
gen tx=_n+1+t
tsset tx,d
gen gap=tx-t
format %9.0g t
format %9.0g tx
gen gapn = tx-t

Posted in Uncategorized | Tagged: | Leave a Comment »

Latex trick

Posted by changjx on October 29, 2009

\newcommand{\argmax}{\operatornamewithlimits{argmax}}
\def\[#1\]{\begin{align}#1\end{align}}
\newcommand{\ve}{\varepsilon}
\newcommand{\nf}{\textrm}
\newcommand{\argmax}{\operatornamewithlimits{argmax}}

Posted in Uncategorized | Tagged: | Leave a Comment »

do stata in mata

Posted by changjx on October 26, 2009

mata: mata clear
sysuse auto, clear
reg price weight trunk
mata
b=st_matrix(“e(b)”)’
b
end
******************************
mata: mata clear
sysuse auto, clear
gen constant=1
capt ssc inst tomata
di in red _rc
tomata price constant weight trunk
mata
X=(weight, trunk, constant)
y=price
b=invsym(X’*X)*X’*y
b
//for comparison:
mata stata reg price weight trunk
end
*****************stata mata interact**********************
mata: mata clear
mata:
void underbar(string scalar mat, real scalar thresh)
{
real scalar i
i = rows(st_matrix(mat))
st_matrix(“B”, (st_matrix(mat) :< thresh) :/ st_matrix(mat) :* (J(i,i,
1) – I(i)) )
}
end
sysuse auto,clear
qui reg price trunk weight length turn
mat list e(V)
mata: underbar(“e(V)”, 10)
mat list B
******************************************************

mata: mata clear

sysuse auto, clear

reg price weight trunk

mata

b=st_matrix(“e(b)”)’

b

end

******************************

mata: mata clear

sysuse auto, clear

gen constant=1

capt ssc inst tomata

di in red _rc

tomata price constant weight trunk

mata

X=(weight, trunk, constant)

y=price

b=invsym(X’*X)*X’*y

b

//for comparison:

mata stata reg price weight trunk

end

*****************stata mata interact**********************

mata: mata clear

mata:

void underbar(string scalar mat, real scalar thresh)

{

real scalar i

i = rows(st_matrix(mat))

st_matrix(“B”, (st_matrix(mat) :< thresh) :/ st_matrix(mat) :* (J(i,i,

1) – I(i)) )

}

end

sysuse auto,clear

qui reg price trunk weight length turn

mat list e(V)

mata: underbar(“e(V)”, 10)

mat list B

******************************************************

Posted in Econ, Stata | Leave a Comment »

latex bold symbol

Posted by changjx on October 26, 2009

use \pmb instead of \boldsymbol

Posted in Uncategorized | Tagged: | Leave a Comment »