比如现在Stata中有一个矩阵A和一个数值型标量a,那么B=st_matrix("A")和b=st_scalar("a")就会在mata中生成对应的矩阵和数值型标量。如果这个标量是字符型的,则需要命令st_strscalar(),如果仍然使用st_numscalar()的话,该标量进入mata之后,就变成一个0行0阶的矩阵。 最后再看看宏如何进入mata: st_local()...
生成对角矩阵及逆矩阵 mat b=(2,3,5) mat D=diag(b) mat list D mat f=inv(D) mat list f 特征根与特征向量 mat X=(2,2,-2\2,5,-4\-2,-4,5) mata V=st_matrix("X") mata p=. mata l=. mata eigensystem(V,p,l) mata p mata l发布...
x=st_data(.,"x*") c=correlation(x) n=rows(c) b=strofreal(lowertriangle(c)+uppertriangle(st_matrix("r(Rho)")),"%9.3f") p=st_matrix("r(P)") for (i=2; i<=n; i++) { for (j=1; j<=i-1; j++) { p[i,j]=2*ttail(rows(x)-2,abs(c[i,...
例如,假设一个二次型趋势: mata st_matrix("deltaquad",0.024 *((-4::3):-(-1)):^2)pretrends,numpre(3)b(beta)v(sigma)deltatrue(deltaquad)coefplot image.png 详细的统计量展示可以用: matlistr(results)returnlist
在Stata 15 中,可通过命令 spmatrix 来方便地定义、导入、导出空间权重矩阵,以及进行相应的操作。比如, .spmatrix create contiguity W 此命令将根据相邻关系(contiguity)来创建空间权重矩阵。如想使用距离的倒数来定义空间权重矩阵,可输入以下命令: .spmatrix create idistance W ...
比如目前Stata中有一台矩阵A和一台数值型标量a,那么B=st_matrix("A")和b=st_scalar("a")就会在mata中生成对应的矩阵和数值型标量。如果这个标量是字符型的,则需要命令st_strscalar(),如果仍然使用st_numscalar()的话,该标量进入mata之后,就变成一台0行0阶的矩阵。
st: matrix in Stata FromAna Venâncio <avenancio@gmail.com> Tostatalist@hsphsun2.harvard.edu Subjectst: matrix in Stata DateSat, 16 Feb 2013 20:39:02 +0000
spearmanx*n matax=st_data(.,"x")c=correlation(x)n=rows(c)b=strofreal(lowertriangle(c)+uppertriangle(st_matrix("r(Rho)")),"%9.3f")p=st_matrix("r(P)")for (i=2; i<=n; i++) {for (j=1; j<=i-1; j++) {p[i,j]=2ttail(rows(x)-2,abs(c[i,j]/sqrt((1-c[i,...
b=strofreal(lowertriangle©+uppertriangle(st_matrix(“r(Rho)”)),“%9.3f”) p=st_matrix(“r§”) for (i=2; i<=n; i++) { for (j=1; j<=i-1; j++) { p[i,j]=2ttail(rows(x)-2,abs(c[i,j]/sqrt((1-c[i,j]^2)/(rows(x)-2))) b...
每组随机抽取一个年份 keep Year save matchyear.dta, replace mkmat Year,matrix(sampleyear) ...